regenerate all machine-generated coefficient files

--> now includes stuff to store interpolation coeffs in COEFF(...) experessions --> will be used to compute Jacobian coeffs git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@38 df1f8a13-aa1d-4dd4-9681-27ded5b42416
author: jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> 2002-05-14 15:03:52 +0000
committer: jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> 2002-05-14 15:03:52 +0000
commit: 2177a5948dc567e37964df7f281d3a3ee11eb8c6 (patch)
tree: 9e63d67ade8cadd69cb3bd3423bc8a461d8ed699 /src/GeneralizedPolynomial-Uniform/2d.coeffs
parent: 918fdb59c1bc57ae64c8f26db86dd50f316d557e (diff)
40 files changed, 1550 insertions, 1226 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c
index 221b33b..31d7188 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c
@@ -5,14 +5,14 @@ fp t7,
    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;
+      t7 = RATIONAL(1.0,2.0);
+      t6 = RATIONAL(-1.0,2.0);
+      t5 = RATIONAL(1.0,4.0);
+      t4 = t7*x;
+      t3 = t7*y;
+      t2 = t6*x;
+      t1 = t6*y;
+      coeff_I_0_0 = RATIONAL(3.0,4.0)+t1+t2;
+      coeff_I_p1_0 = t5+t1+t4;
+      coeff_I_0_p1 = t2+t5+t3;
+      coeff_I_p1_p1 = t3+RATIONAL(-1.0,4.0)+t4;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c
index 3d89f01..8ac23a5 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c
@@ -1,8 +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;
+      t9 = RATIONAL(1.0,2.0);
+      t8 = RATIONAL(-1.0,2.0);
+      coeff_dx_0_0 = t8;
+      coeff_dx_p1_0 = t9;
+      coeff_dx_0_p1 = t8;
+      coeff_dx_p1_p1 = t9;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c
index 42819a8..107293d 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c
@@ -1,8 +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;
+      t11 = RATIONAL(1.0,2.0);
+      t10 = RATIONAL(-1.0,2.0);
+      coeff_dy_0_0 = t10;
+      coeff_dy_p1_0 = t10;
+      coeff_dy_0_p1 = t11;
+      coeff_dy_p1_p1 = t11;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c
index eae89ba..8a310e0 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c
@@ -1,53 +1,49 @@
-fp t27,
-   t23,
-   t17,
+fp t26,
+   t27,
+   t18,
    t28,
-   t19,
+   t21,
    t34,
    t33,
    t22,
    t15,
-   t26,
+   t23,
    t32,
-   t21,
+   t17,
    t31,
-   t24,
-   t16,
    t30,
-   t14,
    t29,
+   t25,
    t20,
-   t18,
-   t13,
-   t12;
+   t19,
+   t16,
+   t14;
+      t26 = RATIONAL(1.0,6.0);
       t27 = x*x;
-      t23 = RATIONAL(1.0,6.0);
-      t17 = t23*t27;
+      t18 = t26*t27;
       t28 = y*y;
-      t19 = t23*t28;
-      t34 = t17+t19+RATIONAL(-1.0,9.0);
+      t21 = t26*t28;
+      t34 = t18+t21+RATIONAL(-1.0,9.0);
       t33 = x*y;
       t22 = RATIONAL(-1.0,3.0);
       t15 = t22*t28;
-      t26 = RATIONAL(2.0,9.0);
-      t32 = t15+t17+t26;
-      t21 = t22*t27;
-      t31 = t19+t21+t26;
-      t24 = RATIONAL(-1.0,6.0);
-      t16 = t24*y;
-      t30 = t16+t34;
-      t14 = t23*y;
-      t29 = t14+t34;
-      t20 = t23*x;
-      t18 = t24*x;
-      t13 = RATIONAL(-1.0,4.0)*t33;
-      t12 = RATIONAL(1.0,4.0)*t33;
-      coeff_I_m1_m1 = t18+t12+t30;
-      coeff_I_0_m1 = t16+t31;
-      coeff_I_p1_m1 = t20+t13+t30;
-      coeff_I_m1_0 = t18+t32;
-      coeff_I_0_0 = t15+t21+RATIONAL(5.0,9.0);
-      coeff_I_p1_0 = t20+t32;
-      coeff_I_m1_p1 = t18+t13+t29;
+      t23 = RATIONAL(2.0,9.0);
+      t32 = t15+t18+t23;
+      t17 = t22*t27;
+      t31 = t17+t21+t23;
+      t30 = RATIONAL(1.0,4.0)*t33+t34;
+      t29 = RATIONAL(-1.0,4.0)*t33+t34;
+      t25 = RATIONAL(-1.0,6.0);
+      t20 = t25*y;
+      t19 = t25*x;
+      t16 = t26*x;
+      t14 = t26*y;
+      coeff_I_m1_m1 = t20+t19+t30;
+      coeff_I_0_m1 = t20+t31;
+      coeff_I_p1_m1 = t20+t16+t29;
+      coeff_I_m1_0 = t19+t32;
+      coeff_I_0_0 = t15+RATIONAL(5.0,9.0)+t17;
+      coeff_I_p1_0 = t16+t32;
+      coeff_I_m1_p1 = t19+t14+t29;
       coeff_I_0_p1 = t14+t31;
-      coeff_I_p1_p1 = t12+t20+t29;
+      coeff_I_p1_p1 = t16+t14+t30;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c
index 795310a..66a8054 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c
@@ -1,21 +1,21 @@
-fp t35,
+fp t36,
    t42,
    t41,
    t40,
    t39,
-   t36;
-      t35 = RATIONAL(1.0,3.0)*x;
-      t42 = t35+RATIONAL(-1.0,4.0)*y;
-      t41 = t35+RATIONAL(1.0,4.0)*y;
-      t40 = RATIONAL(-1.0,6.0);
-      t39 = RATIONAL(1.0,6.0);
-      t36 = RATIONAL(-2.0,3.0)*x;
-      coeff_dx_m1_m1 = t40+t41;
-      coeff_dx_0_m1 = t36;
-      coeff_dx_p1_m1 = t39+t42;
-      coeff_dx_m1_0 = t40+t35;
-      coeff_dx_0_0 = t36;
-      coeff_dx_p1_0 = t35+t39;
-      coeff_dx_m1_p1 = t40+t42;
-      coeff_dx_0_p1 = t36;
-      coeff_dx_p1_p1 = t39+t41;
+   t35;
+      t36 = RATIONAL(1.0,3.0)*x;
+      t42 = t36+RATIONAL(-1.0,4.0)*y;
+      t41 = t36+RATIONAL(1.0,4.0)*y;
+      t40 = RATIONAL(1.0,6.0);
+      t39 = RATIONAL(-1.0,6.0);
+      t35 = RATIONAL(-2.0,3.0)*x;
+      coeff_dx_m1_m1 = t39+t41;
+      coeff_dx_0_m1 = t35;
+      coeff_dx_p1_m1 = t40+t42;
+      coeff_dx_m1_0 = t39+t36;
+      coeff_dx_0_0 = t35;
+      coeff_dx_p1_0 = t36+t40;
+      coeff_dx_m1_p1 = t39+t42;
+      coeff_dx_0_p1 = t35;
+      coeff_dx_p1_p1 = t40+t41;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c
index c795e1d..299cc38 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c
@@ -2,14 +2,14 @@ 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;
+      t54 = RATIONAL(-1.0,4.0);
+      t53 = RATIONAL(1.0,4.0);
+      coeff_dxy_m1_m1 = t53;
       coeff_dxy_0_m1 = t55;
-      coeff_dxy_p1_m1 = t53;
+      coeff_dxy_p1_m1 = t54;
       coeff_dxy_m1_0 = t55;
       coeff_dxy_0_0 = t55;
       coeff_dxy_p1_0 = t55;
-      coeff_dxy_m1_p1 = t53;
+      coeff_dxy_m1_p1 = t54;
       coeff_dxy_0_p1 = t55;
-      coeff_dxy_p1_p1 = t54;
+      coeff_dxy_p1_p1 = t53;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c
index c4f9522..ae6d9d9 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c
@@ -1,21 +1,21 @@
-fp t43,
+fp t45,
    t50,
    t49,
-   t46,
-   t45,
+   t48,
+   t47,
    t44;
-      t43 = RATIONAL(1.0,3.0)*y;
-      t50 = RATIONAL(1.0,6.0)+t43;
-      t49 = RATIONAL(-1.0,6.0)+t43;
-      t46 = RATIONAL(-1.0,4.0)*x;
-      t45 = RATIONAL(-2.0,3.0)*y;
-      t44 = RATIONAL(1.0,4.0)*x;
-      coeff_dy_m1_m1 = t44+t49;
-      coeff_dy_0_m1 = t49;
-      coeff_dy_p1_m1 = t46+t49;
-      coeff_dy_m1_0 = t45;
-      coeff_dy_0_0 = t45;
-      coeff_dy_p1_0 = t45;
-      coeff_dy_m1_p1 = t46+t50;
-      coeff_dy_0_p1 = t50;
-      coeff_dy_p1_p1 = t44+t50;
+      t45 = RATIONAL(1.0,3.0)*y;
+      t50 = RATIONAL(-1.0,4.0)*x+t45;
+      t49 = t45+RATIONAL(1.0,4.0)*x;
+      t48 = RATIONAL(1.0,6.0);
+      t47 = RATIONAL(-1.0,6.0);
+      t44 = RATIONAL(-2.0,3.0)*y;
+      coeff_dy_m1_m1 = t47+t49;
+      coeff_dy_0_m1 = t47+t45;
+      coeff_dy_p1_m1 = t47+t50;
+      coeff_dy_m1_0 = t44;
+      coeff_dy_0_0 = t44;
+      coeff_dy_p1_0 = t44;
+      coeff_dy_m1_p1 = t48+t50;
+      coeff_dy_0_p1 = t45+t48;
+      coeff_dy_p1_p1 = t48+t49;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c
index 5a2fe9b..f7614d5 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c
@@ -1,70 +1,72 @@
-fp t159,
+fp t161,
+   t131,
+   t160,
    t132,
+   t159,
+   t130,
+   t129,
+   t87,
+   t124,
+   t71,
    t158,
-   t131,
+   t126,
+   t83,
+   t123,
+   t67,
    t157,
-   t125,
-   t63,
-   t118,
-   t130,
-   t76,
+   t128,
+   t133,
+   t84,
+   t66,
    t156,
-   t126,
-   t82,
    t122,
-   t70,
+   t65,
    t155,
-   t133,
-   t78,
-   t67,
+   t119,
+   t79,
    t154,
+   t70,
+   t76,
    t153,
-   t128,
-   t84,
-   t124,
-   t66,
+   t127,
+   t63,
    t152,
-   t151,
-   t83,
+   t82,
    t62,
+   t151,
    t150,
-   t87,
    t69,
    t149,
+   t85,
    t148,
-   t121,
-   t120,
-   t68,
+   t86,
    t147,
-   t65,
    t146,
-   t129,
-   t86,
    t145,
    t144,
-   t85,
    t143,
-   t123,
    t142,
+   t125,
    t141,
    t140,
    t139,
-   t105,
    t138,
+   t105,
    t137,
    t136,
-   t119,
+   t121,
+   t120,
+   t118,
    t117,
    t116,
-   t114,
    t113,
    t112,
-   t111,
+   t110,
    t108,
    t107,
    t106,
    t104,
-   t103,
+   t102,
    t101,
    t99,
    t98,
@@ -78,115 +80,121 @@ fp t159,
    t90,
    t89,
    t88,
+   t81,
    t80,
-   t79,
+   t78,
+   t77,
    t75,
    t74,
    t73,
    t72;
-      t159 = x*y;
-      t132 = y*y;
-      t158 = x*t132;
-      t131 = x*x;
-      t157 = y*t131;
-      t125 = RATIONAL(-3.0,40.0);
-      t63 = t125*t157;
-      t118 = RATIONAL(1.0,24.0);
-      t130 = t132*y;
-      t76 = t118*t130;
-      t156 = t63+t76;
+      t161 = x*y;
+      t131 = y*y;
+      t160 = x*t131;
+      t132 = x*x;
+      t159 = y*t132;
+      t130 = x*t132;
+      t129 = RATIONAL(1.0,8.0);
+      t87 = t129*t130;
+      t124 = RATIONAL(-1.0,40.0);
+      t71 = t124*t160;
+      t158 = t87+t71;
       t126 = RATIONAL(-1.0,24.0);
-      t82 = t126*t130;
-      t122 = RATIONAL(3.0,40.0);
-      t70 = t122*t157;
-      t155 = t82+t70;
-      t133 = x*t131;
-      t78 = t118*t133;
-      t67 = t122*t158;
-      t154 = t78+t67;
-      t153 = t70+t76;
+      t83 = t126*t130;
+      t123 = RATIONAL(3.0,40.0);
+      t67 = t123*t160;
+      t157 = t83+t67;
       t128 = RATIONAL(-1.0,8.0);
-      t84 = t128*t130;
-      t124 = RATIONAL(-1.0,40.0);
-      t66 = t124*t157;
-      t152 = t84+t66;
-      t151 = t82+t63;
-      t83 = t126*t133;
-      t62 = t125*t158;
-      t150 = t83+t62;
-      t87 = t128*t133;
-      t69 = t124*t158;
+      t133 = y*t131;
+      t84 = t128*t133;
+      t66 = t124*t159;
+      t156 = t84+t66;
+      t122 = RATIONAL(1.0,40.0);
+      t65 = t122*t159;
+      t155 = t84+t65;
+      t119 = RATIONAL(1.0,24.0);
+      t79 = t119*t130;
+      t154 = t67+t79;
+      t70 = t123*t159;
+      t76 = t119*t133;
+      t153 = t70+t76;
+      t127 = RATIONAL(-3.0,40.0);
+      t63 = t127*t160;
+      t152 = t63+t79;
+      t82 = t126*t133;
+      t62 = t127*t159;
+      t151 = t82+t62;
+      t150 = t83+t63;
+      t69 = t122*t160;
       t149 = t87+t69;
-      t148 = t83+t67;
-      t121 = RATIONAL(9.0,80.0);
-      t120 = RATIONAL(1.0,40.0);
-      t68 = t120*t158;
-      t147 = t121*t132+t68;
-      t65 = t120*t157;
-      t146 = t121*t131+t65;
-      t129 = RATIONAL(1.0,8.0);
-      t86 = t129*t133;
+      t85 = t129*t133;
+      t148 = t85+t65;
+      t86 = t128*t130;
+      t147 = t86+t71;
+      t146 = t62+t76;
       t145 = t86+t69;
-      t144 = t78+t62;
-      t85 = t129*t130;
+      t144 = t82+t70;
       t143 = t85+t66;
-      t123 = RATIONAL(-1.0,50.0);
-      t142 = RATIONAL(11.0,50.0)+t123*t159;
-      t141 = RATIONAL(-2.0,25.0)+RATIONAL(2.0,25.0)*t159;
-      t140 = RATIONAL(-13.0,100.0)*t159+t123;
-      t139 = RATIONAL(37.0,200.0)+RATIONAL(-1.0,100.0)*t159;
+      t142 = RATIONAL(-2.0,25.0)+RATIONAL(2.0,25.0)*t161;
+      t125 = RATIONAL(-1.0,50.0);
+      t141 = RATIONAL(-13.0,100.0)*t161+t125;
+      t140 = RATIONAL(13.0,100.0)+RATIONAL(7.0,100.0)*t161;
+      t139 = RATIONAL(11.0,50.0)+t125*t161;
+      t138 = RATIONAL(37.0,200.0)+RATIONAL(-1.0,100.0)*t161;
       t105 = RATIONAL(-17.0,200.0);
-      t138 = t105+RATIONAL(-9.0,100.0)*t159;
-      t137 = RATIONAL(13.0,100.0)+RATIONAL(7.0,100.0)*t159;
-      t136 = t131+t132;
-      t119 = RATIONAL(7.0,80.0);
-      t117 = RATIONAL(3.0,80.0);
-      t116 = RATIONAL(1.0,80.0);
-      t114 = RATIONAL(-1.0,80.0);
+      t137 = RATIONAL(-9.0,100.0)*t161+t105;
+      t136 = t132+t131;
+      t121 = RATIONAL(9.0,80.0);
+      t120 = RATIONAL(7.0,80.0);
+      t118 = RATIONAL(3.0,80.0);
+      t117 = RATIONAL(1.0,80.0);
+      t116 = RATIONAL(-1.0,80.0);
       t113 = RATIONAL(13.0,80.0);
-      t112 = RATIONAL(11.0,80.0);
-      t111 = RATIONAL(-23.0,80.0);
+      t112 = RATIONAL(-23.0,80.0);
+      t110 = RATIONAL(11.0,80.0);
       t108 = RATIONAL(-19.0,80.0);
       t107 = RATIONAL(-17.0,80.0);
       t106 = RATIONAL(-31.0,400.0);
       t104 = RATIONAL(-57.0,400.0);
-      t103 = RATIONAL(-63.0,400.0);
+      t102 = RATIONAL(-63.0,400.0);
       t101 = RATIONAL(-49.0,400.0);
       t99 = RATIONAL(69.0,400.0);
       t98 = RATIONAL(111.0,400.0);
       t97 = RATIONAL(89.0,1200.0);
-      t96 = RATIONAL(117.0,400.0);
-      t95 = RATIONAL(43.0,1200.0);
+      t96 = RATIONAL(43.0,1200.0);
+      t95 = RATIONAL(117.0,400.0);
       t94 = RATIONAL(103.0,400.0);
       t93 = RATIONAL(37.0,1200.0);
