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

6 files changed, 97 insertions, 91 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c
deleted file mode 100644
index b189aa6..0000000
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c
+++ /dev/null
@@ -1,38 +0,0 @@
-fp t103,
-   t107,
-   t106,
-   t102,
-   t113,
-   t112,
-   t96,
-   t101,
-   t100,
-   t111,
-   t105,
-   t99,
-   t98,
-   t97;
-      t103 = x*x;
-      t107 = x*t103;
-      t106 = t107*t107;
-      t102 = t103*t103;
-      t113 = RATIONAL(-1.0,144.0)*t102+RATIONAL(1.0,720.0)*t106+RATIONAL(1.0,
-180.0)*t103;
-      t112 = RATIONAL(-1.0,120.0)*t106+RATIONAL(-3.0,40.0)*t103+RATIONAL(1.0,
-12.0)*t102;
-      t96 = RATIONAL(-13.0,48.0);
-      t101 = RATIONAL(3.0,4.0);
-      t100 = RATIONAL(1.0,48.0);
-      t111 = t96*t102+t100*t106+t101*t103;
-      t105 = t107*t103;
-      t99 = RATIONAL(-1.0,60.0);
-      t98 = RATIONAL(1.0,60.0);
-      t97 = RATIONAL(-1.0,48.0);
-      coeff_I_m3 = t99*x+t100*t107+RATIONAL(-1.0,240.0)*t105+t113;
-      coeff_I_m2 = RATIONAL(3.0,20.0)*x+RATIONAL(-1.0,6.0)*t107+t98*t105+t112;
-      coeff_I_m1 = RATIONAL(-3.0,4.0)*x+RATIONAL(13.0,48.0)*t107+t97*t105+t111;
-      coeff_I_0 = RATIONAL(7.0,18.0)*t102+RATIONAL(-49.0,36.0)*t103+RATIONAL(
-1.0,1.0)+RATIONAL(-1.0,36.0)*t106;
-      coeff_I_p1 = t96*t107+t101*x+t100*t105+t111;
-      coeff_I_p2 = RATIONAL(1.0,6.0)*t107+RATIONAL(-3.0,20.0)*x+t99*t105+t112;
-      coeff_I_p3 = t97*t107+t98*x+RATIONAL(1.0,240.0)*t105+t113;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c
deleted file mode 100644
index ab1b33f..0000000
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c
+++ /dev/null
@@ -1,31 +0,0 @@
-fp t123,
-   t125,
-   t126,
-   t129,
-   t127,
-   t133,
-   t132,
-   t131;
-      t123 = RATIONAL(-3.0,20.0);
-      t125 = x*x;
-      t126 = x*t125;
-      t129 = t125*t125;
-      t127 = x*t129;
-      t133 = RATIONAL(-1.0,20.0)*t127+t123*x+RATIONAL(1.0,3.0)*t126;
-      t132 = RATIONAL(1.0,90.0)*x+RATIONAL(-1.0,36.0)*t126+RATIONAL(1.0,120.0)*
-t127;
-      t131 = RATIONAL(-13.0,12.0)*t126+RATIONAL(3.0,2.0)*x+RATIONAL(1.0,8.0)*
-t127;
-      coeff_dx_m3 = RATIONAL(-1.0,60.0)+RATIONAL(1.0,16.0)*t125+RATIONAL(-1.0,
-48.0)*t129+t132;
-      coeff_dx_m2 = RATIONAL(3.0,20.0)+RATIONAL(1.0,12.0)*t129+RATIONAL(-1.0,
-2.0)*t125+t133;
-      coeff_dx_m1 = RATIONAL(-3.0,4.0)+RATIONAL(13.0,16.0)*t125+RATIONAL(-5.0,
-48.0)*t129+t131;
-      coeff_dx_0 = RATIONAL(-1.0,6.0)*t127+RATIONAL(14.0,9.0)*t126+RATIONAL(
--49.0,18.0)*x;
-      coeff_dx_p1 = RATIONAL(5.0,48.0)*t129+RATIONAL(3.0,4.0)+RATIONAL(-13.0,
-16.0)*t125+t131;
-      coeff_dx_p2 = t123+RATIONAL(1.0,2.0)*t125+RATIONAL(-1.0,12.0)*t129+t133;
-      coeff_dx_p3 = RATIONAL(1.0,60.0)+RATIONAL(1.0,48.0)*t129+RATIONAL(-1.0,
-16.0)*t125+t132;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c
deleted file mode 100644
index 46979da..0000000
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c
+++ /dev/null
@@ -1,22 +0,0 @@
-fp t144,
-   t143,
-   t141,
-   t150,
-   t149,
-   t148,
-   t145;
-      t144 = x*x;
-      t143 = t144*t144;
-      t141 = RATIONAL(-1.0,12.0);
-      t150 = RATIONAL(1.0,90.0)+RATIONAL(1.0,24.0)*t143+t141*t144;
-      t149 = RATIONAL(-3.0,20.0)+RATIONAL(-1.0,4.0)*t143+t144;
-      t148 = RATIONAL(3.0,2.0)+RATIONAL(-13.0,4.0)*t144+RATIONAL(5.0,8.0)*t143;
-      t145 = x*t144;
-      coeff_dxx_m3 = RATIONAL(1.0,8.0)*x+t141*t145+t150;
-      coeff_dxx_m2 = RATIONAL(1.0,3.0)*t145-x+t149;
-      coeff_dxx_m1 = RATIONAL(-5.0,12.0)*t145+RATIONAL(13.0,8.0)*x+t148;
-      coeff_dxx_0 = RATIONAL(14.0,3.0)*t144+RATIONAL(-5.0,6.0)*t143+RATIONAL(
--49.0,18.0);
-      coeff_dxx_p1 = RATIONAL(5.0,12.0)*t145+RATIONAL(-13.0,8.0)*x+t148;
-      coeff_dxx_p2 = RATIONAL(-1.0,3.0)*t145+x+t149;
-      coeff_dxx_p3 = RATIONAL(-1.0,8.0)*x+RATIONAL(1.0,12.0)*t145+t150;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-I.compute.c
new file mode 100644
index 0000000..28b9679
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-I.compute.c
@@ -0,0 +1,41 @@
+fp t103;
+fp t106;
+fp t107;
+fp t105;
+fp t113;
+fp t112;
+fp t100;
+fp t96;
+fp t101;
+fp t111;
+fp t104;
+fp t99;
+fp t98;
+fp t97;
+      t103 = x*x;
+      t106 = t103*t103;
+      t107 = x*t103;
+      t105 = t107*t107;
+      t113 = RATIONAL(1.0,720.0)*t105+RATIONAL(1.0,180.0)*t103+RATIONAL(-1.0,
+144.0)*t106;
+      t112 = RATIONAL(-1.0,120.0)*t105+RATIONAL(-3.0,40.0)*t103+RATIONAL(1.0,
+12.0)*t106;
+      t100 = RATIONAL(1.0,48.0);
+      t96 = RATIONAL(-13.0,48.0);
+      t101 = RATIONAL(3.0,4.0);
+      t111 = t100*t105+t96*t106+t101*t103;
+      t104 = t107*t103;
+      t99 = RATIONAL(-1.0,60.0);
+      t98 = RATIONAL(1.0,60.0);
+      t97 = RATIONAL(-1.0,48.0);
+      coeffs_I->coeff_m3 = t99*x+t100*t107+RATIONAL(-1.0,240.0)*t104+t113;
+      coeffs_I->coeff_m2 = RATIONAL(3.0,20.0)*x+RATIONAL(-1.0,6.0)*t107+t98*
+t104+t112;
+      coeffs_I->coeff_m1 = RATIONAL(-3.0,4.0)*x+RATIONAL(13.0,48.0)*t107+t97*
+t104+t111;
+      coeffs_I->coeff_0 = RATIONAL(7.0,18.0)*t106+RATIONAL(-49.0,36.0)*t103+
+RATIONAL(1.0,1.0)+RATIONAL(-1.0,36.0)*t105;
+      coeffs_I->coeff_p1 = t96*t107+t101*x+t100*t104+t111;
+      coeffs_I->coeff_p2 = RATIONAL(1.0,6.0)*t107+RATIONAL(-3.0,20.0)*x+t99*
+t104+t112;
+      coeffs_I->coeff_p3 = t97*t107+t98*x+RATIONAL(1.0,240.0)*t104+t113;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-dx.compute.c
new file mode 100644
index 0000000..d4d1004
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-dx.compute.c
@@ -0,0 +1,32 @@
+fp t123;
+fp t125;
+fp t128;
+fp t126;
+fp t133;
+fp t132;
+fp t131;
+fp t127;
+      t123 = RATIONAL(-3.0,20.0);
+      t125 = x*x;
+      t128 = x*t125;
+      t126 = t128*t125;
+      t133 = t123*x+RATIONAL(-1.0,20.0)*t126+RATIONAL(1.0,3.0)*t128;
+      t132 = RATIONAL(1.0,120.0)*t126+RATIONAL(1.0,90.0)*x+RATIONAL(-1.0,36.0)*
+t128;
+      t131 = RATIONAL(3.0,2.0)*x+RATIONAL(-13.0,12.0)*t128+RATIONAL(1.0,8.0)*
+t126;
+      t127 = t125*t125;
+      coeffs_dx->coeff_m3 = RATIONAL(-1.0,48.0)*t127+RATIONAL(-1.0,60.0)+
+RATIONAL(1.0,16.0)*t125+t132;
+      coeffs_dx->coeff_m2 = RATIONAL(-1.0,2.0)*t125+RATIONAL(3.0,20.0)+RATIONAL
+(1.0,12.0)*t127+t133;
+      coeffs_dx->coeff_m1 = RATIONAL(-5.0,48.0)*t127+RATIONAL(13.0,16.0)*t125+
+RATIONAL(-3.0,4.0)+t131;
+      coeffs_dx->coeff_0 = RATIONAL(-1.0,6.0)*t126+RATIONAL(14.0,9.0)*t128+
+RATIONAL(-49.0,18.0)*x;
+      coeffs_dx->coeff_p1 = RATIONAL(5.0,48.0)*t127+RATIONAL(3.0,4.0)+RATIONAL(
+-13.0,16.0)*t125+t131;
+      coeffs_dx->coeff_p2 = RATIONAL(1.0,2.0)*t125+RATIONAL(-1.0,12.0)*t127+
+t123+t133;
+      coeffs_dx->coeff_p3 = RATIONAL(-1.0,16.0)*t125+RATIONAL(1.0,60.0)+
+RATIONAL(1.0,48.0)*t127+t132;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-dxx.compute.c
new file mode 100644
index 0000000..02aa677
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-dxx.compute.c
@@ -0,0 +1,24 @@
+fp t144;
+fp t145;
+fp t150;
+fp t149;
+fp t141;
+fp t148;
+fp t143;
+      t144 = x*x;
+      t145 = t144*t144;
+      t150 = RATIONAL(-13.0,4.0)*t144+RATIONAL(3.0,2.0)+RATIONAL(5.0,8.0)*t145;
+      t149 = t144+RATIONAL(-3.0,20.0)+RATIONAL(-1.0,4.0)*t145;
+      t141 = RATIONAL(-1.0,12.0);
+      t148 = RATIONAL(1.0,90.0)+RATIONAL(1.0,24.0)*t145+t141*t144;
+      t143 = x*t144;
+      coeffs_dxx->coeff_m3 = RATIONAL(1.0,8.0)*x+t141*t143+t148;
+      coeffs_dxx->coeff_m2 = -x+RATIONAL(1.0,3.0)*t143+t149;
+      coeffs_dxx->coeff_m1 = RATIONAL(-5.0,12.0)*t143+RATIONAL(13.0,8.0)*x+t150
+;
+      coeffs_dxx->coeff_0 = RATIONAL(14.0,3.0)*t144+RATIONAL(-5.0,6.0)*t145+
+RATIONAL(-49.0,18.0);
+      coeffs_dxx->coeff_p1 = RATIONAL(5.0,12.0)*t143+RATIONAL(-13.0,8.0)*x+t150
+;
+      coeffs_dxx->coeff_p2 = x+RATIONAL(-1.0,3.0)*t143+t149;
+      coeffs_dxx->coeff_p3 = RATIONAL(-1.0,8.0)*x+RATIONAL(1.0,12.0)*t143+t148;