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

21 files changed, 1566 insertions, 0 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c
new file mode 100644
index 0000000..221b33b
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c
@@ -0,0 +1,18 @@
+fp t7,
+   t6,
+   t5,
+   t4,
+   t3,
+   t2,
+   t1;
+      t7 = RATIONAL(-1.0,2.0);
+      t6 = RATIONAL(1.0,4.0);
+      t5 = RATIONAL(1.0,2.0);
+      t4 = t7*y;
+      t3 = t7*x;
+      t2 = t5*x;
+      t1 = t5*y;
+      coeff_I_0_0 = RATIONAL(3.0,4.0)+t4+t3;
+      coeff_I_p1_0 = t6+t4+t2;
+      coeff_I_0_p1 = t3+t6+t1;
+      coeff_I_p1_p1 = t1+RATIONAL(-1.0,4.0)+t2;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c
new file mode 100644
index 0000000..3d89f01
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c
@@ -0,0 +1,8 @@
+fp t9,
+   t8;
+      t9 = RATIONAL(-1.0,2.0);
+      t8 = RATIONAL(1.0,2.0);
+      coeff_dx_0_0 = t9;
+      coeff_dx_p1_0 = t8;
+      coeff_dx_0_p1 = t9;
+      coeff_dx_p1_p1 = t8;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c
new file mode 100644
index 0000000..42819a8
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c
@@ -0,0 +1,8 @@
+fp t11,
+   t10;
+      t11 = RATIONAL(-1.0,2.0);
+      t10 = RATIONAL(1.0,2.0);
+      coeff_dy_0_0 = t11;
+      coeff_dy_p1_0 = t11;
+      coeff_dy_0_p1 = t10;
+      coeff_dy_p1_p1 = t10;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c
new file mode 100644
index 0000000..34f402b
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c
@@ -0,0 +1,49 @@
+fp t25,
+   t28,
+   t20,
+   t27,
+   t18,
+   t34,
+   t33,
+   t23,
+   t22,
+   t17,
+   t32,
+   t15,
+   t31,
+   t30,
+   t29,
+   t26,
+   t21,
+   t19,
+   t16,
+   t14;
+      t25 = RATIONAL(1.0,6.0);
+      t28 = y*y;
+      t20 = t25*t28;
+      t27 = x*x;
+      t18 = t25*t27;
+      t34 = t20+RATIONAL(-1.0,9.0)+t18;
+      t33 = x*y;
+      t23 = RATIONAL(2.0,9.0);
+      t22 = RATIONAL(-1.0,3.0);
+      t17 = t22*t27;
+      t32 = t20+t23+t17;
+      t15 = t22*t28;
+      t31 = t15+t23+t18;
+      t30 = RATIONAL(-1.0,4.0)*t33+t34;
+      t29 = RATIONAL(1.0,4.0)*t33+t34;
+      t26 = RATIONAL(-1.0,6.0);
+      t21 = t26*y;
+      t19 = t26*x;
+      t16 = t25*x;
+      t14 = t25*y;
+      coeff_I_m1_m1 = t21+t19+t29;
+      coeff_I_0_m1 = t21+t32;
+      coeff_I_p1_m1 = t21+t16+t30;
+      coeff_I_m1_0 = t19+t31;
+      coeff_I_0_0 = t15+t17+RATIONAL(5.0,9.0);
+      coeff_I_p1_0 = t16+t31;
+      coeff_I_m1_p1 = t14+t19+t30;
+      coeff_I_0_p1 = t14+t32;
+      coeff_I_p1_p1 = t14+t16+t29;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c
new file mode 100644
index 0000000..22a3439
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c
@@ -0,0 +1,21 @@
+fp t37,
+   t42,
+   t41,
+   t40,
+   t39,
+   t35;
+      t37 = RATIONAL(1.0,3.0)*x;
+      t42 = t37+RATIONAL(1.0,4.0)*y;
+      t41 = RATIONAL(-1.0,4.0)*y+t37;
+      t40 = RATIONAL(-1.0,6.0);
+      t39 = RATIONAL(1.0,6.0);
+      t35 = RATIONAL(-2.0,3.0)*x;
+      coeff_dx_m1_m1 = t40+t42;
+      coeff_dx_0_m1 = t35;
+      coeff_dx_p1_m1 = t39+t41;
+      coeff_dx_m1_0 = t40+t37;
+      coeff_dx_0_0 = t35;
+      coeff_dx_p1_0 = t39+t37;
+      coeff_dx_m1_p1 = t40+t41;
+      coeff_dx_0_p1 = t35;
+      coeff_dx_p1_p1 = t39+t42;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c
new file mode 100644
index 0000000..5f6427a
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c
@@ -0,0 +1,13 @@
+fp t52,
+   t51;
+      t52 = RATIONAL(1.0,3.0);
+      t51 = RATIONAL(-2.0,3.0);
+      coeff_dxx_m1_m1 = t52;
+      coeff_dxx_0_m1 = t51;
+      coeff_dxx_p1_m1 = t52;
+      coeff_dxx_m1_0 = t52;
+      coeff_dxx_0_0 = t51;
+      coeff_dxx_p1_0 = t52;
+      coeff_dxx_m1_p1 = t52;
+      coeff_dxx_0_p1 = t51;
+      coeff_dxx_p1_p1 = t52;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c
new file mode 100644
index 0000000..c795e1d
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c
@@ -0,0 +1,15 @@
+fp t55,
+   t54,
+   t53;
+      t55 = RATIONAL(0.0,1.0);
+      t54 = RATIONAL(1.0,4.0);
+      t53 = RATIONAL(-1.0,4.0);
+      coeff_dxy_m1_m1 = t54;
+      coeff_dxy_0_m1 = t55;
+      coeff_dxy_p1_m1 = t53;
+      coeff_dxy_m1_0 = t55;
+      coeff_dxy_0_0 = t55;
+      coeff_dxy_p1_0 = t55;
+      coeff_dxy_m1_p1 = t53;
+      coeff_dxy_0_p1 = t55;
+      coeff_dxy_p1_p1 = t54;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c
new file mode 100644
index 0000000..4afa754
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c
@@ -0,0 +1,21 @@
+fp t46,
+   t50,
+   t49,
+   t48,
+   t47,
+   t43;
+      t46 = RATIONAL(1.0,3.0)*y;
+      t50 = RATIONAL(-1.0,4.0)*x+t46;
+      t49 = RATIONAL(1.0,4.0)*x+t46;
+      t48 = RATIONAL(-1.0,6.0);
+      t47 = RATIONAL(1.0,6.0);
+      t43 = RATIONAL(-2.0,3.0)*y;
+      coeff_dy_m1_m1 = t48+t49;
+      coeff_dy_0_m1 = t46+t48;
+      coeff_dy_p1_m1 = t48+t50;
+      coeff_dy_m1_0 = t43;
+      coeff_dy_0_0 = t43;
+      coeff_dy_p1_0 = t43;
+      coeff_dy_m1_p1 = t47+t50;
+      coeff_dy_0_p1 = t47+t46;
+      coeff_dy_p1_p1 = t47+t49;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c
new file mode 100644
index 0000000..a643108
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c
@@ -0,0 +1,13 @@
+fp t57,
+   t56;
+      t57 = RATIONAL(1.0,3.0);
+      t56 = RATIONAL(-2.0,3.0);
+      coeff_dyy_m1_m1 = t57;
+      coeff_dyy_0_m1 = t57;
+      coeff_dyy_p1_m1 = t57;
+      coeff_dyy_m1_0 = t56;
+      coeff_dyy_0_0 = t56;
+      coeff_dyy_p1_0 = t56;
+      coeff_dyy_m1_p1 = t57;
+      coeff_dyy_0_p1 = t57;
+      coeff_dyy_p1_p1 = t57;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c
new file mode 100644
index 0000000..55cc52e
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c
@@ -0,0 +1,196 @@
+fp t160,
+   t131,
+   t159,
+   t130,
+   t158,
+   t124,
+   t68,
+   t133,
+   t126,
+   t82,
+   t157,
+   t71,
+   t132,
+   t83,
+   t156,
+   t121,
+   t70,
+   t128,
+   t87,
+   t155,
+   t67,
+   t129,
+   t86,
+   t154,
+   t85,
+   t153,
+   t117,
+   t77,
+   t152,
+   t127,
+   t62,
+   t151,
+   t63,
+   t150,
+   t125,
+   t64,
+   t149,
+   t79,
+   t148,
+   t147,
+   t146,
+   t69,
+   t110,
+   t145,
+   t84,
+   t144,
+   t143,
+   t123,
+   t142,
+   t104,
+   t141,
+   t140,
+   t139,
+   t138,
+   t137,
+   t136,
+   t122,
+   t120,
+   t119,
+   t116,
+   t114,
+   t113,
+   t111,
+   t108,
+   t107,
+   t105,
+   t103,
+   t102,
+   t101,
+   t99,
+   t98,
+   t97,
+   t96,
+   t95,
+   t94,
+   t93,
+   t92,
+   t91,
+   t90,
+   t89,
+   t88,
+   t81,
+   t80,
+   t78,
+   t76,
+   t75,
+   t74,
+   t73;
+      t160 = x*y;
+      t131 = y*y;
+      t159 = x*t131;
+      t130 = x*x;
+      t158 = t130*y;
+      t124 = RATIONAL(3.0,40.0);
+      t68 = t124*t158;
+      t133 = t131*y;
+      t126 = RATIONAL(-1.0,24.0);
+      t82 = t126*t133;
+      t157 = t68+t82;
+      t71 = t124*t159;
+      t132 = t130*x;
+      t83 = t126*t132;
+      t156 = t71+t83;
+      t121 = RATIONAL(1.0,40.0);
+      t70 = t121*t158;
+      t128 = RATIONAL(-1.0,8.0);
+      t87 = t128*t133;
+      t155 = t70+t87;
+      t67 = t121*t159;
+      t129 = RATIONAL(1.0,8.0);
+      t86 = t129*t132;
+      t154 = t67+t86;
+      t85 = t128*t132;
+      t153 = t67+t85;
+      t117 = RATIONAL(1.0,24.0);
+      t77 = t117*t132;
+      t152 = t71+t77;
+      t127 = RATIONAL(-3.0,40.0);
+      t62 = t127*t158;
+      t151 = t62+t82;
+      t63 = t127*t159;
+      t150 = t63+t77;
+      t125 = RATIONAL(-1.0,40.0);
+      t64 = t125*t159;
+      t149 = t64+t86;
+      t79 = t117*t133;
+      t148 = t68+t79;
+      t147 = t63+t83;
+      t146 = t62+t79;
+      t69 = t125*t158;
+      t110 = RATIONAL(11.0,80.0);
+      t145 = t69+t110*t130;
+      t84 = t129*t133;
+      t144 = t70+t84;
+      t143 = t64+t85;
+      t123 = RATIONAL(-1.0,50.0);
+      t142 = RATIONAL(-13.0,100.0)*t160+t123;
+      t104 = RATIONAL(-17.0,200.0);
+      t141 = RATIONAL(-9.0,100.0)*t160+t104;
+      t140 = t123*t160+RATIONAL(11.0,50.0);
+      t139 = RATIONAL(-1.0,100.0)*t160+RATIONAL(37.0,200.0);
+      t138 = RATIONAL(-2.0,25.0)+RATIONAL(2.0,25.0)*t160;
+      t137 = RATIONAL(7.0,100.0)*t160+RATIONAL(13.0,100.0);
+      t136 = t130+t131;
+      t122 = RATIONAL(9.0,80.0);
+      t120 = RATIONAL(7.0,80.0);
+      t119 = RATIONAL(-1.0,80.0);
+      t116 = RATIONAL(3.0,80.0);
+      t114 = RATIONAL(1.0,80.0);
+      t113 = RATIONAL(13.0,80.0);
+      t111 = RATIONAL(-23.0,80.0);
+      t108 = RATIONAL(-19.0,80.0);
+      t107 = RATIONAL(-17.0,80.0);
+      t105 = RATIONAL(-31.0,400.0);
+      t103 = RATIONAL(69.0,400.0);
+      t102 = RATIONAL(-57.0,400.0);
+      t101 = RATIONAL(-63.0,400.0);
+      t99 = RATIONAL(-49.0,400.0);
+      t98 = RATIONAL(43.0,1200.0);
+      t97 = RATIONAL(103.0,400.0);
+      t96 = RATIONAL(111.0,400.0);
+      t95 = RATIONAL(37.0,1200.0);
+      t94 = RATIONAL(89.0,1200.0);
+      t93 = RATIONAL(117.0,400.0);
+      t92 = RATIONAL(-109.0,1200.0);
+      t91 = RATIONAL(-223.0,1200.0);
+      t90 = RATIONAL(-157.0,1200.0);
+      t89 = RATIONAL(-131.0,1200.0);
+      t88 = RATIONAL(-149.0,1200.0);
+      t81 = t122*t131;
+      t80 = t120*t130;
+      t78 = t120*t131;
+      t76 = t122*t130;
+      t75 = t113*t131;
+      t74 = t113*t130;
+      t73 = t110*t131;
+      coeff_I_m1_m1 = t75+RATIONAL(-23.0,200.0)+t74+t92*y+(t92+RATIONAL(6.0,
+25.0)*y)*x+t147+t151;
+      coeff_I_0_m1 = t111*t130+t105*x+t73+t91*y+t140+t149+t157;
+      coeff_I_p1_m1 = t81+t96*x+t90*y+t80+t142+t153+t157;
+      coeff_I_p2_m1 = t78+t94*y+t89*x+t116*t130+t141+t151+t152;
+      coeff_I_m1_0 = t105*y+t91*x+t111*t131+t84+t140+t145+t156;
+      coeff_I_0_0 = RATIONAL(93.0,200.0)+t102*y+t136*RATIONAL(-21.0,80.0)+(
+RATIONAL(-1.0,25.0)*y+t102)*x+t144+t154;
+      coeff_I_p1_0 = t108*t131+t101*y+t76+t93*x+t139+t143+t144;
+      coeff_I_p2_0 = t107*t131+t114*t130+t99*y+t84+t69+t98*x+t137+t150;
+      coeff_I_m1_p1 = t90*x+t96*y+t78+t76+t142+t155+t156;
+      coeff_I_0_p1 = t69+t101*x+t108*t130+t87+t93*y+t81+t139+t154;
+      coeff_I_p1_p1 = t97*y+t73+t104+t87+(RATIONAL(3.0,50.0)*y+t97)*x+t143+t145
+;
+      coeff_I_p2_p1 = t95*x+t119*t130+t75+t103*y+t138+t150+t155;
+      coeff_I_m1_p2 = t116*t131+t94*x+t89*y+t80+t141+t147+t148;
+      coeff_I_0_p2 = t99*x+t98*y+t114*t131+t107*t130+t137+t146+t149;
+      coeff_I_p1_p2 = t103*x+t119*t131+t95*y+t74+t138+t146+t153;
+      coeff_I_p2_p2 = RATIONAL(7.0,200.0)+t88*y+t136*RATIONAL(-3.0,80.0)+(t88+
+RATIONAL(-3.0,50.0)*y)*x+t148+t152;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c
new file mode 100644
index 0000000..4650764
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c
@@ -0,0 +1,95 @@
+fp t201,
+   t164,
+   t200,
+   t163,
+   t199,
+   t183,
+   t188,
+   t173,
+   t187,
+   t179,
+   t198,
+   t161,
+   t197,
+   t162,
+   t196,
+   t185,
+   t176,
+   t195,
+   t181,
+   t194,
+   t193,
+   t186,
+   t177,
+   t182,
+   t192,
+   t180,
+   t191,
+   t184,
+   t172,
+   t190,
+   t189,
+   t178,
+   t175,
+   t170,
+   t168,
+   t166,
+   t165;
+      t201 = x*y;
+      t164 = RATIONAL(-3.0,20.0)*t201;
+      t200 = RATIONAL(13.0,40.0)*x+t164;
+      t163 = RATIONAL(3.0,20.0)*t201;
+      t199 = RATIONAL(7.0,40.0)*x+t163;
+      t183 = RATIONAL(1.0,40.0);
+      t188 = y*y;
+      t173 = t183*t188;
+      t187 = x*x;
+      t179 = RATIONAL(-3.0,8.0)*t187;
+      t198 = t173+t179;
+      t161 = RATIONAL(1.0,20.0)*t201;
+      t197 = RATIONAL(9.0,40.0)*x+t161;
+      t162 = RATIONAL(-1.0,20.0)*t201;
+      t196 = RATIONAL(11.0,40.0)*x+t162;
+      t185 = RATIONAL(-1.0,40.0);
+      t176 = t185*t188;
+      t195 = t176+t179;
+      t181 = RATIONAL(3.0,8.0)*t187;
+      t194 = t173+t181;
+      t193 = t176+t181;
+      t186 = RATIONAL(-3.0,40.0);
+      t177 = t186*t188;
+      t182 = RATIONAL(-1.0,8.0)*t187;
+      t192 = t177+t182;
+      t180 = RATIONAL(1.0,8.0)*t187;
+      t191 = t177+t180;
+      t184 = RATIONAL(3.0,40.0);
+      t172 = t184*t188;
+      t190 = t172+t180;
+      t189 = t172+t182;
+      t178 = RATIONAL(2.0,25.0)*y;
+      t175 = RATIONAL(-1.0,50.0)*y;
+      t170 = RATIONAL(7.0,100.0)*y;
+      t168 = RATIONAL(-9.0,100.0)*y;
+      t166 = RATIONAL(-1.0,100.0)*y;
+      t165 = RATIONAL(-13.0,100.0)*y;
+      coeff_dx_m1_m1 = RATIONAL(-109.0,1200.0)+RATIONAL(6.0,25.0)*y+t192+t200;
+      coeff_dx_0_m1 = RATIONAL(-31.0,400.0)+RATIONAL(-23.0,40.0)*x+t175+t163+
+t193;
+      coeff_dx_p1_m1 = RATIONAL(111.0,400.0)+t165+t198+t199;
+      coeff_dx_p2_m1 = t168+t164+RATIONAL(-131.0,1200.0)+t184*x+t190;
+      coeff_dx_m1_0 = RATIONAL(-223.0,1200.0)+t175+t189+t196;
+      coeff_dx_0_0 = RATIONAL(-57.0,400.0)+RATIONAL(-21.0,40.0)*x+RATIONAL(-1.0
+,25.0)*y+t161+t194;
+      coeff_dx_p1_0 = RATIONAL(117.0,400.0)+t166+t195+t197;
+      coeff_dx_p2_0 = RATIONAL(43.0,1200.0)+t183*x+t162+t170+t191;
+      coeff_dx_m1_p1 = RATIONAL(-157.0,1200.0)+t165+t189+t197;
+      coeff_dx_0_p1 = RATIONAL(-63.0,400.0)+t166+t162+RATIONAL(-19.0,40.0)*x+
+t194;
+      coeff_dx_p1_p1 = RATIONAL(3.0,50.0)*y+RATIONAL(103.0,400.0)+t195+t196;
+      coeff_dx_p2_p1 = t185*x+t161+RATIONAL(37.0,1200.0)+t178+t191;
+      coeff_dx_m1_p2 = t168+RATIONAL(89.0,1200.0)+t192+t199;
+      coeff_dx_0_p2 = t164+t170+RATIONAL(-17.0,40.0)*x+RATIONAL(-49.0,400.0)+
+t193;
+      coeff_dx_p1_p2 = RATIONAL(69.0,400.0)+t178+t198+t200;
+      coeff_dx_p2_p2 = t163+RATIONAL(-149.0,1200.0)+t186*x+RATIONAL(-3.0,50.0)*
+y+t190;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c
new file mode 100644
index 0000000..395eceb
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c
@@ -0,0 +1,40 @@
+fp t243,
+   t258,
+   t245,
+   t257,
+   t246,
+   t256,
+   t244,
+   t255,
+   t250,
+   t249,
+   t248,
+   t247;
+      t243 = RATIONAL(1.0,20.0)*y;
+      t258 = t243+RATIONAL(9.0,40.0);
+      t245 = RATIONAL(3.0,20.0)*y;
+      t257 = t245+RATIONAL(7.0,40.0);
+      t246 = RATIONAL(-3.0,20.0)*y;
+      t256 = t246+RATIONAL(13.0,40.0);
+      t244 = RATIONAL(-1.0,20.0)*y;
+      t255 = t244+RATIONAL(11.0,40.0);
+      t250 = RATIONAL(-1.0,4.0)*x;
+      t249 = RATIONAL(3.0,4.0)*x;
+      t248 = RATIONAL(-3.0,4.0)*x;
+      t247 = RATIONAL(1.0,4.0)*x;
+      coeff_dxx_m1_m1 = t250+t256;
+      coeff_dxx_0_m1 = RATIONAL(-23.0,40.0)+t245+t249;
+      coeff_dxx_p1_m1 = t248+t257;
+      coeff_dxx_p2_m1 = t246+t247+RATIONAL(3.0,40.0);
+      coeff_dxx_m1_0 = t250+t255;
+      coeff_dxx_0_0 = t243+RATIONAL(-21.0,40.0)+t249;
+      coeff_dxx_p1_0 = t248+t258;
+      coeff_dxx_p2_0 = RATIONAL(1.0,40.0)+t244+t247;
+      coeff_dxx_m1_p1 = t250+t258;
+      coeff_dxx_0_p1 = t244+RATIONAL(-19.0,40.0)+t249;
+      coeff_dxx_p1_p1 = t248+t255;
+      coeff_dxx_p2_p1 = t243+t247+RATIONAL(-1.0,40.0);
+      coeff_dxx_m1_p2 = t250+t257;
+      coeff_dxx_0_p2 = RATIONAL(-17.0,40.0)+t246+t249;
+      coeff_dxx_p1_p2 = t248+t256;
+      coeff_dxx_p2_p2 = RATIONAL(-3.0,40.0)+t247+t245;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c
new file mode 100644
index 0000000..b48e5bc
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c
@@ -0,0 +1,52 @@
+fp t276,
+   t275,
+   t274,
+   t273,
+   t272,
+   t271,
+   t270,
+   t269,
+   t268,
+   t267,
+   t266,
+   t265,
+   t264,
+   t263,
+   t262,
+   t261,
+   t260,
+   t259;
+      t276 = RATIONAL(-1.0,50.0);
+      t275 = RATIONAL(-3.0,20.0);
+      t274 = RATIONAL(3.0,20.0);
+      t273 = RATIONAL(2.0,25.0);
+      t272 = RATIONAL(-1.0,20.0);
+      t271 = RATIONAL(1.0,20.0);
+      t270 = RATIONAL(-9.0,100.0);
+      t269 = RATIONAL(7.0,100.0);
+      t268 = RATIONAL(-1.0,100.0);
+      t267 = RATIONAL(-13.0,100.0);
+      t266 = t275*y;
+      t265 = t274*y;
+      t264 = t275*x;
+      t263 = t274*x;
+      t262 = t271*x;
+      t261 = t272*x;
+      t260 = t272*y;
+      t259 = t271*y;
+      coeff_dxy_m1_m1 = t266+t264+RATIONAL(6.0,25.0);
+      coeff_dxy_0_m1 = t263+t260+t276;
+      coeff_dxy_p1_m1 = t263+t259+t267;
+      coeff_dxy_p2_m1 = t270+t264+t265;
+      coeff_dxy_m1_0 = t276+t261+t265;
+      coeff_dxy_0_0 = RATIONAL(-1.0,25.0)+t259+t262;
+      coeff_dxy_p1_0 = t268+t262+t260;
+      coeff_dxy_p2_0 = t261+t269+t266;
+      coeff_dxy_m1_p1 = t267+t262+t265;
+      coeff_dxy_0_p1 = t261+t268+t259;
+      coeff_dxy_p1_p1 = RATIONAL(3.0,50.0)+t261+t260;
+      coeff_dxy_p2_p1 = t266+t262+t273;
+      coeff_dxy_m1_p2 = t266+t263+t270;
+      coeff_dxy_0_p2 = t269+t264+t260;
+      coeff_dxy_p1_p2 = t264+t259+t273;
+      coeff_dxy_p2_p2 = t263+t265+RATIONAL(-3.0,50.0);
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c
new file mode 100644
index 0000000..678f775
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c
@@ -0,0 +1,95 @@
+fp t242,
+   t228,
+   t226,
+   t212,
+   t229,
+   t220,
+   t241,
+   t224,
+   t213,
+   t240,
+   t202,
+   t239,
+   t223,
+   t238,
+   t204,
+   t237,
+   t236,
+   t227,
+   t217,
+   t221,
+   t235,
+   t205,
+   t234,
+   t203,
+   t233,
+   t225,
+   t218,
+   t232,
+   t222,
+   t231,
+   t230,
+   t219,
+   t214,
+   t209,
+   t208,
+   t207,
+   t206;
+      t242 = x*y;
+      t228 = x*x;
+      t226 = RATIONAL(-1.0,40.0);
+      t212 = t226*t228;
+      t229 = y*y;
+      t220 = RATIONAL(3.0,8.0)*t229;
+      t241 = t212+t220;
+      t224 = RATIONAL(1.0,40.0);
+      t213 = t224*t228;
+      t240 = t213+t220;
+      t202 = RATIONAL(3.0,20.0)*t242;
+      t239 = t202+RATIONAL(7.0,40.0)*y;
+      t223 = RATIONAL(-3.0,8.0)*t229;
+      t238 = t213+t223;
+      t204 = RATIONAL(1.0,20.0)*t242;
+      t237 = t204+RATIONAL(9.0,40.0)*y;
+      t236 = t212+t223;
+      t227 = RATIONAL(-3.0,40.0);
+      t217 = t227*t228;
+      t221 = RATIONAL(-1.0,8.0)*t229;
+      t235 = t217+t221;
+      t205 = RATIONAL(-3.0,20.0)*t242;
+      t234 = t205+RATIONAL(13.0,40.0)*y;
+      t203 = RATIONAL(-1.0,20.0)*t242;
+      t233 = t203+RATIONAL(11.0,40.0)*y;
+      t225 = RATIONAL(3.0,40.0);
+      t218 = t225*t228;
+      t232 = t218+t221;
+      t222 = RATIONAL(1.0,8.0)*t229;
+      t231 = t218+t222;
+      t230 = t217+t222;
+      t219 = RATIONAL(2.0,25.0)*x;
+      t214 = RATIONAL(-1.0,50.0)*x;
+      t209 = RATIONAL(-1.0,100.0)*x;
+      t208 = RATIONAL(-9.0,100.0)*x;
+      t207 = RATIONAL(7.0,100.0)*x;
+      t206 = RATIONAL(-13.0,100.0)*x;
+      coeff_dy_m1_m1 = RATIONAL(-109.0,1200.0)+RATIONAL(6.0,25.0)*x+t234+t235;
+      coeff_dy_0_m1 = t214+RATIONAL(-223.0,1200.0)+t232+t233;
+      coeff_dy_p1_m1 = t206+RATIONAL(-157.0,1200.0)+t232+t237;
+      coeff_dy_p2_m1 = t208+RATIONAL(89.0,1200.0)+t235+t239;
+      coeff_dy_m1_0 = RATIONAL(-23.0,40.0)*y+t214+RATIONAL(-31.0,400.0)+t202+
+t241;
+      coeff_dy_0_0 = RATIONAL(-21.0,40.0)*y+RATIONAL(-57.0,400.0)+t204+RATIONAL
+(-1.0,25.0)*x+t240;
+      coeff_dy_p1_0 = t209+RATIONAL(-63.0,400.0)+t203+RATIONAL(-19.0,40.0)*y+
+t240;
+      coeff_dy_p2_0 = t207+RATIONAL(-49.0,400.0)+RATIONAL(-17.0,40.0)*y+t205+
+t241;
+      coeff_dy_m1_p1 = RATIONAL(111.0,400.0)+t206+t238+t239;
+      coeff_dy_0_p1 = RATIONAL(117.0,400.0)+t209+t236+t237;
+      coeff_dy_p1_p1 = RATIONAL(103.0,400.0)+RATIONAL(3.0,50.0)*x+t233+t236;
+      coeff_dy_p2_p1 = t219+RATIONAL(69.0,400.0)+t234+t238;
+      coeff_dy_m1_p2 = t225*y+t208+RATIONAL(-131.0,1200.0)+t205+t231;
+      coeff_dy_0_p2 = t207+t203+t224*y+RATIONAL(43.0,1200.0)+t230;
+      coeff_dy_p1_p2 = t219+t226*y+RATIONAL(37.0,1200.0)+t204+t230;
+      coeff_dy_p2_p2 = t202+RATIONAL(-3.0,50.0)*x+RATIONAL(-149.0,1200.0)+t227*
+y+t231;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c
new file mode 100644
index 0000000..dea29f0
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c
@@ -0,0 +1,40 @@
+fp t279,
+   t292,
+   t280,
+   t291,
+   t278,
+   t290,
+   t277,
+   t289,
+   t284,
+   t283,
+   t282,
+   t281;
+      t279 = RATIONAL(-1.0,20.0)*x;
+      t292 = t279+RATIONAL(11.0,40.0);
+      t280 = RATIONAL(-3.0,20.0)*x;
+      t291 = RATIONAL(13.0,40.0)+t280;
+      t278 = RATIONAL(1.0,20.0)*x;
+      t290 = t278+RATIONAL(9.0,40.0);
+      t277 = RATIONAL(3.0,20.0)*x;
+      t289 = t277+RATIONAL(7.0,40.0);
+      t284 = RATIONAL(-1.0,4.0)*y;
+      t283 = RATIONAL(1.0,4.0)*y;
+      t282 = RATIONAL(3.0,4.0)*y;
+      t281 = RATIONAL(-3.0,4.0)*y;
+      coeff_dyy_m1_m1 = t284+t291;
+      coeff_dyy_0_m1 = t284+t292;
+      coeff_dyy_p1_m1 = t284+t290;
+      coeff_dyy_p2_m1 = t284+t289;
+      coeff_dyy_m1_0 = t282+RATIONAL(-23.0,40.0)+t277;
+      coeff_dyy_0_0 = RATIONAL(-21.0,40.0)+t278+t282;
+      coeff_dyy_p1_0 = RATIONAL(-19.0,40.0)+t282+t279;
+      coeff_dyy_p2_0 = t280+t282+RATIONAL(-17.0,40.0);
+      coeff_dyy_m1_p1 = t281+t289;
+      coeff_dyy_0_p1 = t281+t290;
+      coeff_dyy_p1_p1 = t281+t292;
+      coeff_dyy_p2_p1 = t281+t291;
+      coeff_dyy_m1_p2 = t280+t283+RATIONAL(3.0,40.0);
+      coeff_dyy_0_p2 = t283+RATIONAL(1.0,40.0)+t279;
+      coeff_dyy_p1_p2 = t278+RATIONAL(-1.0,40.0)+t283;
+      coeff_dyy_p2_p2 = t277+t283+RATIONAL(-3.0,40.0);
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
new file mode 100644
index 0000000..49655f0
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
@@ -0,0 +1,239 @@
+fp t413,
+   t418,
+   t466,
+   t403,
+   t417,
+   t358,
+   t414,
+   t373,
+   t438,
+   t444,
+   t465,
+   t405,
+   t364,
+   t367,
+   t439,
+   t464,
+   t436,
+   t463,
+   t410,
+   t359,
+   t442,
+   t409,
+   t318,
+   t392,
+   t462,
+   t411,
+   t360,
+   t441,
+   t400,
+   t298,
+   t461,
+   t366,
+   t443,
+   t395,
+   t401,
+   t460,
+   t440,
+   t296,
+   t399,
+   t459,
+   t372,
+   t458,
+   t333,
+   t457,
+   t396,
+   t412,
+   t327,
+   t456,
+   t407,
+   t331,
+   t391,
+   t455,
+   t402,
+   t313,
+   t397,
+   t454,
+   t325,
+   t453,
+   t324,
+   t452,
+   t311,
+   t451,
+   t450,
+   t387,
+   t398,
+   t449,
+   t415,
+   t371,
+   t416,
+   t357,
+   t383,
+   t382,
+   t448,
+   t317,
+   t356,
+   t447,
+   t362,
+   t446,
+   t445,
+   t437,
+   t435,
+   t393,
+   t404,
+   t369,
+   t376,
+   t434,
+   t365,
+   t384,
+   t386,
+   t433,
+   t432,
+   t431,
+   t430,
+   t429,
+   t428,
+   t427,
+   t406,
+   t426,
+   t425,
+   t408,
+   t424,
+   t423,
+   t394,
+   t389,
+   t388,
+   t385;
+      t413 = x*x;
+      t418 = y*y;
+      t466 = t418+t413;
+      t403 = RATIONAL(-1.0,60.0);
+      t417 = y*t418;
+      t358 = t403*t417;
+      t414 = t413*x;
+      t373 = t403*t414;
+      t438 = y*t373+x*t358;
+      t444 = x*y;
+      t465 = t438+RATIONAL(37.0,300.0)*t444;
+      t405 = RATIONAL(1.0,60.0);
+      t364 = t405*t414;
+      t367 = t405*t417;
+      t439 = y*t364+x*t367;
+      t464 = t439+RATIONAL(-37.0,300.0)*t444;
+      t436 = t413*t418;
+      t463 = RATIONAL(-99.0,1225.0)+RATIONAL(-1.0,98.0)*t436;
+      t410 = RATIONAL(-1.0,30.0);
+      t359 = t410*t417;
+      t442 = y*t413;
+      t409 = RATIONAL(1.0,70.0);
+      t318 = t409*t442;
+      t392 = RATIONAL(11.0,105.0);
+      t462 = t359+t318+t392*y;
+      t411 = RATIONAL(1.0,30.0);
+      t360 = t411*t417;
+      t441 = y*t414;
+      t400 = RATIONAL(-1.0,120.0);
+      t298 = RATIONAL(-21.0,200.0)*t444;
+      t461 = x*t360+t298+t400*t441;
+      t366 = t410*t414;
+      t443 = x*t418;
+      t395 = RATIONAL(31.0,210.0);
+      t401 = RATIONAL(-1.0,140.0);
+      t460 = t395*x+t366+t401*t443;
+      t440 = x*t417;
+      t296 = RATIONAL(21.0,200.0)*t444;
+      t399 = RATIONAL(1.0,120.0);
+      t459 = t296+y*t366+t399*t440;
+      t372 = t411*t414;
+      t458 = t400*t440+y*t372+t298;
+      t333 = t409*t443;
+      t457 = t366+t333+t392*x;
+      t396 = RATIONAL(-11.0,105.0);
+      t412 = RATIONAL(-1.0,70.0);
+      t327 = t412*t443;
+      t456 = t396*x+t372+t327;
+      t407 = RATIONAL(1.0,35.0);
+      t331 = t407*t442;
+      t391 = RATIONAL(-31.0,420.0);
+      t455 = t331+t391*y+t367;
+      t402 = RATIONAL(-1.0,35.0);
+      t313 = t402*t442;
+      t397 = RATIONAL(31.0,420.0);
+      t454 = t358+t397*y+t313;
+      t325 = t402*t443;
+      t453 = t397*x+t325+t373;
+      t324 = t407*t443;
+      t452 = t364+t324+t391*x;
+      t311 = t412*t442;
+      t451 = t360+t311+t396*y;
+      t450 = x*t359+t296+t399*t441;
+      t387 = RATIONAL(-31.0,210.0);
+      t398 = RATIONAL(1.0,140.0);
+      t449 = t372+t387*x+t398*t443;
+      t415 = t413*t413;
+      t371 = t410*t415;
+      t416 = t418*t418;
+      t357 = t399*t416;
+      t383 = RATIONAL(71.0,5880.0);
+      t382 = RATIONAL(113.0,735.0);
+      t448 = t371+t382*t413+t383*t418+t357+t463;
+      t317 = RATIONAL(1.0,49.0)*t436;
+      t356 = t399*t415;
+      t447 = t317+RATIONAL(51.0,1225.0)+t356+t357+t466*RATIONAL(-289.0,5880.0);
+      t362 = t410*t416;
+      t446 = t382*t418+t356+t383*t413+t362+t463;
+      t445 = t362+RATIONAL(1.0,196.0)*t436+t371+RATIONAL(-24.0,1225.0)+t466*
+RATIONAL(181.0,1470.0);
+      t437 = RATIONAL(-1.0,49.0)*t436+RATIONAL(96.0,1225.0);
+      t435 = RATIONAL(1.0,98.0)*t436+RATIONAL(246.0,1225.0);
+      t393 = RATIONAL(-41.0,196.0);
+      t404 = RATIONAL(1.0,20.0);
+      t369 = t404*t415;
+      t376 = RATIONAL(191.0,5880.0);
+      t434 = t376*t418+t369+t357+t393*t413+t437;
+      t365 = t404*t416;
+      t384 = RATIONAL(-53.0,196.0);
+      t386 = RATIONAL(83.0,735.0);
+      t433 = t365+t371+t384*t418+t386*t413+t435;
+      t432 = t365+t376*t413+t356+t393*t418+t437;
+      t431 = t369+t384*t413+t362+t386*t418+t435;
+      t430 = RATIONAL(-11.0,150.0)*t444+t439+t447;
+      t429 = t359+t395*y+t401*t442+t445;
+      t428 = t360+t387*y+t398*t442+t445;
+      t427 = RATIONAL(11.0,150.0)*t444+t438+t447;
+      t406 = RATIONAL(1.0,84.0);
+      t426 = t367+t311+t406*y+t448;
+      t425 = t364+t406*x+t327+t446;
+      t408 = RATIONAL(-1.0,84.0);
+      t424 = t358+t318+t408*y+t448;
+      t423 = t333+t408*x+t373+t446;
+      t394 = RATIONAL(-17.0,420.0);
+      t389 = RATIONAL(17.0,105.0);
+      t388 = RATIONAL(17.0,420.0);
+      t385 = RATIONAL(-17.0,105.0);
+      coeff_I_m2_m2 = t430+t453+t454;
+      coeff_I_m1_m2 = t424+t456+t459;
+      coeff_I_0_m2 = t331+t394*y+t358+t434;
+      coeff_I_p1_m2 = t424+t457+t458;
+      coeff_I_p2_m2 = t427+t452+t454;
+      coeff_I_m2_m1 = t423+t450+t451;
+      coeff_I_m1_m1 = t428+t449+t465;
+      coeff_I_0_m1 = t360+t385*y+t318+t431;
+      coeff_I_p1_m1 = t428+t460+t464;
+      coeff_I_p2_m1 = t425+t451+t461;
+      coeff_I_m2_0 = t373+t394*x+t324+t432;
+      coeff_I_m1_0 = t385*x+t333+t372+t433;
+      coeff_I_0_0 = t317+t365+t369+RATIONAL(541.0,1225.0)+t466*RATIONAL(-57.0,
+196.0);
+      coeff_I_p1_0 = t389*x+t327+t366+t433;
+      coeff_I_p2_0 = t325+t364+t388*x+t432;
+      coeff_I_m2_p1 = t423+t461+t462;
+      coeff_I_m1_p1 = t429+t449+t464;
+      coeff_I_0_p1 = t359+t389*y+t311+t431;
+      coeff_I_p1_p1 = t429+t460+t465;
+      coeff_I_p2_p1 = t425+t450+t462;
+      coeff_I_m2_p2 = t427+t453+t455;
+      coeff_I_m1_p2 = t426+t456+t458;
+      coeff_I_0_p2 = t313+t388*y+t367+t434;
+      coeff_I_p1_p2 = t426+t457+t459;
+      coeff_I_p2_p2 = t430+t452+t455;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c
new file mode 100644
index 0000000..83ba345
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c
@@ -0,0 +1,147 @@
+fp t531,
+   t534,
+   t535,
+   t513,
+   t570,
+   t508,
+   t569,
+   t537,
+   t536,
+   t502,
+   t552,
+   t568,
+   t512,
+   t555,
+   t567,
+   t556,
+   t566,
+   t509,
+   t533,
+   t510,
+   t565,
+   t530,
+   t506,
+   t507,
+   t564,
+   t563,
+   t471,
+   t562,
+   t561,
+   t523,
+   t511,
+   t560,
+   t528,
+   t504,
+   t559,
+   t558,
+   t557,
+   t478,
+   t554,
+   t553,
+   t516,
+   t476,
+   t551,
+   t483,
+   t550,
+   t549,
+   t548,
+   t547,
+   t546,
+   t545,
+   t544,
+   t543,
+   t542,
+   t541,
+   t540,
+   t484,
+   t477,
+   t475,
+   t473,
+   t469,
+   t467;
+      t531 = RATIONAL(1.0,30.0);
+      t534 = x*x;
+      t535 = t534*x;
+      t513 = t531*t535;
+      t570 = t513+RATIONAL(71.0,2940.0)*x;
+      t508 = RATIONAL(-2.0,15.0)*t535;
+      t569 = t508+RATIONAL(226.0,735.0)*x;
+      t537 = y*y;
+      t536 = y*t537;
+      t502 = RATIONAL(1.0,20.0)*t534;
+      t552 = RATIONAL(1.0,60.0)*t536+y*t502;
+      t568 = t552+RATIONAL(-37.0,300.0)*y;
+      t512 = RATIONAL(-1.0,20.0)*t534;
+      t555 = y*t512+RATIONAL(-1.0,60.0)*t536;
+      t567 = t555+RATIONAL(37.0,300.0)*y;
+      t556 = x*t537;
+      t566 = t508+RATIONAL(181.0,735.0)*x+RATIONAL(1.0,98.0)*t556;
+      t509 = RATIONAL(1.0,10.0)*t534;
+      t533 = RATIONAL(-1.0,70.0);
+      t510 = t533*t537;
+      t565 = t509+t510+RATIONAL(-11.0,105.0);
+      t530 = RATIONAL(1.0,70.0);
+      t506 = t530*t537;
+      t507 = RATIONAL(-1.0,10.0)*t534;
+      t564 = t506+t507+RATIONAL(11.0,105.0);
+      t563 = t510+RATIONAL(1.0,84.0)+t502;
+      t471 = RATIONAL(2.0,49.0)*t556;
+      t562 = t513+t471+RATIONAL(-289.0,2940.0)*x;
+      t561 = t506+t512+RATIONAL(-1.0,84.0);
+      t523 = RATIONAL(-1.0,35.0);
+      t511 = t523*t537;
+      t560 = t511+t512+RATIONAL(31.0,420.0);
+      t528 = RATIONAL(1.0,35.0);
+      t504 = t528*t537;
+      t559 = RATIONAL(-31.0,420.0)+t502+t504;
+      t558 = x*y;
+      t557 = y*t534;
+      t478 = RATIONAL(-1.0,49.0)*t556;
+      t554 = RATIONAL(-21.0,200.0)*y+t478;
+      t553 = RATIONAL(21.0,200.0)*y+t478;
+      t516 = RATIONAL(1.0,5.0)*t535;
+      t476 = RATIONAL(-2.0,49.0)*t556;
+      t551 = t516+t476+RATIONAL(-41.0,98.0)*x;
+      t483 = RATIONAL(1.0,49.0)*t556;
+      t550 = t516+t483+RATIONAL(-53.0,98.0)*x;
+      t549 = t508+RATIONAL(166.0,735.0)*x+t483;
+      t548 = t513+t476+RATIONAL(191.0,2940.0)*x;
+      t547 = t509+RATIONAL(-31.0,210.0)+RATIONAL(1.0,140.0)*t537+t566;
+      t546 = RATIONAL(11.0,150.0)*y+t555+t562;
+      t545 = y*t509+RATIONAL(-1.0,120.0)*t536+t554+t569;
+      t544 = t507+RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t537+t566;
+      t543 = y*t507+RATIONAL(1.0,120.0)*t536+t553+t569;
+      t542 = t531*t536+RATIONAL(-1.0,40.0)*t557+t554+t570;
+      t541 = RATIONAL(-11.0,150.0)*y+t552+t562;
+      t540 = RATIONAL(1.0,40.0)*t557+RATIONAL(-1.0,30.0)*t536+t553+t570;
+      t484 = t530*t558;
+      t477 = t528*t558;
+      t475 = RATIONAL(2.0,35.0)*t558;
+      t473 = t533*t558;
+      t469 = RATIONAL(-2.0,35.0)*t558;
+      t467 = t523*t558;
+      coeff_dx_m2_m2 = t469+t541+t560;
+      coeff_dx_m1_m2 = t477+t543+t565;
+      coeff_dx_0_m2 = t475+t551;
+      coeff_dx_p1_m2 = t477+t545+t564;
+      coeff_dx_p2_m2 = t469+t546+t559;
+      coeff_dx_m2_m1 = t467+t540+t561;
+      coeff_dx_m1_m1 = t484+t547+t567;
+      coeff_dx_0_m1 = t477+t550;
+      coeff_dx_p1_m1 = t484+t544+t568;
+      coeff_dx_p2_m1 = t467+t542+t563;
+      coeff_dx_m2_0 = t504+t512+RATIONAL(-17.0,420.0)+t548;
+      coeff_dx_m1_0 = RATIONAL(-17.0,105.0)+t509+t506+t549;
+      coeff_dx_0_0 = RATIONAL(-57.0,98.0)*x+t516+t471;
+      coeff_dx_p1_0 = t510+RATIONAL(17.0,105.0)+t507+t549;
+      coeff_dx_p2_0 = t511+RATIONAL(17.0,420.0)+t502+t548;
+      coeff_dx_m2_p1 = t477+t542+t561;
+      coeff_dx_m1_p1 = t473+t547+t568;
+      coeff_dx_0_p1 = t467+t550;
+      coeff_dx_p1_p1 = t473+t544+t567;
+      coeff_dx_p2_p1 = t477+t540+t563;
+      coeff_dx_m2_p2 = t475+t546+t560;
+      coeff_dx_m1_p2 = t467+t545+t565;
+      coeff_dx_0_p2 = t469+t551;
+      coeff_dx_p1_p2 = t467+t543+t564;
+      coeff_dx_p2_p2 = t475+t541+t559;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c
new file mode 100644
index 0000000..385d8a7
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c
@@ -0,0 +1,111 @@
+fp t712,
+   t692,
+   t707,
+   t711,
+   t693,
+   t729,
+   t687,
+   t728,
+   t697,
+   t727,
+   t726,
+   t725,
+   t686,
+   t694,
+   t695,
+   t698,
+   t689,
+   t724,
+   t681,
+   t723,
+   t696,
+   t722,
+   t721,
+   t688,
+   t720,
+   t691,
+   t719,
+   t685,
+   t718,
+   t717,
+   t716,
+   t690,
+   t715,
+   t714,
+   t713,
+   t684,
+   t683,
+   t680,
+   t679,
+   t678,
+   t677,
+   t676,
+   t675;
+      t712 = y*y;
+      t692 = RATIONAL(2.0,49.0)*t712;
+      t707 = RATIONAL(1.0,10.0);
+      t711 = x*x;
+      t693 = t707*t711;
+      t729 = t692+t693+RATIONAL(-289.0,2940.0);
+      t687 = RATIONAL(-1.0,49.0)*t712;
+      t728 = t693+RATIONAL(71.0,2940.0)+t687;
+      t697 = RATIONAL(-2.0,5.0)*t711;
+      t727 = t697+RATIONAL(226.0,735.0)+t687;
+      t726 = t697+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t712;
+      t725 = x*y;
+      t686 = t707*x;
+      t694 = RATIONAL(-1.0,10.0)*x;
+      t695 = RATIONAL(-1.0,5.0)*x;
+      t698 = RATIONAL(1.0,5.0)*x;
+      t689 = RATIONAL(-2.0,49.0)*t712;
+      t724 = t693+RATIONAL(191.0,2940.0)+t689;
+      t681 = RATIONAL(1.0,49.0)*t712;
+      t723 = t697+RATIONAL(166.0,735.0)+t681;
+      t696 = RATIONAL(3.0,5.0)*t711;
+      t722 = t696+RATIONAL(-53.0,98.0)+t681;
+      t721 = t696+RATIONAL(-41.0,98.0)+t689;
+      t688 = RATIONAL(2.0,35.0)*y;
+      t720 = t688+t729;
+      t691 = RATIONAL(-2.0,35.0)*y;
+      t719 = t691+t729;
+      t685 = RATIONAL(-1.0,35.0)*y;
+      t718 = t685+t727;
+      t717 = t686+t728;
+      t716 = t695+t726;
+      t690 = RATIONAL(1.0,35.0)*y;
+      t715 = t690+t727;
+      t714 = t698+t726;
+      t713 = t694+t728;
+      t684 = RATIONAL(1.0,70.0)*y;
+      t683 = RATIONAL(-1.0,70.0)*y;
+      t680 = y*t695;
+      t679 = y*t698;
+      t678 = RATIONAL(-1.0,20.0)*t725;
+      t677 = y*t686;
+      t676 = y*t694;
+      t675 = RATIONAL(1.0,20.0)*t725;
+      coeff_dxx_m2_m2 = t677+t694+t719;
+      coeff_dxx_m1_m2 = t680+t698+t715;
+      coeff_dxx_0_m2 = t688+t721;
+      coeff_dxx_p1_m2 = t679+t695+t715;
+      coeff_dxx_p2_m2 = t686+t676+t719;
+      coeff_dxx_m2_m1 = t685+t675+t713;
+      coeff_dxx_m1_m1 = t676+t684+t714;
+      coeff_dxx_0_m1 = t690+t722;
+      coeff_dxx_p1_m1 = t684+t677+t716;
+      coeff_dxx_p2_m1 = t685+t678+t717;
+      coeff_dxx_m2_0 = t694+t724;
+      coeff_dxx_m1_0 = t698+t723;
+      coeff_dxx_0_0 = t692+RATIONAL(-57.0,98.0)+t696;
+      coeff_dxx_p1_0 = t695+t723;
+      coeff_dxx_p2_0 = t686+t724;
+      coeff_dxx_m2_p1 = t690+t678+t713;
+      coeff_dxx_m1_p1 = t677+t683+t714;
+      coeff_dxx_0_p1 = t685+t722;
+      coeff_dxx_p1_p1 = t683+t676+t716;
+      coeff_dxx_p2_p1 = t690+t675+t717;
+      coeff_dxx_m2_p2 = t676+t694+t720;
+      coeff_dxx_m1_p2 = t698+t679+t718;
+      coeff_dxx_0_p2 = t691+t721;
+      coeff_dxx_p1_p2 = t695+t680+t718;
+      coeff_dxx_p2_p2 = t686+t677+t720;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c
new file mode 100644
index 0000000..9388906
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c
@@ -0,0 +1,119 @@
+fp t778,
+   t777,
+   t796,
+   t766,
+   t735,
+   t772,
+   t755,
+   t795,
+   t751,
+   t748,
+   t794,
+   t793,
+   t792,
+   t791,
+   t733,
+   t790,
+   t789,
+   t788,
+   t787,
+   t765,
+   t770,
+   t786,
+   t734,
+   t785,
+   t767,
+   t771,
+   t784,
+   t783,
+   t782,
+   t730,
+   t781,
+   t780,
+   t779,
+   t776,
+   t775,
+   t773,
+   t769,
+   t756,
+   t754,
+   t750,
+   t749,
+   t746,
+   t745,
+   t743,
+   t738,
+   t732,
+   t731;
+      t778 = y*y;
+      t777 = x*x;
+      t796 = t778+t777;
+      t766 = RATIONAL(-1.0,35.0);
+      t735 = t766*x;
+      t772 = RATIONAL(1.0,35.0);
+      t755 = t772*y;
+      t795 = t735+t755;
+      t751 = t766*y;
+      t748 = t772*x;
+      t794 = t751+t748;
+      t793 = t735+t751;
+      t792 = t755+t748;
+      t791 = x*y;
+      t733 = RATIONAL(-2.0,49.0)*t791;
+      t790 = t733+RATIONAL(-21.0,200.0);
+      t789 = t733+RATIONAL(21.0,200.0);
+      t788 = t796*RATIONAL(-1.0,20.0);
+      t787 = t796*RATIONAL(1.0,20.0);
+      t765 = RATIONAL(-1.0,40.0);
+      t770 = RATIONAL(1.0,10.0);
+      t786 = t765*t777+t770*t778+t790;
+      t734 = RATIONAL(4.0,49.0)*t791;
+      t785 = t734+RATIONAL(11.0,150.0)+t788;
+      t767 = RATIONAL(1.0,40.0);
+      t771 = RATIONAL(-1.0,10.0);
+      t784 = t771*t778+t767*t777+t789;
+      t783 = t771*t777+t767*t778+t789;
+      t782 = t734+RATIONAL(-11.0,150.0)+t787;
+      t730 = RATIONAL(1.0,49.0)*t791;
+      t781 = RATIONAL(37.0,300.0)+t730+t788;
+      t780 = RATIONAL(-37.0,300.0)+t730+t787;
+      t779 = t765*t778+t770*t777+t790;
+      t776 = RATIONAL(-2.0,35.0);
+      t775 = RATIONAL(-1.0,70.0);
+      t773 = RATIONAL(1.0,70.0);
+      t769 = RATIONAL(2.0,35.0);
+      t756 = t776*y;
+      t754 = t769*y;
+      t750 = t776*x;
+      t749 = t775*x;
+      t746 = t773*y;
+      t745 = t775*y;
+      t743 = t769*x;
+      t738 = t773*x;
+      t732 = RATIONAL(-4.0,49.0)*t791;
+      t731 = RATIONAL(2.0,49.0)*t791;
+      coeff_dxy_m2_m2 = t756+t750+t782;
+      coeff_dxy_m1_m2 = t783+t794;
+      coeff_dxy_0_m2 = t743+t732;
+      coeff_dxy_p1_m2 = t779+t792;
+      coeff_dxy_p2_m2 = t750+t754+t785;
+      coeff_dxy_m2_m1 = t784+t795;
+      coeff_dxy_m1_m1 = t738+t746+t781;
+      coeff_dxy_0_m1 = t748+t731;
+      coeff_dxy_p1_m1 = t745+t738+t780;
+      coeff_dxy_p2_m1 = t786+t793;
+      coeff_dxy_m2_0 = t754+t732;
+      coeff_dxy_m1_0 = t755+t731;
+      coeff_dxy_0_0 = t734;
+      coeff_dxy_p1_0 = t751+t731;
+      coeff_dxy_p2_0 = t732+t756;
+      coeff_dxy_m2_p1 = t786+t792;
+      coeff_dxy_m1_p1 = t749+t746+t780;
+      coeff_dxy_0_p1 = t735+t731;
+      coeff_dxy_p1_p1 = t749+t745+t781;
+      coeff_dxy_p2_p1 = t784+t794;
+      coeff_dxy_m2_p2 = t743+t756+t785;
+      coeff_dxy_m1_p2 = t779+t793;
+      coeff_dxy_0_p2 = t750+t732;
+      coeff_dxy_p1_p2 = t783+t795;
+      coeff_dxy_p2_p2 = t754+t743+t782;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c
new file mode 100644
index 0000000..03c68c3
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c
@@ -0,0 +1,151 @@
+fp t641,
+   t611,
+   t638,
+   t639,
+   t657,
+   t674,
+   t608,
+   t656,
+   t673,
+   t661,
+   t589,
+   t672,
+   t671,
+   t635,
+   t640,
+   t609,
+   t670,
+   t574,
+   t669,
+   t614,
+   t598,
+   t668,
+   t632,
+   t615,
+   t667,
+   t634,
+   t619,
+   t666,
+   t617,
+   t596,
+   t665,
+   t627,
+   t606,
+   t664,
+   t637,
+   t616,
+   t663,
+   t662,
+   t660,
+   t618,
+   t659,
+   t658,
+   t586,
+   t655,
+   t578,
+   t654,
+   t620,
+   t653,
+   t652,
+   t651,
+   t650,
+   t649,
+   t648,
+   t647,
+   t646,
+   t645,
+   t644,
+   t588,
+   t583,
+   t582,
+   t579,
+   t576,
+   t572;
+      t641 = y*y;
+      t611 = RATIONAL(1.0,20.0)*t641;
+      t638 = x*x;
+      t639 = t638*x;
+      t657 = RATIONAL(1.0,60.0)*t639+x*t611;
+      t674 = t657+RATIONAL(-37.0,300.0)*x;
+      t608 = RATIONAL(-1.0,20.0)*t641;
+      t656 = RATIONAL(-1.0,60.0)*t639+x*t608;
+      t673 = t656+RATIONAL(37.0,300.0)*x;
+      t661 = y*t638;
+      t589 = RATIONAL(-1.0,49.0)*t661;
+      t672 = t589+RATIONAL(226.0,735.0)*y;
+      t671 = RATIONAL(1.0,98.0)*t661+RATIONAL(181.0,735.0)*y;
+      t635 = RATIONAL(1.0,30.0);
+      t640 = y*t641;
+      t609 = t635*t640;
+      t670 = t589+t609+RATIONAL(71.0,2940.0)*y;
+      t574 = RATIONAL(2.0,49.0)*t661;
+      t669 = RATIONAL(-289.0,2940.0)*y+t574+t609;
+      t614 = RATIONAL(-1.0,10.0)*t641;
+      t598 = RATIONAL(21.0,200.0)*x;
+      t668 = x*t614+RATIONAL(1.0,120.0)*t639+t598;
+      t632 = RATIONAL(1.0,35.0);
+      t615 = t632*t638;
+      t667 = t611+RATIONAL(-31.0,420.0)+t615;
+      t634 = RATIONAL(1.0,70.0);
+      t619 = t634*t638;
+      t666 = t619+RATIONAL(-1.0,84.0)+t608;
+      t617 = RATIONAL(1.0,10.0)*t641;
+      t596 = RATIONAL(-21.0,200.0)*x;
+      t665 = RATIONAL(-1.0,120.0)*t639+x*t617+t596;
+      t627 = RATIONAL(-1.0,35.0);
+      t606 = t627*t638;
+      t664 = RATIONAL(31.0,420.0)+t606+t608;
+      t637 = RATIONAL(-1.0,70.0);
+      t616 = t637*t638;
+      t663 = t611+t616+RATIONAL(1.0,84.0);
+      t662 = x*y;
+      t660 = x*t641;
+      t618 = RATIONAL(-2.0,15.0)*t640;
+      t659 = t617+t618;
+      t658 = t618+t614;
+      t586 = RATIONAL(-2.0,49.0)*t661;
+      t655 = t586+t609+RATIONAL(191.0,2940.0)*y;
+      t578 = RATIONAL(1.0,49.0)*t661;
+      t654 = t618+RATIONAL(166.0,735.0)*y+t578;
+      t620 = RATIONAL(1.0,5.0)*t640;
+      t653 = RATIONAL(-53.0,98.0)*y+t620+t578;
+      t652 = t620+RATIONAL(-41.0,98.0)*y+t586;
+      t651 = RATIONAL(1.0,140.0)*t638+RATIONAL(-31.0,210.0)+t659+t671;
+      t650 = t619+RATIONAL(11.0,105.0)+t658+t672;
+      t649 = RATIONAL(-11.0,150.0)*x+t657+t669;
+      t648 = RATIONAL(11.0,150.0)*x+t656+t669;
+      t647 = t616+RATIONAL(-11.0,105.0)+t659+t672;
+      t646 = t635*t639+RATIONAL(-1.0,40.0)*t660+t596+t670;
+      t645 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t638+t658+t671;
+      t644 = RATIONAL(-1.0,30.0)*t639+RATIONAL(1.0,40.0)*t660+t598+t670;
+      t588 = t627*t662;
+      t583 = t632*t662;
+      t582 = RATIONAL(2.0,35.0)*t662;
+      t579 = t634*t662;
+      t576 = t637*t662;
+      t572 = RATIONAL(-2.0,35.0)*t662;
+      coeff_dy_m2_m2 = t572+t649+t664;
+      coeff_dy_m1_m2 = t588+t644+t666;
+      coeff_dy_0_m2 = RATIONAL(-17.0,420.0)+t615+t608+t655;
+      coeff_dy_p1_m2 = t583+t646+t666;
+      coeff_dy_p2_m2 = t582+t648+t664;
+      coeff_dy_m2_m1 = t583+t647+t668;
+      coeff_dy_m1_m1 = t579+t651+t673;
+      coeff_dy_0_m1 = RATIONAL(-17.0,105.0)+t617+t619+t654;
+      coeff_dy_p1_m1 = t576+t651+t674;
+      coeff_dy_p2_m1 = t588+t647+t665;
+      coeff_dy_m2_0 = t582+t652;
+      coeff_dy_m1_0 = t583+t653;
+      coeff_dy_0_0 = t574+RATIONAL(-57.0,98.0)*y+t620;
+      coeff_dy_p1_0 = t588+t653;
+      coeff_dy_p2_0 = t572+t652;
+      coeff_dy_m2_p1 = t583+t650+t665;
+      coeff_dy_m1_p1 = t579+t645+t674;
+      coeff_dy_0_p1 = RATIONAL(17.0,105.0)+t616+t614+t654;
+      coeff_dy_p1_p1 = t576+t645+t673;
+      coeff_dy_p2_p1 = t588+t650+t668;
+      coeff_dy_m2_p2 = t572+t648+t667;
+      coeff_dy_m1_p2 = t588+t646+t663;
+      coeff_dy_0_p2 = t611+RATIONAL(17.0,420.0)+t606+t655;
+      coeff_dy_p1_p2 = t583+t644+t663;
+      coeff_dy_p2_p2 = t582+t649+t667;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c
new file mode 100644
index 0000000..cc0ccef
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c
@@ -0,0 +1,115 @@
+fp t833,
+   t808,
+   t853,
+   t807,
+   t852,
+   t834,
+   t818,
+   t851,
+   t850,
+   t849,
+   t829,
+   t815,
+   t817,
+   t806,
+   t820,
+   t804,
+   t848,
+   t847,
+   t812,
+   t846,
+   t819,
+   t845,
+   t803,
+   t844,
+   t843,
+   t798,
+   t842,
+   t841,
+   t840,
+   t839,
+   t838,
+   t837,
+   t800,
+   t836,
+   t835,
+   t816,
+   t814,
+   t811,
+   t810,
+   t809,
+   t805,
+   t802,
+   t801,
+   t799,
+   t797;
+      t833 = x*x;
+      t808 = RATIONAL(-1.0,49.0)*t833;
+      t853 = t808+RATIONAL(71.0,2940.0);
+      t807 = RATIONAL(2.0,49.0)*t833;
+      t852 = t807+RATIONAL(-289.0,2940.0);
+      t834 = y*y;
+      t818 = RATIONAL(-2.0,5.0)*t834;
+      t851 = t818+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t833;
+      t850 = t818+t808+RATIONAL(226.0,735.0);
+      t849 = x*y;
+      t829 = RATIONAL(1.0,10.0);
+      t815 = t829*y;
+      t817 = RATIONAL(-1.0,5.0)*y;
+      t806 = RATIONAL(-1.0,10.0)*y;
+      t820 = RATIONAL(1.0,5.0)*y;
+      t804 = t829*t834;
+      t848 = t815+t804;
+      t847 = t806+t804;
+      t812 = RATIONAL(1.0,49.0)*t833;
+      t846 = t818+RATIONAL(166.0,735.0)+t812;
+      t819 = RATIONAL(3.0,5.0)*t834;
+      t845 = t819+RATIONAL(-53.0,98.0)+t812;
+      t803 = RATIONAL(-2.0,49.0)*t833;
+      t844 = t803+RATIONAL(191.0,2940.0)+t804;
+      t843 = t803+t819+RATIONAL(-41.0,98.0);
+      t798 = x*t815;
+      t842 = t798+t851;
+      t841 = t848+t853;
+      t840 = t847+t853;
+      t839 = t847+t852;
+      t838 = t820+t850;
+      t837 = t817+t850;
+      t800 = x*t806;
+      t836 = t800+t851;
+      t835 = t848+t852;
+      t816 = RATIONAL(2.0,35.0)*x;
+      t814 = RATIONAL(1.0,70.0)*x;
+      t811 = RATIONAL(-1.0,70.0)*x;
+      t810 = RATIONAL(1.0,35.0)*x;
+      t809 = RATIONAL(-1.0,35.0)*x;
+      t805 = RATIONAL(-2.0,35.0)*x;
+      t802 = x*t817;
+      t801 = x*t820;
+      t799 = RATIONAL(-1.0,20.0)*t849;
+      t797 = RATIONAL(1.0,20.0)*t849;
+      coeff_dyy_m2_m2 = t798+t805+t839;
+      coeff_dyy_m1_m2 = t797+t809+t840;
+      coeff_dyy_0_m2 = t806+t844;
+      coeff_dyy_p1_m2 = t810+t799+t840;
+      coeff_dyy_p2_m2 = t800+t816+t839;
+      coeff_dyy_m2_m1 = t810+t802+t838;
+      coeff_dyy_m1_m1 = t820+t814+t836;
+      coeff_dyy_0_m1 = t820+t846;
+      coeff_dyy_p1_m1 = t811+t820+t842;
+      coeff_dyy_p2_m1 = t809+t801+t838;
+      coeff_dyy_m2_0 = t816+t843;
+      coeff_dyy_m1_0 = t810+t845;
+      coeff_dyy_0_0 = t819+RATIONAL(-57.0,98.0)+t807;
+      coeff_dyy_p1_0 = t809+t845;
+      coeff_dyy_p2_0 = t805+t843;
+      coeff_dyy_m2_p1 = t810+t801+t837;
+      coeff_dyy_m1_p1 = t817+t814+t842;
+      coeff_dyy_0_p1 = t817+t846;
+      coeff_dyy_p1_p1 = t811+t817+t836;
+      coeff_dyy_p2_p1 = t809+t802+t837;
+      coeff_dyy_m2_p2 = t805+t800+t835;
+      coeff_dyy_m1_p2 = t809+t799+t841;
+      coeff_dyy_0_p2 = t815+t844;
+      coeff_dyy_p1_p2 = t797+t810+t841;
+      coeff_dyy_p2_p2 = t816+t798+t835;