-      t92 = RATIONAL(-149.0,1200.0);
-      t91 = RATIONAL(-109.0,1200.0);
-      t90 = RATIONAL(-223.0,1200.0);
-      t89 = RATIONAL(-157.0,1200.0);
-      t88 = RATIONAL(-131.0,1200.0);
-      t80 = t119*t131;
-      t79 = t119*t132;
+      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 = t121*t131;
+      t80 = t120*t132;
+      t78 = t120*t131;
+      t77 = t121*t132;
       t75 = t113*t132;
       t74 = t113*t131;
-      t73 = t112*t132;
-      t72 = t112*t131;
-      coeff_I_m1_m1 = RATIONAL(-23.0,200.0)+t74+t75+t91*y+(RATIONAL(6.0,25.0)*y
-+t91)*x+t150+t151;
-      coeff_I_0_m1 = t73+t90*y+t111*t131+t106*x+t142+t145+t155;
-      coeff_I_p1_m1 = t87+t80+t98*x+t89*y+t140+t147+t155;
-      coeff_I_p2_m1 = t117*t131+t88*x+t79+t97*y+t138+t151+t154;
-      coeff_I_m1_0 = t111*t132+t106*y+t72+t90*x+t142+t143+t148;
-      coeff_I_0_0 = t85+RATIONAL(93.0,200.0)+t86+t104*y+t68+t65+t136*RATIONAL(
--21.0,80.0)+(RATIONAL(-1.0,25.0)*y+t104)*x;
-      coeff_I_p1_0 = t85+t108*t132+t96*x+t103*y+t139+t146+t149;
-      coeff_I_p2_0 = t95*x+t116*t131+t107*t132+t101*y+t137+t143+t144;
-      coeff_I_m1_p1 = t89*x+t84+t79+t98*y+t140+t146+t148;
-      coeff_I_0_p1 = t86+t96*y+t103*x+t108*t131+t139+t147+t152;
-      coeff_I_p1_p1 = t105+t73+t94*y+t72+(t94+RATIONAL(3.0,50.0)*y)*x+t149+t152
+      t73 = t110*t131;
+      t72 = t110*t132;
+      coeff_I_m1_m1 = t74+t92*y+t75+RATIONAL(-23.0,200.0)+(RATIONAL(6.0,25.0)*y
++t92)*x+t150+t151;
+      coeff_I_0_m1 = t91*y+t112*t132+t73+t106*x+t139+t144+t158;
+      coeff_I_p1_m1 = t81+t90*y+t98*x+t80+t141+t144+t145;
+      coeff_I_p2_m1 = t78+t89*x+t97*y+t118*t132+t137+t151+t154;
+      coeff_I_m1_0 = t91*x+t72+t112*t131+t106*y+t139+t143+t157;
+      coeff_I_0_0 = t104*y+RATIONAL(93.0,200.0)+t136*RATIONAL(-21.0,80.0)+(t104
++RATIONAL(-1.0,25.0)*y)*x+t148+t149;
+      coeff_I_p1_0 = t102*y+t108*t131+t77+t95*x+t138+t147+t148;
+      coeff_I_p2_0 = t117*t132+t96*x+t101*y+t107*t131+t140+t143+t152;
+      coeff_I_m1_p1 = t98*y+t90*x+t77+t78+t141+t155+t157;
+      coeff_I_0_p1 = t95*y+t108*t132+t102*x+t81+t138+t149+t156;
+      coeff_I_p1_p1 = t94*y+t73+t105+t72+(RATIONAL(3.0,50.0)*y+t94)*x+t147+t156
 ;
-      coeff_I_p2_p1 = t99*y+t84+t93*x+t65+t75+t114*t131+t141+t144;
-      coeff_I_m1_p2 = t97*x+t117*t132+t80+t88*y+t138+t150+t153;
-      coeff_I_0_p2 = t95*y+t107*t131+t101*x+t116*t132+t137+t145+t156;
-      coeff_I_p1_p2 = t93*y+t68+t87+t114*t132+t74+t99*x+t141+t156;
-      coeff_I_p2_p2 = t92*y+RATIONAL(7.0,200.0)+t136*RATIONAL(-3.0,80.0)+(
-RATIONAL(-3.0,50.0)*y+t92)*x+t153+t154;
+      coeff_I_p2_p1 = t116*t132+t74+t99*y+t93*x+t142+t152+t155;
+      coeff_I_m1_p2 = t89*y+t97*x+t80+t118*t131+t137+t150+t153;
+      coeff_I_0_p2 = t101*x+t107*t132+t117*t131+t96*y+t140+t146+t158;
+      coeff_I_p1_p2 = t93*y+t116*t131+t75+t99*x+t142+t145+t146;
+      coeff_I_p2_p2 = RATIONAL(7.0,200.0)+t88*y+t136*RATIONAL(-3.0,80.0)+(
+RATIONAL(-3.0,50.0)*y+t88)*x+t153+t154;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c
index f6d13a5..300fa79 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c
@@ -1,95 +1,95 @@
-fp t200,
+fp t202,
    t184,
-   t187,
-   t170,
-   t186,
+   t188,
    t178,
+   t189,
+   t183,
+   t201,
+   t187,
+   t173,
+   t182,
+   t200,
+   t185,
+   t176,
    t199,
-   t181,
+   t186,
+   t177,
    t198,
-   t180,
-   t183,
-   t175,
+   t181,
    t197,
-   t185,
-   t172,
+   t164,
    t196,
-   t182,
-   t176,
+   t163,
    t195,
-   t161,
    t194,
-   t179,
+   t165,
    t193,
-   t163,
+   t180,
    t192,
    t162,
    t191,
    t190,
-   t189,
-   t160,
-   t188,
-   t173,
+   t174,
+   t172,
    t171,
    t169,
-   t167,
-   t166,
-   t164;
-      t200 = x*y;
-      t184 = RATIONAL(-1.0,40.0);
-      t187 = y*y;
-      t170 = t184*t187;
-      t186 = x*x;
-      t178 = RATIONAL(3.0,8.0)*t186;
-      t199 = t170+t178;
-      t181 = RATIONAL(-3.0,8.0)*t186;
-      t198 = t170+t181;
-      t180 = RATIONAL(1.0,8.0)*t186;
-      t183 = RATIONAL(3.0,40.0);
-      t175 = t183*t187;
-      t197 = t180+t175;
-      t185 = RATIONAL(-3.0,40.0);
-      t172 = t185*t187;
-      t196 = t172+t180;
-      t182 = RATIONAL(1.0,40.0);
-      t176 = t182*t187;
-      t195 = t176+t181;
-      t161 = RATIONAL(-3.0,20.0)*t200;
-      t194 = t161+RATIONAL(13.0,40.0)*x;
-      t179 = RATIONAL(-1.0,8.0)*t186;
-      t193 = t172+t179;
-      t163 = RATIONAL(-1.0,20.0)*t200;
-      t192 = t163+RATIONAL(11.0,40.0)*x;
-      t162 = RATIONAL(1.0,20.0)*t200;
-      t191 = t162+RATIONAL(9.0,40.0)*x;
-      t190 = t175+t179;
-      t189 = t178+t176;
-      t160 = RATIONAL(3.0,20.0)*t200;
-      t188 = t160+RATIONAL(7.0,40.0)*x;
-      t173 = RATIONAL(2.0,25.0)*y;
-      t171 = RATIONAL(-1.0,50.0)*y;
-      t169 = RATIONAL(-9.0,100.0)*y;
-      t167 = RATIONAL(-1.0,100.0)*y;
-      t166 = RATIONAL(7.0,100.0)*y;
-      t164 = RATIONAL(-13.0,100.0)*y;
-      coeff_dx_m1_m1 = RATIONAL(-109.0,1200.0)+RATIONAL(6.0,25.0)*y+t193+t194;
-      coeff_dx_0_m1 = RATIONAL(-31.0,400.0)+t171+t160+RATIONAL(-23.0,40.0)*x+
-t199;
-      coeff_dx_p1_m1 = t164+RATIONAL(111.0,400.0)+t188+t195;
-      coeff_dx_p2_m1 = t161+t183*x+RATIONAL(-131.0,1200.0)+t169+t197;
-      coeff_dx_m1_0 = t171+RATIONAL(-223.0,1200.0)+t190+t192;
-      coeff_dx_0_0 = RATIONAL(-1.0,25.0)*y+RATIONAL(-57.0,400.0)+RATIONAL(-21.0
-,40.0)*x+t162+t189;
-      coeff_dx_p1_0 = RATIONAL(117.0,400.0)+t167+t191+t198;
-      coeff_dx_p2_0 = t163+t166+t182*x+RATIONAL(43.0,1200.0)+t196;
-      coeff_dx_m1_p1 = RATIONAL(-157.0,1200.0)+t164+t190+t191;
-      coeff_dx_0_p1 = t163+RATIONAL(-19.0,40.0)*x+t167+RATIONAL(-63.0,400.0)+
-t189;
-      coeff_dx_p1_p1 = RATIONAL(103.0,400.0)+RATIONAL(3.0,50.0)*y+t192+t198;
-      coeff_dx_p2_p1 = t173+RATIONAL(37.0,1200.0)+t162+t184*x+t196;
-      coeff_dx_m1_p2 = RATIONAL(89.0,1200.0)+t169+t188+t193;
-      coeff_dx_0_p2 = RATIONAL(-49.0,400.0)+RATIONAL(-17.0,40.0)*x+t166+t161+
-t199;
-      coeff_dx_p1_p2 = RATIONAL(69.0,400.0)+t173+t194+t195;
-      coeff_dx_p2_p2 = t160+RATIONAL(-149.0,1200.0)+RATIONAL(-3.0,50.0)*y+t185*
-x+t197;
+   t168,
+   t166;
+      t202 = x*y;
+      t184 = RATIONAL(1.0,40.0);
+      t188 = y*y;
+      t178 = t184*t188;
+      t189 = x*x;
+      t183 = RATIONAL(-3.0,8.0)*t189;
+      t201 = t178+t183;
+      t187 = RATIONAL(-3.0,40.0);
+      t173 = t187*t188;
+      t182 = RATIONAL(1.0,8.0)*t189;
+      t200 = t173+t182;
+      t185 = RATIONAL(3.0,40.0);
+      t176 = t185*t188;
+      t199 = t176+t182;
+      t186 = RATIONAL(-1.0,40.0);
+      t177 = t186*t188;
+      t198 = t177+t183;
+      t181 = RATIONAL(-1.0,8.0)*t189;
+      t197 = t173+t181;
+      t164 = RATIONAL(1.0,20.0)*t202;
+      t196 = t164+RATIONAL(9.0,40.0)*x;
+      t163 = RATIONAL(-3.0,20.0)*t202;
+      t195 = t163+RATIONAL(13.0,40.0)*x;
+      t194 = t176+t181;
+      t165 = RATIONAL(-1.0,20.0)*t202;
+      t193 = t165+RATIONAL(11.0,40.0)*x;
+      t180 = RATIONAL(3.0,8.0)*t189;
+      t192 = t177+t180;
+      t162 = RATIONAL(3.0,20.0)*t202;
+      t191 = t162+RATIONAL(7.0,40.0)*x;
+      t190 = t178+t180;
+      t174 = RATIONAL(2.0,25.0)*y;
+      t172 = RATIONAL(-1.0,50.0)*y;
+      t171 = RATIONAL(-9.0,100.0)*y;
+      t169 = RATIONAL(7.0,100.0)*y;
+      t168 = RATIONAL(-1.0,100.0)*y;
+      t166 = RATIONAL(-13.0,100.0)*y;
+      coeff_dx_m1_m1 = RATIONAL(6.0,25.0)*y+RATIONAL(-109.0,1200.0)+t195+t197;
+      coeff_dx_0_m1 = RATIONAL(-23.0,40.0)*x+t172+t162+RATIONAL(-31.0,400.0)+
+t192;
+      coeff_dx_p1_m1 = t166+RATIONAL(111.0,400.0)+t191+t201;
+      coeff_dx_p2_m1 = t171+t185*x+t163+RATIONAL(-131.0,1200.0)+t199;
+      coeff_dx_m1_0 = RATIONAL(-223.0,1200.0)+t172+t193+t194;
+      coeff_dx_0_0 = RATIONAL(-21.0,40.0)*x+t164+RATIONAL(-1.0,25.0)*y+RATIONAL
+(-57.0,400.0)+t190;
+      coeff_dx_p1_0 = t168+RATIONAL(117.0,400.0)+t196+t198;
+      coeff_dx_p2_0 = t165+t184*x+t169+RATIONAL(43.0,1200.0)+t200;
+      coeff_dx_m1_p1 = t166+RATIONAL(-157.0,1200.0)+t194+t196;
+      coeff_dx_0_p1 = t165+RATIONAL(-63.0,400.0)+RATIONAL(-19.0,40.0)*x+t168+
+t190;
+      coeff_dx_p1_p1 = RATIONAL(103.0,400.0)+RATIONAL(3.0,50.0)*y+t193+t198;
+      coeff_dx_p2_p1 = t164+t186*x+RATIONAL(37.0,1200.0)+t174+t200;
+      coeff_dx_m1_p2 = t171+RATIONAL(89.0,1200.0)+t191+t197;
+      coeff_dx_0_p2 = t163+RATIONAL(-49.0,400.0)+RATIONAL(-17.0,40.0)*x+t169+
+t192;
+      coeff_dx_p1_p2 = RATIONAL(69.0,400.0)+t174+t195+t201;
+      coeff_dx_p2_p2 = RATIONAL(-149.0,1200.0)+t162+t187*x+RATIONAL(-3.0,50.0)*
+y+t199;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c
index b113dfe..f28076b 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c
@@ -1,40 +1,40 @@
-fp t244,
+fp t247,
+   t259,
+   t244,
+   t258,
+   t245,
    t257,
-   t243,
+   t246,
    t256,
-   t242,
-   t255,
-   t245,
-   t254,
+   t251,
+   t250,
    t249,
-   t248,
-   t247,
-   t246;
-      t244 = RATIONAL(3.0,20.0)*y;
-      t257 = t244+RATIONAL(7.0,40.0);
-      t243 = RATIONAL(1.0,20.0)*y;
-      t256 = RATIONAL(9.0,40.0)+t243;
-      t242 = RATIONAL(-1.0,20.0)*y;
-      t255 = t242+RATIONAL(11.0,40.0);
-      t245 = RATIONAL(-3.0,20.0)*y;
-      t254 = t245+RATIONAL(13.0,40.0);
-      t249 = RATIONAL(-1.0,4.0)*x;
-      t248 = RATIONAL(3.0,4.0)*x;
-      t247 = RATIONAL(-3.0,4.0)*x;
-      t246 = RATIONAL(1.0,4.0)*x;
-      coeff_dxx_m1_m1 = t249+t254;
-      coeff_dxx_0_m1 = RATIONAL(-23.0,40.0)+t248+t244;
-      coeff_dxx_p1_m1 = t247+t257;
-      coeff_dxx_p2_m1 = t246+RATIONAL(3.0,40.0)+t245;
-      coeff_dxx_m1_0 = t249+t255;
-      coeff_dxx_0_0 = t248+t243+RATIONAL(-21.0,40.0);
-      coeff_dxx_p1_0 = t247+t256;
-      coeff_dxx_p2_0 = t242+RATIONAL(1.0,40.0)+t246;
-      coeff_dxx_m1_p1 = t249+t256;
-      coeff_dxx_0_p1 = t242+t248+RATIONAL(-19.0,40.0);
-      coeff_dxx_p1_p1 = t247+t255;
-      coeff_dxx_p2_p1 = RATIONAL(-1.0,40.0)+t243+t246;
-      coeff_dxx_m1_p2 = t249+t257;
-      coeff_dxx_0_p2 = RATIONAL(-17.0,40.0)+t248+t245;
-      coeff_dxx_p1_p2 = t247+t254;
-      coeff_dxx_p2_p2 = t244+RATIONAL(-3.0,40.0)+t246;
+   t248;
+      t247 = RATIONAL(-3.0,20.0)*y;
+      t259 = RATIONAL(13.0,40.0)+t247;
+      t244 = RATIONAL(-1.0,20.0)*y;
+      t258 = RATIONAL(11.0,40.0)+t244;
+      t245 = RATIONAL(1.0,20.0)*y;
+      t257 = RATIONAL(9.0,40.0)+t245;
+      t246 = RATIONAL(3.0,20.0)*y;
+      t256 = RATIONAL(7.0,40.0)+t246;
+      t251 = RATIONAL(-1.0,4.0)*x;
+      t250 = RATIONAL(-3.0,4.0)*x;
+      t249 = RATIONAL(3.0,4.0)*x;
+      t248 = RATIONAL(1.0,4.0)*x;
+      coeff_dxx_m1_m1 = t251+t259;
+      coeff_dxx_0_m1 = RATIONAL(-23.0,40.0)+t246+t249;
+      coeff_dxx_p1_m1 = t250+t256;
+      coeff_dxx_p2_m1 = t248+RATIONAL(3.0,40.0)+t247;
+      coeff_dxx_m1_0 = t251+t258;
+      coeff_dxx_0_0 = RATIONAL(-21.0,40.0)+t245+t249;
+      coeff_dxx_p1_0 = t250+t257;
+      coeff_dxx_p2_0 = t248+RATIONAL(1.0,40.0)+t244;
+      coeff_dxx_m1_p1 = t251+t257;
+      coeff_dxx_0_p1 = t244+t249+RATIONAL(-19.0,40.0);
+      coeff_dxx_p1_p1 = t250+t258;
+      coeff_dxx_p2_p1 = RATIONAL(-1.0,40.0)+t248+t245;
+      coeff_dxx_m1_p2 = t251+t256;
+      coeff_dxx_0_p2 = RATIONAL(-17.0,40.0)+t249+t247;
+      coeff_dxx_p1_p2 = t250+t259;
+      coeff_dxx_p2_p2 = t248+t246+RATIONAL(-3.0,40.0);
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c
index f0c2888..f68960f 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c
@@ -1,4 +1,6 @@
-fp t275,
+fp t277,
+   t276,
+   t275,
    t274,
    t273,
    t272,
