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

1 files changed, 17 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;