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/3d.coeffs/3d.cube.order2.smooth0/coeffs-dy.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/3d.coeffs/3d.cube.order2.smooth0/coeffs-dy.compute.c')
-rw-r--r-- | src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeffs-dy.compute.c | 59 |
1 files changed, 59 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeffs-dy.compute.c new file mode 100644 index 0000000..e3774d8 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeffs-dy.compute.c @@ -0,0 +1,59 @@ +fp t110; +fp t122; +fp t121; +fp t113; +fp t108; +fp t120; +fp t112; +fp t106; +fp t119; +fp t105; +fp t118; +fp t117; +fp t116; +fp t107; +fp t115; +fp t109; + t110 = RATIONAL(1.0,9.0)*y; + t122 = RATIONAL(-1.0,18.0)+t110; + t121 = t110+RATIONAL(1.0,18.0); + t113 = RATIONAL(1.0,12.0); + t108 = t113*x; + t120 = t108+t122; + t112 = RATIONAL(-1.0,12.0); + t106 = t112*x; + t119 = t106+t122; + t105 = t112*z; + t118 = t105+t122; + t117 = t106+t121; + t116 = t108+t121; + t107 = t113*z; + t115 = t107+t121; + t109 = RATIONAL(-2.0,9.0)*y; + coeffs_dy->coeff_m1_m1_m1 = t107+t120; + coeffs_dy->coeff_0_m1_m1 = t107+t122; + coeffs_dy->coeff_p1_m1_m1 = t107+t119; + coeffs_dy->coeff_m1_0_m1 = t109; + coeffs_dy->coeff_0_0_m1 = t109; + coeffs_dy->coeff_p1_0_m1 = t109; + coeffs_dy->coeff_m1_p1_m1 = t105+t117; + coeffs_dy->coeff_0_p1_m1 = t105+t121; + coeffs_dy->coeff_p1_p1_m1 = t105+t116; + coeffs_dy->coeff_m1_m1_0 = t120; + coeffs_dy->coeff_0_m1_0 = t122; + coeffs_dy->coeff_p1_m1_0 = t119; + coeffs_dy->coeff_m1_0_0 = t109; + coeffs_dy->coeff_0_0_0 = t109; + coeffs_dy->coeff_p1_0_0 = t109; + coeffs_dy->coeff_m1_p1_0 = t117; + coeffs_dy->coeff_0_p1_0 = t121; + coeffs_dy->coeff_p1_p1_0 = t116; + coeffs_dy->coeff_m1_m1_p1 = t108+t118; + coeffs_dy->coeff_0_m1_p1 = t118; + coeffs_dy->coeff_p1_m1_p1 = t106+t118; + coeffs_dy->coeff_m1_0_p1 = t109; + coeffs_dy->coeff_0_0_p1 = t109; + coeffs_dy->coeff_p1_0_p1 = t109; + coeffs_dy->coeff_m1_p1_p1 = t106+t115; + coeffs_dy->coeff_0_p1_p1 = t115; + coeffs_dy->coeff_p1_p1_p1 = t108+t115; |