@@ -13,40 +15,38 @@ fp t275,
    t263,
    t262,
    t261,
-   t260,
-   t259,
-   t258;
-      t275 = RATIONAL(-1.0,20.0);
-      t274 = RATIONAL(1.0,20.0);
-      t273 = RATIONAL(-1.0,50.0);
-      t272 = RATIONAL(-3.0,20.0);
-      t271 = RATIONAL(3.0,20.0);
-      t270 = RATIONAL(2.0,25.0);
-      t269 = RATIONAL(-9.0,100.0);
-      t268 = RATIONAL(7.0,100.0);
-      t267 = RATIONAL(-1.0,100.0);
-      t266 = RATIONAL(-13.0,100.0);
-      t265 = t272*y;
-      t264 = t271*y;
-      t263 = t274*y;
-      t262 = t275*y;
-      t261 = t271*x;
-      t260 = t272*x;
-      t259 = t274*x;
-      t258 = t275*x;
-      coeff_dxy_m1_m1 = t260+t265+RATIONAL(6.0,25.0);
-      coeff_dxy_0_m1 = t273+t262+t261;
-      coeff_dxy_p1_m1 = t266+t261+t263;
-      coeff_dxy_p2_m1 = t260+t269+t264;
-      coeff_dxy_m1_0 = t264+t273+t258;
-      coeff_dxy_0_0 = t263+RATIONAL(-1.0,25.0)+t259;
-      coeff_dxy_p1_0 = t262+t267+t259;
-      coeff_dxy_p2_0 = t258+t265+t268;
-      coeff_dxy_m1_p1 = t264+t259+t266;
-      coeff_dxy_0_p1 = t258+t263+t267;
-      coeff_dxy_p1_p1 = t262+t258+RATIONAL(3.0,50.0);
-      coeff_dxy_p2_p1 = t265+t259+t270;
-      coeff_dxy_m1_p2 = t265+t261+t269;
-      coeff_dxy_0_p2 = t268+t262+t260;
-      coeff_dxy_p1_p2 = t263+t270+t260;
-      coeff_dxy_p2_p2 = RATIONAL(-3.0,50.0)+t261+t264;
+   t260;
+      t277 = RATIONAL(-1.0,20.0);
+      t276 = RATIONAL(1.0,20.0);
+      t275 = RATIONAL(-1.0,50.0);
+      t274 = RATIONAL(2.0,25.0);
+      t273 = RATIONAL(-3.0,20.0);
+      t272 = RATIONAL(3.0,20.0);
+      t271 = RATIONAL(-9.0,100.0);
+      t270 = RATIONAL(-1.0,100.0);
+      t269 = RATIONAL(7.0,100.0);
+      t268 = RATIONAL(-13.0,100.0);
+      t267 = t273*y;
+      t266 = t272*y;
+      t265 = t276*y;
+      t264 = t277*y;
+      t263 = t273*x;
+      t262 = t272*x;
+      t261 = t277*x;
+      t260 = t276*x;
+      coeff_dxy_m1_m1 = t267+t263+RATIONAL(6.0,25.0);
+      coeff_dxy_0_m1 = t275+t264+t262;
+      coeff_dxy_p1_m1 = t262+t265+t268;
+      coeff_dxy_p2_m1 = t263+t266+t271;
+      coeff_dxy_m1_0 = t261+t275+t266;
+      coeff_dxy_0_0 = t265+RATIONAL(-1.0,25.0)+t260;
+      coeff_dxy_p1_0 = t260+t270+t264;
+      coeff_dxy_p2_0 = t267+t261+t269;
+      coeff_dxy_m1_p1 = t260+t268+t266;
+      coeff_dxy_0_p1 = t270+t261+t265;
+      coeff_dxy_p1_p1 = t261+RATIONAL(3.0,50.0)+t264;
+      coeff_dxy_p2_p1 = t274+t267+t260;
+      coeff_dxy_m1_p2 = t262+t267+t271;
+      coeff_dxy_0_p2 = t264+t263+t269;
+      coeff_dxy_p1_p2 = t274+t263+t265;
+      coeff_dxy_p2_p2 = t262+t266+RATIONAL(-3.0,50.0);
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c
index 017eadf..c063405 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c
@@ -1,95 +1,95 @@
-fp t241,
-   t227,
-   t224,
-   t213,
+fp t243,
+   t204,
+   t242,
    t228,
+   t230,
+   t215,
+   t229,
    t221,
+   t241,
+   t203,
    t240,
-   t202,
+   t226,
+   t213,
    t239,
-   t223,
-   t215,
-   t222,
+   t206,
    t238,
-   t201,
+   t205,
    t237,
    t225,
-   t216,
+   t217,
+   t223,
    t236,
-   t204,
+   t224,
    t235,
-   t203,
    t234,
-   t226,
-   t214,
+   t227,
+   t218,
    t233,
-   t219,
+   t222,
    t232,
    t231,
-   t220,
-   t230,
-   t229,
-   t217,
+   t216,
+   t214,
    t212,
+   t211,
    t210,
-   t209,
-   t206,
-   t205;
-      t241 = x*y;
-      t227 = x*x;
-      t224 = RATIONAL(3.0,40.0);
-      t213 = t224*t227;
-      t228 = y*y;
-      t221 = RATIONAL(-1.0,8.0)*t228;
-      t240 = t213+t221;
-      t202 = RATIONAL(3.0,20.0)*t241;
-      t239 = t202+RATIONAL(7.0,40.0)*y;
-      t223 = RATIONAL(1.0,40.0);
-      t215 = t223*t227;
-      t222 = RATIONAL(-3.0,8.0)*t228;
-      t238 = t215+t222;
-      t201 = RATIONAL(1.0,20.0)*t241;
-      t237 = t201+RATIONAL(9.0,40.0)*y;
-      t225 = RATIONAL(-1.0,40.0);
-      t216 = t225*t227;
-      t236 = t216+t222;
-      t204 = RATIONAL(-3.0,20.0)*t241;
-      t235 = t204+RATIONAL(13.0,40.0)*y;
-      t203 = RATIONAL(-1.0,20.0)*t241;
-      t234 = t203+RATIONAL(11.0,40.0)*y;
-      t226 = RATIONAL(-3.0,40.0);
-      t214 = t226*t227;
-      t233 = t214+t221;
-      t219 = RATIONAL(3.0,8.0)*t228;
-      t232 = t215+t219;
-      t231 = t216+t219;
-      t220 = RATIONAL(1.0,8.0)*t228;
-      t230 = t213+t220;
-      t229 = t214+t220;
-      t217 = RATIONAL(2.0,25.0)*x;
-      t212 = RATIONAL(-1.0,50.0)*x;
-      t210 = RATIONAL(7.0,100.0)*x;
-      t209 = RATIONAL(-9.0,100.0)*x;
-      t206 = RATIONAL(-1.0,100.0)*x;
-      t205 = RATIONAL(-13.0,100.0)*x;
-      coeff_dy_m1_m1 = RATIONAL(6.0,25.0)*x+RATIONAL(-109.0,1200.0)+t233+t235;
-      coeff_dy_0_m1 = RATIONAL(-223.0,1200.0)+t212+t234+t240;
-      coeff_dy_p1_m1 = t205+RATIONAL(-157.0,1200.0)+t237+t240;
-      coeff_dy_p2_m1 = t209+RATIONAL(89.0,1200.0)+t233+t239;
-      coeff_dy_m1_0 = t202+RATIONAL(-23.0,40.0)*y+RATIONAL(-31.0,400.0)+t212+
+   t207;
+      t243 = x*y;
+      t204 = RATIONAL(-3.0,20.0)*t243;
+      t242 = t204+RATIONAL(13.0,40.0)*y;
+      t228 = RATIONAL(-3.0,40.0);
+      t230 = x*x;
+      t215 = t228*t230;
+      t229 = y*y;
+      t221 = RATIONAL(-1.0,8.0)*t229;
+      t241 = t215+t221;
+      t203 = RATIONAL(-1.0,20.0)*t243;
+      t240 = t203+RATIONAL(11.0,40.0)*y;
+      t226 = RATIONAL(3.0,40.0);
+      t213 = t226*t230;
+      t239 = t213+t221;
+      t206 = RATIONAL(1.0,20.0)*t243;
+      t238 = t206+RATIONAL(9.0,40.0)*y;
+      t205 = RATIONAL(3.0,20.0)*t243;
+      t237 = t205+RATIONAL(7.0,40.0)*y;
+      t225 = RATIONAL(1.0,40.0);
+      t217 = t225*t230;
+      t223 = RATIONAL(-3.0,8.0)*t229;
+      t236 = t217+t223;
+      t224 = RATIONAL(1.0,8.0)*t229;
+      t235 = t215+t224;
+      t234 = t213+t224;
+      t227 = RATIONAL(-1.0,40.0);
+      t218 = t227*t230;
+      t233 = t218+t223;
+      t222 = RATIONAL(3.0,8.0)*t229;
+      t232 = t217+t222;
+      t231 = t218+t222;
+      t216 = RATIONAL(2.0,25.0)*x;
+      t214 = RATIONAL(-1.0,50.0)*x;
+      t212 = RATIONAL(7.0,100.0)*x;
+      t211 = RATIONAL(-9.0,100.0)*x;
+      t210 = RATIONAL(-1.0,100.0)*x;
+      t207 = RATIONAL(-13.0,100.0)*x;
+      coeff_dy_m1_m1 = RATIONAL(-109.0,1200.0)+RATIONAL(6.0,25.0)*x+t241+t242;
+      coeff_dy_0_m1 = RATIONAL(-223.0,1200.0)+t214+t239+t240;
+      coeff_dy_p1_m1 = t207+RATIONAL(-157.0,1200.0)+t238+t239;
+      coeff_dy_p2_m1 = RATIONAL(89.0,1200.0)+t211+t237+t241;
+      coeff_dy_m1_0 = t205+RATIONAL(-31.0,400.0)+t214+RATIONAL(-23.0,40.0)*y+
 t231;
-      coeff_dy_0_0 = t201+RATIONAL(-1.0,25.0)*x+RATIONAL(-57.0,400.0)+RATIONAL(
--21.0,40.0)*y+t232;
-      coeff_dy_p1_0 = t206+RATIONAL(-19.0,40.0)*y+RATIONAL(-63.0,400.0)+t203+
+      coeff_dy_0_0 = t206+RATIONAL(-57.0,400.0)+RATIONAL(-21.0,40.0)*y+RATIONAL
+(-1.0,25.0)*x+t232;
+      coeff_dy_p1_0 = RATIONAL(-19.0,40.0)*y+t203+RATIONAL(-63.0,400.0)+t210+
 t232;
-      coeff_dy_p2_0 = t204+RATIONAL(-49.0,400.0)+RATIONAL(-17.0,40.0)*y+t210+
+      coeff_dy_p2_0 = t212+RATIONAL(-17.0,40.0)*y+RATIONAL(-49.0,400.0)+t204+
 t231;
-      coeff_dy_m1_p1 = t205+RATIONAL(111.0,400.0)+t238+t239;
-      coeff_dy_0_p1 = t206+RATIONAL(117.0,400.0)+t236+t237;
-      coeff_dy_p1_p1 = RATIONAL(3.0,50.0)*x+RATIONAL(103.0,400.0)+t234+t236;
-      coeff_dy_p2_p1 = t217+RATIONAL(69.0,400.0)+t235+t238;
-      coeff_dy_m1_p2 = RATIONAL(-131.0,1200.0)+t224*y+t209+t204+t230;
-      coeff_dy_0_p2 = t223*y+t210+t203+RATIONAL(43.0,1200.0)+t229;
-      coeff_dy_p1_p2 = RATIONAL(37.0,1200.0)+t201+t225*y+t217+t229;
-      coeff_dy_p2_p2 = t226*y+RATIONAL(-3.0,50.0)*x+t202+RATIONAL(-149.0,1200.0
-)+t230;
+      coeff_dy_m1_p1 = t207+RATIONAL(111.0,400.0)+t236+t237;
+      coeff_dy_0_p1 = t210+RATIONAL(117.0,400.0)+t233+t238;
+      coeff_dy_p1_p1 = RATIONAL(3.0,50.0)*x+RATIONAL(103.0,400.0)+t233+t240;
+      coeff_dy_p2_p1 = t216+RATIONAL(69.0,400.0)+t236+t242;
+      coeff_dy_m1_p2 = t226*y+t211+RATIONAL(-131.0,1200.0)+t204+t234;
+      coeff_dy_0_p2 = t203+t225*y+t212+RATIONAL(43.0,1200.0)+t235;
+      coeff_dy_p1_p2 = t206+t216+RATIONAL(37.0,1200.0)+t227*y+t235;
+      coeff_dy_p2_p2 = t228*y+RATIONAL(-149.0,1200.0)+RATIONAL(-3.0,50.0)*x+
+t205+t234;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c
index 822ddd4..ac7c0fa 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c
@@ -1,40 +1,40 @@
-fp t279,
+fp t278,
+   t293,
+   t279,
+   t292,
+   t281,
    t291,
-   t276,
+   t280,
    t290,
-   t277,
-   t289,
-   t278,
-   t288,
+   t285,
+   t284,
    t283,
-   t282,
-   t281,
-   t280;
-      t279 = RATIONAL(-3.0,20.0)*x;
-      t291 = t279+RATIONAL(13.0,40.0);
-      t276 = RATIONAL(-1.0,20.0)*x;
-      t290 = RATIONAL(11.0,40.0)+t276;
-      t277 = RATIONAL(1.0,20.0)*x;
-      t289 = RATIONAL(9.0,40.0)+t277;
+   t282;
       t278 = RATIONAL(3.0,20.0)*x;
-      t288 = RATIONAL(7.0,40.0)+t278;
-      t283 = RATIONAL(-1.0,4.0)*y;
-      t282 = RATIONAL(3.0,4.0)*y;
-      t281 = RATIONAL(-3.0,4.0)*y;
-      t280 = RATIONAL(1.0,4.0)*y;
-      coeff_dyy_m1_m1 = t283+t291;
-      coeff_dyy_0_m1 = t283+t290;
-      coeff_dyy_p1_m1 = t283+t289;
-      coeff_dyy_p2_m1 = t283+t288;
-      coeff_dyy_m1_0 = RATIONAL(-23.0,40.0)+t282+t278;
-      coeff_dyy_0_0 = RATIONAL(-21.0,40.0)+t277+t282;
-      coeff_dyy_p1_0 = t276+t282+RATIONAL(-19.0,40.0);
-      coeff_dyy_p2_0 = t282+t279+RATIONAL(-17.0,40.0);
-      coeff_dyy_m1_p1 = t281+t288;
-      coeff_dyy_0_p1 = t281+t289;
-      coeff_dyy_p1_p1 = t281+t290;
-      coeff_dyy_p2_p1 = t281+t291;
-      coeff_dyy_m1_p2 = t280+t279+RATIONAL(3.0,40.0);
-      coeff_dyy_0_p2 = t276+t280+RATIONAL(1.0,40.0);
-      coeff_dyy_p1_p2 = t280+t277+RATIONAL(-1.0,40.0);
-      coeff_dyy_p2_p2 = RATIONAL(-3.0,40.0)+t280+t278;
+      t293 = t278+RATIONAL(7.0,40.0);
+      t279 = RATIONAL(1.0,20.0)*x;
+      t292 = t279+RATIONAL(9.0,40.0);
+      t281 = RATIONAL(-3.0,20.0)*x;
+      t291 = t281+RATIONAL(13.0,40.0);
+      t280 = RATIONAL(-1.0,20.0)*x;
+      t290 = t280+RATIONAL(11.0,40.0);
+      t285 = RATIONAL(3.0,4.0)*y;
+      t284 = RATIONAL(-3.0,4.0)*y;
+      t283 = RATIONAL(1.0,4.0)*y;
+      t282 = RATIONAL(-1.0,4.0)*y;
+      coeff_dyy_m1_m1 = t282+t291;
+      coeff_dyy_0_m1 = t282+t290;
+      coeff_dyy_p1_m1 = t282+t292;
+      coeff_dyy_p2_m1 = t282+t293;
+      coeff_dyy_m1_0 = t278+t285+RATIONAL(-23.0,40.0);
+      coeff_dyy_0_0 = RATIONAL(-21.0,40.0)+t285+t279;
+      coeff_dyy_p1_0 = t280+RATIONAL(-19.0,40.0)+t285;
+      coeff_dyy_p2_0 = t285+t281+RATIONAL(-17.0,40.0);
+      coeff_dyy_m1_p1 = t284+t293;
+      coeff_dyy_0_p1 = t284+t292;
+      coeff_dyy_p1_p1 = t284+t290;
+      coeff_dyy_p2_p1 = t284+t291;
+      coeff_dyy_m1_p2 = RATIONAL(3.0,40.0)+t283+t281;
+      coeff_dyy_0_p2 = RATIONAL(1.0,40.0)+t280+t283;
+      coeff_dyy_p1_p2 = RATIONAL(-1.0,40.0)+t279+t283;
+      coeff_dyy_p2_p2 = RATIONAL(-3.0,40.0)+t278+t283;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
index e30de22..8b9827c 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
@@ -1,251 +1,239 @@
-fp t412,
-   t413,
-   t471,
-   t449,
-   t404,
-   t417,
-   t359,
-   t414,
-   t358,
-   t444,
-   t470,
-   t469,
-   t406,
-   t367,
-   t362,
-   t437,
-   t468,
+fp t418,
+   t419,
    t467,
-   t443,
+   t436,
    t466,
-   t383,
-   t409,
-   t371,
-   t446,
-   t402,
-   t310,
+   t415,
+   t410,
+   t365,
+   t414,
+   t366,
+   t440,
+   t445,
    t465,
-   t384,
-   t405,
-   t313,
+   t412,
+   t369,
+   t372,
+   t439,
    t464,
