This is a major reworking of this interpolator.

The interpolation code in template.c is now split up into subfunctions
(in common/, Lagrange/, and Hermite/), so it should now compile in
finite cpu/memory even with optimization.


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

12 files changed, 516 insertions, 518 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c
deleted file mode 100644
index a022b8a..0000000
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c
+++ /dev/null
@@ -1,196 +0,0 @@
-fp t160,
-   t131,
-   t159,
-   t133,
-   t158,
-   t117,
-   t66,
-   t132,
-   t128,
-   t87,
-   t157,
-   t129,
-   t84,
-   t156,
-   t124,
-   t68,
-   t130,
-   t116,
-   t77,
-   t155,
-   t126,
-   t70,
-   t125,
-   t83,
-   t154,
-   t71,
-   t153,
-   t108,
-   t122,
-   t69,
-   t152,
-   t151,
-   t63,
-   t150,
-   t67,
-   t86,
-   t149,
-   t65,
-   t148,
-   t85,
-   t147,
-   t79,
-   t146,
-   t82,
-   t145,
-   t144,
-   t143,
-   t142,
-   t123,
-   t141,
-   t101,
-   t140,
-   t139,
-   t138,
-   t137,
-   t136,
-   t127,
-   t121,
-   t119,
-   t118,
-   t115,
-   t113,
-   t112,
-   t111,
-   t109,
-   t106,
-   t105,
-   t103,
-   t102,
-   t100,
-   t98,
-   t97,
-   t96,
-   t95,
-   t94,
-   t93,
-   t92,
-   t91,
-   t90,
-   t89,
-   t88,
-   t81,
-   t80,
-   t78,
-   t76,
-   t75,
-   t74,
-   t72;
-      t160 = x*y;
-      t131 = x*x;
-      t159 = t131*y;
-      t133 = y*y;
-      t158 = x*t133;
-      t117 = RATIONAL(1.0,40.0);
-      t66 = t117*t159;
-      t132 = t133*y;
-      t128 = RATIONAL(-1.0,8.0);
-      t87 = t128*t132;
-      t157 = t66+t87;
-      t129 = RATIONAL(1.0,8.0);
-      t84 = t129*t132;
-      t156 = t66+t84;
-      t124 = RATIONAL(3.0,40.0);
-      t68 = t124*t158;
-      t130 = t131*x;
-      t116 = RATIONAL(1.0,24.0);
-      t77 = t116*t130;
-      t155 = t68+t77;
-      t126 = RATIONAL(-3.0,40.0);
-      t70 = t126*t159;
-      t125 = RATIONAL(-1.0,24.0);
-      t83 = t125*t132;
-      t154 = t70+t83;
-      t71 = t124*t159;
-      t153 = t71+t83;
-      t108 = RATIONAL(11.0,80.0);
-      t122 = RATIONAL(-1.0,40.0);
-      t69 = t122*t158;
-      t152 = t69+t108*t133;
-      t151 = RATIONAL(-1.0,100.0)*t160+RATIONAL(37.0,200.0);
-      t63 = t122*t159;
-      t150 = t63+t87;
-      t67 = t117*t158;
-      t86 = t129*t130;
-      t149 = t67+t86;
-      t65 = t126*t158;
-      t148 = t65+t77;
-      t85 = t128*t130;
-      t147 = t67+t85;
-      t79 = t116*t132;
-      t146 = t70+t79;
-      t82 = t125*t130;
-      t145 = t68+t82;
-      t144 = t63+t84;
-      t143 = t65+t82;
-      t142 = t71+t79;
-      t123 = RATIONAL(-1.0,50.0);
-      t141 = t123+RATIONAL(-13.0,100.0)*t160;
-      t101 = RATIONAL(-17.0,200.0);
-      t140 = RATIONAL(-9.0,100.0)*t160+t101;
-      t139 = RATIONAL(11.0,50.0)+t123*t160;
-      t138 = RATIONAL(2.0,25.0)*t160+RATIONAL(-2.0,25.0);
-      t137 = RATIONAL(7.0,100.0)*t160+RATIONAL(13.0,100.0);
-      t136 = t133+t131;
-      t127 = RATIONAL(1.0,80.0);
-      t121 = RATIONAL(9.0,80.0);
-      t119 = RATIONAL(-1.0,80.0);
-      t118 = RATIONAL(7.0,80.0);
-      t115 = RATIONAL(3.0,80.0);
-      t113 = RATIONAL(-19.0,80.0);
-      t112 = RATIONAL(-17.0,80.0);
-      t111 = RATIONAL(13.0,80.0);
-      t109 = RATIONAL(-23.0,80.0);
-      t106 = RATIONAL(-63.0,400.0);
-      t105 = RATIONAL(-49.0,400.0);
-      t103 = RATIONAL(-31.0,400.0);
-      t102 = RATIONAL(69.0,400.0);
-      t100 = RATIONAL(-57.0,400.0);
-      t98 = RATIONAL(117.0,400.0);
-      t97 = RATIONAL(43.0,1200.0);
-      t96 = RATIONAL(103.0,400.0);
-      t95 = RATIONAL(37.0,1200.0);
-      t94 = RATIONAL(111.0,400.0);
-      t93 = RATIONAL(89.0,1200.0);
-      t92 = RATIONAL(-109.0,1200.0);
-      t91 = RATIONAL(-223.0,1200.0);
-      t90 = RATIONAL(-157.0,1200.0);
-      t89 = RATIONAL(-149.0,1200.0);
-      t88 = RATIONAL(-131.0,1200.0);
-      t81 = t121*t133;
-      t80 = t118*t131;
-      t78 = t118*t133;
-      t76 = t121*t131;
-      t75 = t111*t131;
-      t74 = t111*t133;
-      t72 = t108*t131;
-      coeff_I_m1_m1 = t75+RATIONAL(-23.0,200.0)+t74+t92*y+(t92+RATIONAL(6.0,
-25.0)*y)*x+t143+t154;
-      coeff_I_0_m1 = t109*t131+t91*y+t86+t103*x+t139+t152+t153;
-      coeff_I_p1_m1 = t81+t94*x+t90*y+t80+t141+t147+t153;
-      coeff_I_p2_m1 = t78+t93*y+t88*x+t115*t131+t140+t154+t155;
-      coeff_I_m1_0 = t72+t109*t133+t91*x+t103*y+t139+t144+t145;
-      coeff_I_0_0 = RATIONAL(93.0,200.0)+t100*y+t136*RATIONAL(-21.0,80.0)+(
-RATIONAL(-1.0,25.0)*y+t100)*x+t149+t156;
-      coeff_I_p1_0 = t69+t76+t98*x+t85+t113*t133+t106*y+t151+t156;
-      coeff_I_p2_0 = t97*x+t127*t131+t105*y+t112*t133+t137+t144+t148;
-      coeff_I_m1_p1 = t76+t90*x+t78+t94*y+t141+t145+t157;
-      coeff_I_0_p1 = t113*t131+t81+t106*x+t98*y+t149+t150+t151;
-      coeff_I_p1_p1 = t72+t101+t96*y+t85+(t96+RATIONAL(3.0,50.0)*y)*x+t150+t152
-;
-      coeff_I_p2_p1 = t102*y+t74+t119*t131+t95*x+t138+t148+t157;
-      coeff_I_m1_p2 = t80+t88*y+t93*x+t115*t133+t140+t142+t143;
-      coeff_I_0_p2 = t127*t133+t86+t97*y+t105*x+t112*t131+t69+t137+t146;
-      coeff_I_p1_p2 = t95*y+t75+t119*t133+t102*x+t138+t146+t147;
-      coeff_I_p2_p2 = t89*y+RATIONAL(7.0,200.0)+t136*RATIONAL(-3.0,80.0)+(t89+
-RATIONAL(-3.0,50.0)*y)*x+t142+t155;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c
deleted file mode 100644
index 5bfa6f4..0000000
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c
+++ /dev/null
@@ -1,95 +0,0 @@
-fp t201,
-   t161,
-   t200,
-   t184,
-   t188,
-   t174,
-   t187,
-   t180,
-   t199,
-   t164,
-   t198,
-   t183,
-   t175,
-   t197,
-   t162,
-   t196,
-   t181,
-   t195,
-   t194,
-   t163,
-   t193,
-   t186,
-   t176,
-   t182,
-   t192,
-   t185,
-   t172,
-   t191,
-   t179,
-   t190,
-   t189,
-   t178,
-   t177,
-   t169,
-   t167,
-   t166,
-   t165;
-      t201 = x*y;
-      t161 = RATIONAL(1.0,20.0)*t201;
-      t200 = t161+RATIONAL(9.0,40.0)*x;
-      t184 = RATIONAL(-1.0,40.0);
-      t188 = y*y;
-      t174 = t184*t188;
-      t187 = x*x;
-      t180 = RATIONAL(-3.0,8.0)*t187;
-      t199 = t174+t180;
-      t164 = RATIONAL(-3.0,20.0)*t201;
-      t198 = t164+RATIONAL(13.0,40.0)*x;
-      t183 = RATIONAL(1.0,40.0);
-      t175 = t183*t188;
-      t197 = t175+t180;
-      t162 = RATIONAL(-1.0,20.0)*t201;
-      t196 = t162+RATIONAL(11.0,40.0)*x;
-      t181 = RATIONAL(3.0,8.0)*t187;
-      t195 = t174+t181;
-      t194 = t175+t181;
-      t163 = RATIONAL(3.0,20.0)*t201;
-      t193 = RATIONAL(7.0,40.0)*x+t163;
-      t186 = RATIONAL(-3.0,40.0);
-      t176 = t186*t188;
-      t182 = RATIONAL(-1.0,8.0)*t187;
-      t192 = t176+t182;
-      t185 = RATIONAL(3.0,40.0);
-      t172 = t185*t188;
-      t191 = t172+t182;
-      t179 = RATIONAL(1.0,8.0)*t187;
-      t190 = t172+t179;
-      t189 = t176+t179;
-      t178 = RATIONAL(2.0,25.0)*y;
-      t177 = 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;
-      t165 = RATIONAL(-13.0,100.0)*y;
-      coeff_dx_m1_m1 = RATIONAL(-109.0,1200.0)+RATIONAL(6.0,25.0)*y+t192+t198;
-      coeff_dx_0_m1 = RATIONAL(-31.0,400.0)+t177+RATIONAL(-23.0,40.0)*x+t163+
-t195;
-      coeff_dx_p1_m1 = RATIONAL(111.0,400.0)+t165+t193+t197;
-      coeff_dx_p2_m1 = t169+RATIONAL(-131.0,1200.0)+t164+t185*x+t190;
-      coeff_dx_m1_0 = RATIONAL(-223.0,1200.0)+t177+t191+t196;
-      coeff_dx_0_0 = RATIONAL(-1.0,25.0)*y+t161+RATIONAL(-21.0,40.0)*x+RATIONAL
-(-57.0,400.0)+t194;
-      coeff_dx_p1_0 = t167+RATIONAL(117.0,400.0)+t199+t200;
-      coeff_dx_p2_0 = t183*x+t162+RATIONAL(43.0,1200.0)+t166+t189;
-      coeff_dx_m1_p1 = RATIONAL(-157.0,1200.0)+t165+t191+t200;
-      coeff_dx_0_p1 = RATIONAL(-19.0,40.0)*x+RATIONAL(-63.0,400.0)+t162+t167+
-t194;
-      coeff_dx_p1_p1 = RATIONAL(3.0,50.0)*y+RATIONAL(103.0,400.0)+t196+t199;
-      coeff_dx_p2_p1 = t161+t184*x+RATIONAL(37.0,1200.0)+t178+t189;
-      coeff_dx_m1_p2 = RATIONAL(89.0,1200.0)+t169+t192+t193;
-      coeff_dx_0_p2 = RATIONAL(-17.0,40.0)*x+t166+RATIONAL(-49.0,400.0)+t164+
-t195;
-      coeff_dx_p1_p2 = t178+RATIONAL(69.0,400.0)+t197+t198;
-      coeff_dx_p2_p2 = t163+RATIONAL(-149.0,1200.0)+t186*x+RATIONAL(-3.0,50.0)*
-y+t190;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c
deleted file mode 100644
index 29b5494..0000000
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c
+++ /dev/null
@@ -1,40 +0,0 @@
-fp t244,
-   t258,
-   t246,
-   t257,
-   t245,
-   t256,
-   t243,
-   t255,
-   t250,
-   t249,
-   t248,
-   t247;
-      t244 = RATIONAL(-3.0,20.0)*y;
-      t258 = t244+RATIONAL(13.0,40.0);
-      t246 = RATIONAL(-1.0,20.0)*y;
-      t257 = RATIONAL(11.0,40.0)+t246;
-      t245 = RATIONAL(1.0,20.0)*y;
-      t256 = RATIONAL(9.0,40.0)+t245;
-      t243 = RATIONAL(3.0,20.0)*y;
-      t255 = t243+RATIONAL(7.0,40.0);
-      t250 = RATIONAL(1.0,4.0)*x;
-      t249 = RATIONAL(-1.0,4.0)*x;
-      t248 = RATIONAL(3.0,4.0)*x;
-      t247 = RATIONAL(-3.0,4.0)*x;
-      coeff_dxx_m1_m1 = t249+t258;
-      coeff_dxx_0_m1 = t248+RATIONAL(-23.0,40.0)+t243;
-      coeff_dxx_p1_m1 = t247+t255;
-      coeff_dxx_p2_m1 = RATIONAL(3.0,40.0)+t250+t244;
-      coeff_dxx_m1_0 = t249+t257;
-      coeff_dxx_0_0 = RATIONAL(-21.0,40.0)+t245+t248;
-      coeff_dxx_p1_0 = t247+t256;
-      coeff_dxx_p2_0 = RATIONAL(1.0,40.0)+t246+t250;
-      coeff_dxx_m1_p1 = t249+t256;
-      coeff_dxx_0_p1 = t248+RATIONAL(-19.0,40.0)+t246;
-      coeff_dxx_p1_p1 = t247+t257;
-      coeff_dxx_p2_p1 = RATIONAL(-1.0,40.0)+t250+t245;
-      coeff_dxx_m1_p2 = t249+t255;
-      coeff_dxx_0_p2 = t244+RATIONAL(-17.0,40.0)+t248;
-      coeff_dxx_p1_p2 = t247+t258;
-      coeff_dxx_p2_p2 = RATIONAL(-3.0,40.0)+t250+t243;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c
deleted file mode 100644
index 20e2106..0000000
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c
+++ /dev/null
@@ -1,52 +0,0 @@
-fp t276,
-   t275,
-   t274,
-   t273,
-   t272,
-   t271,
-   t270,
-   t269,
-   t268,
-   t267,
-   t266,
-   t265,
-   t264,
-   t263,
-   t262,
-   t261,
-   t260,
-   t259;
-      t276 = RATIONAL(-1.0,50.0);
-      t275 = RATIONAL(2.0,25.0);
-      t274 = RATIONAL(-3.0,20.0);
-      t273 = RATIONAL(3.0,20.0);
-      t272 = RATIONAL(-1.0,20.0);
-      t271 = RATIONAL(1.0,20.0);
-      t270 = RATIONAL(-9.0,100.0);
-      t269 = RATIONAL(-1.0,100.0);
-      t268 = RATIONAL(7.0,100.0);
-      t267 = RATIONAL(-13.0,100.0);
-      t266 = t272*y;
-      t265 = t271*y;
-      t264 = t273*x;
-      t263 = t274*x;
-      t262 = t271*x;
-      t261 = t272*x;
-      t260 = t274*y;
-      t259 = t273*y;
-      coeff_dxy_m1_m1 = t260+RATIONAL(6.0,25.0)+t263;
-      coeff_dxy_0_m1 = t276+t264+t266;
-      coeff_dxy_p1_m1 = t264+t267+t265;
-      coeff_dxy_p2_m1 = t263+t270+t259;
-      coeff_dxy_m1_0 = t276+t259+t261;
-      coeff_dxy_0_0 = t262+RATIONAL(-1.0,25.0)+t265;
-      coeff_dxy_p1_0 = t269+t262+t266;
-      coeff_dxy_p2_0 = t260+t268+t261;
-      coeff_dxy_m1_p1 = t262+t267+t259;
-      coeff_dxy_0_p1 = t269+t265+t261;
-      coeff_dxy_p1_p1 = RATIONAL(3.0,50.0)+t266+t261;
-      coeff_dxy_p2_p1 = t260+t275+t262;
-      coeff_dxy_m1_p2 = t264+t260+t270;
-      coeff_dxy_0_p2 = t268+t263+t266;
-      coeff_dxy_p1_p2 = t263+t275+t265;
-      coeff_dxy_p2_p2 = t259+RATIONAL(-3.0,50.0)+t264;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c
deleted file mode 100644
index c766157..0000000
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c
+++ /dev/null
@@ -1,95 +0,0 @@
-fp t242,
-   t225,
-   t228,
-   t212,
-   t229,
-   t220,
-   t241,
-   t203,
-   t240,
-   t204,
-   t239,
-   t205,
-   t238,
-   t224,
-   t213,
-   t237,
-   t202,
-   t236,
-   t226,
-   t215,
-   t222,
-   t235,
-   t227,
-   t218,
-   t234,
-   t221,
-   t233,
-   t232,
-   t223,
-   t231,
-   t230,
-   t219,
-   t217,
-   t209,
-   t208,
-   t207,
-   t206;
-      t242 = x*y;
-      t225 = RATIONAL(-1.0,40.0);
-      t228 = x*x;
-      t212 = t225*t228;
-      t229 = y*y;
-      t220 = RATIONAL(-3.0,8.0)*t229;
-      t241 = t212+t220;
-      t203 = RATIONAL(1.0,20.0)*t242;
-      t240 = t203+RATIONAL(9.0,40.0)*y;
-      t204 = RATIONAL(-1.0,20.0)*t242;
-      t239 = RATIONAL(11.0,40.0)*y+t204;
-      t205 = RATIONAL(-3.0,20.0)*t242;
-      t238 = t205+RATIONAL(13.0,40.0)*y;
-      t224 = RATIONAL(1.0,40.0);
-      t213 = t224*t228;
-      t237 = t213+t220;
-      t202 = RATIONAL(3.0,20.0)*t242;
-      t236 = RATIONAL(7.0,40.0)*y+t202;
-      t226 = RATIONAL(3.0,40.0);
-      t215 = t226*t228;
-      t222 = RATIONAL(-1.0,8.0)*t229;
-      t235 = t215+t222;
-      t227 = RATIONAL(-3.0,40.0);
-      t218 = t227*t228;
-      t234 = t218+t222;
-      t221 = RATIONAL(3.0,8.0)*t229;
-      t233 = t212+t221;
-      t232 = t213+t221;
-      t223 = RATIONAL(1.0,8.0)*t229;
-      t231 = t218+t223;
-      t230 = t215+t223;
-      t219 = RATIONAL(2.0,25.0)*x;
-      t217 = RATIONAL(-1.0,50.0)*x;
-      t209 = RATIONAL(-9.0,100.0)*x;
-      t208 = RATIONAL(7.0,100.0)*x;
-      t207 = RATIONAL(-1.0,100.0)*x;
-      t206 = RATIONAL(-13.0,100.0)*x;
-      coeff_dy_m1_m1 = RATIONAL(6.0,25.0)*x+RATIONAL(-109.0,1200.0)+t234+t238;
-      coeff_dy_0_m1 = RATIONAL(-223.0,1200.0)+t217+t235+t239;
-      coeff_dy_p1_m1 = t206+RATIONAL(-157.0,1200.0)+t235+t240;
-      coeff_dy_p2_m1 = t209+RATIONAL(89.0,1200.0)+t234+t236;
-      coeff_dy_m1_0 = RATIONAL(-23.0,40.0)*y+t202+RATIONAL(-31.0,400.0)+t217+
-t233;
-      coeff_dy_0_0 = RATIONAL(-1.0,25.0)*x+RATIONAL(-57.0,400.0)+t203+RATIONAL(
--21.0,40.0)*y+t232;
-      coeff_dy_p1_0 = RATIONAL(-63.0,400.0)+t207+t204+RATIONAL(-19.0,40.0)*y+
-t232;
-      coeff_dy_p2_0 = t205+RATIONAL(-17.0,40.0)*y+RATIONAL(-49.0,400.0)+t208+
-t233;
-      coeff_dy_m1_p1 = RATIONAL(111.0,400.0)+t206+t236+t237;
-      coeff_dy_0_p1 = RATIONAL(117.0,400.0)+t207+t240+t241;
-      coeff_dy_p1_p1 = RATIONAL(3.0,50.0)*x+RATIONAL(103.0,400.0)+t239+t241;
-      coeff_dy_p2_p1 = t219+RATIONAL(69.0,400.0)+t237+t238;
-      coeff_dy_m1_p2 = t209+RATIONAL(-131.0,1200.0)+t205+t226*y+t230;
-      coeff_dy_0_p2 = t204+t208+t224*y+RATIONAL(43.0,1200.0)+t231;
-      coeff_dy_p1_p2 = t225*y+t219+t203+RATIONAL(37.0,1200.0)+t231;
-      coeff_dy_p2_p2 = RATIONAL(-149.0,1200.0)+t202+RATIONAL(-3.0,50.0)*x+t227*
-y+t230;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c
deleted file mode 100644
index b449df9..0000000
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c
+++ /dev/null
@@ -1,40 +0,0 @@
-fp t277,
-   t292,
-   t279,
-   t291,
-   t280,
-   t290,
-   t278,
-   t289,
-   t284,
-   t283,
-   t282,
-   t281;
-      t277 = RATIONAL(3.0,20.0)*x;
-      t292 = t277+RATIONAL(7.0,40.0);
-      t279 = RATIONAL(-1.0,20.0)*x;
-      t291 = t279+RATIONAL(11.0,40.0);
-      t280 = RATIONAL(-3.0,20.0)*x;
-      t290 = t280+RATIONAL(13.0,40.0);
-      t278 = RATIONAL(1.0,20.0)*x;
-      t289 = t278+RATIONAL(9.0,40.0);
-      t284 = RATIONAL(-1.0,4.0)*y;
-      t283 = RATIONAL(-3.0,4.0)*y;
-      t282 = RATIONAL(1.0,4.0)*y;
-      t281 = RATIONAL(3.0,4.0)*y;
-      coeff_dyy_m1_m1 = t284+t290;
-      coeff_dyy_0_m1 = t284+t291;
-      coeff_dyy_p1_m1 = t284+t289;
-      coeff_dyy_p2_m1 = t284+t292;
-      coeff_dyy_m1_0 = RATIONAL(-23.0,40.0)+t281+t277;
-      coeff_dyy_0_0 = t281+RATIONAL(-21.0,40.0)+t278;
-      coeff_dyy_p1_0 = RATIONAL(-19.0,40.0)+t279+t281;
-      coeff_dyy_p2_0 = t281+RATIONAL(-17.0,40.0)+t280;
-      coeff_dyy_m1_p1 = t283+t292;
-      coeff_dyy_0_p1 = t283+t289;
-      coeff_dyy_p1_p1 = t283+t291;
-      coeff_dyy_p2_p1 = t283+t290;
-      coeff_dyy_m1_p2 = t282+RATIONAL(3.0,40.0)+t280;
-      coeff_dyy_0_p2 = t282+RATIONAL(1.0,40.0)+t279;
-      coeff_dyy_p1_p2 = RATIONAL(-1.0,40.0)+t278+t282;
-      coeff_dyy_p2_p2 = t277+t282+RATIONAL(-3.0,40.0);
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c
new file mode 100644
index 0000000..ac11344
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c
@@ -0,0 +1,190 @@
+fp t158;
+fp t131;
+fp t157;
+fp t132;
+fp t156;
+fp t117;
+fp t63;
+fp t129;
+fp t130;
+fp t84;
+fp t155;
+fp t154;
+fp t128;
+fp t87;
+fp t153;
+fp t126;
+fp t64;
+fp t133;
+fp t116;
+fp t76;
+fp t152;
+fp t67;
+fp t86;
+fp t151;
+fp t123;
+fp t71;
+fp t150;
+fp t124;
+fp t65;
+fp t149;
+fp t125;
+fp t81;
+fp t148;
+fp t85;
+fp t147;
+fp t146;
+fp t112;
+fp t69;
+fp t145;
+fp t110;
+fp t66;
+fp t144;
+fp t143;
+fp t70;
+fp t120;
+fp t142;
+fp t122;
+fp t141;
+fp t140;
+fp t101;
+fp t139;
+fp t138;
+fp t137;
+fp t136;
+fp t127;
+fp t121;
+fp t118;
+fp t115;
+fp t113;
+fp t109;
+fp t107;
+fp t106;
+fp t104;
+fp t103;
+fp t102;
+fp t100;
+fp t98;
+fp t97;
+fp t96;
+fp t95;
+fp t94;
+fp t93;
+fp t92;
+fp t91;
+fp t90;
+fp t89;
+fp t88;
+fp t83;
+fp t82;
+fp t80;
+fp t78;
+fp t77;
+fp t74;
+fp t73;
+      t158 = x*y;
+      t131 = x*x;
+      t157 = t131*y;
+      t132 = y*y;
+      t156 = x*t132;
+      t117 = RATIONAL(1.0,40.0);
+      t63 = t117*t157;
+      t129 = RATIONAL(1.0,8.0);
+      t130 = t132*y;
+      t84 = t129*t130;
+      t155 = t63+t84;
+      t154 = RATIONAL(-1.0,100.0)*t158+RATIONAL(37.0,200.0);
+      t128 = RATIONAL(-1.0,8.0);
+      t87 = t128*t130;
+      t153 = t63+t87;
+      t126 = RATIONAL(-3.0,40.0);
+      t64 = t126*t156;
+      t133 = t131*x;
+      t116 = RATIONAL(1.0,24.0);
+      t76 = t116*t133;
+      t152 = t64+t76;
+      t67 = t117*t156;
+      t86 = t129*t133;
+      t151 = t67+t86;
+      t123 = RATIONAL(-1.0,40.0);
+      t71 = t123*t156;
+      t150 = t71+t86;
+      t124 = RATIONAL(3.0,40.0);
+      t65 = t124*t156;
+      t149 = t65+t76;
+      t125 = RATIONAL(-1.0,24.0);
+      t81 = t125*t133;
+      t148 = t64+t81;
+      t85 = t128*t133;
+      t147 = t71+t85;
+      t146 = t65+t81;
+      t112 = RATIONAL(13.0,80.0);
+      t69 = t126*t157;
+      t145 = t69+t112*t131;
+      t110 = RATIONAL(11.0,80.0);
+      t66 = t123*t157;
+      t144 = t66+t110*t131;
+      t143 = t67+t85;
+      t70 = t124*t157;
+      t120 = RATIONAL(7.0,80.0);
+      t142 = t70+t120*t131;
+      t122 = RATIONAL(-1.0,50.0);
+      t141 = t122+RATIONAL(-13.0,100.0)*t158;
+      t140 = RATIONAL(2.0,25.0)*t158+RATIONAL(-2.0,25.0);
+      t101 = RATIONAL(-17.0,200.0);
+      t139 = RATIONAL(-9.0,100.0)*t158+t101;
+      t138 = RATIONAL(11.0,50.0)+t122*t158;
+      t137 = RATIONAL(7.0,100.0)*t158+RATIONAL(13.0,100.0);
+      t136 = t131+t132;
+      t127 = RATIONAL(1.0,80.0);
+      t121 = RATIONAL(9.0,80.0);
+      t118 = RATIONAL(-1.0,80.0);
+      t115 = RATIONAL(3.0,80.0);
+      t113 = RATIONAL(-17.0,80.0);
+      t109 = RATIONAL(-23.0,80.0);
+      t107 = RATIONAL(-19.0,80.0);
+      t106 = RATIONAL(-63.0,400.0);
+      t104 = RATIONAL(-49.0,400.0);
+      t103 = RATIONAL(-31.0,400.0);
+      t102 = RATIONAL(69.0,400.0);
+      t100 = RATIONAL(-57.0,400.0);
+      t98 = RATIONAL(117.0,400.0);
+      t97 = RATIONAL(43.0,1200.0);
+      t96 = RATIONAL(103.0,400.0);
+      t95 = RATIONAL(111.0,400.0);
+      t94 = RATIONAL(37.0,1200.0);
+      t93 = RATIONAL(89.0,1200.0);
+      t92 = RATIONAL(-109.0,1200.0);
+      t91 = RATIONAL(-223.0,1200.0);
+      t90 = RATIONAL(-157.0,1200.0);
+      t89 = RATIONAL(-131.0,1200.0);
+      t88 = RATIONAL(-149.0,1200.0);
+      t83 = t121*t131;
+      t82 = t125*t130;
+      t80 = t121*t132;
+      t78 = t116*t130;
+      t77 = t120*t132;
+      t74 = t112*t132;
+      t73 = t110*t132;
+      coeffs_I->coeff_m1_m1 = t74+t82+t92*y+RATIONAL(-23.0,200.0)+(t92+RATIONAL
+(6.0,25.0)*y)*x+t145+t148;
+      coeffs_I->coeff_0_m1 = t70+t91*y+t73+t103*x+t82+t109*t131+t138+t150;
+      coeffs_I->coeff_p1_m1 = t80+t82+t95*x+t90*y+t141+t142+t143;
+      coeffs_I->coeff_p2_m1 = t93*y+t89*x+t77+t82+t69+t115*t131+t139+t149;
+      coeffs_I->coeff_m1_0 = t109*t132+t84+t91*x+t103*y+t138+t144+t146;
+      coeffs_I->coeff_0_0 = RATIONAL(93.0,200.0)+t100*y+t136*RATIONAL(-21.0,
+80.0)+(t100+RATIONAL(-1.0,25.0)*y)*x+t151+t155;
+      coeffs_I->coeff_p1_0 = t107*t132+t106*y+t83+t98*x+t147+t154+t155;
+      coeffs_I->coeff_p2_0 = t66+t113*t132+t127*t131+t84+t104*y+t97*x+t137+t152
+;
+      coeffs_I->coeff_m1_p1 = t95*y+t90*x+t77+t83+t141+t146+t153;
+      coeffs_I->coeff_0_p1 = t87+t106*x+t66+t80+t98*y+t107*t131+t151+t154;
+      coeffs_I->coeff_p1_p1 = t87+t73+t101+t96*y+(t96+RATIONAL(3.0,50.0)*y)*x+
+t144+t147;
+      coeffs_I->coeff_p2_p1 = t118*t131+t102*y+t74+t94*x+t140+t152+t153;
+      coeffs_I->coeff_m1_p2 = t93*x+t89*y+t115*t132+t78+t139+t142+t148;
+      coeffs_I->coeff_0_p2 = t97*y+t69+t127*t132+t113*t131+t78+t104*x+t137+t150
+;
+      coeffs_I->coeff_p1_p2 = t102*x+t94*y+t78+t118*t132+t140+t143+t145;
+      coeffs_I->coeff_p2_p2 = t70+t88*y+t78+RATIONAL(7.0,200.0)+t136*RATIONAL(
+-3.0,80.0)+(t88+RATIONAL(-3.0,50.0)*y)*x+t149;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c
new file mode 100644
index 0000000..0ea2d65
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c
@@ -0,0 +1,97 @@
+fp t199;
+fp t185;
+fp t178;
+fp t181;
+fp t186;
+fp t173;
+fp t198;
+fp t159;
+fp t197;
+fp t182;
+fp t175;
+fp t196;
+fp t161;
+fp t195;
+fp t162;
+fp t194;
+fp t160;
+fp t193;
+fp t179;
+fp t184;
+fp t176;
+fp t192;
+fp t183;
+fp t171;
+fp t191;
+fp t180;
+fp t190;
+fp t189;
+fp t177;
+fp t188;
+fp t187;
+fp t174;
+fp t169;
+fp t167;
+fp t165;
+fp t164;
+fp t163;
+      t199 = x*y;
+      t185 = x*x;
+      t178 = RATIONAL(-3.0,8.0)*t185;
+      t181 = RATIONAL(1.0,40.0);
+      t186 = y*y;
+      t173 = t181*t186;
+      t198 = t178+t173;
+      t159 = RATIONAL(-3.0,20.0)*t199;
+      t197 = RATIONAL(13.0,40.0)*x+t159;
+      t182 = RATIONAL(-1.0,40.0);
+      t175 = t182*t186;
+      t196 = t178+t175;
+      t161 = RATIONAL(-1.0,20.0)*t199;
+      t195 = t161+RATIONAL(11.0,40.0)*x;
+      t162 = RATIONAL(3.0,20.0)*t199;
+      t194 = RATIONAL(7.0,40.0)*x+t162;
+      t160 = RATIONAL(1.0,20.0)*t199;
+      t193 = RATIONAL(9.0,40.0)*x+t160;
+      t179 = RATIONAL(-1.0,8.0)*t185;
+      t184 = RATIONAL(-3.0,40.0);
+      t176 = t184*t186;
+      t192 = t179+t176;
+      t183 = RATIONAL(3.0,40.0);
+      t171 = t183*t186;
+      t191 = t179+t171;
+      t180 = RATIONAL(3.0,8.0)*t185;
+      t190 = t180+t175;
+      t189 = t180+t173;
+      t177 = RATIONAL(1.0,8.0)*t185;
+      t188 = t171+t177;
+      t187 = t176+t177;
+      t174 = RATIONAL(-1.0,50.0)*y;
+      t169 = RATIONAL(2.0,25.0)*y;
+      t167 = RATIONAL(-9.0,100.0)*y;
+      t165 = RATIONAL(-1.0,100.0)*y;
+      t164 = RATIONAL(7.0,100.0)*y;
+      t163 = RATIONAL(-13.0,100.0)*y;
+      coeffs_dx->coeff_m1_m1 = RATIONAL(6.0,25.0)*y+RATIONAL(-109.0,1200.0)+
+t192+t197;
+      coeffs_dx->coeff_0_m1 = t162+t174+RATIONAL(-31.0,400.0)+RATIONAL(-23.0,
+40.0)*x+t190;
+      coeffs_dx->coeff_p1_m1 = RATIONAL(111.0,400.0)+t163+t194+t198;
+      coeffs_dx->coeff_p2_m1 = t183*x+RATIONAL(-131.0,1200.0)+t167+t159+t188;
+      coeffs_dx->coeff_m1_0 = RATIONAL(-223.0,1200.0)+t174+t191+t195;
+      coeffs_dx->coeff_0_0 = RATIONAL(-1.0,25.0)*y+t160+RATIONAL(-21.0,40.0)*x+
+RATIONAL(-57.0,400.0)+t189;
+      coeffs_dx->coeff_p1_0 = t165+RATIONAL(117.0,400.0)+t193+t196;
+      coeffs_dx->coeff_p2_0 = t161+t181*x+RATIONAL(43.0,1200.0)+t164+t187;
+      coeffs_dx->coeff_m1_p1 = RATIONAL(-157.0,1200.0)+t163+t191+t193;
+      coeffs_dx->coeff_0_p1 = RATIONAL(-63.0,400.0)+t165+t161+RATIONAL(-19.0,
+40.0)*x+t189;
+      coeffs_dx->coeff_p1_p1 = RATIONAL(103.0,400.0)+RATIONAL(3.0,50.0)*y+t195+
+t196;
+      coeffs_dx->coeff_p2_p1 = t160+t169+RATIONAL(37.0,1200.0)+t182*x+t187;
+      coeffs_dx->coeff_m1_p2 = RATIONAL(89.0,1200.0)+t167+t192+t194;
+      coeffs_dx->coeff_0_p2 = t159+t164+RATIONAL(-17.0,40.0)*x+RATIONAL(-49.0,
+400.0)+t190;
+      coeffs_dx->coeff_p1_p2 = RATIONAL(69.0,400.0)+t169+t197+t198;
+      coeffs_dx->coeff_p2_p2 = RATIONAL(-149.0,1200.0)+t162+RATIONAL(-3.0,50.0)
+*y+t184*x+t188;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c
new file mode 100644
index 0000000..ba139bf
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c
@@ -0,0 +1,40 @@
+fp t244;
+fp t256;
+fp t242;
+fp t255;
+fp t241;
+fp t254;
+fp t243;
+fp t253;
+fp t248;
+fp t247;
+fp t246;
+fp t245;
+      t244 = RATIONAL(-3.0,20.0)*y;
+      t256 = t244+RATIONAL(13.0,40.0);
+      t242 = RATIONAL(-1.0,20.0)*y;
+      t255 = t242+RATIONAL(11.0,40.0);
+      t241 = RATIONAL(1.0,20.0)*y;
+      t254 = t241+RATIONAL(9.0,40.0);
+      t243 = RATIONAL(3.0,20.0)*y;
+      t253 = t243+RATIONAL(7.0,40.0);
+      t248 = RATIONAL(-1.0,4.0)*x;
+      t247 = RATIONAL(1.0,4.0)*x;
+      t246 = RATIONAL(-3.0,4.0)*x;
+      t245 = RATIONAL(3.0,4.0)*x;
+      coeffs_dxx->coeff_m1_m1 = t248+t256;
+      coeffs_dxx->coeff_0_m1 = t245+t243+RATIONAL(-23.0,40.0);
+      coeffs_dxx->coeff_p1_m1 = t246+t253;
+      coeffs_dxx->coeff_p2_m1 = t244+RATIONAL(3.0,40.0)+t247;
+      coeffs_dxx->coeff_m1_0 = t248+t255;
+      coeffs_dxx->coeff_0_0 = RATIONAL(-21.0,40.0)+t241+t245;
+      coeffs_dxx->coeff_p1_0 = t246+t254;
+      coeffs_dxx->coeff_p2_0 = t247+t242+RATIONAL(1.0,40.0);
+      coeffs_dxx->coeff_m1_p1 = t248+t254;
+      coeffs_dxx->coeff_0_p1 = RATIONAL(-19.0,40.0)+t242+t245;
+      coeffs_dxx->coeff_p1_p1 = t246+t255;
+      coeffs_dxx->coeff_p2_p1 = RATIONAL(-1.0,40.0)+t241+t247;
+      coeffs_dxx->coeff_m1_p2 = t248+t253;
+      coeffs_dxx->coeff_0_p2 = t244+t245+RATIONAL(-17.0,40.0);
+      coeffs_dxx->coeff_p1_p2 = t246+t256;
+      coeffs_dxx->coeff_p2_p2 = RATIONAL(-3.0,40.0)+t247+t243;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c
new file mode 100644
index 0000000..e9208a2
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c
@@ -0,0 +1,52 @@
+fp t274;
+fp t273;
+fp t272;
+fp t271;
+fp t270;
+fp t269;
+fp t268;
+fp t267;
+fp t266;
+fp t265;
+fp t264;
+fp t263;
+fp t262;
+fp t261;
+fp t260;
+fp t259;
+fp t258;
+fp t257;
+      t274 = RATIONAL(-3.0,20.0);
+      t273 = RATIONAL(3.0,20.0);
+      t272 = RATIONAL(2.0,25.0);
+      t271 = RATIONAL(-1.0,20.0);
+      t270 = RATIONAL(1.0,20.0);
+      t269 = RATIONAL(-1.0,50.0);
+      t268 = RATIONAL(-1.0,100.0);
+      t267 = RATIONAL(7.0,100.0);
+      t266 = RATIONAL(-9.0,100.0);
+      t265 = RATIONAL(-13.0,100.0);
+      t264 = t274*x;
+      t263 = t273*x;
+      t262 = t271*x;
+      t261 = t274*y;
+      t260 = t270*x;
+      t259 = t273*y;
+      t258 = t271*y;
+      t257 = t270*y;
+      coeffs_dxy->coeff_m1_m1 = t264+RATIONAL(6.0,25.0)+t261;
+      coeffs_dxy->coeff_0_m1 = t263+t269+t258;
+      coeffs_dxy->coeff_p1_m1 = t257+t265+t263;
+      coeffs_dxy->coeff_p2_m1 = t259+t264+t266;
+      coeffs_dxy->coeff_m1_0 = t259+t262+t269;
+      coeffs_dxy->coeff_0_0 = t260+RATIONAL(-1.0,25.0)+t257;
+      coeffs_dxy->coeff_p1_0 = t258+t268+t260;
+      coeffs_dxy->coeff_p2_0 = t261+t262+t267;
+      coeffs_dxy->coeff_m1_p1 = t259+t265+t260;
+      coeffs_dxy->coeff_0_p1 = t262+t257+t268;
+      coeffs_dxy->coeff_p1_p1 = t262+RATIONAL(3.0,50.0)+t258;
+      coeffs_dxy->coeff_p2_p1 = t261+t260+t272;
+      coeffs_dxy->coeff_m1_p2 = t266+t261+t263;
+      coeffs_dxy->coeff_0_p2 = t267+t258+t264;
+      coeffs_dxy->coeff_p1_p2 = t264+t257+t272;
+      coeffs_dxy->coeff_p2_p2 = t263+t259+RATIONAL(-3.0,50.0);
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c
new file mode 100644
index 0000000..c12d9d1
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c
@@ -0,0 +1,97 @@
+fp t240;
+fp t227;
+fp t219;
+fp t226;
+fp t222;
+fp t211;
+fp t239;
+fp t201;
+fp t238;
+fp t223;
+fp t213;
+fp t237;
+fp t200;
+fp t236;
+fp t203;
+fp t235;
+fp t202;
+fp t234;
+fp t224;
+fp t216;
+fp t221;
+fp t233;
+fp t225;
+fp t217;
+fp t232;
+fp t218;
+fp t231;
+fp t230;
+fp t220;
+fp t229;
+fp t228;
+fp t215;
+fp t210;
+fp t207;
+fp t206;
+fp t205;
+fp t204;
+      t240 = x*y;
+      t227 = y*y;
+      t219 = RATIONAL(-3.0,8.0)*t227;
+      t226 = x*x;
+      t222 = RATIONAL(1.0,40.0);
+      t211 = t222*t226;
+      t239 = t219+t211;
+      t201 = RATIONAL(1.0,20.0)*t240;
+      t238 = t201+RATIONAL(9.0,40.0)*y;
+      t223 = RATIONAL(-1.0,40.0);
+      t213 = t223*t226;
+      t237 = t219+t213;
+      t200 = RATIONAL(3.0,20.0)*t240;
+      t236 = t200+RATIONAL(7.0,40.0)*y;
+      t203 = RATIONAL(-3.0,20.0)*t240;
+      t235 = t203+RATIONAL(13.0,40.0)*y;
+      t202 = RATIONAL(-1.0,20.0)*t240;
+      t234 = t202+RATIONAL(11.0,40.0)*y;
+      t224 = RATIONAL(3.0,40.0);
+      t216 = t224*t226;
+      t221 = RATIONAL(-1.0,8.0)*t227;
+      t233 = t216+t221;
+      t225 = RATIONAL(-3.0,40.0);
+      t217 = t225*t226;
+      t232 = t217+t221;
+      t218 = RATIONAL(1.0,8.0)*t227;
+      t231 = t218+t217;
+      t230 = t216+t218;
+      t220 = RATIONAL(3.0,8.0)*t227;
+      t229 = t220+t213;
+      t228 = t220+t211;
+      t215 = RATIONAL(-1.0,50.0)*x;
+      t210 = RATIONAL(2.0,25.0)*x;
+      t207 = RATIONAL(-9.0,100.0)*x;
+      t206 = RATIONAL(-1.0,100.0)*x;
+      t205 = RATIONAL(7.0,100.0)*x;
+      t204 = RATIONAL(-13.0,100.0)*x;
+      coeffs_dy->coeff_m1_m1 = RATIONAL(6.0,25.0)*x+RATIONAL(-109.0,1200.0)+
+t232+t235;
+      coeffs_dy->coeff_0_m1 = t215+RATIONAL(-223.0,1200.0)+t233+t234;
+      coeffs_dy->coeff_p1_m1 = RATIONAL(-157.0,1200.0)+t204+t233+t238;
+      coeffs_dy->coeff_p2_m1 = t207+RATIONAL(89.0,1200.0)+t232+t236;
+      coeffs_dy->coeff_m1_0 = RATIONAL(-23.0,40.0)*y+RATIONAL(-31.0,400.0)+t200
++t215+t229;
+      coeffs_dy->coeff_0_0 = RATIONAL(-21.0,40.0)*y+RATIONAL(-57.0,400.0)+
+RATIONAL(-1.0,25.0)*x+t201+t228;
+      coeffs_dy->coeff_p1_0 = t202+RATIONAL(-19.0,40.0)*y+RATIONAL(-63.0,400.0)
++t206+t228;
+      coeffs_dy->coeff_p2_0 = t203+RATIONAL(-49.0,400.0)+RATIONAL(-17.0,40.0)*y
++t205+t229;
+      coeffs_dy->coeff_m1_p1 = t204+RATIONAL(111.0,400.0)+t236+t239;
+      coeffs_dy->coeff_0_p1 = t206+RATIONAL(117.0,400.0)+t237+t238;
+      coeffs_dy->coeff_p1_p1 = RATIONAL(3.0,50.0)*x+RATIONAL(103.0,400.0)+t234+
+t237;
+      coeffs_dy->coeff_p2_p1 = t210+RATIONAL(69.0,400.0)+t235+t239;
+      coeffs_dy->coeff_m1_p2 = t207+t203+RATIONAL(-131.0,1200.0)+t224*y+t230;
+      coeffs_dy->coeff_0_p2 = t205+t202+t222*y+RATIONAL(43.0,1200.0)+t231;
+      coeffs_dy->coeff_p1_p2 = t201+RATIONAL(37.0,1200.0)+t223*y+t210+t231;
+      coeffs_dy->coeff_p2_p2 = t225*y+t200+RATIONAL(-3.0,50.0)*x+RATIONAL(
+-149.0,1200.0)+t230;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c
new file mode 100644
index 0000000..6bbdde7
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c
@@ -0,0 +1,40 @@
+fp t276;
+fp t290;
+fp t277;
+fp t289;
+fp t278;
+fp t288;
+fp t275;
+fp t287;
+fp t282;
+fp t281;
+fp t280;
+fp t279;
+      t276 = RATIONAL(-1.0,20.0)*x;
+      t290 = t276+RATIONAL(11.0,40.0);
+      t277 = RATIONAL(3.0,20.0)*x;
+      t289 = RATIONAL(7.0,40.0)+t277;
+      t278 = RATIONAL(-3.0,20.0)*x;
+      t288 = t278+RATIONAL(13.0,40.0);
+      t275 = RATIONAL(1.0,20.0)*x;
+      t287 = t275+RATIONAL(9.0,40.0);
+      t282 = RATIONAL(1.0,4.0)*y;
+      t281 = RATIONAL(-1.0,4.0)*y;
+      t280 = RATIONAL(3.0,4.0)*y;
+      t279 = RATIONAL(-3.0,4.0)*y;
+      coeffs_dyy->coeff_m1_m1 = t281+t288;
+      coeffs_dyy->coeff_0_m1 = t281+t290;
+      coeffs_dyy->coeff_p1_m1 = t281+t287;
+      coeffs_dyy->coeff_p2_m1 = t281+t289;
+      coeffs_dyy->coeff_m1_0 = t277+t280+RATIONAL(-23.0,40.0);
+      coeffs_dyy->coeff_0_0 = t280+t275+RATIONAL(-21.0,40.0);
+      coeffs_dyy->coeff_p1_0 = t280+RATIONAL(-19.0,40.0)+t276;
+      coeffs_dyy->coeff_p2_0 = RATIONAL(-17.0,40.0)+t280+t278;
+      coeffs_dyy->coeff_m1_p1 = t279+t289;
+      coeffs_dyy->coeff_0_p1 = t279+t287;
+      coeffs_dyy->coeff_p1_p1 = t279+t290;
+      coeffs_dyy->coeff_p2_p1 = t279+t288;
+      coeffs_dyy->coeff_m1_p2 = t278+t282+RATIONAL(3.0,40.0);
+      coeffs_dyy->coeff_0_p2 = t282+RATIONAL(1.0,40.0)+t276;
+      coeffs_dyy->coeff_p1_p2 = t275+t282+RATIONAL(-1.0,40.0);
+      coeffs_dyy->coeff_p2_p2 = t277+RATIONAL(-3.0,40.0)+t282;