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/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.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/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.compute.c')
-rw-r--r-- | src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.compute.c | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..c30526f --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.compute.c @@ -0,0 +1,24 @@ +fp t25; +fp t23; +fp t29; +fp t21; +fp t22; +fp t28; +fp t24; +fp t20; +fp t19; + t25 = x*x; + t23 = t25*t25; + t29 = RATIONAL(1.0,24.0)*t23+RATIONAL(-1.0,24.0)*t25; + t21 = RATIONAL(-1.0,6.0); + t22 = RATIONAL(2.0,3.0); + t28 = t21*t23+t22*t25; + t24 = x*t25; + t20 = RATIONAL(1.0,12.0); + t19 = RATIONAL(-1.0,12.0); + coeffs_I->coeff_m2 = t20*x+t19*t24+t29; + coeffs_I->coeff_m1 = RATIONAL(-2.0,3.0)*x+RATIONAL(1.0,6.0)*t24+t28; + coeffs_I->coeff_0 = RATIONAL(-5.0,4.0)*t25+RATIONAL(1.0,1.0)+RATIONAL(1.0 +,4.0)*t23; + coeffs_I->coeff_p1 = t21*t24+t22*x+t28; + coeffs_I->coeff_p2 = t20*t24+t19*x+t29; |