+   t442,
+   t406,
+   t316,
    t403,
-   t366,
-   t389,
-   t399,
    t463,
-   t292,
-   t397,
-   t445,
+   t409,
+   t332,
+   t408,
    t462,
-   t396,
-   t400,
+   t441,
+   t325,
+   t407,
+   t370,
+   t397,
    t461,
-   t381,
-   t382,
-   t415,
-   t355,
+   t405,
+   t368,
+   t393,
    t460,
-   t370,
-   t398,
-   t297,
-   t447,
+   t374,
    t459,
-   t411,
-   t329,
-   t401,
+   t402,
+   t297,
+   t443,
    t458,
+   t361,
+   t392,
+   t400,
    t457,
-   t407,
+   t296,
+   t401,
    t456,
-   t369,
+   t404,
+   t328,
+   t395,
    t455,
-   t360,
+   t321,
    t454,
-   t410,
-   t323,
-   t390,
+   t386,
+   t399,
    t453,
-   t385,
+   t411,
+   t331,
+   t394,
    t452,
-   t314,
+   t327,
    t451,
+   t314,
    t450,
-   t448,
-   t442,
-   t441,
    t416,
-   t356,
-   t440,
+   t357,
+   t417,
    t363,
-   t439,
+   t381,
+   t380,
+   t449,
+   t358,
+   t326,
+   t448,
+   t362,
+   t447,
+   t446,
+   t444,
    t438,
-   t436,
-   t322,
+   t437,
+   t377,
+   t413,
+   t359,
+   t390,
    t435,
-   t332,
+   t388,
+   t389,
    t434,
-   t391,
-   t393,
-   t408,
-   t364,
+   t373,
    t433,
-   t368,
    t432,
-   t374,
-   t392,
    t431,
    t430,
    t429,
-   t330,
    t428,
-   t318,
    t427,
    t426,
    t425,
    t424,
-   t423,
-   t422,
-   t395,
-   t394,
-   t387,
-   t386;
-      t412 = x*x;
-      t413 = y*y;
-      t471 = t412+t413;
-      t449 = x*y;
-      t404 = RATIONAL(-1.0,60.0);
-      t417 = y*t413;
-      t359 = t404*t417;
-      t414 = t412*x;
-      t358 = t404*t414;
-      t444 = y*t358+x*t359;
-      t470 = RATIONAL(11.0,150.0)*t449+t444;
-      t469 = RATIONAL(37.0,300.0)*t449+t444;
-      t406 = RATIONAL(1.0,60.0);
-      t367 = t406*t414;
-      t362 = t406*t417;
-      t437 = y*t367+x*t362;
-      t468 = t437+RATIONAL(-37.0,300.0)*t449;
-      t467 = RATIONAL(-11.0,150.0)*t449+t437;
-      t443 = t412*t413;
-      t466 = RATIONAL(-1.0,98.0)*t443+RATIONAL(-99.0,1225.0);
-      t383 = RATIONAL(-11.0,105.0);
-      t409 = RATIONAL(1.0,30.0);
-      t371 = t409*t414;
-      t446 = x*t413;
-      t402 = RATIONAL(-1.0,70.0);
-      t310 = t402*t446;
-      t465 = t371+t383*x+t310;
-      t384 = RATIONAL(31.0,420.0);
-      t405 = RATIONAL(-1.0,35.0);
-      t313 = t405*t446;
-      t464 = t384*x+t358+t313;
-      t403 = RATIONAL(-1.0,30.0);
-      t366 = t403*t414;
-      t389 = RATIONAL(31.0,210.0);
-      t399 = RATIONAL(-1.0,140.0);
-      t463 = t366+t389*x+t399*t446;
-      t292 = RATIONAL(21.0,200.0)*t449;
-      t397 = RATIONAL(1.0,120.0);
-      t445 = x*t417;
-      t462 = t292+t397*t445+y*t366;
-      t396 = RATIONAL(-31.0,210.0);
+   t398,
+   t396,
+   t391,
+   t387;
+      t418 = y*y;
+      t419 = x*x;
+      t467 = t418+t419;
+      t436 = t418*t419;
+      t466 = RATIONAL(-99.0,1225.0)+RATIONAL(-1.0,98.0)*t436;
+      t415 = t419*x;
+      t410 = RATIONAL(-1.0,60.0);
+      t365 = t410*t415;
+      t414 = t418*y;
+      t366 = t410*t414;
+      t440 = x*t366+y*t365;
+      t445 = x*y;
+      t465 = RATIONAL(37.0,300.0)*t445+t440;
+      t412 = RATIONAL(1.0,60.0);
+      t369 = t412*t415;
+      t372 = t412*t414;
+      t439 = x*t372+y*t369;
+      t464 = t439+RATIONAL(-37.0,300.0)*t445;
+      t442 = y*t419;
+      t406 = RATIONAL(-1.0,70.0);
+      t316 = t406*t442;
+      t403 = RATIONAL(1.0,84.0);
+      t463 = t316+t403*y+t372;
+      t409 = RATIONAL(1.0,70.0);
+      t332 = t409*t442;
+      t408 = RATIONAL(-1.0,84.0);
+      t462 = t408*y+t366+t332;
+      t441 = x*t418;
+      t325 = t409*t441;
+      t407 = RATIONAL(-1.0,30.0);
+      t370 = t407*t415;
+      t397 = RATIONAL(11.0,105.0);
+      t461 = t325+t397*x+t370;
+      t405 = RATIONAL(1.0,30.0);
+      t368 = t405*t414;
+      t393 = RATIONAL(-11.0,105.0);
+      t460 = t393*y+t316+t368;
+      t374 = t407*t414;
+      t459 = t374+t397*y+t332;
+      t402 = RATIONAL(-1.0,120.0);
+      t297 = RATIONAL(-21.0,200.0)*t445;
+      t443 = y*t415;
+      t458 = t402*t443+x*t368+t297;
+      t361 = t405*t415;
+      t392 = RATIONAL(-31.0,210.0);
       t400 = RATIONAL(1.0,140.0);
-      t461 = t400*t446+t371+t396*x;
+      t457 = t392*x+t400*t441+t361;
+      t296 = RATIONAL(21.0,200.0)*t445;
+      t401 = RATIONAL(1.0,120.0);
+      t456 = t296+x*t374+t401*t443;
+      t404 = RATIONAL(1.0,35.0);
+      t328 = t404*t442;
+      t395 = RATIONAL(-31.0,420.0);
+      t455 = t328+t372+t395*y;
+      t321 = t404*t441;
+      t454 = t395*x+t369+t321;
+      t386 = RATIONAL(31.0,210.0);
+      t399 = RATIONAL(-1.0,140.0);
+      t453 = t399*t441+t370+t386*x;
+      t411 = RATIONAL(-1.0,35.0);
+      t331 = t411*t441;
+      t394 = RATIONAL(31.0,420.0);
+      t452 = t394*x+t365+t331;
+      t327 = t411*t442;
+      t451 = t394*y+t327+t366;
+      t314 = t406*t441;
+      t450 = t314+t361+t393*x;
+      t416 = t418*t418;
+      t357 = t401*t416;
+      t417 = t419*t419;
+      t363 = t407*t417;
       t381 = RATIONAL(71.0,5880.0);
