diff options
author | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-09-02 10:36:32 +0000 |
---|---|---|
committer | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-09-02 10:36:32 +0000 |
commit | 80df27de110b6e739fda09ba4d9eae9fdbefc342 (patch) | |
tree | 3c17be6653cf3bd67a40fdf79f0821ff1afdd333 /src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c | |
parent | 2248cabada859337c9945616de9a3d612d8f2199 (diff) |
This is a major reworking of this interpolator.
The interpolation code in template.c is now split up into subfunctions
(in common/, Lagrange/, and Hermite/), so it should now compile in
finite cpu/memory even with optimization.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@115 df1f8a13-aa1d-4dd4-9681-27ded5b42416
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c')
-rw-r--r-- | src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c | 52 |
1 files changed, 52 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c new file mode 100644 index 0000000..e9208a2 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c @@ -0,0 +1,52 @@ +fp t274; +fp t273; +fp t272; +fp t271; +fp t270; +fp t269; +fp t268; +fp t267; +fp t266; +fp t265; +fp t264; +fp t263; +fp t262; +fp t261; +fp t260; +fp t259; +fp t258; +fp t257; + t274 = RATIONAL(-3.0,20.0); + t273 = RATIONAL(3.0,20.0); + t272 = RATIONAL(2.0,25.0); + t271 = RATIONAL(-1.0,20.0); + t270 = RATIONAL(1.0,20.0); + t269 = RATIONAL(-1.0,50.0); + t268 = RATIONAL(-1.0,100.0); + t267 = RATIONAL(7.0,100.0); + t266 = RATIONAL(-9.0,100.0); + t265 = RATIONAL(-13.0,100.0); + t264 = t274*x; + t263 = t273*x; + t262 = t271*x; + t261 = t274*y; + t260 = t270*x; + t259 = t273*y; + t258 = t271*y; + t257 = t270*y; + coeffs_dxy->coeff_m1_m1 = t264+RATIONAL(6.0,25.0)+t261; + coeffs_dxy->coeff_0_m1 = t263+t269+t258; + coeffs_dxy->coeff_p1_m1 = t257+t265+t263; + coeffs_dxy->coeff_p2_m1 = t259+t264+t266; + coeffs_dxy->coeff_m1_0 = t259+t262+t269; + coeffs_dxy->coeff_0_0 = t260+RATIONAL(-1.0,25.0)+t257; + coeffs_dxy->coeff_p1_0 = t258+t268+t260; + coeffs_dxy->coeff_p2_0 = t261+t262+t267; + coeffs_dxy->coeff_m1_p1 = t259+t265+t260; + coeffs_dxy->coeff_0_p1 = t262+t257+t268; + coeffs_dxy->coeff_p1_p1 = t262+RATIONAL(3.0,50.0)+t258; + coeffs_dxy->coeff_p2_p1 = t261+t260+t272; + coeffs_dxy->coeff_m1_p2 = t266+t261+t263; + coeffs_dxy->coeff_0_p2 = t267+t258+t264; + coeffs_dxy->coeff_p1_p2 = t264+t257+t272; + coeffs_dxy->coeff_p2_p2 = t263+t259+RATIONAL(-3.0,50.0); |