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@116 df1f8a13-aa1d-4dd4-9681-27ded5b42416

1 files changed, 21 insertions, 0 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-dxx.compute.c
new file mode 100644
index 0000000..cc8dc91
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-dxx.compute.c
@@ -0,0 +1,21 @@
+fp t37;
+fp t36;
+fp t35;
+fp t34;
+fp t33;
+      t37 = x*x;
+      t36 = x*t37;
+      t35 = RATIONAL(-1.0,4.0);
+      t34 = RATIONAL(-25.0,6.0);
+      t33 = t35*x;
+      coeffs_dxx->coeff_m2 = t35+t33+RATIONAL(3.0,2.0)*t37+RATIONAL(-5.0,6.0)*
+t36;
+      coeffs_dxx->coeff_m1 = RATIONAL(13.0,6.0)+t33+RATIONAL(-7.0,1.0)*t37+
+RATIONAL(25.0,6.0)*t36;
+      coeffs_dxx->coeff_0 = RATIONAL(-25.0,3.0)*t36+RATIONAL(13.0,1.0)*t37+t34+
+RATIONAL(5.0,2.0)*x;
+      coeffs_dxx->coeff_p1 = RATIONAL(3.0,1.0)+RATIONAL(-7.0,2.0)*x+RATIONAL(
+-12.0,1.0)*t37+RATIONAL(25.0,3.0)*t36;
+      coeffs_dxx->coeff_p2 = RATIONAL(11.0,2.0)*t37+t34*t36+RATIONAL(7.0,4.0)*x
++RATIONAL(-11.0,12.0);
+      coeffs_dxx->coeff_p3 = RATIONAL(5.0,6.0)*t36+RATIONAL(1.0,6.0)+t33-t37;