add Maple code and C skeletons for Hermite polynomial interpolation

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

12 files changed, 2332 insertions, 0 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-I.compute.c
new file mode 100644
index 0000000..b18719b
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-I.compute.c
@@ -0,0 +1,132 @@
+fp t41,
+   t40,
+   t68,
+   t38,
+   t39,
+   t67,
+   t48,
+   t49,
+   t36,
+   t69,
+   t28,
+   t66,
+   t35,
+   t34,
+   t56,
+   t7,
+   t52,
+   t65,
+   t60,
+   t37,
+   t64,
+   t61,
+   t63,
+   t27,
+   t55,
+   t21,
+   t32,
+   t51,
+   t62,
+   t59,
+   t58,
+   t57,
+   t54,
+   t53,
+   t31,
+   t50,
+   t25,
+   t29,
+   t47,
+   t46,
+   t33,
+   t30,
+   t26,
+   t24,
+   t23,
+   t22,
+   t20,
+   t19,
+   t18,
+   t16,
+   t12,
+   t11,
+   t10,
+   t9,
+   t8,
+   t6,
+   t4;
+      t41 = x*x;
+      t40 = t41*x;
+      t68 = t40+x;
+      t38 = y*y;
+      t39 = y*t38;
+      t67 = y+t39;
+      t48 = t68*t38;
+      t49 = t67*t41;
+      t36 = RATIONAL(1.0,4.0);
+      t69 = x*t36;
+      t28 = RATIONAL(-1.0,4.0);
+      t66 = t28*y;
+      t35 = RATIONAL(-5.0,2.0);
+      t34 = RATIONAL(-3.0,4.0);
+      t56 = t34*t40;
+      t7 = t39*t56;
+      t52 = t41*t38;
+      t65 = t7+t35*t52;
+      t60 = t28*t39;
+      t37 = RATIONAL(-1.0,2.0);
+      t64 = t40*t60+t37*t52;
+      t61 = t39*t40;
+      t63 = RATIONAL(-9.0,4.0)*t61+RATIONAL(-5.0,1.0)*t52;
+      t27 = RATIONAL(3.0,4.0);
+      t55 = t27*t39;
+      t21 = t40*t55;
+      t32 = RATIONAL(5.0,4.0);
+      t51 = t32*t38;
+      t62 = t41*t51+t21;
+      t59 = t36*t40;
+      t58 = t34*t39;
+      t57 = t27*t40;
+      t54 = t41*t39;
+      t53 = t40*t38;
+      t31 = RATIONAL(2.0,1.0);
+      t50 = t21+t31*t52+x*t66;
+      t25 = RATIONAL(-5.0,4.0);
+      t29 = RATIONAL(1.0,2.0);
+      t47 = t25*t38+t29;
+      t46 = t29+t25*t41;
+      t33 = RATIONAL(3.0,2.0);
+      t30 = RATIONAL(-3.0,2.0);
+      t26 = RATIONAL(3.0,1.0);
+      t24 = RATIONAL(-3.0,1.0);
+      t23 = RATIONAL(15.0,4.0);
+      t22 = RATIONAL(-15.0,4.0);
+      t20 = y*t69;
+      t19 = x*t60;
+      t18 = x*t55;
+      t16 = t40*t66;
+      t12 = RATIONAL(9.0,4.0)*t61;
+      t11 = x*t58;
+      t10 = t39*t69;
+      t9 = y*t59;
+      t8 = t39*t59;
+      t6 = y*t56;
+      t4 = y*t57;
+      coeff_I_m1_m1 = t20+t10+t9+t8+t52+(t49+t48)*t37;
+      coeff_I_0_m1 = t38+t6+t33*t53+t65+t67*(t32*t41+t37);
+      coeff_I_p1_m1 = t19+t4-t49+(t30*t40+t29*x)*t38+t50;
+      coeff_I_p2_m1 = t29*t53+t16+t36*t49+t64;
+      coeff_I_m1_0 = t11+t41+t33*t54+t65+t68*(t51+t37);
+      coeff_I_0_0 = t12+RATIONAL(1.0,1.0)+(t35+t22*t39)*t41+(t39+t40)*t33+(t35+
+t22*t40+RATIONAL(25.0,4.0)*t41)*t38;
+      coeff_I_p1_0 = t18+(t31+t26*t39)*t41+(t30+t23*t38)*t40+t47*x+t63;
+      coeff_I_p2_0 = (t37+t58)*t41+t47*t40+t62;
+      coeff_I_m1_p1 = t16+t18+(t30*t39+t29*y)*t41-t48+t50;
+      coeff_I_0_p1 = t4+(t23*t41+t30)*t39+(t26*t40+t31)*t38+t46*y+t63;
+      coeff_I_p1_p1 = t20+t11+t6+t12+(t24*t39+y)*t41+(RATIONAL(4.0,1.0)*t41+x+
+t24*t40)*t38;
+      coeff_I_p2_p1 = t7+t9+t53+(t55+t66-t38)*t41;
+      coeff_I_m1_p2 = t29*t54+t19+t36*t48+t64;
+      coeff_I_0_p2 = t46*t39+(t56+t37)*t38+t62;
+      coeff_I_p1_p2 = t7+t54+t10+(t57+t28*x-t41)*t38;
+      coeff_I_p2_p2 = t8+t36*t52+(t54+t53)*t28;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dx.compute.c
new file mode 100644
index 0000000..a4abe98
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dx.compute.c
@@ -0,0 +1,125 @@
+fp t111,
+   t112,
+   t130,
+   t129,
+   t128,
+   t99,
+   t113,
+   t119,
+   t84,
+   t93,
+   t127,
+   t96,
+   t75,
+   t98,
+   t126,
+   t110,
+   t125,
+   t124,
+   t123,
+   t122,
+   t100,
+   t121,
+   t106,
+   t120,
+   t95,
+   t101,
+   t118,
+   t117,
+   t102,
+   t116,
+   t103,
+   t109,
+   t115,
+   t108,
+   t107,
+   t105,
+   t104,
+   t97,
+   t94,
+   t92,
+   t90,
+   t89,
+   t88,
+   t87,
+   t86,
+   t83,
+   t82,
+   t81,
+   t78,
+   t77,
+   t76,
+   t72;
+      t111 = y*y;
+      t112 = t111*y;
+      t130 = t112+y;
+      t129 = x*t111;
+      t128 = x*t112;
+      t99 = RATIONAL(9.0,4.0);
+      t113 = x*x;
+      t119 = t112*t113;
+      t84 = t99*t119;
+      t93 = RATIONAL(5.0,2.0);
+      t127 = t84+t93*t129;
+      t96 = RATIONAL(-9.0,4.0);
+      t75 = t96*t119;
+      t98 = RATIONAL(-5.0,1.0);
+      t126 = t75+t98*t129;
+      t110 = RATIONAL(-2.0,1.0);
+      t125 = t75+t110*t129;
+      t124 = RATIONAL(-10.0,1.0)*t129+RATIONAL(-27.0,4.0)*t119;
+      t123 = t99*t113;
+      t122 = t111*t113;
+      t100 = RATIONAL(3.0,4.0);
+      t121 = t100*t113;
+      t106 = RATIONAL(-3.0,4.0);
+      t120 = t106*t113;
+      t95 = RATIONAL(4.0,1.0);
+      t101 = RATIONAL(-1.0,4.0);
+      t118 = t84+t101*y+t95*t129;
+      t117 = x*t130;
+      t102 = RATIONAL(1.0,2.0);
+      t116 = RATIONAL(-9.0,2.0)*t113+t102;
+      t103 = RATIONAL(-3.0,2.0);
+      t109 = RATIONAL(-1.0,2.0);
+      t115 = t109+t103*t113;
+      t108 = RATIONAL(1.0,4.0);
+      t107 = RATIONAL(-5.0,2.0);
+      t105 = RATIONAL(3.0,2.0);
+      t104 = RATIONAL(2.0,1.0);
+      t97 = RATIONAL(3.0,1.0);
+      t94 = RATIONAL(-3.0,1.0);
+      t92 = RATIONAL(9.0,2.0);
+      t90 = t100*t112;
+      t89 = t106*t112;
+      t88 = t108*t112;
+      t87 = t108*y;
+      t86 = t101*t112;
+      t83 = y*t123;
+      t82 = t106*t119;
+      t81 = y*t120;
+      t78 = t100*t119;
+      t77 = y*t121;
+      t76 = t96*t113*y;
+      t72 = RATIONAL(27.0,4.0)*t119;
+      coeff_dx_m1_m1 = t88+t87+t78+t77+t115*t111+(t104*t111-t130)*x;
+      coeff_dx_0_m1 = t76+t92*t122+t93*t117+t126;
+      coeff_dx_p1_m1 = t86+t83+t116*t111+t110*t117+t118;
+      coeff_dx_p2_m1 = t82+t105*t122+t81+(-t111+t130*t102)*x;
+      coeff_dx_m1_0 = t89+(RATIONAL(5.0,4.0)+RATIONAL(15.0,4.0)*t113)*t111+(t97
+*t112+t104)*x+t115+t126;
+      coeff_dx_0_0 = t72+(t92+RATIONAL(-45.0,4.0)*t111)*t113+(RATIONAL(25.0,2.0
+)*t111+t98+RATIONAL(-15.0,2.0)*t112)*x;
+      coeff_dx_p1_0 = t90+(RATIONAL(45.0,4.0)*t113+RATIONAL(-5.0,4.0))*t111+(
+t95+RATIONAL(6.0,1.0)*t112)*x+t116+t124;
+      coeff_dx_p2_0 = -x+t103*t128+(RATIONAL(-15.0,4.0)*t111+t105)*t113+t127;
+      coeff_dx_m1_p1 = t94*t122+t90-t111+t81+(y+t94*t112)*x+t118;
+      coeff_dx_0_p1 = t83+RATIONAL(9.0,1.0)*t122+(RATIONAL(15.0,2.0)*t112+t107*
+y)*x+t124;
+      coeff_dx_p1_p1 = t89+t76+RATIONAL(-9.0,1.0)*t122+t111+t87+t72+(RATIONAL(
+8.0,1.0)*t111+RATIONAL(-6.0,1.0)*t112+t104*y)*x;
+      coeff_dx_p2_p1 = t97*t122+t77+(t109*y+t105*t112)*x+t125;
+      coeff_dx_m1_p2 = t86+t128+t82+(t121+t108-x)*t111;
+      coeff_dx_0_p2 = t107*t128+t96*t122+t127;
+      coeff_dx_p1_p2 = t88+t104*t128+(t123+t101)*t111+t125;
+      coeff_dx_p2_p2 = t111*t120+t78+(t109*t112+t102*t111)*x;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dxx.compute.c
new file mode 100644
index 0000000..05e8adc
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dxx.compute.c
@@ -0,0 +1,104 @@
+fp t240,
+   t215,
+   t239,
+   t225,
+   t226,
+   t238,
+   t237,
+   t236,
+   t207,
+   t199,
+   t224,
+   t235,
+   t234,
+   t211,
+   t208,
+   t198,
+   t233,
+   t217,
+   t232,
+   t231,
+   t213,
+   t230,
+   t229,
+   t209,
+   t228,
+   t223,
+   t222,
+   t221,
+   t220,
+   t219,
+   t218,
+   t216,
+   t214,
+   t212,
+   t210,
+   t206,
+   t205,
+   t197,
+   t196,
+   t194,
+   t193,
+   t192,
+   t191;
+      t240 = x*y;
+      t215 = RATIONAL(-3.0,2.0);
+      t239 = t215*x;
+      t225 = y*y;
+      t226 = t225*y;
+      t238 = x*t226;
+      t237 = y+t226;
+      t236 = x*t225;
+      t207 = RATIONAL(-9.0,2.0);
+      t199 = t207*t238;
+      t224 = RATIONAL(-2.0,1.0);
+      t235 = t199+t224*t225;
+      t234 = RATIONAL(-10.0,1.0)*t225+RATIONAL(-27.0,2.0)*t238;
+      t211 = RATIONAL(4.0,1.0);
+      t208 = RATIONAL(9.0,2.0);
+      t198 = t208*t238;
+      t233 = t198+t211*t225;
+      t217 = RATIONAL(3.0,2.0);
+      t232 = t217*t226;
+      t231 = t215*t226;
+      t213 = RATIONAL(-5.0,1.0);
+      t230 = t199+t213*t225;
+      t229 = -t225+x*t231;
+      t209 = RATIONAL(5.0,2.0);
+      t228 = t198+t209*t225;
+      t223 = RATIONAL(6.0,1.0);
+      t222 = RATIONAL(9.0,1.0);
+      t221 = RATIONAL(-6.0,1.0);
+      t220 = RATIONAL(-1.0,2.0);
+      t219 = RATIONAL(-9.0,1.0);
+      t218 = RATIONAL(-5.0,2.0);
+      t216 = RATIONAL(2.0,1.0);
+      t214 = RATIONAL(1.0,2.0);
+      t212 = RATIONAL(3.0,1.0);
+      t210 = RATIONAL(-3.0,1.0);
+      t206 = RATIONAL(-15.0,2.0);
+      t205 = RATIONAL(15.0,2.0);
+      t197 = t208*t240;
+      t196 = y*t239;
+      t194 = t207*t240;
+      t193 = x*t232;
+      t192 = t217*t240;
+      t191 = RATIONAL(27.0,2.0)*t238;
+      coeff_dxx_m1_m1 = t192+t193+(t216+t210*x)*t225-t237;
+      coeff_dxx_0_m1 = t222*t236+t194+t237*t209+t230;
+      coeff_dxx_p1_m1 = t219*t236+t197+t237*t224+t233;
+      coeff_dxx_p2_m1 = t196+t212*t236+t237*t214+t229;
+      coeff_dxx_m1_0 = t216+t212*t226+(t210+t205*t225)*x+t230;
+      coeff_dxx_0_0 = RATIONAL(25.0,2.0)*t225+t191+t213+t206*t226+(t222+
+RATIONAL(-45.0,2.0)*t225)*x;
+      coeff_dxx_p1_0 = t223*t226+t211+(t219+RATIONAL(45.0,2.0)*t225)*x+t234;
+      coeff_dxx_p2_0 = t231+RATIONAL(-1.0,1.0)+(t206*t225+t212)*x+t228;
+      coeff_dxx_m1_p1 = y+t221*t236+t196+t210*t226+t233;
+      coeff_dxx_0_p1 = t218*y+t205*t226+t197+RATIONAL(18.0,1.0)*t236+t234;
+      coeff_dxx_p1_p1 = t191+t194+t221*t226+t216*y+(RATIONAL(-18.0,1.0)*x+
+RATIONAL(8.0,1.0))*t225;
+      coeff_dxx_p2_p1 = t223*t236+t232+t192+t220*y+t235;
+      coeff_dxx_m1_p2 = t217*t236+t226+t229;
+      coeff_dxx_0_p2 = t218*t226+t207*t236+t228;
+      coeff_dxx_p1_p2 = t208*t236+t216*t226+t235;
+      coeff_dxx_p2_p2 = t193+t220*t226+(t214+t239)*t225;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dxy.compute.c
new file mode 100644
index 0000000..6df54ad
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dxy.compute.c
@@ -0,0 +1,120 @@
+fp t300,
+   t289,
+   t299,
+   t290,
+   t298,
+   t294,
+   t241,
+   t297,
+   t296,
+   t242,
+   t295,
+   t293,
+   t268,
+   t288,
+   t292,
+   t291,
+   t287,
+   t286,
+   t285,
+   t284,
+   t283,
+   t282,
+   t281,
+   t280,
+   t279,
+   t278,
+   t277,
+   t276,
+   t275,
+   t274,
+   t273,
+   t272,
+   t271,
+   t270,
+   t269,
+   t266,
+   t265,
+   t264,
+   t263,
+   t262,
+   t261,
+   t260,
+   t259,
+   t258,
+   t257,
+   t256,
+   t255,
+   t254,
+   t253,
+   t249,
+   t246;
+      t300 = x*y;
+      t289 = y*y;
+      t299 = x*t289;
+      t290 = x*x;
+      t298 = t290*y;
+      t294 = t290*t289;
+      t241 = RATIONAL(-27.0,4.0)*t294;
+      t297 = t241+RATIONAL(-10.0,1.0)*t300;
+      t296 = RATIONAL(-81.0,4.0)*t294+RATIONAL(-20.0,1.0)*t300;
+      t242 = RATIONAL(27.0,4.0)*t294;
+      t295 = t242+RATIONAL(5.0,1.0)*t300;
+      t293 = t241+RATIONAL(-4.0,1.0)*t300;
+      t268 = RATIONAL(-9.0,4.0);
+      t288 = RATIONAL(-2.0,1.0);
+      t292 = t288*t300+t268*t294;
+      t291 = RATIONAL(-1.0,4.0)+t242+RATIONAL(8.0,1.0)*t300;
+      t287 = RATIONAL(6.0,1.0);
+      t286 = RATIONAL(9.0,1.0);
+      t285 = RATIONAL(-6.0,1.0);
+      t284 = RATIONAL(-1.0,2.0);
+      t283 = RATIONAL(-9.0,1.0);
+      t282 = RATIONAL(-5.0,2.0);
+      t281 = RATIONAL(3.0,2.0);
+      t280 = RATIONAL(2.0,1.0);
+      t279 = RATIONAL(-3.0,2.0);
+      t278 = RATIONAL(1.0,2.0);
+      t277 = RATIONAL(3.0,1.0);
+      t276 = RATIONAL(-3.0,1.0);
+      t275 = RATIONAL(9.0,4.0);
+      t274 = RATIONAL(3.0,4.0);
+      t273 = RATIONAL(1.0,4.0);
+      t272 = RATIONAL(5.0,2.0);
+      t271 = RATIONAL(9.0,2.0);
+      t270 = RATIONAL(-9.0,2.0);
+      t269 = RATIONAL(-3.0,4.0);
+      t266 = RATIONAL(-15.0,2.0);
+      t265 = RATIONAL(15.0,2.0);
+      t264 = RATIONAL(-45.0,2.0);
+      t263 = RATIONAL(45.0,2.0);
+      t262 = RATIONAL(18.0,1.0);
+      t261 = RATIONAL(-18.0,1.0);
+      t260 = t268*t289;
+      t259 = t275*t289;
+      t258 = t274*t289;
+      t257 = t274*t290;
+      t256 = t275*t290;
+      t255 = t269*t289;
+      t254 = t268*t290;
+      t253 = t269*t290;
+      t249 = t275*t294;
+      t246 = RATIONAL(81.0,4.0)*t294;
+      coeff_dxy_m1_m1 = t276*t298+t257-y+t273+t258+t249+(RATIONAL(4.0,1.0)*y
+-1.0+t276*t289)*x;
+      coeff_dxy_0_m1 = t254+t286*t298+(t272+t265*t289)*x+t297;
+      coeff_dxy_p1_m1 = t283*t298+t255+t256+y+(t288+t285*t289)*x+t291;
+      coeff_dxy_p2_m1 = t253+t277*t298+(t281*t289+t278)*x+t292;
+      coeff_dxy_m1_0 = t286*t299+t260+(t265*t290+t272)*y+t297;
+      coeff_dxy_0_0 = t246+t264*t298+(t264*t289+RATIONAL(25.0,1.0)*y)*x;
+      coeff_dxy_p1_0 = t262*t299+t259+(t282+t263*t290)*y+t296;
+      coeff_dxy_p2_0 = t270*t299+t266*t298+t295;
+      coeff_dxy_m1_p1 = t259+t253+x+t283*t299+(t288+t285*t290)*y+t291;
+      coeff_dxy_0_p1 = t262*t298+t256+(t263*t289+t282)*x+t296;
+      coeff_dxy_p1_p1 = t254+t246+t260+t273+(t280+t261*t290)*y+(t261*t289+t280+
+RATIONAL(16.0,1.0)*y)*x;
+      coeff_dxy_p2_p1 = t287*t298+t257+(t271*t289+t284)*x+t293;
+      coeff_dxy_m1_p2 = t255+t277*t299+(t278+t281*t290)*y+t292;
+      coeff_dxy_0_p2 = t266*t299+t270*t298+t295;
+      coeff_dxy_p1_p2 = t287*t299+t258+(t284+t271*t290)*y+t293;
+      coeff_dxy_p2_p2 = t249+t279*t298+(t279*t289+y)*x;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dy.compute.c
new file mode 100644
index 0000000..229c17a
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dy.compute.c
@@ -0,0 +1,124 @@
+fp t174,
+   t173,
+   t189,
+   t172,
+   t188,
+   t187,
+   t186,
+   t160,
+   t185,
+   t157,
+   t182,
+   t140,
+   t159,
+   t184,
+   t183,
+   t171,
+   t181,
+   t180,
+   t161,
+   t179,
+   t145,
+   t154,
+   t178,
+   t177,
+   t156,
+   t162,
+   t176,
+   t170,
+   t169,
+   t168,
+   t167,
+   t166,
+   t165,
+   t164,
+   t163,
+   t158,
+   t155,
+   t153,
+   t152,
+   t151,
+   t150,
+   t148,
+   t147,
+   t146,
+   t144,
+   t143,
+   t139,
+   t138,
+   t137,
+   t135,
+   t133;
+      t174 = x*x;
+      t173 = t174*x;
+      t189 = t173*y;
+      t172 = y*y;
+      t188 = x*t172;
+      t187 = t174*y;
+      t186 = x+t173;
+      t160 = RATIONAL(9.0,4.0);
+      t185 = t160*t172;
+      t157 = RATIONAL(-9.0,4.0);
+      t182 = t172*t173;
+      t140 = t157*t182;
+      t159 = RATIONAL(-5.0,1.0);
+      t184 = t140+t159*t187;
+      t183 = RATIONAL(-27.0,4.0)*t182+RATIONAL(-10.0,1.0)*t187;
+      t171 = RATIONAL(-2.0,1.0);
+      t181 = t171*t187+t140;
+      t180 = t172*t174;
+      t161 = RATIONAL(3.0,4.0);
+      t179 = t161*t172;
+      t145 = t160*t182;
+      t154 = RATIONAL(5.0,2.0);
+      t178 = t154*t187+t145;
+      t177 = y*t186;
+      t156 = RATIONAL(4.0,1.0);
+      t162 = RATIONAL(-1.0,4.0);
+      t176 = t156*t187+t145+t162*x;
+      t170 = RATIONAL(-1.0,2.0);
+      t169 = RATIONAL(1.0,4.0);
+      t168 = RATIONAL(-5.0,2.0);
+      t167 = RATIONAL(-3.0,4.0);
+      t166 = RATIONAL(3.0,2.0);
+      t165 = RATIONAL(2.0,1.0);
+      t164 = RATIONAL(-3.0,2.0);
+      t163 = RATIONAL(1.0,2.0);
+      t158 = RATIONAL(3.0,1.0);
+      t155 = RATIONAL(-3.0,1.0);
+      t153 = RATIONAL(9.0,2.0);
+      t152 = RATIONAL(-9.0,2.0);
+      t151 = t167*t173;
+      t150 = t161*t173;
+      t148 = t162*t173;
+      t147 = t169*x;
+      t146 = t169*t173;
+      t144 = x*t179;
+      t143 = t167*t188;
+      t139 = t167*t182;
+      t138 = t157*t188;
+      t137 = x*t185;
+      t135 = t173*t179;
+      t133 = RATIONAL(27.0,4.0)*t182;
+      coeff_dy_m1_m1 = t144+t135+t147+t146+(t164*t172+t170)*t174+(t165*t174-
+t186)*y;
+      coeff_dy_0_m1 = RATIONAL(5.0,4.0)*t174+t151+t170+(t164+RATIONAL(15.0,4.0)
+*t174)*t172+(t158*t173+t165)*y+t184;
+      coeff_dy_p1_m1 = t155*t180+t143+t150-t174+(t155*t173+x)*y+t176;
+      coeff_dy_p2_m1 = t148+t189+t139+(t179+t169-y)*t174;
+      coeff_dy_m1_0 = t153*t180+t138+t154*t177+t184;
+      coeff_dy_0_0 = t133+(RATIONAL(-45.0,4.0)*t174+t153)*t172+(t159+RATIONAL(
+-15.0,2.0)*t173+RATIONAL(25.0,2.0)*t174)*y;
+      coeff_dy_p1_0 = t137+RATIONAL(9.0,1.0)*t180+(RATIONAL(15.0,2.0)*t173+t168
+*x)*y+t183;
+      coeff_dy_p2_0 = t168*t189+t157*t180+t178;
+      coeff_dy_m1_p1 = t148+t137+(t152*t172+t163)*t174+t171*t177+t176;
+      coeff_dy_0_p1 = t163+t150+RATIONAL(-5.0,4.0)*t174+(t152+RATIONAL(45.0,4.0
+)*t174)*t172+(RATIONAL(6.0,1.0)*t173+t156)*y+t183;
+      coeff_dy_p1_p1 = t151+RATIONAL(-9.0,1.0)*t180+t133+t147+t174+t138+(
+RATIONAL(-6.0,1.0)*t173+RATIONAL(8.0,1.0)*t174+t165*x)*y;
+      coeff_dy_p2_p1 = t146+t165*t189+(t185+t162)*t174+t181;
+      coeff_dy_m1_p2 = t166*t180+t139+t143+(-t174+t186*t163)*y;
+      coeff_dy_0_p2 = t164*t189-y+(t166+RATIONAL(-15.0,4.0)*t174)*t172+t178;
+      coeff_dy_p1_p2 = t158*t180+t144+(t170*x+t166*t173)*y+t181;
+      coeff_dy_p2_p2 = t135+t167*t180+(t170*t173+t163*t174)*y;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dyy.compute.c
new file mode 100644
index 0000000..c1a4e13
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-dyy.compute.c
@@ -0,0 +1,106 @@
+fp t352,
+   t337,
+   t336,
+   t351,
+   t350,
+   t325,
+   t349,
+   t318,
+   t348,
+   t347,
+   t346,
+   t335,
+   t305,
+   t345,
+   t344,
+   t327,
+   t343,
+   t322,
+   t319,
+   t306,
+   t342,
+   t320,
+   t341,
+   t334,
+   t340,
+   t339,
+   t333,
+   t332,
+   t331,
+   t330,
+   t329,
+   t328,
+   t326,
+   t324,
+   t323,
+   t321,
+   t317,
+   t316,
+   t310,
+   t308,
+   t307,
+   t304,
+   t303,
+   t302;
+      t352 = x*y;
+      t337 = x*x;
+      t336 = t337*x;
+      t351 = x+t336;
+      t350 = y*t337;
+      t325 = RATIONAL(-3.0,2.0);
+      t349 = t325*y;
+      t318 = RATIONAL(-9.0,2.0);
+      t348 = y*t318;
+      t347 = y*t336;
+      t346 = RATIONAL(-10.0,1.0)*t337+RATIONAL(-27.0,2.0)*t347;
+      t335 = RATIONAL(-2.0,1.0);
+      t305 = t318*t347;
+      t345 = t335*t337+t305;
+      t344 = t325*t336;
+      t327 = RATIONAL(3.0,2.0);
+      t343 = t327*t336;
+      t322 = RATIONAL(4.0,1.0);
+      t319 = RATIONAL(9.0,2.0);
+      t306 = t319*t347;
+      t342 = t322*t337+t306;
+      t320 = RATIONAL(5.0,2.0);
+      t341 = t320*t337+t306;
+      t334 = RATIONAL(-5.0,1.0);
+      t340 = t305+t334*t337;
+      t339 = y*t344-t337;
+      t333 = RATIONAL(6.0,1.0);
+      t332 = RATIONAL(9.0,1.0);
+      t331 = RATIONAL(-6.0,1.0);
+      t330 = RATIONAL(-1.0,2.0);
+      t329 = RATIONAL(-9.0,1.0);
+      t328 = RATIONAL(-5.0,2.0);
+      t326 = RATIONAL(2.0,1.0);
+      t324 = RATIONAL(1.0,2.0);
+      t323 = RATIONAL(3.0,1.0);
+      t321 = RATIONAL(-3.0,1.0);
+      t317 = RATIONAL(-15.0,2.0);
+      t316 = RATIONAL(15.0,2.0);
+      t310 = x*t349;
+      t308 = x*t348;
+      t307 = t319*t352;
+      t304 = y*t343;
+      t303 = t327*t352;
+      t302 = RATIONAL(27.0,2.0)*t347;
+      coeff_dyy_m1_m1 = t303+t304+(t321*y+t326)*t337-t351;
+      coeff_dyy_0_m1 = t326+t323*t336+(t321+t316*t337)*y+t340;
+      coeff_dyy_p1_m1 = x+t310+t331*t350+t321*t336+t342;
+      coeff_dyy_p2_m1 = t327*t350+t336+t339;
+      coeff_dyy_m1_0 = t308+t332*t350+t351*t320+t340;
+      coeff_dyy_0_0 = t317*t336+t302+t334+RATIONAL(25.0,2.0)*t337+(t332+
+RATIONAL(-45.0,2.0)*t337)*y;
+      coeff_dyy_p1_0 = t328*x+t316*t336+RATIONAL(18.0,1.0)*t350+t307+t346;
+      coeff_dyy_p2_0 = t337*t348+t328*t336+t341;
+      coeff_dyy_m1_p1 = t329*t350+t307+t351*t335+t342;
+      coeff_dyy_0_p1 = t333*t336+t322+(RATIONAL(45.0,2.0)*t337+t329)*y+t346;
+      coeff_dyy_p1_p1 = t302+t331*t336+t308+t326*x+(RATIONAL(8.0,1.0)+RATIONAL(
+-18.0,1.0)*y)*t337;
+      coeff_dyy_p2_p1 = t319*t350+t326*t336+t345;
+      coeff_dyy_m1_p2 = t323*t350+t310+t351*t324+t339;
+      coeff_dyy_0_p2 = t344+RATIONAL(-1.0,1.0)+(t323+t317*t337)*y+t341;
+      coeff_dyy_p1_p2 = t303+t330*x+t333*t350+t343+t345;
+      coeff_dyy_p2_p2 = t304+t330*t336+(t349+t324)*t337;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-I.compute.c
new file mode 100644
index 0000000..c3036a7
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-I.compute.c
@@ -0,0 +1,341 @@
+fp t470,
+   t471,
+   t519,
+   t518,
+   t473,
+   t472,
+   t517,
+   t516,
+   t449,
+   t515,
+   t425,
+   t514,
+   t414,
+   t513,
+   t512,
+   t511,
+   t411,
+   t510,
+   t504,
+   t506,
+   t509,
+   t410,
+   t494,
+   t360,
+   t435,
+   t508,
+   t396,
+   t507,
+   t426,
+   t486,
+   t505,
+   t462,
+   t503,
+   t451,
+   t488,
+   t394,
+   t448,
+   t502,
+   t416,
+   t501,
+   t500,
+   t420,
+   t489,
+   t491,
+   t408,
+   t499,
+   t424,
+   t492,
+   t369,
+   t498,
+   t497,
+   t496,
+   t495,
+   t421,
+   t493,
+   t490,
+   t406,
+   t487,
+   t372,
+   t485,
+   t457,
+   t484,
+   t483,
+   t422,
+   t482,
+   t392,
+   t481,
+   t480,
+   t479,
+   t478,
+   t434,
+   t436,
+   t477,
+   t469,
+   t468,
+   t467,
+   t466,
+   t465,
+   t464,
+   t463,
+   t461,
+   t460,
+   t459,
+   t458,
+   t456,
+   t455,
+   t454,
+   t453,
+   t452,
+   t450,
+   t447,
+   t446,
+   t445,
+   t444,
+   t443,
+   t442,
+   t441,
+   t440,
+   t439,
+   t438,
+   t437,
+   t433,
+   t432,
+   t431,
+   t430,
+   t429,
+   t428,
+   t427,
+   t423,
+   t419,
+   t418,
+   t417,
+   t415,
+   t413,
+   t412,
+   t409,
+   t407,
+   t405,
+   t402,
+   t401,
+   t400,
+   t398,
+   t397,
+   t393,
+   t390,
+   t387,
+   t384,
+   t383,
+   t382,
+   t378,
+   t377,
+   t374,
+   t373,
+   t368,
+   t365,
+   t363,
+   t362,
+   t361,
+   t359,
+   t358,
+   t357,
+   t356,
+   t355,
+   t353;
+      t470 = y*y;
+      t471 = y*t470;
+      t519 = t471+y;
+      t518 = x*y;
+      t473 = x*x;
+      t472 = t473*x;
+      t517 = t472+x;
+      t516 = t472*y;
+      t449 = RATIONAL(1.0,18.0);
+      t515 = t449*y;
+      t425 = RATIONAL(-1.0,18.0);
+      t514 = t425*y;
+      t414 = RATIONAL(-1.0,144.0);
+      t513 = t414*x;
+      t512 = t449*x;
+      t511 = x*t471;
+      t411 = RATIONAL(1.0,144.0);
+      t510 = x*t411;
+      t504 = t470*t473;
+      t506 = t472*t471;
+      t509 = RATIONAL(-16.0,9.0)*t506+RATIONAL(-35.0,9.0)*t504;
+      t410 = RATIONAL(-7.0,144.0);
+      t494 = t410*t471;
+      t360 = t472*t494;
+      t435 = RATIONAL(-1.0,24.0);
+      t508 = t435*t504+t360;
+      t396 = RATIONAL(-7.0,9.0)*t506;
+      t507 = t396+RATIONAL(-35.0,12.0)*t504;
+      t426 = RATIONAL(-1.0,72.0);
+      t486 = t414*t472;
+      t505 = t426*t504+t471*t486;
+      t462 = RATIONAL(-1.0,9.0);
+      t503 = t462*t471;
+      t451 = RATIONAL(1.0,9.0);
+      t488 = t451*t472;
+      t394 = t471*t488;
+      t448 = RATIONAL(5.0,36.0);
+      t502 = t448*t504+t394;
+      t416 = RATIONAL(7.0,144.0);
+      t501 = t416*t471;
+      t500 = t451*t471;
+      t420 = RATIONAL(-7.0,36.0);
+      t489 = t420*t473;
+      t491 = t462*t472;
+      t408 = t471*t491;
+      t499 = t408+t470*t489;
+      t424 = RATIONAL(5.0,48.0);
+      t492 = t416*t472;
+      t369 = t471*t492;
+      t498 = t369+t424*t504;
+      t497 = t473*t471;
+      t496 = t472*t470;
+      t495 = t410*t472;
+      t421 = RATIONAL(7.0,18.0);
+      t493 = t421*t504+t394;
+      t490 = t472*t411;
+      t406 = RATIONAL(7.0,9.0)*t506;
+      t487 = RATIONAL(7.0,6.0)*t504+t406;
+      t372 = x*t514;
+      t485 = t360+t372+RATIONAL(-5.0,24.0)*t504;
+      t457 = RATIONAL(-5.0,6.0);
+      t484 = t457*t504+t396+t372;
+      t483 = RATIONAL(-4.0,9.0)*t518+RATIONAL(25.0,12.0)*t504+t406;
+      t422 = RATIONAL(1.0,12.0);
+      t482 = t422*t504+y*t513+t369;
+      t392 = y*t512;
+      t481 = RATIONAL(-5.0,18.0)*t504+t392+t408;
+      t480 = t392+RATIONAL(-49.0,144.0)*t506+RATIONAL(-5.0,8.0)*t504;
+      t479 = t517*t470;
+      t478 = t519*t473;
+      t434 = RATIONAL(-1.0,12.0);
+      t436 = RATIONAL(7.0,36.0);
+      t477 = t434+t436*t473;
+      t469 = RATIONAL(5.0,4.0);
+      t468 = RATIONAL(-2.0,9.0);
+      t467 = RATIONAL(-8.0,9.0);
+      t466 = RATIONAL(2.0,9.0);
+      t465 = RATIONAL(-2.0,3.0);
+      t464 = RATIONAL(-1.0,3.0);
+      t463 = RATIONAL(-5.0,3.0);
+      t461 = RATIONAL(5.0,6.0);
+      t460 = RATIONAL(8.0,9.0);
+      t459 = RATIONAL(1.0,3.0);
+      t458 = RATIONAL(-1.0,2.0);
+      t456 = RATIONAL(4.0,3.0);
+      t455 = RATIONAL(-7.0,3.0);
+      t454 = RATIONAL(-4.0,3.0);
+      t453 = RATIONAL(2.0,3.0);
+      t452 = RATIONAL(5.0,3.0);
+      t450 = RATIONAL(-1.0,6.0);
+      t447 = RATIONAL(-10.0,9.0);
+      t446 = RATIONAL(-7.0,18.0);
+      t445 = RATIONAL(7.0,24.0);
+      t444 = RATIONAL(-5.0,48.0);
+      t443 = RATIONAL(-5.0,36.0);
+      t442 = RATIONAL(-7.0,72.0);
+      t441 = RATIONAL(-28.0,9.0);
+      t440 = RATIONAL(20.0,9.0);
+      t439 = RATIONAL(-14.0,9.0);
+      t438 = RATIONAL(1.0,72.0);
+      t437 = RATIONAL(28.0,9.0);
+      t433 = RATIONAL(7.0,12.0);
+      t432 = RATIONAL(-7.0,12.0);
+      t431 = RATIONAL(14.0,9.0);
+      t430 = RATIONAL(10.0,9.0);
+      t429 = RATIONAL(-20.0,9.0);
+      t428 = RATIONAL(-7.0,24.0);
+      t427 = RATIONAL(1.0,24.0);
+      t423 = RATIONAL(7.0,72.0);
+      t419 = RATIONAL(-35.0,36.0);
+      t418 = RATIONAL(35.0,48.0);
+      t417 = RATIONAL(35.0,36.0);
+      t415 = RATIONAL(-49.0,36.0);
+      t413 = RATIONAL(-35.0,48.0);
+      t412 = RATIONAL(49.0,36.0);
+      t409 = y*t491;
+      t407 = t467*t516;
+      t405 = t460*t516;
+      t402 = x*t500;
+      t401 = t467*t511;
+      t400 = t460*t511;
+      t398 = x*t503;
+      t397 = RATIONAL(4.0,9.0)*t518;
+      t393 = y*t488;
+      t390 = t449*t511;
+      t387 = t446*t511;
+      t384 = t446*t516;
+      t383 = t472*t515;
+      t382 = RATIONAL(16.0,9.0)*t506;
+      t378 = t472*t514;
+      t377 = t421*t511;
+      t374 = t421*t516;
+      t373 = t425*t511;
+      t368 = y*t492;
+      t365 = t414*t511;
+      t363 = x*t501;
+      t362 = y*t486;
+      t361 = x*t494;
+      t359 = y*t495;
+      t358 = t471*t490;
+      t357 = y*t490;
+      t356 = t471*t510;
+      t355 = y*t510;
+      t353 = RATIONAL(49.0,144.0)*t506;
+      coeff_I_m2_m2 = t355+t356+t357+t358+RATIONAL(1.0,36.0)*t504+(t478+t479)*
+t426;
+      coeff_I_m1_m2 = t359+t373+t424*t478+(t423*t472+t451*x)*t470+t485;
+      coeff_I_0_m2 = t393+(t450+t468*t472)*t470+t493+t519*(t422+t489);
+      coeff_I_p1_m2 = t390+t409+t448*t478+(t466*t472+t462*x)*t470+t481;
+      coeff_I_p2_m2 = t365+t368+t435*t478+(t442*t472+t438*x)*t470+t482;
+      coeff_I_p3_m2 = t362+t438*t496+t411*t478+t505;
+      coeff_I_m2_m1 = t361+t378+(t423*t471+t451*y)*t473+t424*t479+t485;
+      coeff_I_m1_m1 = t377+t397+t353+t374+(t457*y+t413*t471)*t473+(t457*x+
+RATIONAL(25.0,16.0)*t473+t413*t472)*t470;
+      coeff_I_0_m1 = t407+(t412*t473+t432)*t471+(t431*t473+t465)*y+(t452*t472+
+t469)*t470+t507;
+      coeff_I_p1_m1 = t405+t387+(t419*t471+t447*y)*t473+(t461*x+t463*t472)*t470
++t483;
+      coeff_I_p2_m1 = t363+t384+(t445*t471+t459*y)*t473+(t418*t472+t444*x)*t470
++t480;
+      coeff_I_p3_m1 = t444*t496+t383+(t514+t494)*t473+t498;
+      coeff_I_m2_0 = t402+t517*t422+(t450+t468*t471)*t473+t420*t479+t493;
+      coeff_I_m1_0 = t432*t472+t401+t465*x+(t452*t471+t469)*t473+(t412*t472+
+t431*x)*t470+t507;
+      coeff_I_0_0 = t455*t473+t456*t472+RATIONAL(1.0,1.0)+t382+(t456+t441*t473)
+*t471+(RATIONAL(49.0,9.0)*t473+t455+t441*t472)*t470;
+      coeff_I_p1_0 = t453*x+t400+t454*t472+(t440*t471+t452)*t473+(t439*x+t437*
+t472)*t470+t509;
+      coeff_I_p2_0 = t434*x+t398+t433*t472+(t465*t471+t458)*t473+(t415*t472+
+t436*x)*t470+t487;
+      coeff_I_p3_0 = (t422+t500)*t473+(t434+t436*t470)*t472+t499;
+      coeff_I_m2_p1 = t398+t383+(t466*t471+t462*y)*t473+t448*t479+t481;
+      coeff_I_m1_p1 = t384+t400+(t461*y+t463*t471)*t473+(t447*x+t419*t472)*t470
++t483;
+      coeff_I_0_p1 = t405+(t454+t437*t473)*t471+(t453+t439*t473)*y+(t452+t440*
+t472)*t470+t509;
+      coeff_I_p1_p1 = t382+t407+t397+t401+(t430*y+t429*t471)*t473+(t430*x+t429*
+t472+RATIONAL(25.0,9.0)*t473)*t470;
+      coeff_I_p2_p1 = t374+t402+(t464*y+t453*t471)*t473+(t417*t472+t443*x)*t470
++t484;
+      coeff_I_p3_p1 = t443*t496+t378+(t515+t503)*t473+t502;
+      coeff_I_m2_p2 = t363+t362+(t442*t471+t438*y)*t473+t435*t479+t482;
+      coeff_I_m1_p2 = t387+t368+(t418*t471+t444*y)*t473+(t459*x+t445*t472)*t470
++t480;
+      coeff_I_0_p2 = t409+(t415*t473+t433)*t471+t477*y+(t458+t465*t472)*t470+
+t487;
+      coeff_I_p1_p2 = t377+t393+(t417*t471+t443*y)*t473+(t464*x+t453*t472)*t470
++t484;
+      coeff_I_p2_p2 = t359+t355+t353+t361+(t428*t471+t427*y)*t473+(t427*x+t428*
+t472+RATIONAL(1.0,4.0)*t473)*t470;
+      coeff_I_p3_p2 = t427*t496+t357+(t414*y+t501)*t473+t508;
+      coeff_I_m2_p3 = t365+t438*t497+t411*t479+t505;
+      coeff_I_m1_p3 = t390+t444*t497+(t425*x+t495)*t470+t498;
+      coeff_I_0_p3 = t477*t471+(t488+t422)*t470+t499;
+      coeff_I_p1_p3 = t373+t443*t497+(t512+t491)*t470+t502;
+      coeff_I_p2_p3 = t427*t497+t356+(t513+t492)*t470+t508;
+      coeff_I_p3_p3 = t358+t411*t504+(t496+t497)*t414;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dx.compute.c
new file mode 100644
index 0000000..034e400
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dx.compute.c
@@ -0,0 +1,275 @@
+fp t620,
+   t649,
+   t618,
+   t648,
+   t619,
+   t647,
+   t646,
+   t598,
+   t639,
+   t535,
+   t645,
+   t593,
+   t644,
+   t610,
+   t640,
+   t553,
+   t583,
+   t643,
+   t642,
+   t641,
+   t638,
+   t637,
+   t636,
+   t587,
+   t566,
+   t546,
+   t635,
+   t634,
+   t584,
+   t633,
+   t602,
+   t614,
+   t631,
+   t556,
+   t632,
+   t605,
+   t545,
+   t630,
+   t582,
+   t530,
+   t629,
+   t628,
+   t611,
+   t627,
+   t604,
+   t568,
+   t626,
+   t591,
+   t570,
+   t586,
+   t625,
+   t606,
+   t578,
+   t624,
+   t613,
+   t623,
+   t622,
+   t617,
+   t616,
+   t615,
+   t612,
+   t609,
+   t608,
+   t607,
+   t603,
+   t601,
+   t600,
+   t599,
+   t597,
+   t596,
+   t595,
+   t594,
+   t592,
+   t590,
+   t589,
+   t585,
+   t579,
+   t577,
+   t576,
+   t575,
+   t573,
+   t572,
+   t571,
+   t569,
+   t567,
+   t565,
+   t564,
+   t563,
+   t561,
+   t560,
+   t559,
+   t555,
+   t554,
+   t551,
+   t548,
+   t547,
+   t544,
+   t543,
+   t541,
+   t537,
+   t536,
+   t531,
+   t529,
+   t528,
+   t524,
+   t521;
+      t620 = x*x;
+      t649 = y*t620;
+      t618 = y*y;
+      t648 = x*t618;
+      t619 = t618*y;
+      t647 = x*t619;
+      t646 = t619+y;
+      t598 = RATIONAL(-1.0,12.0);
+      t639 = RATIONAL(-7.0,48.0)*t620;
+      t535 = t619*t639;
+      t645 = t535+t598*t648;
+      t593 = RATIONAL(1.0,48.0);
+      t644 = t593*t620;
+      t610 = RATIONAL(1.0,3.0);
+      t640 = t610*t620;
+      t553 = t619*t640;
+      t583 = RATIONAL(5.0,18.0);
+      t643 = t553+t583*t648;
+      t642 = t620*t619;
+      t641 = t618*t620;
+      t638 = RATIONAL(-16.0,3.0)*t642+RATIONAL(-70.0,9.0)*t648;
+      t637 = RATIONAL(7.0,9.0)*t648+t553;
+      t636 = RATIONAL(-1.0,48.0)*t620;
+      t587 = RATIONAL(7.0,18.0);
+      t566 = t587*t619;
+      t546 = RATIONAL(-7.0,3.0)*t642;
+      t635 = t546+RATIONAL(-35.0,6.0)*t648;
+      t634 = RATIONAL(7.0,48.0)*t620;
+      t584 = RATIONAL(-1.0,36.0);
+      t633 = t619*t636+t584*t648;
+      t602 = RATIONAL(-7.0,18.0);
+      t614 = RATIONAL(-1.0,3.0);
+      t631 = t614*t620;
+      t556 = t619*t631;
+      t632 = t556+t602*t648;
+      t605 = RATIONAL(7.0,3.0);
+      t545 = t605*t642;
+      t630 = t605*t648+t545;
+      t582 = RATIONAL(5.0,24.0);
+      t530 = t619*t634;
+      t629 = t582*t648+t530;
+      t628 = x*t646;
+      t611 = RATIONAL(-4.0,9.0);
+      t627 = RATIONAL(25.0,6.0)*t648+t545+t611*y;
+      t604 = RATIONAL(1.0,18.0);
+      t568 = t604*y;
+      t626 = t556+RATIONAL(-5.0,9.0)*t648+t568;
+      t591 = RATIONAL(-1.0,18.0);
+      t570 = t591*y;
+      t586 = RATIONAL(-5.0,12.0);
+      t625 = t570+t535+t586*t648;
+      t606 = RATIONAL(1.0,6.0);
+      t578 = RATIONAL(-1.0,144.0);
+      t624 = t578*y+t530+t606*t648;
+      t613 = RATIONAL(-5.0,3.0);
+      t623 = t613*t648+t546+t570;
+      t622 = RATIONAL(-49.0,48.0)*t642+RATIONAL(-5.0,4.0)*t648+t568;
+      t617 = RATIONAL(-2.0,9.0);
+      t616 = RATIONAL(2.0,9.0);
+      t615 = RATIONAL(-2.0,3.0);
+      t612 = RATIONAL(-1.0,9.0);
+      t609 = RATIONAL(2.0,3.0);
+      t608 = RATIONAL(4.0,9.0);
+      t607 = RATIONAL(1.0,9.0);
+      t603 = RATIONAL(-5.0,18.0);
+      t601 = RATIONAL(1.0,72.0);
+      t600 = RATIONAL(7.0,36.0);
+      t599 = RATIONAL(-1.0,24.0);
+      t597 = RATIONAL(7.0,12.0);
+      t596 = RATIONAL(-7.0,12.0);
+      t595 = RATIONAL(1.0,24.0);
+      t594 = RATIONAL(-1.0,72.0);
+      t592 = RATIONAL(1.0,36.0);
+      t590 = RATIONAL(-5.0,24.0);
+      t589 = RATIONAL(1.0,12.0);
+      t585 = RATIONAL(-7.0,36.0);
+      t579 = RATIONAL(10.0,3.0);
+      t577 = RATIONAL(1.0,144.0);
+      t576 = t612*t619;
+      t575 = t608*y;
+      t573 = t607*t619;
+      t572 = RATIONAL(-8.0,9.0)*t619;
+      t571 = RATIONAL(8.0,9.0)*t619;
+      t569 = t591*t619;
+      t567 = t604*t619;
+      t565 = t602*t619;
+      t564 = t577*t619;
+      t563 = t577*y;
+      t561 = t578*t619;
+      t560 = RATIONAL(-7.0,144.0)*t619;
+      t559 = RATIONAL(7.0,144.0)*t619;
+      t555 = y*t631;
+      t554 = y*t640;
+      t551 = RATIONAL(-1.0,6.0)*t649;
+      t548 = RATIONAL(7.0,6.0)*t649;
+      t547 = RATIONAL(-8.0,3.0)*t649;
+      t544 = RATIONAL(8.0,3.0)*t649;
+      t543 = RATIONAL(-7.0,6.0)*t649;
+      t541 = t606*t649;
+      t537 = y*t636;
+      t536 = y*t639;
+      t531 = y*t634;
+      t529 = t593*t642;
+      t528 = y*t644;
+      t524 = RATIONAL(16.0,3.0)*t642;
+      t521 = RATIONAL(49.0,48.0)*t642;
+      coeff_dx_m2_m2 = t529+t564+t528+t563+(t599*t620+t594)*t618+(t604*t618+
+t646*t584)*x;
+      coeff_dx_m1_m2 = t569+t536+(t607+RATIONAL(7.0,24.0)*t620)*t618+t582*t628+
+t625;
+      coeff_dx_0_m2 = t615*t641+t554+t602*t628+t637;
+      coeff_dx_p1_m2 = t555+t567+(t612+t609*t620)*t618+t583*t628+t626;
+      coeff_dx_p2_m2 = t561+t531+(t601+RATIONAL(-7.0,24.0)*t620)*t618+t598*t628
++t624;
+      coeff_dx_p3_m2 = t537+t595*t641+t601*t628+t633;
+      coeff_dx_m2_m1 = t551+t560+(RATIONAL(5.0,48.0)+RATIONAL(5.0,16.0)*t620)*
+t618+(t600*t619+t616*y)*x+t625;
+      coeff_dx_m1_m1 = t548+t575+t566+t521+(RATIONAL(-35.0,16.0)*t620+RATIONAL(
+-5.0,6.0))*t618+(t613*y+RATIONAL(25.0,8.0)*t618+RATIONAL(-35.0,24.0)*t619)*x;
+      coeff_dx_0_m1 = t547+RATIONAL(5.0,1.0)*t641+(RATIONAL(28.0,9.0)*y+
+RATIONAL(49.0,18.0)*t619)*x+t635;
+      coeff_dx_p1_m1 = t565+t544+(RATIONAL(-5.0,1.0)*t620+RATIONAL(5.0,6.0))*
+t618+(RATIONAL(-35.0,18.0)*t619+RATIONAL(-20.0,9.0)*y)*x+t627;
+      coeff_dx_p2_m1 = t543+t559+(RATIONAL(35.0,16.0)*t620+RATIONAL(-5.0,48.0))
+*t618+(t597*t619+t609*y)*x+t622;
+      coeff_dx_p3_m1 = t541+RATIONAL(-5.0,16.0)*t641+(RATIONAL(-7.0,72.0)*t619+
+t612*y)*x+t629;
+      coeff_dx_m2_0 = t589+t573+RATIONAL(1.0,4.0)*t620+(t585+t596*t620)*t618+(
+t614+t611*t619)*x+t637;
+      coeff_dx_m1_0 = RATIONAL(-7.0,4.0)*t620+t615+t572+(RATIONAL(14.0,9.0)+
+RATIONAL(49.0,12.0)*t620)*t618+(t579*t619+RATIONAL(5.0,2.0))*x+t635;
+      coeff_dx_0_0 = t524+(RATIONAL(-28.0,3.0)*t618+RATIONAL(4.0,1.0))*t620+(
+RATIONAL(-14.0,3.0)+RATIONAL(98.0,9.0)*t618+RATIONAL(-56.0,9.0)*t619)*x;
+      coeff_dx_p1_0 = RATIONAL(-4.0,1.0)*t620+t571+t609+(RATIONAL(28.0,3.0)*
+t620+RATIONAL(-14.0,9.0))*t618+(t579+RATIONAL(40.0,9.0)*t619)*x+t638;
+      coeff_dx_p2_0 = t598+RATIONAL(-4.0,3.0)*t647-x+t576+RATIONAL(7.0,4.0)*
+t620+(RATIONAL(-49.0,12.0)*t620+t600)*t618+t630;
+      coeff_dx_p3_0 = (t597*t618+RATIONAL(-1.0,4.0))*t620+(t616*t619+t606)*x+
+t632;
+      coeff_dx_m2_p1 = t576+t541+(RATIONAL(5.0,12.0)*t620+RATIONAL(5.0,36.0))*
+t618+(t617*y+t608*t619)*x+t626;
+      coeff_dx_m1_p1 = t571+t543+(RATIONAL(-10.0,9.0)+RATIONAL(-35.0,12.0)*t620
+)*t618+(RATIONAL(-10.0,3.0)*t619+RATIONAL(5.0,3.0)*y)*x+t627;
+      coeff_dx_0_p1 = t544+RATIONAL(20.0,3.0)*t641+(RATIONAL(56.0,9.0)*t619+
+RATIONAL(-28.0,9.0)*y)*x+t638;
+      coeff_dx_p1_p1 = t575+t547+t572+t524+(RATIONAL(-20.0,3.0)*t620+RATIONAL(
+10.0,9.0))*t618+(RATIONAL(50.0,9.0)*t618+RATIONAL(20.0,9.0)*y+RATIONAL(-40.0,
+9.0)*t619)*x;
+      coeff_dx_p2_p1 = t548+t573+(RATIONAL(-5.0,36.0)+RATIONAL(35.0,12.0)*t620)
+*t618+(RATIONAL(4.0,3.0)*t619+t615*y)*x+t623;
+      coeff_dx_p3_p1 = t551+t586*t641+(t607*y+t617*t619)*x+t643;
+      coeff_dx_m2_p2 = t537+t559+(RATIONAL(-1.0,8.0)*t620+t599)*t618+(t592*y+
+t585*t619)*x+t624;
+      coeff_dx_m1_p2 = t565+t531+(RATIONAL(7.0,8.0)*t620+t610)*t618+(RATIONAL(
+35.0,24.0)*t619+t590*y)*x+t622;
+      coeff_dx_0_p2 = t555+RATIONAL(-2.0,1.0)*t641+(RATIONAL(-49.0,18.0)*t619+
+t587*y)*x+t630;
+      coeff_dx_p1_p2 = t554+t566+(RATIONAL(2.0,1.0)*t620+t614)*t618+(t603*y+
+RATIONAL(35.0,18.0)*t619)*x+t623;
+      coeff_dx_p2_p2 = t563+t521+t536+t560+(RATIONAL(-7.0,8.0)*t620+t595)*t618+
+(t596*t619+t589*y+RATIONAL(1.0,2.0)*t618)*x;
+      coeff_dx_p3_p2 = t528+RATIONAL(1.0,8.0)*t641+(t594*y+RATIONAL(7.0,72.0)*
+t619)*x+t645;
+      coeff_dx_m2_p3 = t592*t647+t561+(t644+t577)*t618+t633;
+      coeff_dx_m1_p3 = t590*t647+t567+(t639+t591)*t618+t629;
+      coeff_dx_0_p3 = x*t566+t618*t640+t632;
+      coeff_dx_p1_p3 = t603*t647+t569+(t604+t631)*t618+t643;
+      coeff_dx_p2_p3 = t564+t589*t647+(t634+t578)*t618+t645;
+      coeff_dx_p3_p3 = t529+t618*t636+(t594*t619+t601*t618)*x;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dxx.compute.c
new file mode 100644
index 0000000..6ee2c84
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dxx.compute.c
@@ -0,0 +1,218 @@
+fp t878,
+   t835,
+   t877,
+   t851,
+   t876,
+   t842,
+   t875,
+   t852,
+   t874,
+   t837,
+   t873,
+   t848,
+   t872,
+   t830,
+   t871,
+   t831,
+   t870,
+   t869,
+   t841,
+   t794,
+   t868,
+   t799,
+   t867,
+   t822,
+   t866,
+   t789,
+   t840,
+   t865,
+   t864,
+   t824,
+   t863,
+   t805,
+   t796,
+   t862,
+   t861,
+   t783,
+   t860,
+   t819,
+   t791,
+   t859,
+   t838,
+   t858,
+   t823,
+   t857,
+   t834,
+   t856,
+   t846,
+   t855,
+   t854,
+   t850,
+   t849,
+   t847,
+   t845,
+   t844,
+   t843,
+   t839,
+   t836,
+   t833,
+   t832,
+   t829,
+   t828,
+   t827,
+   t826,
+   t825,
+   t820,
+   t803,
+   t802,
+   t801,
+   t800,
+   t798,
+   t797,
+   t795,
+   t793,
+   t790,
+   t788,
+   t787,
+   t786,
+   t785,
+   t784,
+   t781;
+      t878 = x*y;
+      t835 = RATIONAL(-1.0,24.0);
+      t877 = t835*x;
+      t851 = y*y;
+      t876 = x*t851;
+      t842 = RATIONAL(2.0,3.0);
+      t875 = t842*y;
+      t852 = t851*y;
+      t874 = t852+y;
+      t837 = RATIONAL(7.0,24.0);
+      t873 = x*t837;
+      t848 = RATIONAL(-2.0,3.0);
+      t872 = x*t848;
+      t830 = RATIONAL(1.0,24.0);
+      t871 = x*t830;
+      t831 = RATIONAL(-7.0,24.0);
+      t870 = x*t831;
+      t869 = x*t852;
+      t841 = RATIONAL(1.0,6.0);
+      t794 = t837*t869;
+      t868 = t794+t841*t851;
+      t799 = t842*t869;
+      t867 = t799+RATIONAL(7.0,9.0)*t851;
+      t822 = RATIONAL(5.0,24.0);
+      t866 = t794+t822*t851;
+      t789 = RATIONAL(14.0,3.0)*t869;
+      t840 = RATIONAL(7.0,3.0);
+      t865 = t789+t840*t851;
+      t864 = RATIONAL(-49.0,24.0)*t869+RATIONAL(-5.0,4.0)*t851;
+      t824 = RATIONAL(-1.0,36.0);
+      t863 = t835*t869+t824*t851;
+      t805 = RATIONAL(-35.0,6.0)*t851;
+      t796 = t848*t869;
+      t862 = t796+RATIONAL(-5.0,9.0)*t851;
+      t861 = RATIONAL(-32.0,3.0)*t869+RATIONAL(-70.0,9.0)*t851;
+      t783 = t831*t869;
+      t860 = t783+RATIONAL(-5.0,12.0)*t851;
+      t819 = RATIONAL(-14.0,3.0);
+      t791 = t819*t869;
+      t859 = t791+t805;
+      t838 = RATIONAL(-7.0,18.0);
+      t858 = t796+t838*t851;
+      t823 = RATIONAL(5.0,18.0);
+      t857 = t799+t823*t851;
+      t834 = RATIONAL(-1.0,12.0);
+      t856 = t783+t834*t851;
+      t846 = RATIONAL(-5.0,3.0);
+      t855 = t791+t846*t851;
+      t854 = t789+RATIONAL(25.0,6.0)*t851;
+      t850 = RATIONAL(-2.0,9.0);
+      t849 = RATIONAL(2.0,9.0);
+      t847 = RATIONAL(-1.0,3.0);
+      t845 = RATIONAL(1.0,2.0);
+      t844 = RATIONAL(4.0,3.0);
+      t843 = RATIONAL(-4.0,3.0);
+      t839 = RATIONAL(-5.0,18.0);
+      t836 = RATIONAL(1.0,72.0);
+      t833 = RATIONAL(7.0,12.0);
+      t832 = RATIONAL(-7.0,12.0);
+      t829 = RATIONAL(-1.0,72.0);
+      t828 = RATIONAL(1.0,36.0);
+      t827 = RATIONAL(-5.0,24.0);
+      t826 = RATIONAL(1.0,12.0);
+      t825 = RATIONAL(7.0,18.0);
+      t820 = RATIONAL(10.0,3.0);
+      t803 = t847*t878;
+      t802 = t840*t878;
+      t801 = RATIONAL(-7.0,3.0)*t878;
+      t800 = RATIONAL(1.0,3.0)*t878;
+      t798 = x*t875;
+      t797 = y*t872;
+      t795 = y*t873;
+      t793 = y*t877;
+      t790 = RATIONAL(-16.0,3.0)*t878;
+      t788 = RATIONAL(16.0,3.0)*t878;
+      t787 = RATIONAL(32.0,3.0)*t869;
+      t786 = t830*t869;
+      t785 = y*t870;
+      t784 = y*t871;
+      t781 = RATIONAL(49.0,24.0)*t869;
+      coeff_dxx_m2_m2 = t786+t784+t874*t824+(RATIONAL(1.0,18.0)+t834*x)*t851;
+      coeff_dxx_m1_m2 = t833*t876+t785+t874*t822+t860;
+      coeff_dxx_0_m2 = t798+t843*t876+t874*t838+t867;
+      coeff_dxx_p1_m2 = t844*t876+t797+t874*t823+t862;
+      coeff_dxx_p2_m2 = t832*t876+t795+t874*t834+t868;
+      coeff_dxx_p3_m2 = t826*t876+t793+t874*t836+t863;
+      coeff_dxx_m2_m1 = t803+RATIONAL(5.0,8.0)*t876+t849*y+RATIONAL(7.0,36.0)*
+t852+t860;
+      coeff_dxx_m1_m1 = t802+t846*y+RATIONAL(-35.0,24.0)*t852+t781+(RATIONAL(
+25.0,8.0)+RATIONAL(-35.0,8.0)*x)*t851;
+      coeff_dxx_0_m1 = RATIONAL(49.0,18.0)*t852+t790+RATIONAL(28.0,9.0)*y+
+RATIONAL(10.0,1.0)*t876+t859;
+      coeff_dxx_p1_m1 = RATIONAL(-35.0,18.0)*t852+RATIONAL(-20.0,9.0)*y+
+RATIONAL(-10.0,1.0)*t876+t788+t854;
+      coeff_dxx_p2_m1 = t801+t833*t852+t875+RATIONAL(35.0,8.0)*t876+t864;
+      coeff_dxx_p3_m1 = t800+RATIONAL(-7.0,72.0)*t852+RATIONAL(-1.0,9.0)*y+
+RATIONAL(-5.0,8.0)*t876+t866;
+      coeff_dxx_m2_0 = t847+RATIONAL(-4.0,9.0)*t852+(t845+RATIONAL(-7.0,6.0)*
+t851)*x+t867;
+      coeff_dxx_m1_0 = t820*t852+RATIONAL(5.0,2.0)+(RATIONAL(49.0,6.0)*t851+
+RATIONAL(-7.0,2.0))*x+t859;
+      coeff_dxx_0_0 = t819+t787+RATIONAL(98.0,9.0)*t851+RATIONAL(-56.0,9.0)*
+t852+(RATIONAL(8.0,1.0)+RATIONAL(-56.0,3.0)*t851)*x;
+      coeff_dxx_p1_0 = t820+RATIONAL(40.0,9.0)*t852+(RATIONAL(-8.0,1.0)+
+RATIONAL(56.0,3.0)*t851)*x+t861;
+      coeff_dxx_p2_0 = RATIONAL(-1.0,1.0)+t843*t852+(RATIONAL(-49.0,6.0)*t851+
+RATIONAL(7.0,2.0))*x+t865;
+      coeff_dxx_p3_0 = t841+t849*t852+(RATIONAL(-1.0,2.0)+RATIONAL(7.0,6.0)*
+t851)*x+t858;
+      coeff_dxx_m2_p1 = t800+RATIONAL(4.0,9.0)*t852+RATIONAL(5.0,6.0)*t876+t850
+*y+t862;
+      coeff_dxx_m1_p1 = t801+RATIONAL(-10.0,3.0)*t852+RATIONAL(5.0,3.0)*y+x*
+t805+t854;
+      coeff_dxx_0_p1 = RATIONAL(-28.0,9.0)*y+RATIONAL(40.0,3.0)*t876+t788+
+RATIONAL(56.0,9.0)*t852+t861;
+      coeff_dxx_p1_p1 = RATIONAL(20.0,9.0)*y+t790+RATIONAL(-40.0,9.0)*t852+t787
++(RATIONAL(-40.0,3.0)*x+RATIONAL(50.0,9.0))*t851;
+      coeff_dxx_p2_p1 = RATIONAL(35.0,6.0)*t876+t802+t844*t852+t848*y+t855;
+      coeff_dxx_p3_p1 = t803+RATIONAL(-5.0,6.0)*t876+RATIONAL(1.0,9.0)*y+t850*
+t852+t857;
+      coeff_dxx_m2_p2 = t793+RATIONAL(-1.0,4.0)*t876+RATIONAL(-7.0,36.0)*t852+
+t828*y+t868;
+      coeff_dxx_m1_p2 = t827*y+t795+RATIONAL(7.0,4.0)*t876+RATIONAL(35.0,24.0)*
+t852+t864;
+      coeff_dxx_0_p2 = RATIONAL(-4.0,1.0)*t876+t825*y+RATIONAL(-49.0,18.0)*t852
++t797+t865;
+      coeff_dxx_p1_p2 = t839*y+t798+RATIONAL(35.0,18.0)*t852+RATIONAL(4.0,1.0)*
+t876+t855;
+      coeff_dxx_p2_p2 = t832*t852+t826*y+t785+t781+(RATIONAL(-7.0,4.0)*x+t845)*
+t851;
+      coeff_dxx_p3_p2 = t784+RATIONAL(1.0,4.0)*t876+RATIONAL(7.0,72.0)*t852+
+t829*y+t856;
+      coeff_dxx_m2_p3 = t851*t871+t828*t852+t863;
+      coeff_dxx_m1_p3 = t851*t870+t827*t852+t866;
+      coeff_dxx_0_p3 = t825*t852+t842*t876+t858;
+      coeff_dxx_p1_p3 = t839*t852+t851*t872+t857;
+      coeff_dxx_p2_p3 = t826*t852+t851*t873+t856;
+      coeff_dxx_p3_p3 = t829*t852+t786+(t877+t836)*t851;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dxy.compute.c
new file mode 100644
index 0000000..84daf78
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dxy.compute.c
@@ -0,0 +1,296 @@
+fp t1020,
+   t1001,
+   t1019,
+   t1002,
+   t1018,
+   t1000,
+   t889,
+   t977,
+   t1017,
+   t898,
+   t1016,
+   t888,
+   t1015,
+   t942,
+   t1014,
+   t1013,
+   t897,
+   t1012,
+   t980,
+   t918,
+   t974,
+   t915,
+   t917,
+   t913,
+   t1011,
+   t1010,
+   t1009,
+   t969,
+   t1008,
+   t986,
+   t1007,
+   t989,
+   t1006,
+   t1005,
+   t1004,
+   t1003,
+   t999,
+   t998,
+   t997,
+   t996,
+   t995,
+   t994,
+   t993,
+   t992,
+   t990,
+   t988,
+   t987,
+   t985,
+   t984,
+   t983,
+   t982,
+   t981,
+   t979,
+   t978,
+   t976,
+   t975,
+   t973,
+   t972,
+   t971,
+   t970,
+   t968,
+   t967,
+   t966,
+   t965,
+   t964,
+   t963,
+   t962,
+   t961,
+   t960,
+   t959,
+   t958,
+   t957,
+   t956,
+   t955,
+   t954,
+   t953,
+   t952,
+   t951,
+   t950,
+   t949,
+   t948,
+   t947,
+   t946,
+   t945,
+   t944,
+   t943,
+   t941,
+   t940,
+   t939,
+   t938,
+   t937,
+   t936,
+   t935,
+   t934,
+   t933,
+   t932,
+   t931,
+   t930,
+   t929,
+   t927,
+   t926,
+   t925,
+   t924,
+   t923,
+   t922,
+   t921,
+   t920,
+   t919,
+   t916,
+   t914,
+   t912,
+   t911,
+   t910,
+   t909,
+   t908,
+   t907,
+   t906,
+   t905,
+   t904,
+   t895,
+   t885,
+   t880;
+      t1020 = x*y;
+      t1001 = y*y;
+      t1019 = x*t1001;
+      t1002 = x*x;
+      t1018 = t1002*y;
+      t1000 = t1002*t1001;
+      t889 = RATIONAL(-7.0,16.0)*t1000;
+      t977 = RATIONAL(-1.0,6.0);
+      t1017 = t889+t977*t1020;
+      t898 = RATIONAL(-7.0,1.0)*t1000;
+      t1016 = RATIONAL(-35.0,3.0)*t1020+t898;
+      t888 = RATIONAL(7.0,16.0)*t1000;
+      t1015 = t888+RATIONAL(5.0,12.0)*t1020;
+      t942 = RATIONAL(-1.0,18.0);
+      t1014 = t942*t1020+RATIONAL(-1.0,16.0)*t1000;
+      t1013 = RATIONAL(-16.0,1.0)*t1000+RATIONAL(-140.0,9.0)*t1020;
+      t897 = RATIONAL(7.0,1.0)*t1000;
+      t1012 = RATIONAL(14.0,3.0)*t1020+t897;
+      t980 = RATIONAL(7.0,6.0);
+      t918 = t980*t1001;
+      t974 = RATIONAL(-7.0,6.0);
+      t915 = t974*t1001;
+      t917 = t980*t1002;
+      t913 = t974*t1002;
+      t1011 = -t1000+RATIONAL(-7.0,9.0)*t1020;
+      t1010 = RATIONAL(14.0,9.0)*t1020+t1000;
+      t1009 = t1000+RATIONAL(5.0,9.0)*t1020;
+      t969 = RATIONAL(1.0,18.0);
+      t1008 = RATIONAL(-5.0,2.0)*t1020+t969+RATIONAL(-49.0,16.0)*t1000;
+      t986 = RATIONAL(-5.0,6.0);
+      t1007 = t889+t942+t986*t1020;
+      t989 = RATIONAL(1.0,3.0);
+      t1006 = RATIONAL(-1.0,144.0)+t888+t989*t1020;
+      t1005 = RATIONAL(25.0,3.0)*t1020+RATIONAL(-4.0,9.0)+t897;
+      t1004 = t942+RATIONAL(-10.0,3.0)*t1020+t898;
+      t1003 = -t1000+RATIONAL(-10.0,9.0)*t1020+t969;
+      t999 = RATIONAL(-2.0,9.0);
+      t998 = RATIONAL(1.0,4.0);
+      t997 = RATIONAL(2.0,9.0);
+      t996 = RATIONAL(-2.0,3.0);
+      t995 = RATIONAL(-1.0,3.0);
+      t994 = RATIONAL(-5.0,3.0);
+      t993 = RATIONAL(-1.0,9.0);
+      t992 = RATIONAL(5.0,6.0);
+      t990 = RATIONAL(7.0,4.0);
+      t988 = RATIONAL(5.0,8.0);
+      t987 = RATIONAL(-1.0,4.0);
+      t985 = RATIONAL(4.0,3.0);
+      t984 = RATIONAL(-5.0,8.0);
+      t983 = RATIONAL(-4.0,3.0);
+      t982 = RATIONAL(2.0,3.0);
+      t981 = RATIONAL(4.0,9.0);
+      t979 = RATIONAL(5.0,3.0);
+      t978 = RATIONAL(1.0,9.0);
+      t976 = RATIONAL(1.0,6.0);
+      t975 = RATIONAL(-8.0,3.0);
+      t973 = RATIONAL(8.0,3.0);
+      t972 = RATIONAL(-7.0,4.0);
+      t971 = RATIONAL(4.0,1.0);
+      t970 = RATIONAL(-4.0,1.0);
+      t968 = RATIONAL(10.0,1.0);
+      t967 = RATIONAL(-49.0,6.0);
+      t966 = RATIONAL(-10.0,1.0);
+      t965 = RATIONAL(-5.0,18.0);
+      t964 = RATIONAL(-7.0,18.0);
+      t963 = RATIONAL(7.0,24.0);
+      t962 = RATIONAL(40.0,3.0);
+      t961 = RATIONAL(35.0,8.0);
+      t960 = RATIONAL(35.0,6.0);
+      t959 = RATIONAL(-40.0,3.0);
+      t958 = RATIONAL(-28.0,9.0);
+      t957 = RATIONAL(20.0,9.0);
+      t956 = RATIONAL(1.0,72.0);
+      t955 = RATIONAL(28.0,9.0);
+      t954 = RATIONAL(-1.0,24.0);
+      t953 = RATIONAL(-1.0,12.0);
+      t952 = RATIONAL(7.0,12.0);
+      t951 = RATIONAL(49.0,6.0);
+      t950 = RATIONAL(-7.0,12.0);
+      t949 = RATIONAL(-35.0,8.0);
+      t948 = RATIONAL(-20.0,9.0);
+      t947 = RATIONAL(-56.0,3.0);
+      t946 = RATIONAL(-7.0,24.0);
+      t945 = RATIONAL(1.0,24.0);
+      t944 = RATIONAL(1.0,48.0);
+      t943 = RATIONAL(-1.0,72.0);
+      t941 = RATIONAL(1.0,36.0);
+      t940 = RATIONAL(-7.0,48.0);
+      t939 = RATIONAL(-5.0,24.0);
+      t938 = RATIONAL(1.0,12.0);
+      t937 = RATIONAL(7.0,18.0);
+      t936 = RATIONAL(-1.0,36.0);
+      t935 = RATIONAL(-1.0,48.0);
+      t934 = RATIONAL(5.0,18.0);
+      t933 = RATIONAL(5.0,24.0);
+      t932 = RATIONAL(7.0,48.0);
+      t931 = RATIONAL(56.0,3.0);
+      t930 = RATIONAL(-35.0,6.0);
+      t929 = RATIONAL(1.0,144.0);
+      t927 = t989*t1001;
+      t926 = t975*t1001;
+      t925 = t973*t1001;
+      t924 = t995*t1001;
+      t923 = t977*t1001;
+      t922 = t995*t1002;
+      t921 = t976*t1001;
+      t920 = t989*t1002;
+      t919 = t977*t1002;
+      t916 = t975*t1002;
+      t914 = t973*t1002;
+      t912 = t976*t1002;
+      t911 = t944*t1002;
+      t910 = t944*t1001;
+      t909 = t940*t1002;
+      t908 = t935*t1001;
+      t907 = t932*t1002;
+      t906 = t940*t1001;
+      t905 = t935*t1002;
+      t904 = t932*t1001;
+      t895 = RATIONAL(16.0,1.0)*t1000;
+      t885 = RATIONAL(1.0,16.0)*t1000;
+      t880 = RATIONAL(49.0,16.0)*t1000;
+      coeff_dxy_m2_m2 = t885+t929+t911+t910+(t953*t1002+t936)*y+(t978*y+t953*
+t1001+t936)*x;
+      coeff_dxy_m1_m2 = t923+t909+(t952*t1002+t997)*y+(t933+t988*t1001)*x+t1007
+;
+      coeff_dxy_0_m2 = t983*t1018+t920+(t915+t964)*x+t1010;
+      coeff_dxy_p1_m2 = t921+t922+(t999+t985*t1002)*y+(t934+t992*t1001)*x+t1003
+;
+      coeff_dxy_p2_m2 = t907+t908+(t950*t1002+t941)*y+(t953+t987*t1001)*x+t1006
+;
+      coeff_dxy_p3_m2 = t938*t1018+t905+(t956+t945*t1001)*x+t1014;
+      coeff_dxy_m2_m1 = t919+t906+(t988*t1002+t933)*y+(t952*t1001+t997)*x+t1007
+;
+      coeff_dxy_m1_m1 = t981+t880+t918+t917+(t949*t1002+t994)*y+(RATIONAL(25.0,
+4.0)*y+t994+t949*t1001)*x;
+      coeff_dxy_0_m1 = t916+t968*t1018+(t951*t1001+t955)*x+t1016;
+      coeff_dxy_p1_m1 = t915+t914+(t979+t966*t1002)*y+(t930*t1001+t948)*x+t1005
+;
+      coeff_dxy_p2_m1 = t913+t904+(t939+t961*t1002)*y+(t990*t1001+t982)*x+t1008
+;
+      coeff_dxy_p3_m1 = t984*t1018+t912+(t993+t946*t1001)*x+t1015;
+      coeff_dxy_m2_0 = t927+t983*t1019+(t913+t964)*y+t1010;
+      coeff_dxy_m1_0 = t968*t1019+t926+(t951*t1002+t955)*y+t1016;
+      coeff_dxy_0_0 = t947*t1018+t895+(t947*t1001+RATIONAL(196.0,9.0)*y)*x;
+      coeff_dxy_p1_0 = t962*t1019+t925+(t931*t1002+t958)*y+t1013;
+      coeff_dxy_p2_0 = t924+t970*t1019+(t937+t967*t1002)*y+t1012;
+      coeff_dxy_p3_0 = y*t917+t982*t1019+t1011;
+      coeff_dxy_m2_p1 = t924+t912+(t992*t1002+t934)*y+(t999+t985*t1001)*x+t1003
+;
+      coeff_dxy_m1_p1 = t913+t925+(t948+t930*t1002)*y+(t966*t1001+t979)*x+t1005
+;
+      coeff_dxy_0_p1 = t962*t1018+t914+(t931*t1001+t958)*x+t1013;
+      coeff_dxy_p1_p1 = t895+t916+t926+t981+(t957+t959*t1002)*y+(t957+RATIONAL(
+100.0,9.0)*y+t959*t1001)*x;
+      coeff_dxy_p2_p1 = t917+t927+(t960*t1002+t965)*y+(t971*t1001+t996)*x+t1004
+;
+      coeff_dxy_p3_p1 = t986*t1018+t919+(t996*t1001+t978)*x+t1009;
+      coeff_dxy_m2_p2 = t904+t905+(t987*t1002+t953)*y+(t950*t1001+t941)*x+t1006
+;
+      coeff_dxy_m1_p2 = t907+t915+(t982+t990*t1002)*y+(t939+t961*t1001)*x+t1008
+;
+      coeff_dxy_0_p2 = t970*t1018+t922+(t967*t1001+t937)*x+t1012;
+      coeff_dxy_p1_p2 = t918+t920+(t971*t1002+t996)*y+(t965+t960*t1001)*x+t1004
+;
+      coeff_dxy_p2_p2 = t906+t929+t909+t880+(t938+t972*t1002)*y+(t972*t1001+y+
+t938)*x;
+      coeff_dxy_p3_p2 = t911+t998*t1018+(t963*t1001+t943)*x+t1017;
+      coeff_dxy_m2_p3 = t908+t938*t1019+(t956+t945*t1002)*y+t1014;
+      coeff_dxy_m1_p3 = t984*t1019+t921+(t993+t946*t1002)*y+t1015;
+      coeff_dxy_0_p3 = x*t918+t982*t1018+t1011;
+      coeff_dxy_p1_p3 = t923+t986*t1019+(t996*t1002+t978)*y+t1009;
+      coeff_dxy_p2_p3 = t910+t998*t1019+(t963*t1002+t943)*y+t1017;
+      coeff_dxy_p3_p3 = t885+t954*t1018+(t954*t1001+t941*y)*x;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dy.compute.c
new file mode 100644
index 0000000..31e4695
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dy.compute.c
@@ -0,0 +1,277 @@
+fp t750,
+   t779,
+   t748,
+   t778,
+   t749,
+   t777,
+   t776,
+   t775,
+   t740,
+   t761,
+   t688,
+   t713,
+   t774,
+   t711,
+   t773,
+   t744,
+   t759,
+   t687,
+   t732,
+   t772,
+   t718,
+   t771,
+   t735,
+   t677,
+   t770,
+   t680,
+   t769,
+   t717,
+   t698,
+   t768,
+   t767,
+   t670,
+   t728,
+   t766,
+   t765,
+   t710,
+   t714,
+   t764,
+   t763,
+   t668,
+   t712,
+   t762,
+   t760,
+   t758,
+   t741,
+   t757,
+   t734,
+   t695,
+   t756,
+   t755,
+   t736,
+   t708,
+   t754,
+   t743,
+   t721,
+   t700,
+   t753,
+   t752,
+   t747,
+   t746,
+   t745,
+   t742,
+   t739,
+   t738,
+   t737,
+   t733,
+   t731,
+   t730,
+   t729,
+   t727,
+   t726,
+   t725,
+   t724,
+   t722,
+   t720,
+   t719,
+   t715,
+   t709,
+   t707,
+   t706,
+   t705,
+   t703,
+   t702,
+   t701,
+   t699,
+   t697,
+   t696,
+   t694,
+   t693,
+   t692,
+   t691,
+   t690,
+   t686,
+   t683,
+   t681,
+   t679,
+   t676,
+   t675,
+   t673,
+   t672,
+   t669,
+   t666,
+   t663,
+   t658,
+   t656,
+   t655,
+   t651;
+      t750 = x*x;
+      t779 = y*t750;
+      t748 = y*y;
+      t778 = x*t748;
+      t749 = t750*x;
+      t777 = t749+x;
+      t776 = y*t749;
+      t775 = t748*t750;
+      t740 = RATIONAL(1.0,3.0);
+      t761 = t748*t749;
+      t688 = t740*t761;
+      t713 = RATIONAL(5.0,18.0);
+      t774 = t688+t713*t779;
+      t711 = RATIONAL(7.0,48.0);
+      t773 = t711*t748;
+      t744 = RATIONAL(-1.0,3.0);
+      t759 = t744*t748;
+      t687 = t749*t759;
+      t732 = RATIONAL(-7.0,18.0);
+      t772 = t687+t732*t779;
+      t718 = RATIONAL(-7.0,48.0);
+      t771 = t718*t748;
+      t735 = RATIONAL(7.0,3.0);
+      t677 = t735*t761;
+      t770 = t735*t779+t677;
+      t680 = RATIONAL(-7.0,3.0)*t761;
+      t769 = RATIONAL(-35.0,6.0)*t779+t680;
+      t717 = RATIONAL(7.0,18.0);
+      t698 = t717*t749;
+      t768 = RATIONAL(-16.0,3.0)*t761+RATIONAL(-70.0,9.0)*t779;
+      t767 = t750*RATIONAL(-5.0,12.0);
+      t670 = t718*t761;
+      t728 = RATIONAL(-1.0,12.0);
+      t766 = t728*t779+t670;
+      t765 = t740*t748;
+      t710 = RATIONAL(-1.0,48.0);
+      t714 = RATIONAL(-1.0,36.0);
+      t764 = t714*t779+t710*t761;
+      t763 = RATIONAL(7.0,9.0)*t779+t688;
+      t668 = t711*t761;
+      t712 = RATIONAL(5.0,24.0);
+      t762 = t712*t779+t668;
+      t760 = RATIONAL(1.0,48.0)*t748;
+      t758 = y*t777;
+      t741 = RATIONAL(-4.0,9.0);
+      t757 = RATIONAL(25.0,6.0)*t779+t741*x+t677;
+      t734 = RATIONAL(1.0,18.0);
+      t695 = t734*x;
+      t756 = t695+t687+RATIONAL(-5.0,9.0)*t779;
+      t755 = t695+RATIONAL(-49.0,48.0)*t761+RATIONAL(-5.0,4.0)*t779;
+      t736 = RATIONAL(1.0,6.0);
+      t708 = RATIONAL(-1.0,144.0);
+      t754 = t668+t708*x+t736*t779;
+      t743 = RATIONAL(-5.0,3.0);
+      t721 = RATIONAL(-1.0,18.0);
+      t700 = t721*x;
+      t753 = t700+t680+t743*t779;
+      t752 = t700+y*t767+t670;
+      t747 = RATIONAL(-2.0,9.0);
+      t746 = RATIONAL(2.0,9.0);
+      t745 = RATIONAL(-2.0,3.0);
+      t742 = RATIONAL(-1.0,9.0);
+      t739 = RATIONAL(2.0,3.0);
+      t738 = RATIONAL(4.0,9.0);
+      t737 = RATIONAL(1.0,9.0);
+      t733 = RATIONAL(-5.0,18.0);
+      t731 = RATIONAL(1.0,72.0);
+      t730 = RATIONAL(7.0,36.0);
+      t729 = RATIONAL(-1.0,24.0);
+      t727 = RATIONAL(7.0,12.0);
+      t726 = RATIONAL(-7.0,12.0);
+      t725 = RATIONAL(1.0,24.0);
+      t724 = RATIONAL(-1.0,72.0);
+      t722 = RATIONAL(1.0,36.0);
+      t720 = RATIONAL(-5.0,24.0);
+      t719 = RATIONAL(1.0,12.0);
+      t715 = RATIONAL(-7.0,36.0);
+      t709 = RATIONAL(10.0,3.0);
+      t707 = RATIONAL(1.0,144.0);
+      t706 = t738*x;
+      t705 = RATIONAL(-8.0,9.0)*t749;
+      t703 = RATIONAL(8.0,9.0)*t749;
+      t702 = t737*t749;
+      t701 = t742*t749;
+      t699 = t721*t749;
+      t697 = t732*t749;
+      t696 = t734*t749;
+      t694 = RATIONAL(7.0,144.0)*t749;
+      t693 = RATIONAL(-7.0,144.0)*t749;
+      t692 = t708*t749;
+      t691 = t707*x;
+      t690 = t707*t749;
+      t686 = t736*t778;
+      t683 = RATIONAL(7.0,6.0)*t778;
+      t681 = RATIONAL(-1.0,6.0)*t778;
+      t679 = RATIONAL(-7.0,6.0)*t778;
+      t676 = x*t765;
+      t675 = RATIONAL(-8.0,3.0)*t778;
+      t673 = RATIONAL(8.0,3.0)*t778;
+      t672 = x*t759;
+      t669 = t710*t778;
+      t666 = x*t760;
+      t663 = t749*t760;
+      t658 = x*t773;
+      t656 = x*t771;
+      t655 = RATIONAL(16.0,3.0)*t761;
+      t651 = RATIONAL(49.0,48.0)*t761;
+      coeff_dy_m2_m2 = t690+t666+t663+t691+(t724+t729*t748)*t750+(t734*t750+
+t777*t714)*y;
+      coeff_dy_m1_m2 = t681+t693+(RATIONAL(5.0,48.0)+RATIONAL(5.0,16.0)*t748)*
+t750+(t730*t749+t746*x)*y+t752;
+      coeff_dy_0_m2 = t719+t715*t750+t702+(t726*t750+RATIONAL(1.0,4.0))*t748+(
+t741*t749+t744)*y+t763;
+      coeff_dy_p1_m2 = t686+t701+(RATIONAL(5.0,36.0)+RATIONAL(5.0,12.0)*t748)*
+t750+(t738*t749+t747*x)*y+t756;
+      coeff_dy_p2_m2 = t669+t694+(t729+RATIONAL(-1.0,8.0)*t748)*t750+(t722*x+
+t715*t749)*y+t754;
+      coeff_dy_p3_m2 = t692+t722*t776+(t707+t760)*t750+t764;
+      coeff_dy_m2_m1 = t656+t699+(t737+RATIONAL(7.0,24.0)*t748)*t750+t712*t758+
+t752;
+      coeff_dy_m1_m1 = t698+t683+t706+t651+(RATIONAL(-5.0,6.0)+RATIONAL(-35.0,
+16.0)*t748)*t750+(RATIONAL(25.0,8.0)*t750+RATIONAL(-35.0,24.0)*t749+t743*x)*y;
+      coeff_dy_0_m1 = t745+RATIONAL(14.0,9.0)*t750+t705+(RATIONAL(49.0,12.0)*
+t750+RATIONAL(-7.0,4.0))*t748+(RATIONAL(5.0,2.0)+t709*t749)*y+t769;
+      coeff_dy_p1_m1 = t703+t679+(RATIONAL(-10.0,9.0)+RATIONAL(-35.0,12.0)*t748
+)*t750+(RATIONAL(-10.0,3.0)*t749+RATIONAL(5.0,3.0)*x)*y+t757;
+      coeff_dy_p2_m1 = t697+t658+(t740+RATIONAL(7.0,8.0)*t748)*t750+(RATIONAL(
+35.0,24.0)*t749+t720*x)*y+t755;
+      coeff_dy_p3_m1 = t720*t776+t696+(t721+t771)*t750+t762;
+      coeff_dy_m2_0 = t745*t775+t676+t732*t758+t763;
+      coeff_dy_m1_0 = RATIONAL(5.0,1.0)*t775+t675+(RATIONAL(28.0,9.0)*x+
+RATIONAL(49.0,18.0)*t749)*y+t769;
+      coeff_dy_0_0 = t655+(RATIONAL(-28.0,3.0)*t750+RATIONAL(4.0,1.0))*t748+(
+RATIONAL(-56.0,9.0)*t749+RATIONAL(-14.0,3.0)+RATIONAL(98.0,9.0)*t750)*y;
+      coeff_dy_p1_0 = t673+RATIONAL(20.0,3.0)*t775+(RATIONAL(-28.0,9.0)*x+
+RATIONAL(56.0,9.0)*t749)*y+t768;
+      coeff_dy_p2_0 = t672+RATIONAL(-2.0,1.0)*t775+(t717*x+RATIONAL(-49.0,18.0)
+*t749)*y+t770;
+      coeff_dy_p3_0 = t750*t765+y*t698+t772;
+      coeff_dy_m2_p1 = t672+t696+(t739*t748+t742)*t750+t713*t758+t756;
+      coeff_dy_m1_p1 = t697+t673+(RATIONAL(-5.0,1.0)*t748+RATIONAL(5.0,6.0))*
+t750+(RATIONAL(-35.0,18.0)*t749+RATIONAL(-20.0,9.0)*x)*y+t757;
+      coeff_dy_0_p1 = t739+t703+RATIONAL(-14.0,9.0)*t750+(RATIONAL(28.0,3.0)*
+t750+RATIONAL(-4.0,1.0))*t748+(t709+RATIONAL(40.0,9.0)*t749)*y+t768;
+      coeff_dy_p1_p1 = t675+t705+t655+t706+(RATIONAL(-20.0,3.0)*t748+RATIONAL(
+10.0,9.0))*t750+(RATIONAL(50.0,9.0)*t750+RATIONAL(20.0,9.0)*x+RATIONAL(-40.0,
+9.0)*t749)*y;
+      coeff_dy_p2_p1 = t676+t698+(RATIONAL(2.0,1.0)*t748+t744)*t750+(RATIONAL(
+35.0,18.0)*t749+t733*x)*y+t753;
+      coeff_dy_p3_p1 = t699+t733*t776+(t759+t734)*t750+t774;
+      coeff_dy_m2_p2 = t658+t692+(t731+RATIONAL(-7.0,24.0)*t748)*t750+t728*t758
++t754;
+      coeff_dy_m1_p2 = t679+t694+(RATIONAL(35.0,16.0)*t748+RATIONAL(-5.0,48.0))
+*t750+(t739*x+t727*t749)*y+t755;
+      coeff_dy_0_p2 = RATIONAL(-4.0,3.0)*t776+t701+t728+t730*t750-y+(RATIONAL(
+-49.0,12.0)*t750+RATIONAL(7.0,4.0))*t748+t770;
+      coeff_dy_p1_p2 = t702+t683+(RATIONAL(35.0,12.0)*t748+RATIONAL(-5.0,36.0))
+*t750+(t745*x+RATIONAL(4.0,3.0)*t749)*y+t753;
+      coeff_dy_p2_p2 = t691+t656+t693+t651+(t725+RATIONAL(-7.0,8.0)*t748)*t750+
+(RATIONAL(1.0,2.0)*t750+t719*x+t726*t749)*y;
+      coeff_dy_p3_p2 = t690+t719*t776+(t773+t708)*t750+t766;
+      coeff_dy_m2_p3 = t725*t775+t669+t731*t758+t764;
+      coeff_dy_m1_p3 = t686+RATIONAL(-5.0,16.0)*t775+(RATIONAL(-7.0,72.0)*t749+
+t742*x)*y+t762;
+      coeff_dy_0_p3 = (RATIONAL(-1.0,4.0)+t727*t750)*t748+(t736+t746*t749)*y+
+t772;
+      coeff_dy_p1_p3 = t681+t748*t767+(t737*x+t747*t749)*y+t774;
+      coeff_dy_p2_p3 = RATIONAL(1.0,8.0)*t775+t666+(RATIONAL(7.0,72.0)*t749+
+t724*x)*y+t766;
+      coeff_dy_p3_p3 = t710*t775+t663+(t731*t750+t724*t749)*y;
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dyy.compute.c
new file mode 100644
index 0000000..56fd101
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeff-dyy.compute.c
@@ -0,0 +1,214 @@
+fp t1119,
+   t1093,
+   t1092,
+   t1118,
+   t1117,
+   t1072,
+   t1116,
+   t1115,
+   t1114,
+   t1113,
+   t1089,
+   t1112,
+   t1084,
+   t1111,
+   t1110,
+   t1109,
+   t1032,
+   t1075,
+   t1108,
+   t1041,
+   t1107,
+   t1081,
+   t1033,
+   t1106,
+   t1064,
+   t1105,
+   t1104,
+   t1042,
+   t1079,
+   t1103,
+   t1102,
+   t1046,
+   t1065,
+   t1029,
+   t1101,
+   t1026,
+   t1082,
+   t1100,
+   t1087,
+   t1099,
+   t1098,
+   t1062,
+   t1097,
+   t1061,
+   t1096,
+   t1095,
+   t1091,
+   t1090,
+   t1088,
+   t1086,
+   t1085,
+   t1083,
+   t1080,
+   t1077,
+   t1074,
+   t1073,
+   t1070,
+   t1069,
+   t1068,
+   t1067,
+   t1066,
+   t1063,
+   t1044,
+   t1043,
+   t1040,
+   t1039,
+   t1038,
+   t1037,
+   t1036,
+   t1034,
+   t1031,
+   t1030,
+   t1028,
+   t1025,
+   t1024,
+   t1023,
+   t1022;
+      t1119 = x*y;
+      t1093 = x*x;
+      t1092 = x*t1093;
+      t1118 = t1092+x;
+      t1117 = y*t1093;
+      t1072 = RATIONAL(-7.0,24.0);
+      t1116 = y*t1072;
+      t1115 = y*t1092;
+      t1114 = RATIONAL(-1.0,24.0)*y;
+      t1113 = y*RATIONAL(7.0,24.0);
+      t1089 = RATIONAL(-2.0,3.0);
+      t1112 = t1089*x;
+      t1084 = RATIONAL(2.0,3.0);
+      t1111 = y*t1084;
+      t1110 = y*RATIONAL(1.0,24.0);
+      t1109 = RATIONAL(-70.0,9.0)*t1093+RATIONAL(-32.0,3.0)*t1115;
+      t1032 = t1072*t1115;
+      t1075 = RATIONAL(-1.0,12.0);
+      t1108 = t1032+t1075*t1093;
+      t1041 = t1092*t1111;
+      t1107 = RATIONAL(7.0,9.0)*t1093+t1041;
+      t1081 = RATIONAL(1.0,6.0);
+      t1033 = t1092*t1113;
+      t1106 = t1033+t1081*t1093;
+      t1064 = RATIONAL(-1.0,36.0);
+      t1105 = t1064*t1093+t1092*t1114;
+      t1104 = t1032+RATIONAL(-5.0,12.0)*t1093;
+      t1042 = t1089*t1115;
+      t1079 = RATIONAL(-7.0,18.0);
+      t1103 = t1079*t1093+t1042;
+      t1102 = RATIONAL(-5.0,9.0)*t1093+t1042;
+      t1046 = RATIONAL(-35.0,6.0)*t1093;
+      t1065 = RATIONAL(-14.0,3.0);
+      t1029 = t1065*t1115;
+      t1101 = t1046+t1029;
+      t1026 = RATIONAL(14.0,3.0)*t1115;
+      t1082 = RATIONAL(7.0,3.0);
+      t1100 = t1026+t1082*t1093;
+      t1087 = RATIONAL(-5.0,3.0);
+      t1099 = t1087*t1093+t1029;
+      t1098 = RATIONAL(-5.0,4.0)*t1093+RATIONAL(-49.0,24.0)*t1115;
+      t1062 = RATIONAL(5.0,18.0);
+      t1097 = t1062*t1093+t1041;
+      t1061 = RATIONAL(5.0,24.0);
+      t1096 = t1061*t1093+t1033;
+      t1095 = t1026+RATIONAL(25.0,6.0)*t1093;
+      t1091 = RATIONAL(-2.0,9.0);
+      t1090 = RATIONAL(2.0,9.0);
+      t1088 = RATIONAL(-1.0,3.0);
+      t1086 = RATIONAL(4.0,3.0);
+      t1085 = RATIONAL(-4.0,3.0);
+      t1083 = RATIONAL(1.0,2.0);
+      t1080 = RATIONAL(-5.0,18.0);
+      t1077 = RATIONAL(1.0,72.0);
+      t1074 = RATIONAL(7.0,12.0);
+      t1073 = RATIONAL(-7.0,12.0);
+      t1070 = RATIONAL(-1.0,72.0);
+      t1069 = RATIONAL(1.0,36.0);
+      t1068 = RATIONAL(-5.0,24.0);
+      t1067 = RATIONAL(1.0,12.0);
+      t1066 = RATIONAL(7.0,18.0);
+      t1063 = RATIONAL(10.0,3.0);
+      t1044 = y*t1112;
+      t1043 = RATIONAL(1.0,3.0)*t1119;
+      t1040 = t1088*t1119;
+      t1039 = t1082*t1119;
+      t1038 = RATIONAL(-7.0,3.0)*t1119;
+      t1037 = x*t1111;
+      t1036 = RATIONAL(16.0,3.0)*t1119;
+      t1034 = x*t1114;
+      t1031 = t1092*t1110;
+      t1030 = x*t1110;
+      t1028 = x*t1116;
+      t1025 = x*t1113;
+      t1024 = RATIONAL(-16.0,3.0)*t1119;
+      t1023 = RATIONAL(32.0,3.0)*t1115;
+      t1022 = RATIONAL(49.0,24.0)*t1115;
+      coeff_dyy_m2_m2 = t1031+t1030+(t1075*y+RATIONAL(1.0,18.0))*t1093+t1118*
+t1064;
+      coeff_dyy_m1_m2 = RATIONAL(7.0,36.0)*t1092+t1090*x+RATIONAL(5.0,8.0)*
+t1117+t1040+t1104;
+      coeff_dyy_0_m2 = t1088+RATIONAL(-4.0,9.0)*t1092+(t1083+RATIONAL(-7.0,6.0)
+*t1093)*y+t1107;
+      coeff_dyy_p1_m2 = RATIONAL(4.0,9.0)*t1092+t1043+RATIONAL(5.0,6.0)*t1117+
+t1091*x+t1102;
+      coeff_dyy_p2_m2 = t1069*x+RATIONAL(-1.0,4.0)*t1117+t1034+RATIONAL(-7.0,
+36.0)*t1092+t1106;
+      coeff_dyy_p3_m2 = t1069*t1092+t1093*t1110+t1105;
+      coeff_dyy_m2_m1 = t1074*t1117+t1028+t1118*t1061+t1104;
+      coeff_dyy_m1_m1 = t1022+t1039+RATIONAL(-35.0,24.0)*t1092+t1087*x+(
+RATIONAL(25.0,8.0)+RATIONAL(-35.0,8.0)*y)*t1093;
+      coeff_dyy_0_m1 = t1063*t1092+RATIONAL(5.0,2.0)+(RATIONAL(-7.0,2.0)+
+RATIONAL(49.0,6.0)*t1093)*y+t1101;
+      coeff_dyy_p1_m1 = RATIONAL(-10.0,3.0)*t1092+y*t1046+RATIONAL(5.0,3.0)*x+
+t1038+t1095;
+      coeff_dyy_p2_m1 = RATIONAL(7.0,4.0)*t1117+t1068*x+t1025+RATIONAL(35.0,
+24.0)*t1092+t1098;
+      coeff_dyy_p3_m1 = t1093*t1116+t1068*t1092+t1096;
+      coeff_dyy_m2_0 = t1037+t1085*t1117+t1118*t1079+t1107;
+      coeff_dyy_m1_0 = RATIONAL(49.0,18.0)*t1092+RATIONAL(10.0,1.0)*t1117+t1024
++RATIONAL(28.0,9.0)*x+t1101;
+      coeff_dyy_0_0 = t1023+t1065+RATIONAL(-56.0,9.0)*t1092+RATIONAL(98.0,9.0)*
+t1093+(RATIONAL(8.0,1.0)+RATIONAL(-56.0,3.0)*t1093)*y;
+      coeff_dyy_p1_0 = RATIONAL(-28.0,9.0)*x+RATIONAL(40.0,3.0)*t1117+t1036+
+RATIONAL(56.0,9.0)*t1092+t1109;
+      coeff_dyy_p2_0 = t1044+RATIONAL(-49.0,18.0)*t1092+t1066*x+RATIONAL(-4.0,
+1.0)*t1117+t1100;
+      coeff_dyy_p3_0 = t1093*t1111+t1066*t1092+t1103;
+      coeff_dyy_m2_p1 = t1044+t1086*t1117+t1118*t1062+t1102;
+      coeff_dyy_m1_p1 = RATIONAL(-35.0,18.0)*t1092+RATIONAL(-20.0,9.0)*x+t1036+
+RATIONAL(-10.0,1.0)*t1117+t1095;
+      coeff_dyy_0_p1 = RATIONAL(40.0,9.0)*t1092+t1063+(RATIONAL(56.0,3.0)*t1093
++RATIONAL(-8.0,1.0))*y+t1109;
+      coeff_dyy_p1_p1 = t1023+RATIONAL(20.0,9.0)*x+t1024+RATIONAL(-40.0,9.0)*
+t1092+(RATIONAL(50.0,9.0)+RATIONAL(-40.0,3.0)*y)*t1093;
+      coeff_dyy_p2_p1 = RATIONAL(4.0,1.0)*t1117+t1080*x+t1037+RATIONAL(35.0,
+18.0)*t1092+t1099;
+      coeff_dyy_p3_p1 = t1080*t1092+t1089*t1117+t1097;
+      coeff_dyy_m2_p2 = t1073*t1117+t1025+t1118*t1075+t1106;
+      coeff_dyy_m1_p2 = t1084*x+t1074*t1092+RATIONAL(35.0,8.0)*t1117+t1038+
+t1098;
+      coeff_dyy_0_p2 = t1085*t1092+RATIONAL(-1.0,1.0)+(RATIONAL(7.0,2.0)+
+RATIONAL(-49.0,6.0)*t1093)*y+t1100;
+      coeff_dyy_p1_p2 = RATIONAL(35.0,6.0)*t1117+t1112+t1086*t1092+t1039+t1099;
+      coeff_dyy_p2_p2 = t1022+t1073*t1092+t1028+t1067*x+(RATIONAL(-7.0,4.0)*y+
+t1083)*t1093;
+      coeff_dyy_p3_p2 = t1093*t1113+t1067*t1092+t1108;
+      coeff_dyy_m2_p3 = t1034+t1067*t1117+t1118*t1077+t1105;
+      coeff_dyy_m1_p3 = RATIONAL(-5.0,8.0)*t1117+t1043+RATIONAL(-7.0,72.0)*
+t1092+RATIONAL(-1.0,9.0)*x+t1096;
+      coeff_dyy_0_p3 = t1090*t1092+t1081+(RATIONAL(7.0,6.0)*t1093+RATIONAL(-1.0
+,2.0))*y+t1103;
+      coeff_dyy_p1_p3 = t1040+t1091*t1092+RATIONAL(1.0,9.0)*x+RATIONAL(-5.0,6.0
+)*t1117+t1097;
+      coeff_dyy_p2_p3 = t1030+RATIONAL(7.0,72.0)*t1092+RATIONAL(1.0,4.0)*t1117+
+t1070*x+t1108;
+      coeff_dyy_p3_p3 = t1070*t1092+t1031+(t1077+t1114)*t1093;