-      t382 = RATIONAL(113.0,735.0);
-      t415 = t412*t412;
-      t355 = t397*t415;
-      t460 = t381*t412+t382*t413+t355;
-      t370 = t409*t417;
-      t398 = RATIONAL(-1.0,120.0);
-      t297 = RATIONAL(-21.0,200.0)*t449;
-      t447 = y*t414;
-      t459 = t297+t398*t447+x*t370;
-      t411 = RATIONAL(1.0,70.0);
-      t329 = t411*t446;
-      t401 = RATIONAL(-1.0,84.0);
-      t458 = t329+t401*x+t358;
-      t457 = y*t371+t398*t445+t297;
-      t407 = RATIONAL(1.0,84.0);
-      t456 = t407*x+t367+t310;
-      t369 = t403*t415;
-      t455 = t369+t382*t412+t381*t413;
-      t360 = t403*t417;
-      t454 = x*t360+t292+t397*t447;
-      t410 = RATIONAL(1.0,35.0);
-      t323 = t410*t446;
-      t390 = RATIONAL(-31.0,420.0);
-      t453 = t323+t367+t390*x;
-      t385 = RATIONAL(11.0,105.0);
-      t452 = t366+t329+t385*x;
-      t314 = RATIONAL(1.0,49.0)*t443;
-      t451 = RATIONAL(51.0,1225.0)+t355+t314+t471*RATIONAL(-289.0,5880.0);
-      t450 = RATIONAL(1.0,196.0)*t443+t369+RATIONAL(-24.0,1225.0)+t471*RATIONAL
-(181.0,1470.0);
-      t448 = y*t412;
-      t442 = RATIONAL(246.0,1225.0)+RATIONAL(1.0,98.0)*t443;
-      t441 = RATIONAL(-1.0,49.0)*t443+RATIONAL(96.0,1225.0);
-      t416 = t413*t413;
-      t356 = t397*t416;
-      t440 = t356+t362;
-      t363 = t403*t416;
-      t439 = t360+t363;
-      t438 = t370+t363;
-      t436 = t356+t359;
-      t322 = t402*t448;
-      t435 = t322+t466;
-      t332 = t411*t448;
-      t434 = t332+t466;
-      t391 = RATIONAL(-53.0,196.0);
-      t393 = RATIONAL(83.0,735.0);
-      t408 = RATIONAL(1.0,20.0);
-      t364 = t408*t416;
-      t433 = t369+t393*t412+t391*t413+t364+t442;
-      t368 = t408*t415;
-      t432 = t368+t393*t413+t391*t412+t363+t442;
-      t374 = RATIONAL(191.0,5880.0);
-      t392 = RATIONAL(-41.0,196.0);
-      t431 = t374*t413+t368+t392*t412+t356+t441;
-      t430 = t392*t413+t355+t374*t412+t364+t441;
-      t429 = t400*t448+t396*y+t438+t450;
-      t330 = t410*t448;
-      t428 = t330+t390*y+t440+t451;
-      t318 = t405*t448;
-      t427 = t318+t384*y+t436+t451;
-      t426 = t407*y+t435+t440+t455;
-      t425 = t385*y+t434+t439+t460;
-      t424 = t383*y+t435+t438+t460;
-      t423 = t401*y+t434+t436+t455;
-      t422 = t389*y+t399*t448+t439+t450;
-      t395 = RATIONAL(-17.0,420.0);
-      t394 = RATIONAL(-17.0,105.0);
-      t387 = RATIONAL(17.0,105.0);
-      t386 = RATIONAL(17.0,420.0);
-      coeff_I_m2_m2 = t427+t464+t467;
-      coeff_I_m1_m2 = t423+t462+t465;
-      coeff_I_0_m2 = t395*y+t330+t359+t431;
-      coeff_I_p1_m2 = t423+t452+t457;
-      coeff_I_p2_m2 = t427+t453+t470;
-      coeff_I_m2_m1 = t424+t454+t458;
-      coeff_I_m1_m1 = t429+t461+t469;
-      coeff_I_0_m1 = t332+t394*y+t370+t432;
-      coeff_I_p1_m1 = t429+t463+t468;
-      coeff_I_p2_m1 = t424+t456+t459;
-      coeff_I_m2_0 = t395*x+t358+t323+t430;
-      coeff_I_m1_0 = t394*x+t371+t329+t433;
-      coeff_I_0_0 = t314+t368+t364+RATIONAL(541.0,1225.0)+t471*RATIONAL(-57.0,
+      t380 = RATIONAL(113.0,735.0);
+      t449 = t357+t363+t380*t419+t381*t418+t466;
+      t358 = t401*t417;
+      t326 = RATIONAL(1.0,49.0)*t436;
+      t448 = t357+t358+t326+RATIONAL(51.0,1225.0)+t467*RATIONAL(-289.0,5880.0);
+      t362 = t407*t416;
+      t447 = RATIONAL(1.0,196.0)*t436+RATIONAL(-24.0,1225.0)+t362+t363+t467*
+RATIONAL(181.0,1470.0);
+      t446 = t358+t362+t381*t419+t380*t418+t466;
+      t444 = x*t414;
+      t438 = RATIONAL(246.0,1225.0)+RATIONAL(1.0,98.0)*t436;
+      t437 = RATIONAL(-1.0,49.0)*t436+RATIONAL(96.0,1225.0);
+      t377 = RATIONAL(191.0,5880.0);
+      t413 = RATIONAL(1.0,20.0);
+      t359 = t413*t416;
+      t390 = RATIONAL(-41.0,196.0);
+      t435 = t377*t419+t358+t359+t390*t418+t437;
+      t388 = RATIONAL(-53.0,196.0);
+      t389 = RATIONAL(83.0,735.0);
+      t434 = t359+t363+t389*t419+t388*t418+t438;
+      t373 = t413*t417;
+      t433 = t388*t419+t389*t418+t362+t373+t438;
+      t432 = t390*t419+t357+t377*t418+t373+t437;
+      t431 = t392*y+t400*t442+t368+t447;
+      t430 = t374+t386*y+t399*t442+t447;
+      t429 = t296+t401*t444+y*t370+t449;
+      t428 = t403*x+t314+t369+t446;
+      t427 = t325+t365+t408*x+t446;
+      t426 = y*t361+t402*t444+t297+t449;
+      t425 = RATIONAL(11.0,150.0)*t445+t440+t448;
+      t424 = RATIONAL(-11.0,150.0)*t445+t439+t448;
+      t398 = RATIONAL(17.0,105.0);
+      t396 = RATIONAL(17.0,420.0);
+      t391 = RATIONAL(-17.0,420.0);
+      t387 = RATIONAL(-17.0,105.0);
+      coeff_I_m2_m2 = t424+t451+t452;
+      coeff_I_m1_m2 = t429+t450+t462;
+      coeff_I_0_m2 = t366+t391*y+t328+t432;
+      coeff_I_p1_m2 = t426+t461+t462;
+      coeff_I_p2_m2 = t425+t451+t454;
+      coeff_I_m2_m1 = t427+t456+t460;
+      coeff_I_m1_m1 = t431+t457+t465;
+      coeff_I_0_m1 = t332+t368+t387*y+t433;
+      coeff_I_p1_m1 = t431+t453+t464;
+      coeff_I_p2_m1 = t428+t458+t460;
+      coeff_I_m2_0 = t391*x+t321+t365+t435;
+      coeff_I_m1_0 = t325+t361+t387*x+t434;
+      coeff_I_0_0 = t326+RATIONAL(541.0,1225.0)+t373+t359+t467*RATIONAL(-57.0,
 196.0);
-      coeff_I_p1_0 = t366+t310+t387*x+t433;
-      coeff_I_p2_0 = t386*x+t367+t313+t430;
-      coeff_I_m2_p1 = t425+t458+t459;
-      coeff_I_m1_p1 = t422+t461+t468;
-      coeff_I_0_p1 = t322+t387*y+t360+t432;
-      coeff_I_p1_p1 = t422+t463+t469;
-      coeff_I_p2_p1 = t425+t454+t456;
-      coeff_I_m2_p2 = t428+t464+t470;
-      coeff_I_m1_p2 = t426+t457+t465;
-      coeff_I_0_p2 = t318+t386*y+t362+t431;
-      coeff_I_p1_p2 = t426+t452+t462;
-      coeff_I_p2_p2 = t428+t453+t467;
+      coeff_I_p1_0 = t314+t398*x+t370+t434;
+      coeff_I_p2_0 = t396*x+t369+t331+t435;
+      coeff_I_m2_p1 = t427+t458+t459;
+      coeff_I_m1_p1 = t430+t457+t464;
+      coeff_I_0_p1 = t316+t374+t398*y+t433;
+      coeff_I_p1_p1 = t430+t453+t465;
+      coeff_I_p2_p1 = t428+t456+t459;
+      coeff_I_m2_p2 = t425+t452+t455;
+      coeff_I_m1_p2 = t426+t450+t463;
+      coeff_I_0_p2 = t327+t396*y+t372+t432;
+      coeff_I_p1_p2 = t429+t461+t463;
+      coeff_I_p2_p2 = t424+t454+t455;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c
index 242cdc8..2faa346 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c
@@ -1,151 +1,155 @@
-fp t536,
-   t539,
-   t541,
-   t515,
-   t575,
-   t509,
-   t574,
-   t540,
-   t542,
-   t495,
-   t510,
-   t573,
-   t530,
-   t512,
-   t514,
+fp t537,
+   t560,
    t572,
-   t562,
-   t497,
+   t535,
+   t538,
+   t504,
+   t554,
    t571,
+   t536,
+   t526,
+   t514,
    t570,
-   t561,
+   t515,
+   t555,
    t569,
+   t494,
    t508,
    t568,
+   t559,
    t567,
+   t493,
    t566,
-   t537,
-   t518,
-   t519,
+   t486,
    t565,
-   t481,
+   t512,
    t564,
+   t530,
+   t516,
    t563,
-   t477,
+   t525,
+   t506,
+   t562,
+   t561,
    t529,
-   t511,
-   t560,
-   t538,
    t507,
-   t559,
+   t473,
    t558,
+   t511,
    t557,
-   t521,
-   t474,
    t556,
-   t475,
-   t555,
-   t554,
+   t527,
+   t510,
    t553,
+   t470,
    t552,
+   t483,
    t551,
+   t517,
    t550,
    t549,
    t548,
    t547,
    t546,
+   t490,
    t545,
-   t487,
+   t544,
+   t543,
+   t542,
+   t541,
    t484,
    t482,
-   t478,
-   t476,
-   t473;
-      t536 = RATIONAL(1.0,30.0);
-      t539 = x*x;
-      t541 = t539*x;
-      t515 = t536*t541;
-      t575 = RATIONAL(71.0,2940.0)*x+t515;
-      t509 = RATIONAL(-2.0,15.0)*t541;
-      t574 = RATIONAL(226.0,735.0)*x+t509;
-      t540 = y*y;
-      t542 = t540*y;
-      t495 = RATIONAL(-21.0,200.0)*y;
-      t510 = RATIONAL(1.0,10.0)*t539;
-      t573 = y*t510+RATIONAL(-1.0,120.0)*t542+t495;
+   t481,
+   t475,
+   t471,
+   t468;
+      t537 = y*y;
+      t560 = x*t537;
+      t572 = RATIONAL(1.0,98.0)*t560+RATIONAL(181.0,735.0)*x;
+      t535 = y*t537;
+      t538 = x*x;
+      t504 = RATIONAL(1.0,20.0)*t538;
+      t554 = RATIONAL(1.0,60.0)*t535+y*t504;
+      t571 = RATIONAL(-37.0,300.0)*y+t554;
+      t536 = x*t538;
+      t526 = RATIONAL(1.0,30.0);
+      t514 = t526*t536;
+      t570 = t514+RATIONAL(71.0,2940.0)*x;
+      t515 = RATIONAL(-1.0,20.0)*t538;
+      t555 = y*t515+RATIONAL(-1.0,60.0)*t535;
+      t569 = RATIONAL(37.0,300.0)*y+t555;
+      t494 = RATIONAL(21.0,200.0)*y;
+      t508 = RATIONAL(-1.0,10.0)*t538;
+      t568 = y*t508+t494+RATIONAL(1.0,120.0)*t535;
+      t559 = y*t538;
+      t567 = t494+RATIONAL(-1.0,30.0)*t535+RATIONAL(1.0,40.0)*t559;
+      t493 = RATIONAL(-21.0,200.0)*y;
+      t566 = t493+RATIONAL(-1.0,40.0)*t559+t526*t535;
+      t486 = RATIONAL(2.0,49.0)*t560;
+      t565 = RATIONAL(-289.0,2940.0)*x+t486+t514;
+      t512 = RATIONAL(1.0,10.0)*t538;
+      t564 = t493+RATIONAL(-1.0,120.0)*t535+y*t512;
       t530 = RATIONAL(-1.0,35.0);
-      t512 = t530*t540;
-      t514 = RATIONAL(-1.0,20.0)*t539;
-      t572 = RATIONAL(31.0,420.0)+t512+t514;
-      t562 = y*t539;
-      t497 = RATIONAL(21.0,200.0)*y;
-      t571 = RATIONAL(1.0,40.0)*t562+t497+RATIONAL(-1.0,30.0)*t542;
-      t570 = RATIONAL(-1.0,40.0)*t562+t495+t536*t542;
-      t561 = x*t540;
-      t569 = RATIONAL(1.0,98.0)*t561+RATIONAL(181.0,735.0)*x+t509;
-      t508 = RATIONAL(-1.0,10.0)*t539;
-      t568 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t540+t508;
-      t567 = RATIONAL(-31.0,210.0)+RATIONAL(1.0,140.0)*t540+t510;
-      t566 = y*t508+RATIONAL(1.0,120.0)*t542+t497;
-      t537 = RATIONAL(1.0,35.0);
-      t518 = t537*t540;
-      t519 = RATIONAL(1.0,20.0)*t539;
-      t565 = RATIONAL(-31.0,420.0)+t518+t519;
-      t481 = RATIONAL(2.0,49.0)*t561;
-      t564 = t481+RATIONAL(-289.0,2940.0)*x+t515;
-      t563 = x*y;
-      t477 = RATIONAL(-1.0,49.0)*t561;
-      t529 = RATIONAL(-1.0,70.0);
-      t511 = t529*t540;
-      t560 = t477+t511;
-      t538 = RATIONAL(1.0,70.0);
-      t507 = t538*t540;
-      t559 = t477+t507;
-      t558 = y*t514+RATIONAL(-1.0,60.0)*t542;
-      t557 = y*t519+RATIONAL(1.0,60.0)*t542;
-      t521 = RATIONAL(1.0,5.0)*t541;
-      t474 = RATIONAL(-2.0,49.0)*t561;
-      t556 = t521+t474+RATIONAL(-41.0,98.0)*x;
-      t475 = RATIONAL(1.0,49.0)*t561;
-      t555 = t521+t475+RATIONAL(-53.0,98.0)*x;
-      t554 = t475+RATIONAL(166.0,735.0)*x+t509;
-      t553 = t474+RATIONAL(191.0,2940.0)*x+t515;
-      t552 = RATIONAL(-11.0,105.0)+t510+t560+t574;
-      t551 = RATIONAL(-1.0,84.0)+t514+t559+t575;
-      t550 = RATIONAL(-37.0,300.0)*y+t557+t569;
-      t549 = RATIONAL(1.0,84.0)+t519+t560+t575;
-      t548 = RATIONAL(11.0,105.0)+t508+t559+t574;
-      t547 = RATIONAL(11.0,150.0)*y+t558+t564;
-      t546 = RATIONAL(-11.0,150.0)*y+t557+t564;
-      t545 = RATIONAL(37.0,300.0)*y+t558+t569;
-      t487 = t529*t563;
-      t484 = RATIONAL(-2.0,35.0)*t563;
-      t482 = t530*t563;
-      t478 = t538*t563;
-      t476 = RATIONAL(2.0,35.0)*t563;
-      t473 = t537*t563;
-      coeff_dx_m2_m2 = t484+t546+t572;
-      coeff_dx_m1_m2 = t473+t552+t566;
-      coeff_dx_0_m2 = t476+t556;
-      coeff_dx_p1_m2 = t473+t548+t573;
-      coeff_dx_p2_m2 = t484+t547+t565;
-      coeff_dx_m2_m1 = t482+t551+t571;
-      coeff_dx_m1_m1 = t478+t545+t567;
-      coeff_dx_0_m1 = t473+t555;
-      coeff_dx_p1_m1 = t478+t550+t568;
-      coeff_dx_p2_m1 = t482+t549+t570;
-      coeff_dx_m2_0 = t514+RATIONAL(-17.0,420.0)+t518+t553;
-      coeff_dx_m1_0 = t510+t507+RATIONAL(-17.0,105.0)+t554;
-      coeff_dx_0_0 = t481+RATIONAL(-57.0,98.0)*x+t521;
-      coeff_dx_p1_0 = RATIONAL(17.0,105.0)+t511+t508+t554;
-      coeff_dx_p2_0 = t519+t512+RATIONAL(17.0,420.0)+t553;
-      coeff_dx_m2_p1 = t473+t551+t570;
-      coeff_dx_m1_p1 = t487+t550+t567;
-      coeff_dx_0_p1 = t482+t555;
-      coeff_dx_p1_p1 = t487+t545+t568;
-      coeff_dx_p2_p1 = t473+t549+t571;
-      coeff_dx_m2_p2 = t476+t547+t572;
-      coeff_dx_m1_p2 = t482+t552+t573;
-      coeff_dx_0_p2 = t484+t556;
-      coeff_dx_p1_p2 = t482+t548+t566;
-      coeff_dx_p2_p2 = t476+t546+t565;
+      t516 = t530*t537;
+      t563 = t515+t516+RATIONAL(31.0,420.0);
+      t525 = RATIONAL(1.0,35.0);
+      t506 = t525*t537;
+      t562 = t504+t506+RATIONAL(-31.0,420.0);
+      t561 = x*y;
+      t529 = RATIONAL(1.0,70.0);
+      t507 = t529*t537;
+      t473 = RATIONAL(-1.0,49.0)*t560;
+      t558 = t507+t473;
+      t511 = RATIONAL(-2.0,15.0)*t536;
+      t557 = t508+t511;
+      t556 = t511+t512;
+      t527 = RATIONAL(-1.0,70.0);
+      t510 = t527*t537;
+      t553 = t510+t473;
+      t470 = RATIONAL(1.0,49.0)*t560;
+      t552 = RATIONAL(166.0,735.0)*x+t511+t470;
+      t483 = RATIONAL(-2.0,49.0)*t560;
+      t551 = RATIONAL(191.0,2940.0)*x+t483+t514;
+      t517 = RATIONAL(1.0,5.0)*t536;
+      t550 = RATIONAL(-41.0,98.0)*x+t483+t517;
+      t549 = RATIONAL(-53.0,98.0)*x+t517+t470;
+      t548 = t504+RATIONAL(1.0,84.0)+t553+t570;
+      t547 = RATIONAL(1.0,140.0)*t537+RATIONAL(-31.0,210.0)+t556+t572;
+      t546 = RATIONAL(-1.0,140.0)*t537+RATIONAL(31.0,210.0)+t557+t572;
+      t490 = RATIONAL(226.0,735.0)*x;
+      t545 = t490+RATIONAL(11.0,105.0)+t557+t558;
+      t544 = t515+RATIONAL(-1.0,84.0)+t558+t570;
+      t543 = RATIONAL(11.0,150.0)*y+t555+t565;
+      t542 = RATIONAL(-11.0,150.0)*y+t554+t565;
+      t541 = t490+RATIONAL(-11.0,105.0)+t553+t556;
+      t484 = RATIONAL(-2.0,35.0)*t561;
+      t482 = RATIONAL(2.0,35.0)*t561;
+      t481 = t525*t561;
+      t475 = t530*t561;
+      t471 = t529*t561;
+      t468 = t527*t561;
+      coeff_dx_m2_m2 = t484+t542+t563;
+      coeff_dx_m1_m2 = t481+t541+t568;
+      coeff_dx_0_m2 = t482+t550;
+      coeff_dx_p1_m2 = t481+t545+t564;
+      coeff_dx_p2_m2 = t484+t543+t562;
+      coeff_dx_m2_m1 = t475+t544+t567;
+      coeff_dx_m1_m1 = t471+t547+t569;
+      coeff_dx_0_m1 = t481+t549;
+      coeff_dx_p1_m1 = t471+t546+t571;
+      coeff_dx_p2_m1 = t475+t548+t566;
+      coeff_dx_m2_0 = t515+t506+RATIONAL(-17.0,420.0)+t551;
+      coeff_dx_m1_0 = t512+RATIONAL(-17.0,105.0)+t507+t552;
+      coeff_dx_0_0 = t517+RATIONAL(-57.0,98.0)*x+t486;
+      coeff_dx_p1_0 = t508+RATIONAL(17.0,105.0)+t510+t552;
+      coeff_dx_p2_0 = t504+RATIONAL(17.0,420.0)+t516+t551;
+      coeff_dx_m2_p1 = t481+t544+t566;
+      coeff_dx_m1_p1 = t468+t547+t571;
+      coeff_dx_0_p1 = t475+t549;
+      coeff_dx_p1_p1 = t468+t546+t569;
+      coeff_dx_p2_p1 = t481+t548+t567;
+      coeff_dx_m2_p2 = t482+t543+t563;
+      coeff_dx_m1_p2 = t475+t541+t564;
+      coeff_dx_0_p2 = t484+t550;
+      coeff_dx_p1_p2 = t475+t545+t568;
+      coeff_dx_p2_p2 = t482+t542+t562;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c
index 51c1f78..27d8e04 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c
@@ -1,103 +1,107 @@
-fp t717,
-   t694,
-   t716,
-   t702,
-   t734,
-   t696,
+fp t711,
+   t687,
+   t708,
    t712,
-   t695,
-   t733,
-   t732,
-   t731,
-   t730,
-   t700,
-   t703,
-   t697,
    t691,
-   t693,
    t729,
-   t701,
+   t698,
    t728,
-   t687,
+   t692,
    t727,
    t726,
-   t699,
    t725,
-   t692,
+   t693,
+   t695,
+   t697,
+   t683,
+   t689,
    t724,
-   t698,
+   t685,
+   t696,
    t723,
    t722,
    t721,
    t720,
    t719,
-   t690,
+   t694,
    t718,
-   t685,
+   t717,
+   t716,
+   t686,
+   t715,
+   t714,
+   t713,
+   t690,
+   t688,
    t684,
-   t683,
-   t680;
-      t717 = y*y;
-      t694 = RATIONAL(-1.0,49.0)*t717;
-      t716 = x*x;
-      t702 = RATIONAL(-2.0,5.0)*t716;
-      t734 = t694+t702+RATIONAL(226.0,735.0);
-      t696 = RATIONAL(2.0,49.0)*t717;
-      t712 = RATIONAL(1.0,10.0);
-      t695 = t712*t716;
-      t733 = t696+t695+RATIONAL(-289.0,2940.0);
-      t732 = RATIONAL(1.0,98.0)*t717+t702+RATIONAL(181.0,735.0);
-      t731 = t694+t695+RATIONAL(71.0,2940.0);
-      t730 = x*y;
-      t700 = RATIONAL(-1.0,5.0)*x;
-      t703 = RATIONAL(1.0,5.0)*x;
-      t697 = RATIONAL(-1.0,10.0)*x;
-      t691 = t712*x;
-      t693 = RATIONAL(-2.0,49.0)*t717;
-      t729 = t693+t695+RATIONAL(191.0,2940.0);
-      t701 = RATIONAL(3.0,5.0)*t716;
-      t728 = t693+t701+RATIONAL(-41.0,98.0);
-      t687 = RATIONAL(1.0,49.0)*t717;
-      t727 = t687+t701+RATIONAL(-53.0,98.0);
-      t726 = t687+t702+RATIONAL(166.0,735.0);
-      t699 = RATIONAL(1.0,35.0)*y;
-      t725 = t699+t734;
-      t692 = RATIONAL(2.0,35.0)*y;
-      t724 = t692+t733;
-      t698 = RATIONAL(-2.0,35.0)*y;
-      t723 = t698+t733;
-      t722 = RATIONAL(1.0,20.0)*t730+t731;
-      t721 = RATIONAL(-1.0,70.0)*y+t732;
-      t720 = RATIONAL(1.0,70.0)*y+t732;
-      t719 = RATIONAL(-1.0,20.0)*t730+t731;
-      t690 = RATIONAL(-1.0,35.0)*y;
-      t718 = t690+t734;
-      t685 = y*t703;
-      t684 = y*t700;
-      t683 = y*t697;
-      t680 = y*t691;
-      coeff_dxx_m2_m2 = t697+t680+t723;
-      coeff_dxx_m1_m2 = t684+t703+t725;
-      coeff_dxx_0_m2 = t692+t728;
-      coeff_dxx_p1_m2 = t700+t685+t725;
-      coeff_dxx_p2_m2 = t691+t683+t723;
-      coeff_dxx_m2_m1 = t697+t690+t722;
-      coeff_dxx_m1_m1 = t683+t703+t720;
-      coeff_dxx_0_m1 = t699+t727;
-      coeff_dxx_p1_m1 = t680+t700+t720;
-      coeff_dxx_p2_m1 = t691+t690+t719;
-      coeff_dxx_m2_0 = t697+t729;
-      coeff_dxx_m1_0 = t703+t726;
-      coeff_dxx_0_0 = t701+RATIONAL(-57.0,98.0)+t696;
-      coeff_dxx_p1_0 = t700+t726;
-      coeff_dxx_p2_0 = t691+t729;
-      coeff_dxx_m2_p1 = t699+t697+t719;
-      coeff_dxx_m1_p1 = t680+t703+t721;
-      coeff_dxx_0_p1 = t690+t727;
-      coeff_dxx_p1_p1 = t683+t700+t721;
-      coeff_dxx_p2_p1 = t691+t699+t722;
-      coeff_dxx_m2_p2 = t697+t683+t724;
-      coeff_dxx_m1_p2 = t685+t703+t718;
-      coeff_dxx_0_p2 = t698+t728;
-      coeff_dxx_p1_p2 = t684+t700+t718;
-      coeff_dxx_p2_p2 = t691+t680+t724;
+   t681,
+   t678,
+   t677,
+   t676,
+   t675;
+      t711 = y*y;
+      t687 = RATIONAL(-1.0,49.0)*t711;
+      t708 = RATIONAL(1.0,10.0);
+      t712 = x*x;
+      t691 = t708*t712;
+      t729 = RATIONAL(71.0,2940.0)+t687+t691;
+      t698 = RATIONAL(-2.0,5.0)*t712;
+      t728 = t698+RATIONAL(226.0,735.0)+t687;
+      t692 = RATIONAL(2.0,49.0)*t711;
+      t727 = RATIONAL(-289.0,2940.0)+t691+t692;
+      t726 = t698+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t711;
+      t725 = x*y;
+      t693 = RATIONAL(-1.0,10.0)*x;
+      t695 = RATIONAL(-1.0,5.0)*x;
+      t697 = RATIONAL(1.0,5.0)*x;
+      t683 = t708*x;
+      t689 = RATIONAL(1.0,49.0)*t711;
+      t724 = t698+RATIONAL(166.0,735.0)+t689;
+      t685 = RATIONAL(-2.0,49.0)*t711;
+      t696 = RATIONAL(3.0,5.0)*t712;
+      t723 = t685+RATIONAL(-41.0,98.0)+t696;
+      t722 = RATIONAL(-53.0,98.0)+t689+t696;
+      t721 = t685+RATIONAL(191.0,2940.0)+t691;
+      t720 = t683+t729;
+      t719 = t693+t729;
+      t694 = RATIONAL(-2.0,35.0)*y;
+      t718 = t694+t727;
+      t717 = t697+t726;
+      t716 = y*t695+t728;
+      t686 = RATIONAL(2.0,35.0)*y;
+      t715 = t686+t727;
+      t714 = y*t697+t728;
+      t713 = t695+t726;
+      t690 = RATIONAL(-1.0,70.0)*y;
+      t688 = RATIONAL(1.0,35.0)*y;
+      t684 = RATIONAL(1.0,70.0)*y;
+      t681 = RATIONAL(-1.0,35.0)*y;
+      t678 = y*t683;
+      t677 = y*t693;
+      t676 = RATIONAL(1.0,20.0)*t725;
+      t675 = RATIONAL(-1.0,20.0)*t725;
+      coeff_dxx_m2_m2 = t693+t678+t718;
+      coeff_dxx_m1_m2 = t697+t688+t716;
+      coeff_dxx_0_m2 = t686+t723;
+      coeff_dxx_p1_m2 = t695+t688+t714;
+      coeff_dxx_p2_m2 = t677+t683+t718;
+      coeff_dxx_m2_m1 = t676+t681+t719;
+      coeff_dxx_m1_m1 = t677+t684+t717;
+      coeff_dxx_0_m1 = t688+t722;
+      coeff_dxx_p1_m1 = t684+t678+t713;
+      coeff_dxx_p2_m1 = t675+t681+t720;
+      coeff_dxx_m2_0 = t693+t721;
+      coeff_dxx_m1_0 = t697+t724;
+      coeff_dxx_0_0 = RATIONAL(-57.0,98.0)+t692+t696;
+      coeff_dxx_p1_0 = t695+t724;
+      coeff_dxx_p2_0 = t683+t721;
+      coeff_dxx_m2_p1 = t675+t688+t719;
+      coeff_dxx_m1_p1 = t690+t678+t717;
+      coeff_dxx_0_p1 = t681+t722;
+      coeff_dxx_p1_p1 = t690+t677+t713;
+      coeff_dxx_p2_p1 = t688+t676+t720;
+      coeff_dxx_m2_p2 = t693+t677+t715;
+      coeff_dxx_m1_p2 = t681+t697+t714;
+      coeff_dxx_0_p2 = t694+t723;
+      coeff_dxx_p1_p2 = t695+t681+t716;
+      coeff_dxx_p2_p2 = t683+t678+t715;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c
index 3f68e3b..1494dfa 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c
@@ -1,119 +1,119 @@
-fp t783,
-   t782,
-   t801,
-   t780,
-   t753,
-   t763,
-   t800,
-   t771,
-   t744,
-   t759,
-   t799,
-   t798,
-   t797,
+fp t778,
+   t777,
    t796,
-   t738,
+   t771,
+   t775,
    t795,
    t794,
+   t770,
+   t776,
    t793,
    t792,
-   t773,
-   t777,
    t791,
    t790,
-   t739,
    t789,
-   t736,
+   t732,
    t788,
-   t774,
-   t778,
    t787,
+   t730,
    t786,
+   t731,
    t785,
+   t768,
+   t737,
+   t765,
+   t755,
    t784,
+   t745,
+   t750,
+   t783,
+   t782,
    t781,
-   t776,
-   t772,
-   t770,
-   t762,
-   t761,
-   t760,
+   t780,
+   t779,
+   t774,
+   t769,
+   t767,
+   t766,
    t758,
-   t756,
+   t754,
    t749,
-   t743,
+   t744,
    t742,
-   t737,
-   t735;
-      t783 = y*y;
-      t782 = x*x;
-      t801 = t783+t782;
-      t780 = RATIONAL(1.0,35.0);
-      t753 = t780*x;
-      t763 = t780*y;
-      t800 = t753+t763;
-      t771 = RATIONAL(-1.0,35.0);
-      t744 = t771*x;
-      t759 = t771*y;
-      t799 = t744+t759;
-      t798 = t753+t759;
-      t797 = t744+t763;
-      t796 = x*y;
-      t738 = RATIONAL(-2.0,49.0)*t796;
-      t795 = t738+RATIONAL(-21.0,200.0);
-      t794 = t801*RATIONAL(-1.0,20.0);
-      t793 = t738+RATIONAL(21.0,200.0);
-      t792 = t801*RATIONAL(1.0,20.0);
-      t773 = RATIONAL(1.0,10.0);
-      t777 = RATIONAL(-1.0,40.0);
-      t791 = t777*t782+t773*t783+t795;
-      t790 = t777*t783+t773*t782+t795;
-      t739 = RATIONAL(4.0,49.0)*t796;
-      t789 = t739+RATIONAL(11.0,150.0)+t794;
-      t736 = RATIONAL(1.0,49.0)*t796;
-      t788 = t736+RATIONAL(37.0,300.0)+t794;
-      t774 = RATIONAL(-1.0,10.0);
-      t778 = RATIONAL(1.0,40.0);
-      t787 = t774*t782+t778*t783+t793;
-      t786 = t739+RATIONAL(-11.0,150.0)+t792;
-      t785 = t736+RATIONAL(-37.0,300.0)+t792;
-      t784 = t774*t783+t778*t782+t793;
-      t781 = RATIONAL(1.0,70.0);
-      t776 = RATIONAL(-2.0,35.0);
-      t772 = RATIONAL(2.0,35.0);
-      t770 = RATIONAL(-1.0,70.0);
-      t762 = t776*y;
-      t761 = t770*x;
-      t760 = t772*y;
-      t758 = t781*y;
-      t756 = t776*x;
-      t749 = t772*x;
-      t743 = t781*x;
-      t742 = t770*y;
-      t737 = RATIONAL(-4.0,49.0)*t796;
-      t735 = RATIONAL(2.0,49.0)*t796;
-      coeff_dxy_m2_m2 = t756+t762+t786;
-      coeff_dxy_m1_m2 = t787+t798;
-      coeff_dxy_0_m2 = t737+t749;
-      coeff_dxy_p1_m2 = t790+t800;
-      coeff_dxy_p2_m2 = t760+t756+t789;
-      coeff_dxy_m2_m1 = t784+t797;
-      coeff_dxy_m1_m1 = t758+t743+t788;
-      coeff_dxy_0_m1 = t735+t753;
-      coeff_dxy_p1_m1 = t743+t742+t785;
-      coeff_dxy_p2_m1 = t791+t799;
-      coeff_dxy_m2_0 = t760+t737;
-      coeff_dxy_m1_0 = t763+t735;
-      coeff_dxy_0_0 = t739;
-      coeff_dxy_p1_0 = t759+t735;
-      coeff_dxy_p2_0 = t762+t737;
-      coeff_dxy_m2_p1 = t791+t800;
-      coeff_dxy_m1_p1 = t758+t761+t785;
-      coeff_dxy_0_p1 = t744+t735;
-      coeff_dxy_p1_p1 = t761+t742+t788;
-      coeff_dxy_p2_p1 = t784+t798;
-      coeff_dxy_m2_p2 = t749+t762+t789;
-      coeff_dxy_m1_p2 = t790+t799;
-      coeff_dxy_0_p2 = t756+t737;
-      coeff_dxy_p1_p2 = t787+t797;
-      coeff_dxy_p2_p2 = t760+t749+t786;
+   t739,
+   t738,
+   t736,
+   t734,
+   t733;
+      t778 = x*x;
+      t777 = y*y;
+      t796 = t778+t777;
+      t771 = RATIONAL(1.0,10.0);
+      t775 = RATIONAL(-1.0,40.0);
+      t795 = t771*t777+t775*t778;
+      t794 = t771*t778+t775*t777;
+      t770 = RATIONAL(-1.0,10.0);
+      t776 = RATIONAL(1.0,40.0);
+      t793 = t770*t778+t776*t777;
+      t792 = t770*t777+t776*t778;
+      t791 = x*y;
+      t790 = t796*RATIONAL(-1.0,20.0);
+      t789 = t796*RATIONAL(1.0,20.0);
+      t732 = RATIONAL(-2.0,49.0)*t791;
+      t788 = t732+RATIONAL(21.0,200.0);
+      t787 = t732+RATIONAL(-21.0,200.0);
+      t730 = RATIONAL(1.0,49.0)*t791;
+      t786 = t730+RATIONAL(-37.0,300.0)+t789;
+      t731 = RATIONAL(4.0,49.0)*t791;
+      t785 = t731+RATIONAL(-11.0,150.0)+t789;
+      t768 = RATIONAL(-1.0,35.0);
+      t737 = t768*y;
+      t765 = RATIONAL(1.0,35.0);
+      t755 = t765*x;
+      t784 = t737+t755+t788;
+      t745 = t768*x;
+      t750 = t765*y;
+      t783 = t745+t750+t788;
+      t782 = t737+t745+t787;
+      t781 = t750+t755+t787;
+      t780 = t731+RATIONAL(11.0,150.0)+t790;
+      t779 = t730+RATIONAL(37.0,300.0)+t790;
+      t774 = RATIONAL(-2.0,35.0);
+      t769 = RATIONAL(2.0,35.0);
+      t767 = RATIONAL(1.0,70.0);
+      t766 = RATIONAL(-1.0,70.0);
+      t758 = t774*y;
+      t754 = t766*y;
+      t749 = t769*x;
+      t744 = t769*y;
+      t742 = t767*x;
+      t739 = t767*y;
+      t738 = t766*x;
+      t736 = t774*x;
+      t734 = RATIONAL(2.0,49.0)*t791;
+      t733 = RATIONAL(-4.0,49.0)*t791;
+      coeff_dxy_m2_m2 = t736+t758+t785;
+      coeff_dxy_m1_m2 = t784+t793;
+      coeff_dxy_0_m2 = t733+t749;
+      coeff_dxy_p1_m2 = t781+t794;
+      coeff_dxy_p2_m2 = t736+t744+t780;
+      coeff_dxy_m2_m1 = t783+t792;
+      coeff_dxy_m1_m1 = t742+t739+t779;
+      coeff_dxy_0_m1 = t755+t734;
+      coeff_dxy_p1_m1 = t754+t742+t786;
+      coeff_dxy_p2_m1 = t782+t795;
+      coeff_dxy_m2_0 = t744+t733;
+      coeff_dxy_m1_0 = t750+t734;
+      coeff_dxy_0_0 = t731;
+      coeff_dxy_p1_0 = t737+t734;
+      coeff_dxy_p2_0 = t758+t733;
+      coeff_dxy_m2_p1 = t781+t795;
+      coeff_dxy_m1_p1 = t738+t739+t786;
+      coeff_dxy_0_p1 = t734+t745;
+      coeff_dxy_p1_p1 = t738+t754+t779;
+      coeff_dxy_p2_p1 = t784+t792;
+      coeff_dxy_m2_p2 = t749+t758+t780;
+      coeff_dxy_m1_p2 = t782+t794;
+      coeff_dxy_0_p2 = t736+t733;
+      coeff_dxy_p1_p2 = t783+t793;
+      coeff_dxy_p2_p2 = t744+t749+t785;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c
index a3866a1..ab52832 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c
@@ -1,147 +1,147 @@
-fp t644,
-   t646,
-   t623,
-   t679,
-   t643,
-   t645,
-   t618,
-   t661,
-   t678,
-   t624,
-   t662,
-   t677,
-   t640,
-   t620,
-   t676,
-   t633,
+fp t643,
+   t641,
+   t642,
    t621,
-   t675,
-   t634,
-   t616,
+   t659,
    t674,
-   t641,
-   t613,
+   t608,
+   t658,
    t673,
-   t642,
-   t615,
-   t619,
+   t632,
+   t612,
    t672,
+   t610,
    t671,
-   t666,
+   t599,
+   t609,
    t670,
-   t588,
+   t640,
+   t611,
+   t661,
    t669,
-   t622,
+   t598,
    t668,
+   t631,
+   t618,
+   t580,
    t667,
+   t577,
+   t666,
    t665,
-   t584,
    t664,
+   t634,
+   t616,
    t663,
-   t582,
-   t625,
+   t662,
    t660,
-   t586,
-   t659,
-   t658,
+   t575,
    t657,
+   t583,
+   t622,
    t656,
    t655,
    t654,
    t653,
+   t630,
+   t615,
    t652,
+   t635,
+   t619,
    t651,
    t650,
    t649,
-   t593,
-   t591,
+   t648,
+   t647,
+   t646,
    t590,
-   t583,
-   t578,
-   t577;
-      t644 = y*y;
-      t646 = y*t644;
-      t623 = RATIONAL(-2.0,15.0)*t646;
-      t679 = RATIONAL(226.0,735.0)*y+t623;
+   t589,
+   t588,
+   t584,
+   t576,
+   t573;
       t643 = x*x;
-      t645 = t643*x;
-      t618 = RATIONAL(-1.0,20.0)*t644;
-      t661 = x*t618+RATIONAL(-1.0,60.0)*t645;
-      t678 = t661+RATIONAL(37.0,300.0)*x;
-      t624 = RATIONAL(1.0,20.0)*t644;
-      t662 = x*t624+RATIONAL(1.0,60.0)*t645;
-      t677 = t662+RATIONAL(-37.0,300.0)*x;
-      t640 = RATIONAL(1.0,30.0);
-      t620 = t640*t646;
-      t676 = RATIONAL(71.0,2940.0)*y+t620;
-      t633 = RATIONAL(-1.0,70.0);
-      t621 = t633*t643;
-      t675 = RATIONAL(1.0,84.0)+t621+t624;
-      t634 = RATIONAL(-1.0,35.0);
+      t641 = x*t643;
+      t642 = y*y;
+      t621 = RATIONAL(-1.0,20.0)*t642;
+      t659 = RATIONAL(-1.0,60.0)*t641+x*t621;
+      t674 = RATIONAL(11.0,150.0)*x+t659;
+      t608 = RATIONAL(1.0,20.0)*t642;
+      t658 = x*t608+RATIONAL(1.0,60.0)*t641;
+      t673 = RATIONAL(-11.0,150.0)*x+t658;
+      t632 = RATIONAL(-1.0,70.0);
+      t612 = t632*t643;
+      t672 = t608+t612+RATIONAL(1.0,84.0);
+      t610 = RATIONAL(1.0,10.0)*t642;
+      t671 = RATIONAL(1.0,140.0)*t643+t610+RATIONAL(-31.0,210.0);
+      t599 = RATIONAL(21.0,200.0)*x;
+      t609 = RATIONAL(-1.0,10.0)*t642;
+      t670 = t599+RATIONAL(1.0,120.0)*t641+x*t609;
+      t640 = y*t642;
+      t611 = RATIONAL(-2.0,15.0)*t640;
+      t661 = y*t643;
+      t669 = RATIONAL(1.0,98.0)*t661+t611+RATIONAL(181.0,735.0)*y;
+      t598 = RATIONAL(-21.0,200.0)*x;
+      t668 = t598+x*t610+RATIONAL(-1.0,120.0)*t641;
+      t631 = RATIONAL(1.0,30.0);
+      t618 = t631*t640;
+      t580 = RATIONAL(-1.0,49.0)*t661;
+      t667 = RATIONAL(71.0,2940.0)*y+t618+t580;
+      t577 = RATIONAL(2.0,49.0)*t661;
+      t666 = RATIONAL(-289.0,2940.0)*y+t618+t577;
+      t665 = t611+RATIONAL(226.0,735.0)*y+t580;
+      t664 = RATIONAL(-1.0,140.0)*t643+t609+RATIONAL(31.0,210.0);
+      t634 = RATIONAL(1.0,70.0);
       t616 = t634*t643;
-      t674 = RATIONAL(31.0,420.0)+t616+t618;
-      t641 = RATIONAL(1.0,35.0);
-      t613 = t641*t643;
-      t673 = RATIONAL(-31.0,420.0)+t613+t624;
-      t642 = RATIONAL(1.0,70.0);
-      t615 = t642*t643;
-      t619 = RATIONAL(-1.0,10.0)*t644;
-      t672 = RATIONAL(11.0,105.0)+t615+t619;
-      t671 = RATIONAL(-1.0,84.0)+t615+t618;
-      t666 = y*t643;
-      t670 = RATIONAL(1.0,98.0)*t666+RATIONAL(181.0,735.0)*y+t623;
-      t588 = RATIONAL(2.0,49.0)*t666;
-      t669 = RATIONAL(-289.0,2940.0)*y+t588+t620;
-      t622 = RATIONAL(1.0,10.0)*t644;
-      t668 = RATIONAL(-11.0,105.0)+t621+t622;
-      t667 = x*y;
-      t665 = x*t644;
-      t584 = RATIONAL(-1.0,49.0)*t666;
-      t664 = t584+RATIONAL(-21.0,200.0)*x;
-      t663 = t584+RATIONAL(21.0,200.0)*x;
-      t582 = RATIONAL(-2.0,49.0)*t666;
-      t625 = RATIONAL(1.0,5.0)*t646;
-      t660 = t582+RATIONAL(-41.0,98.0)*y+t625;
-      t586 = RATIONAL(1.0,49.0)*t666;
-      t659 = t586+RATIONAL(-53.0,98.0)*y+t625;
-      t658 = RATIONAL(191.0,2940.0)*y+t582+t620;
-      t657 = t586+RATIONAL(166.0,735.0)*y+t623;
-      t656 = t640*t645+RATIONAL(-1.0,40.0)*t665+t664+t676;
-      t655 = RATIONAL(1.0,40.0)*t665+RATIONAL(-1.0,30.0)*t645+t663+t676;
-      t654 = x*t622+RATIONAL(-1.0,120.0)*t645+t664+t679;
-      t653 = RATIONAL(1.0,120.0)*t645+x*t619+t663+t679;
-      t652 = RATIONAL(-11.0,150.0)*x+t662+t669;
-      t651 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t643+t619+t670;
-      t650 = RATIONAL(-31.0,210.0)+RATIONAL(1.0,140.0)*t643+t622+t670;
-      t649 = RATIONAL(11.0,150.0)*x+t661+t669;
-      t593 = t634*t667;
-      t591 = RATIONAL(-2.0,35.0)*t667;
-      t590 = t642*t667;
-      t583 = t633*t667;
-      t578 = t641*t667;
-      t577 = RATIONAL(2.0,35.0)*t667;
-      coeff_dy_m2_m2 = t591+t652+t674;
-      coeff_dy_m1_m2 = t593+t655+t671;
-      coeff_dy_0_m2 = RATIONAL(-17.0,420.0)+t618+t613+t658;
-      coeff_dy_p1_m2 = t578+t656+t671;
-      coeff_dy_p2_m2 = t577+t649+t674;
-      coeff_dy_m2_m1 = t578+t653+t668;
-      coeff_dy_m1_m1 = t590+t650+t678;
-      coeff_dy_0_m1 = RATIONAL(-17.0,105.0)+t622+t615+t657;
-      coeff_dy_p1_m1 = t583+t650+t677;
-      coeff_dy_p2_m1 = t593+t654+t668;
-      coeff_dy_m2_0 = t577+t660;
-      coeff_dy_m1_0 = t578+t659;
-      coeff_dy_0_0 = RATIONAL(-57.0,98.0)*y+t625+t588;
-      coeff_dy_p1_0 = t593+t659;
-      coeff_dy_p2_0 = t591+t660;
-      coeff_dy_m2_p1 = t578+t654+t672;
-      coeff_dy_m1_p1 = t590+t651+t677;
-      coeff_dy_0_p1 = RATIONAL(17.0,105.0)+t621+t619+t657;
-      coeff_dy_p1_p1 = t583+t651+t678;
-      coeff_dy_p2_p1 = t593+t653+t672;
-      coeff_dy_m2_p2 = t591+t649+t673;
-      coeff_dy_m1_p2 = t593+t656+t675;
-      coeff_dy_0_p2 = t616+t624+RATIONAL(17.0,420.0)+t658;
-      coeff_dy_p1_p2 = t578+t655+t675;
-      coeff_dy_p2_p2 = t577+t652+t673;
+      t663 = t616+t621+RATIONAL(-1.0,84.0);
+      t662 = x*y;
+      t660 = x*t642;
+      t575 = RATIONAL(1.0,49.0)*t661;
+      t657 = t611+RATIONAL(166.0,735.0)*y+t575;
+      t583 = RATIONAL(-2.0,49.0)*t661;
+      t622 = RATIONAL(1.0,5.0)*t640;
+      t656 = t583+RATIONAL(-41.0,98.0)*y+t622;
+      t655 = t583+RATIONAL(191.0,2940.0)*y+t618;
+      t654 = RATIONAL(-53.0,98.0)*y+t622+t575;
+      t653 = t598+RATIONAL(-1.0,40.0)*t660+t631*t641+t667;
+      t630 = RATIONAL(1.0,35.0);
+      t615 = t630*t643;
+      t652 = t608+t615+RATIONAL(-31.0,420.0)+t666;
+      t635 = RATIONAL(-1.0,35.0);
+      t619 = t635*t643;
+      t651 = t619+t621+RATIONAL(31.0,420.0)+t666;
+      t650 = t610+t612+RATIONAL(-11.0,105.0)+t665;
+      t649 = t609+t616+RATIONAL(11.0,105.0)+t665;
+      t648 = RATIONAL(1.0,40.0)*t660+t599+RATIONAL(-1.0,30.0)*t641+t667;
+      t647 = RATIONAL(37.0,300.0)*x+t659+t669;
+      t646 = RATIONAL(-37.0,300.0)*x+t658+t669;
+      t590 = RATIONAL(2.0,35.0)*t662;
+      t589 = t632*t662;
+      t588 = t630*t662;
+      t584 = t635*t662;
+      t576 = t634*t662;
+      t573 = RATIONAL(-2.0,35.0)*t662;
+      coeff_dy_m2_m2 = t573+t651+t673;
+      coeff_dy_m1_m2 = t584+t648+t663;
+      coeff_dy_0_m2 = RATIONAL(-17.0,420.0)+t615+t621+t655;
+      coeff_dy_p1_m2 = t588+t653+t663;
+      coeff_dy_p2_m2 = t590+t651+t674;
+      coeff_dy_m2_m1 = t588+t650+t670;
+      coeff_dy_m1_m1 = t576+t647+t671;
+      coeff_dy_0_m1 = RATIONAL(-17.0,105.0)+t616+t610+t657;
+      coeff_dy_p1_m1 = t589+t646+t671;
+      coeff_dy_p2_m1 = t584+t650+t668;
+      coeff_dy_m2_0 = t590+t656;
+      coeff_dy_m1_0 = t588+t654;
+      coeff_dy_0_0 = t577+t622+RATIONAL(-57.0,98.0)*y;
+      coeff_dy_p1_0 = t584+t654;
+      coeff_dy_p2_0 = t573+t656;
+      coeff_dy_m2_p1 = t588+t649+t668;
+      coeff_dy_m1_p1 = t576+t646+t664;
+      coeff_dy_0_p1 = t609+RATIONAL(17.0,105.0)+t612+t657;
+      coeff_dy_p1_p1 = t589+t647+t664;
+      coeff_dy_p2_p1 = t584+t649+t670;
+      coeff_dy_m2_p2 = t573+t652+t674;
+      coeff_dy_m1_p2 = t584+t653+t672;
+      coeff_dy_0_p2 = t619+t608+RATIONAL(17.0,420.0)+t655;
+      coeff_dy_p1_p2 = t588+t648+t672;
+      coeff_dy_p2_p2 = t590+t652+t673;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c
index f7dc31e..927104b 100644
--- a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c
@@ -1,107 +1,119 @@
-fp t838,
-   t818,
-   t839,
-   t824,
-   t856,
+fp t834,
+   t804,
    t855,
-   t834,
-   t816,
+   t815,
    t854,
-   t820,
    t853,
    t852,
-   t823,
-   t808,
-   t821,
-   t825,
-   t815,
    t851,
-   t810,
+   t816,
+   t819,
+   t830,
+   t806,
+   t818,
+   t833,
+   t811,
    t850,
-   t822,
+   t817,
    t849,
    t848,
-   t802,
    t847,
-   t805,
+   t809,
    t846,
+   t820,
    t845,
+   t812,
    t844,
    t843,
    t842,
    t841,
    t840,
-   t819,
-   t817,
+   t839,
+   t838,
+   t837,
+   t836,
+   t835,
    t814,
    t813,
+   t810,
+   t808,
    t807,
-   t806,
-   t804,
-   t803;
-      t838 = x*x;
-      t818 = RATIONAL(-1.0,49.0)*t838;
-      t839 = y*y;
-      t824 = RATIONAL(-2.0,5.0)*t839;
-      t856 = RATIONAL(226.0,735.0)+t818+t824;
-      t855 = RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t838+t824;
-      t834 = RATIONAL(1.0,10.0);
-      t816 = t834*t839;
-      t854 = RATIONAL(71.0,2940.0)+t816+t818;
-      t820 = RATIONAL(2.0,49.0)*t838;
-      t853 = RATIONAL(-289.0,2940.0)+t816+t820;
-      t852 = x*y;
-      t823 = RATIONAL(1.0,5.0)*y;
-      t808 = t834*y;
-      t821 = RATIONAL(-1.0,10.0)*y;
-      t825 = RATIONAL(-1.0,5.0)*y;
-      t815 = RATIONAL(-2.0,49.0)*t838;
-      t851 = RATIONAL(191.0,2940.0)+t815+t816;
-      t810 = RATIONAL(1.0,49.0)*t838;
-      t850 = RATIONAL(166.0,735.0)+t810+t824;
-      t822 = RATIONAL(3.0,5.0)*t839;
-      t849 = RATIONAL(-53.0,98.0)+t810+t822;
-      t848 = RATIONAL(-41.0,98.0)+t815+t822;
-      t802 = x*t821;
-      t847 = t802+t853;
-      t805 = x*t808;
-      t846 = t805+t853;
-      t845 = RATIONAL(1.0,70.0)*x+t855;
-      t844 = RATIONAL(-1.0,70.0)*x+t855;
-      t843 = t823+t856;
-      t842 = t825+t856;
-      t841 = t821+t854;
-      t840 = t808+t854;
-      t819 = RATIONAL(-2.0,35.0)*x;
-      t817 = RATIONAL(-1.0,35.0)*x;
-      t814 = RATIONAL(1.0,35.0)*x;
-      t813 = RATIONAL(2.0,35.0)*x;
-      t807 = x*t825;
-      t806 = x*t823;
-      t804 = RATIONAL(1.0,20.0)*t852;
-      t803 = RATIONAL(-1.0,20.0)*t852;
-      coeff_dyy_m2_m2 = t821+t819+t846;
-      coeff_dyy_m1_m2 = t817+t804+t841;
-      coeff_dyy_0_m2 = t821+t851;
-      coeff_dyy_p1_m2 = t803+t814+t841;
-      coeff_dyy_p2_m2 = t813+t821+t847;
-      coeff_dyy_m2_m1 = t807+t814+t843;
-      coeff_dyy_m1_m1 = t802+t823+t845;
-      coeff_dyy_0_m1 = t823+t850;
-      coeff_dyy_p1_m1 = t823+t805+t844;
-      coeff_dyy_p2_m1 = t817+t806+t843;
-      coeff_dyy_m2_0 = t813+t848;
-      coeff_dyy_m1_0 = t814+t849;
-      coeff_dyy_0_0 = t820+RATIONAL(-57.0,98.0)+t822;
-      coeff_dyy_p1_0 = t817+t849;
-      coeff_dyy_p2_0 = t819+t848;
-      coeff_dyy_m2_p1 = t806+t814+t842;
-      coeff_dyy_m1_p1 = t825+t805+t845;
-      coeff_dyy_0_p1 = t825+t850;
-      coeff_dyy_p1_p1 = t802+t825+t844;
-      coeff_dyy_p2_p1 = t817+t807+t842;
-      coeff_dyy_m2_p2 = t819+t808+t847;
-      coeff_dyy_m1_p2 = t803+t817+t840;
-      coeff_dyy_0_p2 = t808+t851;
-      coeff_dyy_p1_p2 = t814+t804+t840;
-      coeff_dyy_p2_p2 = t808+t813+t846;
+   t803,
+   t802,
+   t801,
+   t800,
+   t799,
+   t798,
+   t797;
+      t834 = x*x;
+      t804 = RATIONAL(2.0,49.0)*t834;
+      t855 = RATIONAL(-289.0,2940.0)+t804;
+      t815 = RATIONAL(-1.0,49.0)*t834;
+      t854 = t815+RATIONAL(226.0,735.0);
+      t853 = t815+RATIONAL(71.0,2940.0);
+      t852 = RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t834;
+      t851 = x*y;
+      t816 = RATIONAL(-1.0,10.0)*y;
+      t819 = RATIONAL(-1.0,5.0)*y;
+      t830 = RATIONAL(1.0,10.0);
+      t806 = t830*y;
+      t818 = RATIONAL(1.0,5.0)*y;
+      t833 = y*y;
+      t811 = t830*t833;
+      t850 = t806+t811;
+      t817 = RATIONAL(-2.0,5.0)*t833;
+      t849 = t817+t819;
+      t848 = t811+t816;
+      t847 = t818+t817;
+      t809 = RATIONAL(-2.0,49.0)*t834;
+      t846 = t809+t811+RATIONAL(191.0,2940.0);
+      t820 = RATIONAL(3.0,5.0)*t833;
+      t845 = t809+t820+RATIONAL(-41.0,98.0);
+      t812 = RATIONAL(1.0,49.0)*t834;
+      t844 = t812+t820+RATIONAL(-53.0,98.0);
+      t843 = t812+t817+RATIONAL(166.0,735.0);
+      t842 = t848+t855;
+      t841 = t850+t853;
+      t840 = t849+t852;
+      t839 = t848+t853;
+      t838 = t849+t854;
+      t837 = t847+t854;
+      t836 = t847+t852;
+      t835 = t850+t855;
+      t814 = RATIONAL(-1.0,35.0)*x;
+      t813 = RATIONAL(-2.0,35.0)*x;
+      t810 = RATIONAL(-1.0,70.0)*x;
+      t808 = RATIONAL(1.0,35.0)*x;
+      t807 = RATIONAL(2.0,35.0)*x;
+      t803 = RATIONAL(1.0,70.0)*x;
+      t802 = x*t818;
+      t801 = x*t819;
+      t800 = x*t816;
+      t799 = RATIONAL(-1.0,20.0)*t851;
+      t798 = RATIONAL(1.0,20.0)*t851;
+      t797 = x*t806;
+      coeff_dyy_m2_m2 = t813+t797+t842;
+      coeff_dyy_m1_m2 = t814+t798+t839;
+      coeff_dyy_0_m2 = t816+t846;
+      coeff_dyy_p1_m2 = t808+t799+t839;
+      coeff_dyy_p2_m2 = t800+t807+t842;
+      coeff_dyy_m2_m1 = t808+t801+t837;
+      coeff_dyy_m1_m1 = t800+t803+t836;
+      coeff_dyy_0_m1 = t818+t843;
+      coeff_dyy_p1_m1 = t797+t810+t836;
+      coeff_dyy_p2_m1 = t802+t814+t837;
+      coeff_dyy_m2_0 = t807+t845;
+      coeff_dyy_m1_0 = t808+t844;
+      coeff_dyy_0_0 = t820+t804+RATIONAL(-57.0,98.0);
+      coeff_dyy_p1_0 = t814+t844;
+      coeff_dyy_p2_0 = t813+t845;
+      coeff_dyy_m2_p1 = t808+t802+t838;
+      coeff_dyy_m1_p1 = t797+t803+t840;
+      coeff_dyy_0_p1 = t819+t843;
+      coeff_dyy_p1_p1 = t800+t810+t840;
+      coeff_dyy_p2_p1 = t801+t814+t838;
+      coeff_dyy_m2_p2 = t813+t800+t835;
+      coeff_dyy_m1_p2 = t799+t814+t841;
+      coeff_dyy_0_p2 = t806+t846;
+      coeff_dyy_p1_p2 = t808+t798+t841;
+      coeff_dyy_p2_p2 = t807+t797+t835;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size2/coeff-I.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size2/coeff-I.store.c
new file mode 100644
index 0000000..6305377
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size2/coeff-I.store.c
@@ -0,0 +1,4 @@
+COEFF(0,0) = coeff_I_0_0;
+COEFF(1,0) = coeff_I_p1_0;
+COEFF(0,1) = coeff_I_0_p1;
+COEFF(1,1) = coeff_I_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size2/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size2/coeff-dx.store.c
new file mode 100644
index 0000000..858f061
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size2/coeff-dx.store.c
@@ -0,0 +1,4 @@
+COEFF(0,0) = coeff_dx_0_0;
+COEFF(1,0) = coeff_dx_p1_0;
+COEFF(0,1) = coeff_dx_0_p1;
+COEFF(1,1) = coeff_dx_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size2/coeff-dy.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size2/coeff-dy.store.c
new file mode 100644
index 0000000..749d755
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size2/coeff-dy.store.c
@@ -0,0 +1,4 @@
+COEFF(0,0) = coeff_dy_0_0;
+COEFF(1,0) = coeff_dy_p1_0;
+COEFF(0,1) = coeff_dy_0_p1;
+COEFF(1,1) = coeff_dy_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-I.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-I.store.c
new file mode 100644
index 0000000..5169a37
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-I.store.c
@@ -0,0 +1,9 @@
+COEFF(-1,-1) = coeff_I_m1_m1;
+COEFF(0,-1) = coeff_I_0_m1;
+COEFF(1,-1) = coeff_I_p1_m1;
+COEFF(-1,0) = coeff_I_m1_0;
+COEFF(0,0) = coeff_I_0_0;
+COEFF(1,0) = coeff_I_p1_0;
+COEFF(-1,1) = coeff_I_m1_p1;
+COEFF(0,1) = coeff_I_0_p1;
+COEFF(1,1) = coeff_I_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dx.store.c
new file mode 100644
index 0000000..4a0c17d
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dx.store.c
@@ -0,0 +1,9 @@
+COEFF(-1,-1) = coeff_dx_m1_m1;
+COEFF(0,-1) = coeff_dx_0_m1;
+COEFF(1,-1) = coeff_dx_p1_m1;
+COEFF(-1,0) = coeff_dx_m1_0;
+COEFF(0,0) = coeff_dx_0_0;
+COEFF(1,0) = coeff_dx_p1_0;
+COEFF(-1,1) = coeff_dx_m1_p1;
+COEFF(0,1) = coeff_dx_0_p1;
+COEFF(1,1) = coeff_dx_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxx.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxx.store.c
new file mode 100644
index 0000000..b29d025
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxx.store.c
@@ -0,0 +1,9 @@
+COEFF(-1,-1) = coeff_dxx_m1_m1;
+COEFF(0,-1) = coeff_dxx_0_m1;
+COEFF(1,-1) = coeff_dxx_p1_m1;
+COEFF(-1,0) = coeff_dxx_m1_0;
+COEFF(0,0) = coeff_dxx_0_0;
+COEFF(1,0) = coeff_dxx_p1_0;
+COEFF(-1,1) = coeff_dxx_m1_p1;
+COEFF(0,1) = coeff_dxx_0_p1;
+COEFF(1,1) = coeff_dxx_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxy.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxy.store.c
new file mode 100644
index 0000000..e10ffac
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxy.store.c
@@ -0,0 +1,9 @@
+COEFF(-1,-1) = coeff_dxy_m1_m1;
+COEFF(0,-1) = coeff_dxy_0_m1;
+COEFF(1,-1) = coeff_dxy_p1_m1;
+COEFF(-1,0) = coeff_dxy_m1_0;
+COEFF(0,0) = coeff_dxy_0_0;
+COEFF(1,0) = coeff_dxy_p1_0;
+COEFF(-1,1) = coeff_dxy_m1_p1;
+COEFF(0,1) = coeff_dxy_0_p1;
+COEFF(1,1) = coeff_dxy_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dy.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dy.store.c
new file mode 100644
index 0000000..b24b576
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dy.store.c
@@ -0,0 +1,9 @@
+COEFF(-1,-1) = coeff_dy_m1_m1;
+COEFF(0,-1) = coeff_dy_0_m1;
+COEFF(1,-1) = coeff_dy_p1_m1;
+COEFF(-1,0) = coeff_dy_m1_0;
+COEFF(0,0) = coeff_dy_0_0;
+COEFF(1,0) = coeff_dy_p1_0;
+COEFF(-1,1) = coeff_dy_m1_p1;
+COEFF(0,1) = coeff_dy_0_p1;
+COEFF(1,1) = coeff_dy_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dyy.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dyy.store.c
new file mode 100644
index 0000000..6041bdd
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dyy.store.c
@@ -0,0 +1,9 @@
+COEFF(-1,-1) = coeff_dyy_m1_m1;
+COEFF(0,-1) = coeff_dyy_0_m1;
+COEFF(1,-1) = coeff_dyy_p1_m1;
+COEFF(-1,0) = coeff_dyy_m1_0;
+COEFF(0,0) = coeff_dyy_0_0;
+COEFF(1,0) = coeff_dyy_p1_0;
+COEFF(-1,1) = coeff_dyy_m1_p1;
+COEFF(0,1) = coeff_dyy_0_p1;
+COEFF(1,1) = coeff_dyy_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-I.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-I.store.c
new file mode 100644
index 0000000..729b8ea
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-I.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = coeff_I_m1_m1;
+COEFF(0,-1) = coeff_I_0_m1;
+COEFF(1,-1) = coeff_I_p1_m1;
+COEFF(2,-1) = coeff_I_p2_m1;
+COEFF(-1,0) = coeff_I_m1_0;
+COEFF(0,0) = coeff_I_0_0;
+COEFF(1,0) = coeff_I_p1_0;
+COEFF(2,0) = coeff_I_p2_0;
+COEFF(-1,1) = coeff_I_m1_p1;
+COEFF(0,1) = coeff_I_0_p1;
+COEFF(1,1) = coeff_I_p1_p1;
+COEFF(2,1) = coeff_I_p2_p1;
+COEFF(-1,2) = coeff_I_m1_p2;
+COEFF(0,2) = coeff_I_0_p2;
+COEFF(1,2) = coeff_I_p1_p2;
+COEFF(2,2) = coeff_I_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dx.store.c
new file mode 100644
index 0000000..b7952fd
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dx.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = coeff_dx_m1_m1;
+COEFF(0,-1) = coeff_dx_0_m1;
+COEFF(1,-1) = coeff_dx_p1_m1;
+COEFF(2,-1) = coeff_dx_p2_m1;
+COEFF(-1,0) = coeff_dx_m1_0;
+COEFF(0,0) = coeff_dx_0_0;
+COEFF(1,0) = coeff_dx_p1_0;
+COEFF(2,0) = coeff_dx_p2_0;
+COEFF(-1,1) = coeff_dx_m1_p1;
+COEFF(0,1) = coeff_dx_0_p1;
+COEFF(1,1) = coeff_dx_p1_p1;
+COEFF(2,1) = coeff_dx_p2_p1;
+COEFF(-1,2) = coeff_dx_m1_p2;
+COEFF(0,2) = coeff_dx_0_p2;
+COEFF(1,2) = coeff_dx_p1_p2;
+COEFF(2,2) = coeff_dx_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dxx.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dxx.store.c
new file mode 100644
index 0000000..1225578
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dxx.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = coeff_dxx_m1_m1;
+COEFF(0,-1) = coeff_dxx_0_m1;
+COEFF(1,-1) = coeff_dxx_p1_m1;
+COEFF(2,-1) = coeff_dxx_p2_m1;
+COEFF(-1,0) = coeff_dxx_m1_0;
+COEFF(0,0) = coeff_dxx_0_0;
+COEFF(1,0) = coeff_dxx_p1_0;
+COEFF(2,0) = coeff_dxx_p2_0;
+COEFF(-1,1) = coeff_dxx_m1_p1;
+COEFF(0,1) = coeff_dxx_0_p1;
+COEFF(1,1) = coeff_dxx_p1_p1;
+COEFF(2,1) = coeff_dxx_p2_p1;
+COEFF(-1,2) = coeff_dxx_m1_p2;
+COEFF(0,2) = coeff_dxx_0_p2;
+COEFF(1,2) = coeff_dxx_p1_p2;
+COEFF(2,2) = coeff_dxx_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dxy.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dxy.store.c
new file mode 100644
index 0000000..b39219e
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dxy.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = coeff_dxy_m1_m1;
+COEFF(0,-1) = coeff_dxy_0_m1;
+COEFF(1,-1) = coeff_dxy_p1_m1;
+COEFF(2,-1) = coeff_dxy_p2_m1;
+COEFF(-1,0) = coeff_dxy_m1_0;
+COEFF(0,0) = coeff_dxy_0_0;
+COEFF(1,0) = coeff_dxy_p1_0;
+COEFF(2,0) = coeff_dxy_p2_0;
+COEFF(-1,1) = coeff_dxy_m1_p1;
+COEFF(0,1) = coeff_dxy_0_p1;
+COEFF(1,1) = coeff_dxy_p1_p1;
+COEFF(2,1) = coeff_dxy_p2_p1;
+COEFF(-1,2) = coeff_dxy_m1_p2;
+COEFF(0,2) = coeff_dxy_0_p2;
+COEFF(1,2) = coeff_dxy_p1_p2;
+COEFF(2,2) = coeff_dxy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dy.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dy.store.c
new file mode 100644
index 0000000..5528003
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dy.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = coeff_dy_m1_m1;
+COEFF(0,-1) = coeff_dy_0_m1;
+COEFF(1,-1) = coeff_dy_p1_m1;
+COEFF(2,-1) = coeff_dy_p2_m1;
+COEFF(-1,0) = coeff_dy_m1_0;
+COEFF(0,0) = coeff_dy_0_0;
+COEFF(1,0) = coeff_dy_p1_0;
+COEFF(2,0) = coeff_dy_p2_0;
+COEFF(-1,1) = coeff_dy_m1_p1;
+COEFF(0,1) = coeff_dy_0_p1;
+COEFF(1,1) = coeff_dy_p1_p1;
+COEFF(2,1) = coeff_dy_p2_p1;
+COEFF(-1,2) = coeff_dy_m1_p2;
+COEFF(0,2) = coeff_dy_0_p2;
+COEFF(1,2) = coeff_dy_p1_p2;
+COEFF(2,2) = coeff_dy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dyy.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dyy.store.c
new file mode 100644
index 0000000..0d6fdf3
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size4/coeff-dyy.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = coeff_dyy_m1_m1;
+COEFF(0,-1) = coeff_dyy_0_m1;
+COEFF(1,-1) = coeff_dyy_p1_m1;
+COEFF(2,-1) = coeff_dyy_p2_m1;
+COEFF(-1,0) = coeff_dyy_m1_0;
+COEFF(0,0) = coeff_dyy_0_0;
+COEFF(1,0) = coeff_dyy_p1_0;
+COEFF(2,0) = coeff_dyy_p2_0;
+COEFF(-1,1) = coeff_dyy_m1_p1;
+COEFF(0,1) = coeff_dyy_0_p1;
+COEFF(1,1) = coeff_dyy_p1_p1;
+COEFF(2,1) = coeff_dyy_p2_p1;
+COEFF(-1,2) = coeff_dyy_m1_p2;
+COEFF(0,2) = coeff_dyy_0_p2;
+COEFF(1,2) = coeff_dyy_p1_p2;
+COEFF(2,2) = coeff_dyy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-I.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-I.store.c
new file mode 100644
index 0000000..e57f3f8
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-I.store.c
@@ -0,0 +1,25 @@
+COEFF(-2,-2) = coeff_I_m2_m2;
+COEFF(-1,-2) = coeff_I_m1_m2;
+COEFF(0,-2) = coeff_I_0_m2;
+COEFF(1,-2) = coeff_I_p1_m2;
+COEFF(2,-2) = coeff_I_p2_m2;
+COEFF(-2,-1) = coeff_I_m2_m1;
+COEFF(-1,-1) = coeff_I_m1_m1;
+COEFF(0,-1) = coeff_I_0_m1;
+COEFF(1,-1) = coeff_I_p1_m1;
+COEFF(2,-1) = coeff_I_p2_m1;
+COEFF(-2,0) = coeff_I_m2_0;
+COEFF(-1,0) = coeff_I_m1_0;
+COEFF(0,0) = coeff_I_0_0;
+COEFF(1,0) = coeff_I_p1_0;
+COEFF(2,0) = coeff_I_p2_0;
+COEFF(-2,1) = coeff_I_m2_p1;
+COEFF(-1,1) = coeff_I_m1_p1;
+COEFF(0,1) = coeff_I_0_p1;
+COEFF(1,1) = coeff_I_p1_p1;
+COEFF(2,1) = coeff_I_p2_p1;
+COEFF(-2,2) = coeff_I_m2_p2;
+COEFF(-1,2) = coeff_I_m1_p2;
+COEFF(0,2) = coeff_I_0_p2;
+COEFF(1,2) = coeff_I_p1_p2;
+COEFF(2,2) = coeff_I_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dx.store.c
new file mode 100644
index 0000000..f9f7aa9
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dx.store.c
@@ -0,0 +1,25 @@
+COEFF(-2,-2) = coeff_dx_m2_m2;
+COEFF(-1,-2) = coeff_dx_m1_m2;
+COEFF(0,-2) = coeff_dx_0_m2;
+COEFF(1,-2) = coeff_dx_p1_m2;
+COEFF(2,-2) = coeff_dx_p2_m2;
+COEFF(-2,-1) = coeff_dx_m2_m1;
+COEFF(-1,-1) = coeff_dx_m1_m1;
+COEFF(0,-1) = coeff_dx_0_m1;
+COEFF(1,-1) = coeff_dx_p1_m1;
+COEFF(2,-1) = coeff_dx_p2_m1;
+COEFF(-2,0) = coeff_dx_m2_0;
+COEFF(-1,0) = coeff_dx_m1_0;
+COEFF(0,0) = coeff_dx_0_0;
+COEFF(1,0) = coeff_dx_p1_0;
+COEFF(2,0) = coeff_dx_p2_0;
+COEFF(-2,1) = coeff_dx_m2_p1;
+COEFF(-1,1) = coeff_dx_m1_p1;
+COEFF(0,1) = coeff_dx_0_p1;
+COEFF(1,1) = coeff_dx_p1_p1;
+COEFF(2,1) = coeff_dx_p2_p1;
+COEFF(-2,2) = coeff_dx_m2_p2;
+COEFF(-1,2) = coeff_dx_m1_p2;
+COEFF(0,2) = coeff_dx_0_p2;
+COEFF(1,2) = coeff_dx_p1_p2;
+COEFF(2,2) = coeff_dx_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dxx.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dxx.store.c
new file mode 100644
index 0000000..1ee29cf
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dxx.store.c
@@ -0,0 +1,25 @@
+COEFF(-2,-2) = coeff_dxx_m2_m2;
+COEFF(-1,-2) = coeff_dxx_m1_m2;
+COEFF(0,-2) = coeff_dxx_0_m2;
+COEFF(1,-2) = coeff_dxx_p1_m2;
+COEFF(2,-2) = coeff_dxx_p2_m2;
+COEFF(-2,-1) = coeff_dxx_m2_m1;
+COEFF(-1,-1) = coeff_dxx_m1_m1;
+COEFF(0,-1) = coeff_dxx_0_m1;
+COEFF(1,-1) = coeff_dxx_p1_m1;
+COEFF(2,-1) = coeff_dxx_p2_m1;
+COEFF(-2,0) = coeff_dxx_m2_0;
+COEFF(-1,0) = coeff_dxx_m1_0;
+COEFF(0,0) = coeff_dxx_0_0;
+COEFF(1,0) = coeff_dxx_p1_0;
+COEFF(2,0) = coeff_dxx_p2_0;
+COEFF(-2,1) = coeff_dxx_m2_p1;
+COEFF(-1,1) = coeff_dxx_m1_p1;
+COEFF(0,1) = coeff_dxx_0_p1;
+COEFF(1,1) = coeff_dxx_p1_p1;
+COEFF(2,1) = coeff_dxx_p2_p1;
+COEFF(-2,2) = coeff_dxx_m2_p2;
+COEFF(-1,2) = coeff_dxx_m1_p2;
+COEFF(0,2) = coeff_dxx_0_p2;
+COEFF(1,2) = coeff_dxx_p1_p2;
+COEFF(2,2) = coeff_dxx_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dxy.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dxy.store.c
new file mode 100644
index 0000000..e660583
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dxy.store.c
@@ -0,0 +1,25 @@
+COEFF(-2,-2) = coeff_dxy_m2_m2;
+COEFF(-1,-2) = coeff_dxy_m1_m2;
+COEFF(0,-2) = coeff_dxy_0_m2;
+COEFF(1,-2) = coeff_dxy_p1_m2;
+COEFF(2,-2) = coeff_dxy_p2_m2;
+COEFF(-2,-1) = coeff_dxy_m2_m1;
+COEFF(-1,-1) = coeff_dxy_m1_m1;
+COEFF(0,-1) = coeff_dxy_0_m1;
+COEFF(1,-1) = coeff_dxy_p1_m1;
+COEFF(2,-1) = coeff_dxy_p2_m1;
+COEFF(-2,0) = coeff_dxy_m2_0;
+COEFF(-1,0) = coeff_dxy_m1_0;
+COEFF(0,0) = coeff_dxy_0_0;
+COEFF(1,0) = coeff_dxy_p1_0;
+COEFF(2,0) = coeff_dxy_p2_0;
+COEFF(-2,1) = coeff_dxy_m2_p1;
+COEFF(-1,1) = coeff_dxy_m1_p1;
+COEFF(0,1) = coeff_dxy_0_p1;
+COEFF(1,1) = coeff_dxy_p1_p1;
+COEFF(2,1) = coeff_dxy_p2_p1;
+COEFF(-2,2) = coeff_dxy_m2_p2;
+COEFF(-1,2) = coeff_dxy_m1_p2;
+COEFF(0,2) = coeff_dxy_0_p2;
+COEFF(1,2) = coeff_dxy_p1_p2;
+COEFF(2,2) = coeff_dxy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dy.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dy.store.c
new file mode 100644
index 0000000..23682af
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dy.store.c
@@ -0,0 +1,25 @@
+COEFF(-2,-2) = coeff_dy_m2_m2;
+COEFF(-1,-2) = coeff_dy_m1_m2;
+COEFF(0,-2) = coeff_dy_0_m2;
+COEFF(1,-2) = coeff_dy_p1_m2;
+COEFF(2,-2) = coeff_dy_p2_m2;
+COEFF(-2,-1) = coeff_dy_m2_m1;
+COEFF(-1,-1) = coeff_dy_m1_m1;
+COEFF(0,-1) = coeff_dy_0_m1;
+COEFF(1,-1) = coeff_dy_p1_m1;
+COEFF(2,-1) = coeff_dy_p2_m1;
+COEFF(-2,0) = coeff_dy_m2_0;
+COEFF(-1,0) = coeff_dy_m1_0;
+COEFF(0,0) = coeff_dy_0_0;
+COEFF(1,0) = coeff_dy_p1_0;
+COEFF(2,0) = coeff_dy_p2_0;
+COEFF(-2,1) = coeff_dy_m2_p1;
+COEFF(-1,1) = coeff_dy_m1_p1;
+COEFF(0,1) = coeff_dy_0_p1;
+COEFF(1,1) = coeff_dy_p1_p1;
+COEFF(2,1) = coeff_dy_p2_p1;
+COEFF(-2,2) = coeff_dy_m2_p2;
+COEFF(-1,2) = coeff_dy_m1_p2;
+COEFF(0,2) = coeff_dy_0_p2;
+COEFF(1,2) = coeff_dy_p1_p2;
+COEFF(2,2) = coeff_dy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dyy.store.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dyy.store.c
new file mode 100644
index 0000000..9c75614
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size5/coeff-dyy.store.c
@@ -0,0 +1,25 @@
+COEFF(-2,-2) = coeff_dyy_m2_m2;
+COEFF(-1,-2) = coeff_dyy_m1_m2;
+COEFF(0,-2) = coeff_dyy_0_m2;
+COEFF(1,-2) = coeff_dyy_p1_m2;
+COEFF(2,-2) = coeff_dyy_p2_m2;
+COEFF(-2,-1) = coeff_dyy_m2_m1;
+COEFF(-1,-1) = coeff_dyy_m1_m1;
+COEFF(0,-1) = coeff_dyy_0_m1;
+COEFF(1,-1) = coeff_dyy_p1_m1;
+COEFF(2,-1) = coeff_dyy_p2_m1;
+COEFF(-2,0) = coeff_dyy_m2_0;
+COEFF(-1,0) = coeff_dyy_m1_0;
+COEFF(0,0) = coeff_dyy_0_0;
+COEFF(1,0) = coeff_dyy_p1_0;
+COEFF(2,0) = coeff_dyy_p2_0;
+COEFF(-2,1) = coeff_dyy_m2_p1;
+COEFF(-1,1) = coeff_dyy_m1_p1;
+COEFF(0,1) = coeff_dyy_0_p1;
+COEFF(1,1) = coeff_dyy_p1_p1;
+COEFF(2,1) = coeff_dyy_p2_p1;
+COEFF(-2,2) = coeff_dyy_m2_p2;
+COEFF(-1,2) = coeff_dyy_m1_p2;
+COEFF(0,2) = coeff_dyy_0_p2;
+COEFF(1,2) = coeff_dyy_p1_p2;
+COEFF(2,2) = coeff_dyy_p2_p2;
author	jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416>	2002-05-14 15:03:52 +0000
committer	jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416>	2002-05-14 15:03:52 +0000
commit	2177a5948dc567e37964df7f281d3a3ee11eb8c6 (patch)
tree	9e63d67ade8cadd69cb3bd3423bc8a461d8ed699 /src/GeneralizedPolynomial-Uniform/2d.coeffs
parent	918fdb59c1bc57ae64c8f26db86dd50f316d557e (diff)