diff options
| author | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-08-28 11:42:25 +0000 |
|---|---|---|
| committer | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-08-28 11:42:25 +0000 |
| commit | 99f3fe1eceecccb7d9487e8ecdb171c7d6865428 (patch) | |
| tree | 23035f3a6053e32549b2537be7db6dcc26855a87 /src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3 | |
| parent | 288605bf7151fda07125fd0bc0e3571ed8bb5516 (diff) | |
add Maple code and C skeletons for Hermite polynomial interpolation
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@101 df1f8a13-aa1d-4dd4-9681-27ded5b42416
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3')
3 files changed, 34 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeff-I.compute.c new file mode 100644 index 0000000..758efe3 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeff-I.compute.c @@ -0,0 +1,17 @@ +fp t10, + t9, + t8, + t7; + t10 = x*x; + t9 = x*t10; + t8 = RATIONAL(1.0,12.0); + t7 = RATIONAL(-1.0,12.0); + coeff_I_m2 = RATIONAL(-1.0,6.0)*t10+(x+t9)*t8; + coeff_I_m1 = RATIONAL(-2.0,3.0)*x+RATIONAL(5.0,4.0)*t10+RATIONAL(-7.0, +12.0)*t9; + coeff_I_0 = RATIONAL(-7.0,3.0)*t10+RATIONAL(1.0,1.0)+RATIONAL(4.0,3.0)*t9 +; + coeff_I_p1 = RATIONAL(-4.0,3.0)*t9+RATIONAL(2.0,3.0)*x+RATIONAL(5.0,3.0)* +t10; + coeff_I_p2 = RATIONAL(7.0,12.0)*t9+t7*x+RATIONAL(-1.0,2.0)*t10; + coeff_I_p3 = t7*t9+t8*t10; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeff-dx.compute.c new file mode 100644 index 0000000..4300a83 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeff-dx.compute.c @@ -0,0 +1,11 @@ +fp t12; + t12 = x*x; + coeff_dx_m2 = RATIONAL(-1.0,3.0)*x+RATIONAL(1.0,12.0)+RATIONAL(1.0,4.0)* +t12; + coeff_dx_m1 = RATIONAL(-2.0,3.0)+RATIONAL(5.0,2.0)*x+RATIONAL(-7.0,4.0)* +t12; + coeff_dx_0 = RATIONAL(-14.0,3.0)*x+RATIONAL(4.0,1.0)*t12; + coeff_dx_p1 = RATIONAL(-4.0,1.0)*t12+RATIONAL(10.0,3.0)*x+RATIONAL(2.0, +3.0); + coeff_dx_p2 = RATIONAL(7.0,4.0)*t12+RATIONAL(-1.0,12.0)-x; + coeff_dx_p3 = RATIONAL(1.0,6.0)*x+RATIONAL(-1.0,4.0)*t12; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeff-dxx.compute.c new file mode 100644 index 0000000..79fc889 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeff-dxx.compute.c @@ -0,0 +1,6 @@ + coeff_dxx_m2 = RATIONAL(-1.0,3.0)+RATIONAL(1.0,2.0)*x; + coeff_dxx_m1 = RATIONAL(-7.0,2.0)*x+RATIONAL(5.0,2.0); + coeff_dxx_0 = RATIONAL(-14.0,3.0)+RATIONAL(8.0,1.0)*x; + coeff_dxx_p1 = RATIONAL(10.0,3.0)+RATIONAL(-8.0,1.0)*x; + coeff_dxx_p2 = RATIONAL(7.0,2.0)*x+RATIONAL(-1.0,1.0); + coeff_dxx_p3 = RATIONAL(1.0,6.0)+RATIONAL(-1.0,2.0)*x; |