This commit reorganizes the Maple scripts which generate the interpolation

coefficients, and the coefficient files themselves, to properly support
multiple interpolation operators, and adds some skeleton support for
Hermite interpolation.  I will add full support for Hermite in a future
checkin.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@81 df1f8a13-aa1d-4dd4-9681-27ded5b42416

17 files changed, 259 insertions, 0 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-I.compute.c
new file mode 100644
index 0000000..b8ea11e
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-I.compute.c
@@ -0,0 +1,2 @@
+      coeff_I_0 = RATIONAL(1.0,1.0)-x;
+      coeff_I_p1 = x;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-dx.compute.c
new file mode 100644
index 0000000..4f5fb07
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-dx.compute.c
@@ -0,0 +1,2 @@
+      coeff_dx_0 = RATIONAL(-1.0,1.0);
+      coeff_dx_p1 = RATIONAL(1.0,1.0);
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-I.compute.c
new file mode 100644
index 0000000..4d0de20
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-I.compute.c
@@ -0,0 +1,9 @@
+fp t3,
+   t2,
+   t1;
+      t3 = x*x;
+      t2 = RATIONAL(1.0,2.0);
+      t1 = t2*t3;
+      coeff_I_m1 = RATIONAL(-1.0,2.0)*x+t1;
+      coeff_I_0 = RATIONAL(1.0,1.0)-t3;
+      coeff_I_p1 = t2*x+t1;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dx.compute.c
new file mode 100644
index 0000000..e2532c9
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dx.compute.c
@@ -0,0 +1,3 @@
+      coeff_dx_m1 = x+RATIONAL(-1.0,2.0);
+      coeff_dx_0 = RATIONAL(-2.0,1.0)*x;
+      coeff_dx_p1 = RATIONAL(1.0,2.0)+x;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dxx.compute.c
new file mode 100644
index 0000000..4e6fd95
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dxx.compute.c
@@ -0,0 +1,5 @@
+fp t4;
+      t4 = RATIONAL(1.0,1.0);
+      coeff_dxx_m1 = t4;
+      coeff_dxx_0 = RATIONAL(-2.0,1.0);
+      coeff_dxx_p1 = t4;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-I.compute.c
new file mode 100644
index 0000000..2df6cd1
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-I.compute.c
@@ -0,0 +1,16 @@
+fp t10,
+   t9,
+   t8,
+   t7,
+   t6,
+   t5;
+      t10 = x*x;
+      t9 = x*t10;
+      t8 = RATIONAL(-1.0,2.0);
+      t7 = RATIONAL(1.0,2.0);
+      t6 = RATIONAL(-1.0,6.0);
+      t5 = t7*t10;
+      coeff_I_m1 = RATIONAL(-1.0,3.0)*x+t5+t6*t9;
+      coeff_I_0 = RATIONAL(1.0,1.0)+t8*x-t10+t7*t9;
+      coeff_I_p1 = x+t5+t8*t9;
+      coeff_I_p2 = t6*x+RATIONAL(1.0,6.0)*t9;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dx.compute.c
new file mode 100644
index 0000000..4306fc4
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dx.compute.c
@@ -0,0 +1,8 @@
+fp t13,
+   t12;
+      t13 = x*x;
+      t12 = RATIONAL(-1.0,2.0);
+      coeff_dx_m1 = x+t12*t13+RATIONAL(-1.0,3.0);
+      coeff_dx_0 = t12+RATIONAL(-2.0,1.0)*x+RATIONAL(3.0,2.0)*t13;
+      coeff_dx_p1 = RATIONAL(-3.0,2.0)*t13+RATIONAL(1.0,1.0)+x;
+      coeff_dx_p2 = RATIONAL(1.0,2.0)*t13+RATIONAL(-1.0,6.0);
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dxx.compute.c
new file mode 100644
index 0000000..d5e199b
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dxx.compute.c
@@ -0,0 +1,6 @@
+fp t14;
+      t14 = RATIONAL(1.0,1.0);
+      coeff_dxx_m1 = t14-x;
+      coeff_dxx_0 = RATIONAL(-2.0,1.0)+RATIONAL(3.0,1.0)*x;
+      coeff_dxx_p1 = RATIONAL(-3.0,1.0)*x+t14;
+      coeff_dxx_p2 = x;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c
new file mode 100644
index 0000000..05b6c68
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c
@@ -0,0 +1,24 @@
+fp t23,
+   t24,
+   t29,
+   t21,
+   t22,
+   t28,
+   t25,
+   t20,
+   t19;
+      t23 = x*x;
+      t24 = t23*t23;
+      t29 = RATIONAL(-1.0,24.0)*t23+RATIONAL(1.0,24.0)*t24;
+      t21 = RATIONAL(2.0,3.0);
+      t22 = RATIONAL(-1.0,6.0);
+      t28 = t21*t23+t22*t24;
+      t25 = x*t23;
+      t20 = RATIONAL(1.0,12.0);
+      t19 = RATIONAL(-1.0,12.0);
+      coeff_I_m2 = t20*x+t19*t25+t29;
+      coeff_I_m1 = RATIONAL(-2.0,3.0)*x+RATIONAL(1.0,6.0)*t25+t28;
+      coeff_I_0 = RATIONAL(-5.0,4.0)*t23+RATIONAL(1.0,1.0)+RATIONAL(1.0,4.0)*
+t24;
+      coeff_I_p1 = t22*t25+t21*x+t28;
+      coeff_I_p2 = t20*t25+t19*x+t29;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dx.compute.c
new file mode 100644
index 0000000..5ac1469
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dx.compute.c
@@ -0,0 +1,17 @@
+fp t36,
+   t37,
+   t34,
+   t40,
+   t35,
+   t39;
+      t36 = x*x;
+      t37 = x*t36;
+      t34 = RATIONAL(-1.0,12.0);
+      t40 = t34*x+RATIONAL(1.0,6.0)*t37;
+      t35 = RATIONAL(-2.0,3.0);
+      t39 = t35*t37+RATIONAL(4.0,3.0)*x;
+      coeff_dx_m2 = RATIONAL(-1.0,4.0)*t36+RATIONAL(1.0,12.0)+t40;
+      coeff_dx_m1 = t35+RATIONAL(1.0,2.0)*t36+t39;
+      coeff_dx_0 = t37+RATIONAL(-5.0,2.0)*x;
+      coeff_dx_p1 = RATIONAL(2.0,3.0)+RATIONAL(-1.0,2.0)*t36+t39;
+      coeff_dx_p2 = t34+RATIONAL(1.0,4.0)*t36+t40;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dxx.compute.c
new file mode 100644
index 0000000..56d8567
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dxx.compute.c
@@ -0,0 +1,13 @@
+fp t44,
+   t46,
+   t48,
+   t47;
+      t44 = RATIONAL(1.0,2.0);
+      t46 = x*x;
+      t48 = t44*t46+RATIONAL(-1.0,12.0);
+      t47 = RATIONAL(-2.0,1.0)*t46+RATIONAL(4.0,3.0);
+      coeff_dxx_m2 = RATIONAL(-1.0,2.0)*x+t48;
+      coeff_dxx_m1 = x+t47;
+      coeff_dxx_0 = RATIONAL(-5.0,2.0)+RATIONAL(3.0,1.0)*t46;
+      coeff_dxx_p1 = -x+t47;
+      coeff_dxx_p2 = t44*x+t48;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-I.compute.c
new file mode 100644
index 0000000..f32ea37
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-I.compute.c
@@ -0,0 +1,25 @@
+fp t58,
+   t56,
+   t64,
+   t54,
+   t55,
+   t63,
+   t60,
+   t57,
+   t50;
+      t58 = x*x;
+      t56 = t58*t58;
+      t64 = RATIONAL(2.0,3.0)*t58+RATIONAL(-1.0,6.0)*t56;
+      t54 = RATIONAL(1.0,24.0);
+      t55 = RATIONAL(-1.0,24.0);
+      t63 = t54*t56+t55*t58;
+      t60 = t58*x;
+      t57 = t58*t60;
+      t50 = t55*t60;
+      coeff_I_m2 = RATIONAL(1.0,20.0)*x+t50+RATIONAL(-1.0,120.0)*t57+t63;
+      coeff_I_m1 = RATIONAL(-1.0,2.0)*x+t50+t54*t57+t64;
+      coeff_I_0 = RATIONAL(1.0,1.0)+RATIONAL(-1.0,3.0)*x+RATIONAL(-5.0,4.0)*t58
++RATIONAL(5.0,12.0)*t60+RATIONAL(1.0,4.0)*t56+RATIONAL(-1.0,12.0)*t57;
+      coeff_I_p1 = x+RATIONAL(-7.0,12.0)*t60+RATIONAL(1.0,12.0)*t57+t64;
+      coeff_I_p2 = RATIONAL(-1.0,4.0)*x+RATIONAL(7.0,24.0)*t60+t55*t57+t63;
+      coeff_I_p3 = RATIONAL(1.0,30.0)*x+RATIONAL(1.0,120.0)*t57+t50;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dx.compute.c
new file mode 100644
index 0000000..8fb00d3
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dx.compute.c
@@ -0,0 +1,21 @@
+fp t71,
+   t72,
+   t76,
+   t75,
+   t70,
+   t67;
+      t71 = x*x;
+      t72 = x*t71;
+      t76 = RATIONAL(-1.0,12.0)*x+RATIONAL(1.0,6.0)*t72;
+      t75 = RATIONAL(-2.0,3.0)*t72+RATIONAL(4.0,3.0)*x;
+      t70 = t71*t71;
+      t67 = RATIONAL(-1.0,8.0)*t71;
+      coeff_dx_m2 = t67+RATIONAL(-1.0,24.0)*t70+RATIONAL(1.0,20.0)+t76;
+      coeff_dx_m1 = RATIONAL(-1.0,2.0)+t67+RATIONAL(5.0,24.0)*t70+t75;
+      coeff_dx_0 = RATIONAL(5.0,4.0)*t71+t72+RATIONAL(-5.0,12.0)*t70+RATIONAL(
+-5.0,2.0)*x+RATIONAL(-1.0,3.0);
+      coeff_dx_p1 = RATIONAL(5.0,12.0)*t70+RATIONAL(-7.0,4.0)*t71+RATIONAL(1.0,
+1.0)+t75;
+      coeff_dx_p2 = RATIONAL(-1.0,4.0)+RATIONAL(7.0,8.0)*t71+RATIONAL(-5.0,24.0
+)*t70+t76;
+      coeff_dx_p3 = RATIONAL(1.0,24.0)*t70+RATIONAL(1.0,30.0)+t67;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dxx.compute.c
new file mode 100644
index 0000000..4c7a405
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dxx.compute.c
@@ -0,0 +1,17 @@
+fp t82,
+   t86,
+   t85,
+   t83,
+   t79;
+      t82 = x*x;
+      t86 = RATIONAL(1.0,2.0)*t82+RATIONAL(-1.0,12.0);
+      t85 = RATIONAL(-2.0,1.0)*t82+RATIONAL(4.0,3.0);
+      t83 = t82*x;
+      t79 = RATIONAL(-1.0,4.0)*x;
+      coeff_dxx_m2 = t79+RATIONAL(-1.0,6.0)*t83+t86;
+      coeff_dxx_m1 = t79+RATIONAL(5.0,6.0)*t83+t85;
+      coeff_dxx_0 = RATIONAL(-5.0,3.0)*t83+RATIONAL(5.0,2.0)*x+RATIONAL(3.0,1.0
+)*t82+RATIONAL(-5.0,2.0);
+      coeff_dxx_p1 = RATIONAL(5.0,3.0)*t83+RATIONAL(-7.0,2.0)*x+t85;
+      coeff_dxx_p2 = RATIONAL(7.0,4.0)*x+RATIONAL(-5.0,6.0)*t83+t86;
+      coeff_dxx_p3 = t79+RATIONAL(1.0,6.0)*t83;
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
new file mode 100644
index 0000000..b189aa6
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c
@@ -0,0 +1,38 @@
+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
new file mode 100644
index 0000000..ab1b33f
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c
@@ -0,0 +1,31 @@
+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
new file mode 100644
index 0000000..46979da
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c
@@ -0,0 +1,22 @@
+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;