aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs
diff options
context:
space:
mode:
authorjthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416>2002-08-20 16:31:28 +0000
committerjthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416>2002-08-20 16:31:28 +0000
commitd5b49cab5a99a22dd26bd55fcf65659aa7f7e020 (patch)
treee3df72179dc33241852470ee9f86b6d9e67b43e4 /src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs
parentd1aaafaad12880241613d3465da77245ab51bbff (diff)
add $Header:$ lines, switch a few $Id:$ --> $Header:$
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@92 df1f8a13-aa1d-4dd4-9681-27ded5b42416
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs')
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c76
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c30
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c22
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c20
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c22
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c264
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c126
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c64
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c52
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c126
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c56
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c394
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c254
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c180
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c190
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c244
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c176
17 files changed, 1152 insertions, 1144 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c
index 34f402b..52b6d48 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c
@@ -1,49 +1,53 @@
-fp t25,
- t28,
- t20,
- t27,
+fp t27,
+ t26,
t18,
+ t28,
+ t21,
t34,
t33,
t23,
+ t15,
t22,
- t17,
t32,
- t15,
+ t17,
t31,
+ t14,
t30,
+ t25,
+ t20,
t29,
- t26,
- t21,
t19,
t16,
- t14;
- t25 = RATIONAL(1.0,6.0);
- t28 = y*y;
- t20 = t25*t28;
+ t13,
+ t12;
t27 = x*x;
- t18 = t25*t27;
- t34 = t20+RATIONAL(-1.0,9.0)+t18;
+ t26 = RATIONAL(1.0,6.0);
+ t18 = t26*t27;
+ t28 = y*y;
+ t21 = t26*t28;
+ t34 = t18+t21+RATIONAL(-1.0,9.0);
t33 = x*y;
- t23 = RATIONAL(2.0,9.0);
- t22 = RATIONAL(-1.0,3.0);
- t17 = t22*t27;
- t32 = t20+t23+t17;
- t15 = t22*t28;
- t31 = t15+t23+t18;
- t30 = RATIONAL(-1.0,4.0)*t33+t34;
- t29 = RATIONAL(1.0,4.0)*t33+t34;
- t26 = RATIONAL(-1.0,6.0);
- t21 = t26*y;
- t19 = t26*x;
- t16 = t25*x;
- t14 = t25*y;
- coeff_I_m1_m1 = t21+t19+t29;
- coeff_I_0_m1 = t21+t32;
- coeff_I_p1_m1 = t21+t16+t30;
- coeff_I_m1_0 = t19+t31;
- coeff_I_0_0 = t15+t17+RATIONAL(5.0,9.0);
- coeff_I_p1_0 = t16+t31;
- coeff_I_m1_p1 = t14+t19+t30;
- coeff_I_0_p1 = t14+t32;
- coeff_I_p1_p1 = t14+t16+t29;
+ t23 = RATIONAL(-1.0,3.0);
+ t15 = t23*t28;
+ t22 = RATIONAL(2.0,9.0);
+ t32 = t15+t18+t22;
+ t17 = t23*t27;
+ t31 = t17+t21+t22;
+ t14 = t26*y;
+ t30 = t14+t34;
+ t25 = RATIONAL(-1.0,6.0);
+ t20 = t25*y;
+ t29 = t20+t34;
+ t19 = t25*x;
+ t16 = t26*x;
+ t13 = RATIONAL(1.0,4.0)*t33;
+ t12 = RATIONAL(-1.0,4.0)*t33;
+ coeff_I_m1_m1 = t13+t19+t29;
+ coeff_I_0_m1 = t20+t31;
+ coeff_I_p1_m1 = t12+t16+t29;
+ coeff_I_m1_0 = t19+t32;
+ coeff_I_0_0 = t15+RATIONAL(5.0,9.0)+t17;
+ coeff_I_p1_0 = t16+t32;
+ coeff_I_m1_p1 = t19+t12+t30;
+ coeff_I_0_p1 = t14+t31;
+ coeff_I_p1_p1 = t16+t13+t30;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c
index 22a3439..14b1a1c 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c
@@ -3,19 +3,19 @@ fp t37,
t41,
t40,
t39,
- t35;
+ t36;
t37 = RATIONAL(1.0,3.0)*x;
- t42 = t37+RATIONAL(1.0,4.0)*y;
- t41 = RATIONAL(-1.0,4.0)*y+t37;
- t40 = RATIONAL(-1.0,6.0);
- t39 = RATIONAL(1.0,6.0);
- t35 = RATIONAL(-2.0,3.0)*x;
- coeff_dx_m1_m1 = t40+t42;
- coeff_dx_0_m1 = t35;
- coeff_dx_p1_m1 = t39+t41;
- coeff_dx_m1_0 = t40+t37;
- coeff_dx_0_0 = t35;
- coeff_dx_p1_0 = t39+t37;
- coeff_dx_m1_p1 = t40+t41;
- coeff_dx_0_p1 = t35;
- coeff_dx_p1_p1 = t39+t42;
+ t42 = RATIONAL(1.0,4.0)*y+t37;
+ t41 = t37+RATIONAL(-1.0,4.0)*y;
+ t40 = RATIONAL(1.0,6.0);
+ t39 = RATIONAL(-1.0,6.0);
+ t36 = RATIONAL(-2.0,3.0)*x;
+ coeff_dx_m1_m1 = t39+t42;
+ coeff_dx_0_m1 = t36;
+ coeff_dx_p1_m1 = t40+t41;
+ coeff_dx_m1_0 = t39+t37;
+ coeff_dx_0_0 = t36;
+ coeff_dx_p1_0 = t40+t37;
+ coeff_dx_m1_p1 = t39+t41;
+ coeff_dx_0_p1 = t36;
+ coeff_dx_p1_p1 = t40+t42;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c
index 5f6427a..3381569 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c
@@ -1,13 +1,13 @@
fp t52,
t51;
- t52 = RATIONAL(1.0,3.0);
- t51 = RATIONAL(-2.0,3.0);
- coeff_dxx_m1_m1 = t52;
- coeff_dxx_0_m1 = t51;
- coeff_dxx_p1_m1 = t52;
- coeff_dxx_m1_0 = t52;
- coeff_dxx_0_0 = t51;
- coeff_dxx_p1_0 = t52;
- coeff_dxx_m1_p1 = t52;
- coeff_dxx_0_p1 = t51;
- coeff_dxx_p1_p1 = t52;
+ t52 = RATIONAL(-2.0,3.0);
+ t51 = RATIONAL(1.0,3.0);
+ coeff_dxx_m1_m1 = t51;
+ coeff_dxx_0_m1 = t52;
+ coeff_dxx_p1_m1 = t51;
+ coeff_dxx_m1_0 = t51;
+ coeff_dxx_0_0 = t52;
+ coeff_dxx_p1_0 = t51;
+ coeff_dxx_m1_p1 = t51;
+ coeff_dxx_0_p1 = t52;
+ coeff_dxx_p1_p1 = t51;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c
index 4afa754..65e74bd 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c
@@ -5,17 +5,17 @@ fp t46,
t47,
t43;
t46 = RATIONAL(1.0,3.0)*y;
- t50 = RATIONAL(-1.0,4.0)*x+t46;
- t49 = RATIONAL(1.0,4.0)*x+t46;
- t48 = RATIONAL(-1.0,6.0);
- t47 = RATIONAL(1.0,6.0);
+ t50 = RATIONAL(1.0,4.0)*x+t46;
+ t49 = RATIONAL(-1.0,4.0)*x+t46;
+ t48 = RATIONAL(1.0,6.0);
+ t47 = RATIONAL(-1.0,6.0);
t43 = RATIONAL(-2.0,3.0)*y;
- coeff_dy_m1_m1 = t48+t49;
- coeff_dy_0_m1 = t46+t48;
- coeff_dy_p1_m1 = t48+t50;
+ coeff_dy_m1_m1 = t47+t50;
+ coeff_dy_0_m1 = t47+t46;
+ coeff_dy_p1_m1 = t47+t49;
coeff_dy_m1_0 = t43;
coeff_dy_0_0 = t43;
coeff_dy_p1_0 = t43;
- coeff_dy_m1_p1 = t47+t50;
- coeff_dy_0_p1 = t47+t46;
- coeff_dy_p1_p1 = t47+t49;
+ coeff_dy_m1_p1 = t48+t49;
+ coeff_dy_0_p1 = t46+t48;
+ coeff_dy_p1_p1 = t48+t50;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c
index a643108..5d0b896 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c
@@ -1,13 +1,13 @@
fp t57,
t56;
- t57 = RATIONAL(1.0,3.0);
- t56 = RATIONAL(-2.0,3.0);
- coeff_dyy_m1_m1 = t57;
- coeff_dyy_0_m1 = t57;
- coeff_dyy_p1_m1 = t57;
- coeff_dyy_m1_0 = t56;
- coeff_dyy_0_0 = t56;
- coeff_dyy_p1_0 = t56;
- coeff_dyy_m1_p1 = t57;
- coeff_dyy_0_p1 = t57;
- coeff_dyy_p1_p1 = t57;
+ t57 = RATIONAL(-2.0,3.0);
+ t56 = RATIONAL(1.0,3.0);
+ coeff_dyy_m1_m1 = t56;
+ coeff_dyy_0_m1 = t56;
+ coeff_dyy_p1_m1 = t56;
+ coeff_dyy_m1_0 = t57;
+ coeff_dyy_0_0 = t57;
+ coeff_dyy_p1_0 = t57;
+ coeff_dyy_m1_p1 = t56;
+ coeff_dyy_0_p1 = t56;
+ coeff_dyy_p1_p1 = t56;
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
index 55cc52e..a022b8a 100644
--- 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
@@ -1,73 +1,73 @@
fp t160,
t131,
t159,
- t130,
- t158,
- t124,
- t68,
t133,
- t126,
- t82,
- t157,
- t71,
+ t158,
+ t117,
+ t66,
t132,
- t83,
- t156,
- t121,
- t70,
t128,
t87,
- t155,
- t67,
+ t157,
t129,
- t86,
+ t84,
+ t156,
+ t124,
+ t68,
+ t130,
+ t116,
+ t77,
+ t155,
+ t126,
+ t70,
+ t125,
+ t83,
t154,
- t85,
+ t71,
t153,
- t117,
- t77,
+ t108,
+ t122,
+ t69,
t152,
- t127,
- t62,
t151,
t63,
t150,
- t125,
- t64,
+ t67,
+ t86,
t149,
- t79,
+ t65,
t148,
+ t85,
t147,
+ t79,
t146,
- t69,
- t110,
+ t82,
t145,
- t84,
t144,
t143,
- t123,
t142,
- t104,
+ t123,
t141,
+ t101,
t140,
t139,
t138,
t137,
t136,
- t122,
- t120,
+ t127,
+ t121,
t119,
- t116,
- t114,
+ t118,
+ t115,
t113,
+ t112,
t111,
- t108,
- t107,
+ t109,
+ t106,
t105,
t103,
t102,
- t101,
- t99,
+ t100,
t98,
t97,
t96,
@@ -85,112 +85,112 @@ fp t160,
t76,
t75,
t74,
- t73;
+ t72;
t160 = x*y;
- t131 = y*y;
- t159 = x*t131;
- t130 = x*x;
- t158 = t130*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;
- t133 = t131*y;
- t126 = RATIONAL(-1.0,24.0);
- t82 = t126*t133;
- t157 = t68+t82;
+ 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;
- t132 = t130*x;
- t83 = t126*t132;
- t156 = t71+t83;
- t121 = RATIONAL(1.0,40.0);
- t70 = t121*t158;
- t128 = RATIONAL(-1.0,8.0);
- t87 = t128*t133;
- t155 = t70+t87;
- t67 = t121*t159;
- t129 = RATIONAL(1.0,8.0);
- t86 = t129*t132;
- t154 = t67+t86;
- t85 = t128*t132;
- t153 = t67+t85;
- t117 = RATIONAL(1.0,24.0);
- t77 = t117*t132;
- t152 = t71+t77;
- t127 = RATIONAL(-3.0,40.0);
- t62 = t127*t158;
- t151 = t62+t82;
- t63 = t127*t159;
- t150 = t63+t77;
- t125 = RATIONAL(-1.0,40.0);
- t64 = t125*t159;
- t149 = t64+t86;
- t79 = t117*t133;
- t148 = t68+t79;
- t147 = t63+t83;
- t146 = t62+t79;
- t69 = t125*t158;
- t110 = RATIONAL(11.0,80.0);
- t145 = t69+t110*t130;
- t84 = t129*t133;
- t144 = t70+t84;
- t143 = t64+t85;
+ 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);
- t142 = RATIONAL(-13.0,100.0)*t160+t123;
- t104 = RATIONAL(-17.0,200.0);
- t141 = RATIONAL(-9.0,100.0)*t160+t104;
- t140 = t123*t160+RATIONAL(11.0,50.0);
- t139 = RATIONAL(-1.0,100.0)*t160+RATIONAL(37.0,200.0);
- t138 = RATIONAL(-2.0,25.0)+RATIONAL(2.0,25.0)*t160;
+ 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 = t130+t131;
- t122 = RATIONAL(9.0,80.0);
- t120 = RATIONAL(7.0,80.0);
+ t136 = t133+t131;
+ t127 = RATIONAL(1.0,80.0);
+ t121 = RATIONAL(9.0,80.0);
t119 = RATIONAL(-1.0,80.0);
- t116 = RATIONAL(3.0,80.0);
- t114 = RATIONAL(1.0,80.0);
- t113 = RATIONAL(13.0,80.0);
- t111 = RATIONAL(-23.0,80.0);
- t108 = RATIONAL(-19.0,80.0);
- t107 = RATIONAL(-17.0,80.0);
- t105 = RATIONAL(-31.0,400.0);
- t103 = RATIONAL(69.0,400.0);
- t102 = RATIONAL(-57.0,400.0);
- t101 = RATIONAL(-63.0,400.0);
- t99 = RATIONAL(-49.0,400.0);
- t98 = RATIONAL(43.0,1200.0);
- t97 = RATIONAL(103.0,400.0);
- t96 = RATIONAL(111.0,400.0);
+ 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(89.0,1200.0);
- t93 = RATIONAL(117.0,400.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(-131.0,1200.0);
- t88 = RATIONAL(-149.0,1200.0);
- t81 = t122*t131;
- t80 = t120*t130;
- t78 = t120*t131;
- t76 = t122*t130;
- t75 = t113*t131;
- t74 = t113*t130;
- t73 = t110*t131;
+ 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+t147+t151;
- coeff_I_0_m1 = t111*t130+t105*x+t73+t91*y+t140+t149+t157;
- coeff_I_p1_m1 = t81+t96*x+t90*y+t80+t142+t153+t157;
- coeff_I_p2_m1 = t78+t94*y+t89*x+t116*t130+t141+t151+t152;
- coeff_I_m1_0 = t105*y+t91*x+t111*t131+t84+t140+t145+t156;
- coeff_I_0_0 = RATIONAL(93.0,200.0)+t102*y+t136*RATIONAL(-21.0,80.0)+(
-RATIONAL(-1.0,25.0)*y+t102)*x+t144+t154;
- coeff_I_p1_0 = t108*t131+t101*y+t76+t93*x+t139+t143+t144;
- coeff_I_p2_0 = t107*t131+t114*t130+t99*y+t84+t69+t98*x+t137+t150;
- coeff_I_m1_p1 = t90*x+t96*y+t78+t76+t142+t155+t156;
- coeff_I_0_p1 = t69+t101*x+t108*t130+t87+t93*y+t81+t139+t154;
- coeff_I_p1_p1 = t97*y+t73+t104+t87+(RATIONAL(3.0,50.0)*y+t97)*x+t143+t145
+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 = t95*x+t119*t130+t75+t103*y+t138+t150+t155;
- coeff_I_m1_p2 = t116*t131+t94*x+t89*y+t80+t141+t147+t148;
- coeff_I_0_p2 = t99*x+t98*y+t114*t131+t107*t130+t137+t146+t149;
- coeff_I_p1_p2 = t103*x+t119*t131+t95*y+t74+t138+t146+t153;
- coeff_I_p2_p2 = RATIONAL(7.0,200.0)+t88*y+t136*RATIONAL(-3.0,80.0)+(t88+
-RATIONAL(-3.0,50.0)*y)*x+t148+t152;
+ 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
index 4650764..5bfa6f4 100644
--- 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
@@ -1,95 +1,95 @@
fp t201,
- t164,
+ t161,
t200,
- t163,
- t199,
- t183,
+ t184,
t188,
- t173,
+ t174,
t187,
- t179,
+ t180,
+ t199,
+ t164,
t198,
- t161,
+ t183,
+ t175,
t197,
t162,
t196,
- t185,
- t176,
- t195,
t181,
+ t195,
t194,
+ t163,
t193,
t186,
- t177,
+ t176,
t182,
t192,
- t180,
- t191,
- t184,
+ t185,
t172,
+ t191,
+ t179,
t190,
t189,
t178,
- t175,
- t170,
- t168,
+ t177,
+ t169,
+ t167,
t166,
t165;
t201 = x*y;
- t164 = RATIONAL(-3.0,20.0)*t201;
- t200 = RATIONAL(13.0,40.0)*x+t164;
- t163 = RATIONAL(3.0,20.0)*t201;
- t199 = RATIONAL(7.0,40.0)*x+t163;
- t183 = RATIONAL(1.0,40.0);
+ t161 = RATIONAL(1.0,20.0)*t201;
+ t200 = t161+RATIONAL(9.0,40.0)*x;
+ t184 = RATIONAL(-1.0,40.0);
t188 = y*y;
- t173 = t183*t188;
+ t174 = t184*t188;
t187 = x*x;
- t179 = RATIONAL(-3.0,8.0)*t187;
- t198 = t173+t179;
- t161 = RATIONAL(1.0,20.0)*t201;
- t197 = RATIONAL(9.0,40.0)*x+t161;
+ 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 = RATIONAL(11.0,40.0)*x+t162;
- t185 = RATIONAL(-1.0,40.0);
- t176 = t185*t188;
- t195 = t176+t179;
+ t196 = t162+RATIONAL(11.0,40.0)*x;
t181 = RATIONAL(3.0,8.0)*t187;
- t194 = t173+t181;
- t193 = t176+t181;
+ 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);
- t177 = t186*t188;
+ t176 = t186*t188;
t182 = RATIONAL(-1.0,8.0)*t187;
- t192 = t177+t182;
- t180 = RATIONAL(1.0,8.0)*t187;
- t191 = t177+t180;
- t184 = RATIONAL(3.0,40.0);
- t172 = t184*t188;
- t190 = t172+t180;
- t189 = t172+t182;
+ 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;
- t175 = RATIONAL(-1.0,50.0)*y;
- t170 = RATIONAL(7.0,100.0)*y;
- t168 = RATIONAL(-9.0,100.0)*y;
- t166 = RATIONAL(-1.0,100.0)*y;
+ 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+t200;
- coeff_dx_0_m1 = RATIONAL(-31.0,400.0)+RATIONAL(-23.0,40.0)*x+t175+t163+
-t193;
- coeff_dx_p1_m1 = RATIONAL(111.0,400.0)+t165+t198+t199;
- coeff_dx_p2_m1 = t168+t164+RATIONAL(-131.0,1200.0)+t184*x+t190;
- coeff_dx_m1_0 = RATIONAL(-223.0,1200.0)+t175+t189+t196;
- coeff_dx_0_0 = RATIONAL(-57.0,400.0)+RATIONAL(-21.0,40.0)*x+RATIONAL(-1.0
-,25.0)*y+t161+t194;
- coeff_dx_p1_0 = RATIONAL(117.0,400.0)+t166+t195+t197;
- coeff_dx_p2_0 = RATIONAL(43.0,1200.0)+t183*x+t162+t170+t191;
- coeff_dx_m1_p1 = RATIONAL(-157.0,1200.0)+t165+t189+t197;
- coeff_dx_0_p1 = RATIONAL(-63.0,400.0)+t166+t162+RATIONAL(-19.0,40.0)*x+
+ 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)+t195+t196;
- coeff_dx_p2_p1 = t185*x+t161+RATIONAL(37.0,1200.0)+t178+t191;
- coeff_dx_m1_p2 = t168+RATIONAL(89.0,1200.0)+t192+t199;
- coeff_dx_0_p2 = t164+t170+RATIONAL(-17.0,40.0)*x+RATIONAL(-49.0,400.0)+
-t193;
- coeff_dx_p1_p2 = RATIONAL(69.0,400.0)+t178+t198+t200;
+ coeff_dx_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
index 395eceb..29b5494 100644
--- 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
@@ -1,40 +1,40 @@
-fp t243,
+fp t244,
t258,
- t245,
- t257,
t246,
+ t257,
+ t245,
t256,
- t244,
+ t243,
t255,
t250,
t249,
t248,
t247;
- t243 = RATIONAL(1.0,20.0)*y;
- t258 = t243+RATIONAL(9.0,40.0);
- t245 = RATIONAL(3.0,20.0)*y;
- t257 = t245+RATIONAL(7.0,40.0);
- t246 = RATIONAL(-3.0,20.0)*y;
- t256 = t246+RATIONAL(13.0,40.0);
- t244 = RATIONAL(-1.0,20.0)*y;
- t255 = t244+RATIONAL(11.0,40.0);
- t250 = RATIONAL(-1.0,4.0)*x;
- t249 = RATIONAL(3.0,4.0)*x;
- t248 = RATIONAL(-3.0,4.0)*x;
- t247 = RATIONAL(1.0,4.0)*x;
- coeff_dxx_m1_m1 = t250+t256;
- coeff_dxx_0_m1 = RATIONAL(-23.0,40.0)+t245+t249;
- coeff_dxx_p1_m1 = t248+t257;
- coeff_dxx_p2_m1 = t246+t247+RATIONAL(3.0,40.0);
- coeff_dxx_m1_0 = t250+t255;
- coeff_dxx_0_0 = t243+RATIONAL(-21.0,40.0)+t249;
- coeff_dxx_p1_0 = t248+t258;
- coeff_dxx_p2_0 = RATIONAL(1.0,40.0)+t244+t247;
- coeff_dxx_m1_p1 = t250+t258;
- coeff_dxx_0_p1 = t244+RATIONAL(-19.0,40.0)+t249;
- coeff_dxx_p1_p1 = t248+t255;
- coeff_dxx_p2_p1 = t243+t247+RATIONAL(-1.0,40.0);
- coeff_dxx_m1_p2 = t250+t257;
- coeff_dxx_0_p2 = RATIONAL(-17.0,40.0)+t246+t249;
- coeff_dxx_p1_p2 = t248+t256;
- coeff_dxx_p2_p2 = RATIONAL(-3.0,40.0)+t247+t245;
+ 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
index b48e5bc..20e2106 100644
--- 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
@@ -17,36 +17,36 @@ fp t276,
t260,
t259;
t276 = RATIONAL(-1.0,50.0);
- t275 = RATIONAL(-3.0,20.0);
- t274 = RATIONAL(3.0,20.0);
- t273 = RATIONAL(2.0,25.0);
+ 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(7.0,100.0);
- t268 = RATIONAL(-1.0,100.0);
+ t269 = RATIONAL(-1.0,100.0);
+ t268 = RATIONAL(7.0,100.0);
t267 = RATIONAL(-13.0,100.0);
- t266 = t275*y;
- t265 = t274*y;
- t264 = t275*x;
+ t266 = t272*y;
+ t265 = t271*y;
+ t264 = t273*x;
t263 = t274*x;
t262 = t271*x;
t261 = t272*x;
- t260 = t272*y;
- t259 = t271*y;
- coeff_dxy_m1_m1 = t266+t264+RATIONAL(6.0,25.0);
- coeff_dxy_0_m1 = t263+t260+t276;
- coeff_dxy_p1_m1 = t263+t259+t267;
- coeff_dxy_p2_m1 = t270+t264+t265;
- coeff_dxy_m1_0 = t276+t261+t265;
- coeff_dxy_0_0 = RATIONAL(-1.0,25.0)+t259+t262;
- coeff_dxy_p1_0 = t268+t262+t260;
- coeff_dxy_p2_0 = t261+t269+t266;
- coeff_dxy_m1_p1 = t267+t262+t265;
- coeff_dxy_0_p1 = t261+t268+t259;
- coeff_dxy_p1_p1 = RATIONAL(3.0,50.0)+t261+t260;
- coeff_dxy_p2_p1 = t266+t262+t273;
- coeff_dxy_m1_p2 = t266+t263+t270;
- coeff_dxy_0_p2 = t269+t264+t260;
- coeff_dxy_p1_p2 = t264+t259+t273;
- coeff_dxy_p2_p2 = t263+t265+RATIONAL(-3.0,50.0);
+ 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
index 678f775..c766157 100644
--- 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
@@ -1,95 +1,95 @@
fp t242,
+ t225,
t228,
- t226,
t212,
t229,
t220,
t241,
- t224,
- t213,
+ t203,
t240,
- t202,
+ t204,
t239,
- t223,
+ t205,
t238,
- t204,
+ t224,
+ t213,
t237,
+ t202,
t236,
- t227,
- t217,
- t221,
+ t226,
+ t215,
+ t222,
t235,
- t205,
+ t227,
+ t218,
t234,
- t203,
+ t221,
t233,
- t225,
- t218,
t232,
- t222,
+ t223,
t231,
t230,
t219,
- t214,
+ t217,
t209,
t208,
t207,
t206;
t242 = x*y;
+ t225 = RATIONAL(-1.0,40.0);
t228 = x*x;
- t226 = RATIONAL(-1.0,40.0);
- t212 = t226*t228;
+ t212 = t225*t228;
t229 = y*y;
- t220 = RATIONAL(3.0,8.0)*t229;
+ 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;
- t240 = t213+t220;
+ t237 = t213+t220;
t202 = RATIONAL(3.0,20.0)*t242;
- t239 = t202+RATIONAL(7.0,40.0)*y;
- t223 = RATIONAL(-3.0,8.0)*t229;
- t238 = t213+t223;
- t204 = RATIONAL(1.0,20.0)*t242;
- t237 = t204+RATIONAL(9.0,40.0)*y;
- t236 = t212+t223;
+ 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);
- t217 = t227*t228;
- t221 = RATIONAL(-1.0,8.0)*t229;
- t235 = t217+t221;
- t205 = RATIONAL(-3.0,20.0)*t242;
- t234 = t205+RATIONAL(13.0,40.0)*y;
- t203 = RATIONAL(-1.0,20.0)*t242;
- t233 = t203+RATIONAL(11.0,40.0)*y;
- t225 = RATIONAL(3.0,40.0);
- t218 = t225*t228;
- t232 = t218+t221;
- t222 = RATIONAL(1.0,8.0)*t229;
- t231 = t218+t222;
- t230 = t217+t222;
+ 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;
- t214 = RATIONAL(-1.0,50.0)*x;
- t209 = RATIONAL(-1.0,100.0)*x;
- t208 = RATIONAL(-9.0,100.0)*x;
- t207 = RATIONAL(7.0,100.0)*x;
+ 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(-109.0,1200.0)+RATIONAL(6.0,25.0)*x+t234+t235;
- coeff_dy_0_m1 = t214+RATIONAL(-223.0,1200.0)+t232+t233;
- coeff_dy_p1_m1 = t206+RATIONAL(-157.0,1200.0)+t232+t237;
- coeff_dy_p2_m1 = t208+RATIONAL(89.0,1200.0)+t235+t239;
- coeff_dy_m1_0 = RATIONAL(-23.0,40.0)*y+t214+RATIONAL(-31.0,400.0)+t202+
-t241;
- coeff_dy_0_0 = RATIONAL(-21.0,40.0)*y+RATIONAL(-57.0,400.0)+t204+RATIONAL
-(-1.0,25.0)*x+t240;
- coeff_dy_p1_0 = t209+RATIONAL(-63.0,400.0)+t203+RATIONAL(-19.0,40.0)*y+
-t240;
- coeff_dy_p2_0 = t207+RATIONAL(-49.0,400.0)+RATIONAL(-17.0,40.0)*y+t205+
-t241;
- coeff_dy_m1_p1 = RATIONAL(111.0,400.0)+t206+t238+t239;
- coeff_dy_0_p1 = RATIONAL(117.0,400.0)+t209+t236+t237;
- coeff_dy_p1_p1 = RATIONAL(103.0,400.0)+RATIONAL(3.0,50.0)*x+t233+t236;
- coeff_dy_p2_p1 = t219+RATIONAL(69.0,400.0)+t234+t238;
- coeff_dy_m1_p2 = t225*y+t208+RATIONAL(-131.0,1200.0)+t205+t231;
- coeff_dy_0_p2 = t207+t203+t224*y+RATIONAL(43.0,1200.0)+t230;
- coeff_dy_p1_p2 = t219+t226*y+RATIONAL(37.0,1200.0)+t204+t230;
- coeff_dy_p2_p2 = t202+RATIONAL(-3.0,50.0)*x+RATIONAL(-149.0,1200.0)+t227*
-y+t231;
+ 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
index dea29f0..b449df9 100644
--- 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
@@ -1,40 +1,40 @@
-fp t279,
+fp t277,
t292,
- t280,
+ t279,
t291,
- t278,
+ t280,
t290,
- t277,
+ 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;
- t292 = t279+RATIONAL(11.0,40.0);
+ t291 = t279+RATIONAL(11.0,40.0);
t280 = RATIONAL(-3.0,20.0)*x;
- t291 = RATIONAL(13.0,40.0)+t280;
+ t290 = t280+RATIONAL(13.0,40.0);
t278 = RATIONAL(1.0,20.0)*x;
- t290 = t278+RATIONAL(9.0,40.0);
- t277 = RATIONAL(3.0,20.0)*x;
- t289 = t277+RATIONAL(7.0,40.0);
+ t289 = t278+RATIONAL(9.0,40.0);
t284 = RATIONAL(-1.0,4.0)*y;
- t283 = RATIONAL(1.0,4.0)*y;
- t282 = RATIONAL(3.0,4.0)*y;
- t281 = RATIONAL(-3.0,4.0)*y;
- coeff_dyy_m1_m1 = t284+t291;
- coeff_dyy_0_m1 = t284+t292;
- coeff_dyy_p1_m1 = t284+t290;
- coeff_dyy_p2_m1 = t284+t289;
- coeff_dyy_m1_0 = t282+RATIONAL(-23.0,40.0)+t277;
- coeff_dyy_0_0 = RATIONAL(-21.0,40.0)+t278+t282;
- coeff_dyy_p1_0 = RATIONAL(-19.0,40.0)+t282+t279;
- coeff_dyy_p2_0 = t280+t282+RATIONAL(-17.0,40.0);
- coeff_dyy_m1_p1 = t281+t289;
- coeff_dyy_0_p1 = t281+t290;
- coeff_dyy_p1_p1 = t281+t292;
- coeff_dyy_p2_p1 = t281+t291;
- coeff_dyy_m1_p2 = t280+t283+RATIONAL(3.0,40.0);
- coeff_dyy_0_p2 = t283+RATIONAL(1.0,40.0)+t279;
- coeff_dyy_p1_p2 = t278+RATIONAL(-1.0,40.0)+t283;
- coeff_dyy_p2_p2 = t277+t283+RATIONAL(-3.0,40.0);
+ 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.order4.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
index 49655f0..cfe391a 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
@@ -1,239 +1,243 @@
-fp t413,
+fp t417,
t418,
- t466,
- t403,
- t417,
- t358,
+ t468,
t414,
+ t411,
+ t361,
+ t416,
t373,
- t438,
- t444,
- t465,
- t405,
- t364,
+ t435,
+ t446,
+ t467,
+ t410,
t367,
- t439,
+ t368,
+ t437,
+ t466,
+ t465,
+ t441,
t464,
- t436,
t463,
- t410,
- t359,
- t442,
- t409,
- t318,
- t392,
- t462,
- t411,
- t360,
- t441,
+ t395,
+ t445,
t400,
- t298,
- t461,
+ t403,
t366,
- t443,
- t395,
+ t462,
+ t397,
+ t404,
+ t362,
t401,
+ t461,
+ t389,
+ t409,
+ t316,
t460,
- t440,
- t296,
- t399,
+ t402,
+ t312,
+ t384,
t459,
- t372,
+ t412,
+ t443,
+ t407,
+ t309,
t458,
- t333,
+ t295,
+ t399,
+ t442,
t457,
- t396,
- t412,
- t327,
+ t406,
+ t405,
+ t319,
t456,
- t407,
- t331,
- t391,
+ t444,
+ t370,
t455,
- t402,
- t313,
- t397,
+ t315,
t454,
t325,
t453,
- t324,
+ t358,
+ t398,
+ t294,
t452,
- t311,
t451,
+ t413,
+ t364,
t450,
- t387,
- t398,
+ t357,
+ t380,
+ t381,
t449,
+ t332,
t415,
- t371,
- t416,
- t357,
- t383,
- t382,
- t448,
- t317,
t356,
+ t448,
t447,
- t362,
- t446,
- t445,
- t437,
- t435,
- t393,
- t404,
- t369,
- t376,
+ t440,
+ t360,
+ t439,
+ t438,
+ t436,
+ t390,
+ t408,
+ t363,
+ t392,
t434,
- t365,
- t384,
t386,
+ t374,
t433,
+ t372,
t432,
t431,
+ t317,
t430,
t429,
+ t326,
t428,
t427,
- t406,
+ t385,
t426,
t425,
- t408,
+ t388,
t424,
t423,
t394,
- t389,
- t388,
- t385;
- t413 = x*x;
+ t393,
+ t391,
+ t387;
+ t417 = x*x;
t418 = y*y;
- t466 = t418+t413;
- t403 = RATIONAL(-1.0,60.0);
- t417 = y*t418;
- t358 = t403*t417;
- t414 = t413*x;
- t373 = t403*t414;
- t438 = y*t373+x*t358;
- t444 = x*y;
- t465 = t438+RATIONAL(37.0,300.0)*t444;
- t405 = RATIONAL(1.0,60.0);
- t364 = t405*t414;
- t367 = t405*t417;
- t439 = y*t364+x*t367;
- t464 = t439+RATIONAL(-37.0,300.0)*t444;
- t436 = t413*t418;
- t463 = RATIONAL(-99.0,1225.0)+RATIONAL(-1.0,98.0)*t436;
- t410 = RATIONAL(-1.0,30.0);
- t359 = t410*t417;
- t442 = y*t413;
- t409 = RATIONAL(1.0,70.0);
- t318 = t409*t442;
- t392 = RATIONAL(11.0,105.0);
- t462 = t359+t318+t392*y;
- t411 = RATIONAL(1.0,30.0);
- t360 = t411*t417;
- t441 = y*t414;
- t400 = RATIONAL(-1.0,120.0);
- t298 = RATIONAL(-21.0,200.0)*t444;
- t461 = x*t360+t298+t400*t441;
- t366 = t410*t414;
- t443 = x*t418;
- t395 = RATIONAL(31.0,210.0);
+ t468 = t418+t417;
+ t414 = t417*x;
+ t411 = RATIONAL(1.0,60.0);
+ t361 = t411*t414;
+ t416 = y*t418;
+ t373 = t411*t416;
+ t435 = x*t373+y*t361;
+ t446 = x*y;
+ t467 = RATIONAL(-11.0,150.0)*t446+t435;
+ t410 = RATIONAL(-1.0,60.0);
+ t367 = t410*t416;
+ t368 = t410*t414;
+ t437 = x*t367+y*t368;
+ t466 = RATIONAL(37.0,300.0)*t446+t437;
+ t465 = t437+RATIONAL(11.0,150.0)*t446;
+ t441 = t417*t418;
+ t464 = RATIONAL(-1.0,98.0)*t441+RATIONAL(-99.0,1225.0);
+ t463 = RATIONAL(-37.0,300.0)*t446+t435;
+ t395 = RATIONAL(-31.0,210.0);
+ t445 = x*t418;
+ t400 = RATIONAL(1.0,140.0);
+ t403 = RATIONAL(1.0,30.0);
+ t366 = t403*t414;
+ t462 = t400*t445+t395*x+t366;
+ t397 = RATIONAL(31.0,210.0);
+ t404 = RATIONAL(-1.0,30.0);
+ t362 = t404*t414;
t401 = RATIONAL(-1.0,140.0);
- t460 = t395*x+t366+t401*t443;
- t440 = x*t417;
- t296 = RATIONAL(21.0,200.0)*t444;
+ t461 = t362+t401*t445+t397*x;
+ t389 = RATIONAL(31.0,420.0);
+ t409 = RATIONAL(-1.0,35.0);
+ t316 = t409*t445;
+ t460 = t389*x+t316+t368;
+ t402 = RATIONAL(1.0,35.0);
+ t312 = t402*t445;
+ t384 = RATIONAL(-31.0,420.0);
+ t459 = t312+t384*x+t361;
+ t412 = RATIONAL(1.0,84.0);
+ t443 = y*t417;
+ t407 = RATIONAL(-1.0,70.0);
+ t309 = t407*t443;
+ t458 = t309+t412*y+t373;
+ t295 = RATIONAL(21.0,200.0)*t446;
t399 = RATIONAL(1.0,120.0);
- t459 = t296+y*t366+t399*t440;
- t372 = t411*t414;
- t458 = t400*t440+y*t372+t298;
- t333 = t409*t443;
- t457 = t366+t333+t392*x;
- t396 = RATIONAL(-11.0,105.0);
- t412 = RATIONAL(-1.0,70.0);
- t327 = t412*t443;
- t456 = t396*x+t372+t327;
- t407 = RATIONAL(1.0,35.0);
- t331 = t407*t442;
- t391 = RATIONAL(-31.0,420.0);
- t455 = t331+t391*y+t367;
- t402 = RATIONAL(-1.0,35.0);
- t313 = t402*t442;
- t397 = RATIONAL(31.0,420.0);
- t454 = t358+t397*y+t313;
- t325 = t402*t443;
- t453 = t397*x+t325+t373;
- t324 = t407*t443;
- t452 = t364+t324+t391*x;
- t311 = t412*t442;
- t451 = t360+t311+t396*y;
- t450 = x*t359+t296+t399*t441;
- t387 = RATIONAL(-31.0,210.0);
- t398 = RATIONAL(1.0,140.0);
- t449 = t372+t387*x+t398*t443;
- t415 = t413*t413;
- t371 = t410*t415;
- t416 = t418*t418;
- t357 = t399*t416;
- t383 = RATIONAL(71.0,5880.0);
- t382 = RATIONAL(113.0,735.0);
- t448 = t371+t382*t413+t383*t418+t357+t463;
- t317 = RATIONAL(1.0,49.0)*t436;
+ t442 = x*t416;
+ t457 = t295+t399*t442+y*t362;
+ t406 = RATIONAL(-1.0,84.0);
+ t405 = RATIONAL(1.0,70.0);
+ t319 = t405*t443;
+ t456 = t406*y+t319+t367;
+ t444 = t414*y;
+ t370 = t404*t416;
+ t455 = t295+t399*t444+x*t370;
+ t315 = t405*t445;
+ t454 = t406*x+t315+t368;
+ t325 = t407*t445;
+ t453 = t361+t412*x+t325;
+ t358 = t403*t416;
+ t398 = RATIONAL(-1.0,120.0);
+ t294 = RATIONAL(-21.0,200.0)*t446;
+ t452 = t398*t444+t294+x*t358;
+ t451 = y*t366+t398*t442+t294;
+ t413 = t417*t417;
+ t364 = t404*t413;
+ t450 = RATIONAL(-24.0,1225.0)+t364+RATIONAL(1.0,196.0)*t441+t468*RATIONAL
+(181.0,1470.0);
+ t357 = t399*t413;
+ t380 = RATIONAL(71.0,5880.0);
+ t381 = RATIONAL(113.0,735.0);
+ t449 = t357+t380*t417+t381*t418+t464;
+ t332 = RATIONAL(1.0,49.0)*t441;
+ t415 = t418*t418;
t356 = t399*t415;
- t447 = t317+RATIONAL(51.0,1225.0)+t356+t357+t466*RATIONAL(-289.0,5880.0);
- t362 = t410*t416;
- t446 = t382*t418+t356+t383*t413+t362+t463;
- t445 = t362+RATIONAL(1.0,196.0)*t436+t371+RATIONAL(-24.0,1225.0)+t466*
-RATIONAL(181.0,1470.0);
- t437 = RATIONAL(-1.0,49.0)*t436+RATIONAL(96.0,1225.0);
- t435 = RATIONAL(1.0,98.0)*t436+RATIONAL(246.0,1225.0);
- t393 = RATIONAL(-41.0,196.0);
- t404 = RATIONAL(1.0,20.0);
- t369 = t404*t415;
- t376 = RATIONAL(191.0,5880.0);
- t434 = t376*t418+t369+t357+t393*t413+t437;
- t365 = t404*t416;
- t384 = RATIONAL(-53.0,196.0);
- t386 = RATIONAL(83.0,735.0);
- t433 = t365+t371+t384*t418+t386*t413+t435;
- t432 = t365+t376*t413+t356+t393*t418+t437;
- t431 = t369+t384*t413+t362+t386*t418+t435;
- t430 = RATIONAL(-11.0,150.0)*t444+t439+t447;
- t429 = t359+t395*y+t401*t442+t445;
- t428 = t360+t387*y+t398*t442+t445;
- t427 = RATIONAL(11.0,150.0)*t444+t438+t447;
- t406 = RATIONAL(1.0,84.0);
- t426 = t367+t311+t406*y+t448;
- t425 = t364+t406*x+t327+t446;
- t408 = RATIONAL(-1.0,84.0);
- t424 = t358+t318+t408*y+t448;
- t423 = t333+t408*x+t373+t446;
- t394 = RATIONAL(-17.0,420.0);
- t389 = RATIONAL(17.0,105.0);
- t388 = RATIONAL(17.0,420.0);
- t385 = RATIONAL(-17.0,105.0);
- coeff_I_m2_m2 = t430+t453+t454;
- coeff_I_m1_m2 = t424+t456+t459;
- coeff_I_0_m2 = t331+t394*y+t358+t434;
- coeff_I_p1_m2 = t424+t457+t458;
- coeff_I_p2_m2 = t427+t452+t454;
- coeff_I_m2_m1 = t423+t450+t451;
- coeff_I_m1_m1 = t428+t449+t465;
- coeff_I_0_m1 = t360+t385*y+t318+t431;
- coeff_I_p1_m1 = t428+t460+t464;
- coeff_I_p2_m1 = t425+t451+t461;
- coeff_I_m2_0 = t373+t394*x+t324+t432;
- coeff_I_m1_0 = t385*x+t333+t372+t433;
- coeff_I_0_0 = t317+t365+t369+RATIONAL(541.0,1225.0)+t466*RATIONAL(-57.0,
+ t448 = RATIONAL(51.0,1225.0)+t357+t332+t356+t468*RATIONAL(-289.0,5880.0);
+ t447 = t380*t418+t364+t381*t417+t356+t464;
+ t440 = RATIONAL(1.0,98.0)*t441+RATIONAL(246.0,1225.0);
+ t360 = t404*t415;
+ t439 = t360+t370;
+ t438 = t358+t360;
+ t436 = RATIONAL(-1.0,49.0)*t441+RATIONAL(96.0,1225.0);
+ t390 = RATIONAL(83.0,735.0);
+ t408 = RATIONAL(1.0,20.0);
+ t363 = t408*t413;
+ t392 = RATIONAL(-53.0,196.0);
+ t434 = t360+t390*t418+t363+t392*t417+t440;
+ t386 = RATIONAL(-41.0,196.0);
+ t374 = RATIONAL(191.0,5880.0);
+ t433 = t386*t417+t374*t418+t363+t356+t436;
+ t372 = t408*t415;
+ t432 = t357+t374*t417+t386*t418+t372+t436;
+ t431 = t364+t392*t418+t390*t417+t372+t440;
+ t317 = t402*t443;
+ t430 = t384*y+t317+t373+t448;
+ t429 = t395*y+t400*t443+t438+t450;
+ t326 = t409*t443;
+ t428 = t389*y+t367+t326+t448;
+ t427 = t397*y+t401*t443+t439+t450;
+ t385 = RATIONAL(11.0,105.0);
+ t426 = t385*x+t315+t362+t447;
+ t425 = t385*y+t319+t439+t449;
+ t388 = RATIONAL(-11.0,105.0);
+ t424 = t388*y+t309+t438+t449;
+ t423 = t388*x+t366+t325+t447;
+ t394 = RATIONAL(17.0,420.0);
+ t393 = RATIONAL(17.0,105.0);
+ t391 = RATIONAL(-17.0,105.0);
+ t387 = RATIONAL(-17.0,420.0);
+ coeff_I_m2_m2 = t428+t460+t467;
+ coeff_I_m1_m2 = t423+t456+t457;
+ coeff_I_0_m2 = t317+t367+t387*y+t433;
+ coeff_I_p1_m2 = t426+t451+t456;
+ coeff_I_p2_m2 = t428+t459+t465;
+ coeff_I_m2_m1 = t424+t454+t455;
+ coeff_I_m1_m1 = t429+t462+t466;
+ coeff_I_0_m1 = t319+t391*y+t358+t434;
+ coeff_I_p1_m1 = t429+t461+t463;
+ coeff_I_p2_m1 = t424+t452+t453;
+ coeff_I_m2_0 = t368+t312+t387*x+t432;
+ coeff_I_m1_0 = t366+t315+t391*x+t431;
+ coeff_I_0_0 = RATIONAL(541.0,1225.0)+t363+t372+t332+t468*RATIONAL(-57.0,
196.0);
- coeff_I_p1_0 = t389*x+t327+t366+t433;
- coeff_I_p2_0 = t325+t364+t388*x+t432;
- coeff_I_m2_p1 = t423+t461+t462;
- coeff_I_m1_p1 = t429+t449+t464;
- coeff_I_0_p1 = t359+t389*y+t311+t431;
- coeff_I_p1_p1 = t429+t460+t465;
- coeff_I_p2_p1 = t425+t450+t462;
- coeff_I_m2_p2 = t427+t453+t455;
- coeff_I_m1_p2 = t426+t456+t458;
- coeff_I_0_p2 = t313+t388*y+t367+t434;
- coeff_I_p1_p2 = t426+t457+t459;
- coeff_I_p2_p2 = t430+t452+t455;
+ coeff_I_p1_0 = t325+t362+t393*x+t431;
+ coeff_I_p2_0 = t394*x+t316+t361+t432;
+ coeff_I_m2_p1 = t425+t452+t454;
+ coeff_I_m1_p1 = t427+t462+t463;
+ coeff_I_0_p1 = t309+t370+t393*y+t434;
+ coeff_I_p1_p1 = t427+t461+t466;
+ coeff_I_p2_p1 = t425+t453+t455;
+ coeff_I_m2_p2 = t430+t460+t465;
+ coeff_I_m1_p2 = t423+t451+t458;
+ coeff_I_0_p2 = t394*y+t326+t373+t433;
+ coeff_I_p1_p2 = t426+t457+t458;
+ coeff_I_p2_p2 = t430+t459+t467;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c
index 83ba345..e56ee74 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c
@@ -1,48 +1,54 @@
-fp t531,
- t534,
- t535,
- t513,
- t570,
- t508,
- t569,
+fp t538,
+ t507,
+ t539,
t537,
+ t559,
+ t573,
+ t560,
+ t572,
+ t514,
+ t554,
+ t571,
+ t526,
t536,
- t502,
- t552,
+ t516,
+ t570,
+ t531,
+ t517,
+ t569,
+ t481,
t568,
- t512,
- t555,
+ t561,
+ t494,
t567,
- t556,
+ t495,
t566,
- t509,
- t533,
- t510,
- t565,
- t530,
+ t525,
t506,
- t507,
+ t565,
+ t510,
t564,
+ t511,
t563,
- t471,
t562,
- t561,
- t523,
- t511,
- t560,
- t528,
- t504,
- t559,
+ t475,
+ t527,
+ t509,
t558,
+ t513,
t557,
+ t529,
+ t512,
+ t556,
+ t555,
+ t518,
t478,
- t554,
t553,
- t516,
- t476,
+ t479,
+ t552,
t551,
- t483,
t550,
+ t490,
t549,
t548,
t547,
@@ -51,97 +57,99 @@ fp t531,
t544,
t543,
t542,
- t541,
- t540,
- t484,
+ t486,
+ t482,
t477,
- t475,
- t473,
- t469,
- t467;
- t531 = RATIONAL(1.0,30.0);
- t534 = x*x;
- t535 = t534*x;
- t513 = t531*t535;
- t570 = t513+RATIONAL(71.0,2940.0)*x;
- t508 = RATIONAL(-2.0,15.0)*t535;
- t569 = t508+RATIONAL(226.0,735.0)*x;
- t537 = y*y;
- t536 = y*t537;
- t502 = RATIONAL(1.0,20.0)*t534;
- t552 = RATIONAL(1.0,60.0)*t536+y*t502;
- t568 = t552+RATIONAL(-37.0,300.0)*y;
- t512 = RATIONAL(-1.0,20.0)*t534;
- t555 = y*t512+RATIONAL(-1.0,60.0)*t536;
- t567 = t555+RATIONAL(37.0,300.0)*y;
- t556 = x*t537;
- t566 = t508+RATIONAL(181.0,735.0)*x+RATIONAL(1.0,98.0)*t556;
- t509 = RATIONAL(1.0,10.0)*t534;
- t533 = RATIONAL(-1.0,70.0);
- t510 = t533*t537;
- t565 = t509+t510+RATIONAL(-11.0,105.0);
- t530 = RATIONAL(1.0,70.0);
- t506 = t530*t537;
- t507 = RATIONAL(-1.0,10.0)*t534;
- t564 = t506+t507+RATIONAL(11.0,105.0);
- t563 = t510+RATIONAL(1.0,84.0)+t502;
- t471 = RATIONAL(2.0,49.0)*t556;
- t562 = t513+t471+RATIONAL(-289.0,2940.0)*x;
- t561 = t506+t512+RATIONAL(-1.0,84.0);
- t523 = RATIONAL(-1.0,35.0);
- t511 = t523*t537;
- t560 = t511+t512+RATIONAL(31.0,420.0);
- t528 = RATIONAL(1.0,35.0);
- t504 = t528*t537;
- t559 = RATIONAL(-31.0,420.0)+t502+t504;
- t558 = x*y;
- t557 = y*t534;
- t478 = RATIONAL(-1.0,49.0)*t556;
- t554 = RATIONAL(-21.0,200.0)*y+t478;
- t553 = RATIONAL(21.0,200.0)*y+t478;
- t516 = RATIONAL(1.0,5.0)*t535;
- t476 = RATIONAL(-2.0,49.0)*t556;
- t551 = t516+t476+RATIONAL(-41.0,98.0)*x;
- t483 = RATIONAL(1.0,49.0)*t556;
- t550 = t516+t483+RATIONAL(-53.0,98.0)*x;
- t549 = t508+RATIONAL(166.0,735.0)*x+t483;
- t548 = t513+t476+RATIONAL(191.0,2940.0)*x;
- t547 = t509+RATIONAL(-31.0,210.0)+RATIONAL(1.0,140.0)*t537+t566;
- t546 = RATIONAL(11.0,150.0)*y+t555+t562;
- t545 = y*t509+RATIONAL(-1.0,120.0)*t536+t554+t569;
- t544 = t507+RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t537+t566;
- t543 = y*t507+RATIONAL(1.0,120.0)*t536+t553+t569;
- t542 = t531*t536+RATIONAL(-1.0,40.0)*t557+t554+t570;
- t541 = RATIONAL(-11.0,150.0)*y+t552+t562;
- t540 = RATIONAL(1.0,40.0)*t557+RATIONAL(-1.0,30.0)*t536+t553+t570;
- t484 = t530*t558;
- t477 = t528*t558;
- t475 = RATIONAL(2.0,35.0)*t558;
- t473 = t533*t558;
- t469 = RATIONAL(-2.0,35.0)*t558;
- t467 = t523*t558;
- coeff_dx_m2_m2 = t469+t541+t560;
- coeff_dx_m1_m2 = t477+t543+t565;
- coeff_dx_0_m2 = t475+t551;
- coeff_dx_p1_m2 = t477+t545+t564;
- coeff_dx_p2_m2 = t469+t546+t559;
- coeff_dx_m2_m1 = t467+t540+t561;
- coeff_dx_m1_m1 = t484+t547+t567;
- coeff_dx_0_m1 = t477+t550;
- coeff_dx_p1_m1 = t484+t544+t568;
- coeff_dx_p2_m1 = t467+t542+t563;
- coeff_dx_m2_0 = t504+t512+RATIONAL(-17.0,420.0)+t548;
- coeff_dx_m1_0 = RATIONAL(-17.0,105.0)+t509+t506+t549;
- coeff_dx_0_0 = RATIONAL(-57.0,98.0)*x+t516+t471;
- coeff_dx_p1_0 = t510+RATIONAL(17.0,105.0)+t507+t549;
- coeff_dx_p2_0 = t511+RATIONAL(17.0,420.0)+t502+t548;
- coeff_dx_m2_p1 = t477+t542+t561;
- coeff_dx_m1_p1 = t473+t547+t568;
- coeff_dx_0_p1 = t467+t550;
- coeff_dx_p1_p1 = t473+t544+t567;
- coeff_dx_p2_p1 = t477+t540+t563;
- coeff_dx_m2_p2 = t475+t546+t560;
- coeff_dx_m1_p2 = t467+t545+t565;
- coeff_dx_0_p2 = t469+t551;
- coeff_dx_p1_p2 = t467+t543+t564;
- coeff_dx_p2_p2 = t475+t541+t559;
+ t474,
+ t471,
+ t469;
+ t538 = x*x;
+ t507 = RATIONAL(1.0,20.0)*t538;
+ t539 = y*y;
+ t537 = y*t539;
+ t559 = y*t507+RATIONAL(1.0,60.0)*t537;
+ t573 = RATIONAL(-37.0,300.0)*y+t559;
+ t560 = x*t539;
+ t572 = RATIONAL(1.0,98.0)*t560+RATIONAL(181.0,735.0)*x;
+ t514 = RATIONAL(-1.0,20.0)*t538;
+ t554 = y*t514+RATIONAL(-1.0,60.0)*t537;
+ t571 = RATIONAL(37.0,300.0)*y+t554;
+ t526 = RATIONAL(1.0,30.0);
+ t536 = x*t538;
+ t516 = t526*t536;
+ t570 = RATIONAL(71.0,2940.0)*x+t516;
+ t531 = RATIONAL(-1.0,35.0);
+ t517 = t531*t539;
+ t569 = RATIONAL(31.0,420.0)+t514+t517;
+ t481 = RATIONAL(2.0,49.0)*t560;
+ t568 = RATIONAL(-289.0,2940.0)*x+t481+t516;
+ t561 = y*t538;
+ t494 = RATIONAL(-21.0,200.0)*y;
+ t567 = RATIONAL(-1.0,40.0)*t561+t494+t526*t537;
+ t495 = RATIONAL(21.0,200.0)*y;
+ t566 = RATIONAL(1.0,40.0)*t561+t495+RATIONAL(-1.0,30.0)*t537;
+ t525 = RATIONAL(1.0,35.0);
+ t506 = t525*t539;
+ t565 = RATIONAL(-31.0,420.0)+t507+t506;
+ t510 = RATIONAL(-1.0,10.0)*t538;
+ t564 = y*t510+t495+RATIONAL(1.0,120.0)*t537;
+ t511 = RATIONAL(1.0,10.0)*t538;
+ t563 = y*t511+t494+RATIONAL(-1.0,120.0)*t537;
+ t562 = x*y;
+ t475 = RATIONAL(-1.0,49.0)*t560;
+ t527 = RATIONAL(1.0,70.0);
+ t509 = t527*t539;
+ t558 = t475+t509;
+ t513 = RATIONAL(-2.0,15.0)*t536;
+ t557 = t510+t513;
+ t529 = RATIONAL(-1.0,70.0);
+ t512 = t529*t539;
+ t556 = t475+t512;
+ t555 = t511+t513;
+ t518 = RATIONAL(1.0,5.0)*t536;
+ t478 = RATIONAL(-2.0,49.0)*t560;
+ t553 = t518+RATIONAL(-41.0,98.0)*x+t478;
+ t479 = RATIONAL(1.0,49.0)*t560;
+ t552 = t518+RATIONAL(-53.0,98.0)*x+t479;
+ t551 = RATIONAL(191.0,2940.0)*x+t478+t516;
+ t550 = RATIONAL(166.0,735.0)*x+t479+t513;
+ t490 = RATIONAL(226.0,735.0)*x;
+ t549 = RATIONAL(11.0,105.0)+t490+t557+t558;
+ t548 = t490+RATIONAL(-11.0,105.0)+t555+t556;
+ t547 = RATIONAL(-31.0,210.0)+RATIONAL(1.0,140.0)*t539+t555+t572;
+ t546 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t539+t557+t572;
+ t545 = RATIONAL(1.0,84.0)+t507+t556+t570;
+ t544 = RATIONAL(-11.0,150.0)*y+t559+t568;
+ t543 = RATIONAL(-1.0,84.0)+t514+t558+t570;
+ t542 = RATIONAL(11.0,150.0)*y+t554+t568;
+ t486 = t525*t562;
+ t482 = RATIONAL(-2.0,35.0)*t562;
+ t477 = RATIONAL(2.0,35.0)*t562;
+ t474 = t531*t562;
+ t471 = t527*t562;
+ t469 = t529*t562;
+ coeff_dx_m2_m2 = t482+t544+t569;
+ coeff_dx_m1_m2 = t486+t548+t564;
+ coeff_dx_0_m2 = t477+t553;
+ coeff_dx_p1_m2 = t486+t549+t563;
+ coeff_dx_p2_m2 = t482+t542+t565;
+ coeff_dx_m2_m1 = t474+t543+t566;
+ coeff_dx_m1_m1 = t471+t547+t571;
+ coeff_dx_0_m1 = t486+t552;
+ coeff_dx_p1_m1 = t471+t546+t573;
+ coeff_dx_p2_m1 = t474+t545+t567;
+ coeff_dx_m2_0 = t514+RATIONAL(-17.0,420.0)+t506+t551;
+ coeff_dx_m1_0 = t511+RATIONAL(-17.0,105.0)+t509+t550;
+ coeff_dx_0_0 = t481+t518+RATIONAL(-57.0,98.0)*x;
+ coeff_dx_p1_0 = RATIONAL(17.0,105.0)+t510+t512+t550;
+ coeff_dx_p2_0 = t507+RATIONAL(17.0,420.0)+t517+t551;
+ coeff_dx_m2_p1 = t486+t543+t567;
+ coeff_dx_m1_p1 = t469+t547+t573;
+ coeff_dx_0_p1 = t474+t552;
+ coeff_dx_p1_p1 = t469+t546+t571;
+ coeff_dx_p2_p1 = t486+t545+t566;
+ coeff_dx_m2_p2 = t477+t542+t569;
+ coeff_dx_m1_p2 = t474+t548+t563;
+ coeff_dx_0_p2 = t482+t553;
+ coeff_dx_p1_p2 = t474+t549+t564;
+ coeff_dx_p2_p2 = t477+t544+t565;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c
index 385d8a7..3ac1e8c 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c
@@ -1,111 +1,107 @@
-fp t712,
- t692,
- t707,
- t711,
- t693,
- t729,
+fp t713,
+ t683,
+ t708,
+ t712,
t687,
+ t730,
+ t696,
+ t729,
t728,
- t697,
+ t690,
t727,
t726,
- t725,
+ t697,
t686,
- t694,
- t695,
t698,
- t689,
+ t691,
+ t693,
+ t725,
+ t699,
t724,
- t681,
+ t682,
t723,
- t696,
t722,
+ t695,
t721,
- t688,
t720,
- t691,
+ t689,
t719,
- t685,
t718,
t717,
t716,
- t690,
+ t678,
t715,
+ t676,
t714,
- t713,
- t684,
- t683,
+ t694,
+ t685,
+ t681,
t680,
t679,
- t678,
- t677,
- t676,
- t675;
- t712 = y*y;
- t692 = RATIONAL(2.0,49.0)*t712;
- t707 = RATIONAL(1.0,10.0);
- t711 = x*x;
- t693 = t707*t711;
- t729 = t692+t693+RATIONAL(-289.0,2940.0);
- t687 = RATIONAL(-1.0,49.0)*t712;
- t728 = t693+RATIONAL(71.0,2940.0)+t687;
- t697 = RATIONAL(-2.0,5.0)*t711;
- t727 = t697+RATIONAL(226.0,735.0)+t687;
- t726 = t697+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t712;
- t725 = x*y;
- t686 = t707*x;
- t694 = RATIONAL(-1.0,10.0)*x;
- t695 = RATIONAL(-1.0,5.0)*x;
- t698 = RATIONAL(1.0,5.0)*x;
- t689 = RATIONAL(-2.0,49.0)*t712;
- t724 = t693+RATIONAL(191.0,2940.0)+t689;
- t681 = RATIONAL(1.0,49.0)*t712;
- t723 = t697+RATIONAL(166.0,735.0)+t681;
- t696 = RATIONAL(3.0,5.0)*t711;
- t722 = t696+RATIONAL(-53.0,98.0)+t681;
- t721 = t696+RATIONAL(-41.0,98.0)+t689;
- t688 = RATIONAL(2.0,35.0)*y;
- t720 = t688+t729;
- t691 = RATIONAL(-2.0,35.0)*y;
- t719 = t691+t729;
- t685 = RATIONAL(-1.0,35.0)*y;
- t718 = t685+t727;
- t717 = t686+t728;
- t716 = t695+t726;
- t690 = RATIONAL(1.0,35.0)*y;
- t715 = t690+t727;
- t714 = t698+t726;
- t713 = t694+t728;
- t684 = RATIONAL(1.0,70.0)*y;
- t683 = RATIONAL(-1.0,70.0)*y;
- t680 = y*t695;
- t679 = y*t698;
- t678 = RATIONAL(-1.0,20.0)*t725;
- t677 = y*t686;
- t676 = y*t694;
- t675 = RATIONAL(1.0,20.0)*t725;
- coeff_dxx_m2_m2 = t677+t694+t719;
- coeff_dxx_m1_m2 = t680+t698+t715;
- coeff_dxx_0_m2 = t688+t721;
- coeff_dxx_p1_m2 = t679+t695+t715;
- coeff_dxx_p2_m2 = t686+t676+t719;
- coeff_dxx_m2_m1 = t685+t675+t713;
- coeff_dxx_m1_m1 = t676+t684+t714;
- coeff_dxx_0_m1 = t690+t722;
- coeff_dxx_p1_m1 = t684+t677+t716;
- coeff_dxx_p2_m1 = t685+t678+t717;
- coeff_dxx_m2_0 = t694+t724;
- coeff_dxx_m1_0 = t698+t723;
- coeff_dxx_0_0 = t692+RATIONAL(-57.0,98.0)+t696;
- coeff_dxx_p1_0 = t695+t723;
- coeff_dxx_p2_0 = t686+t724;
- coeff_dxx_m2_p1 = t690+t678+t713;
- coeff_dxx_m1_p1 = t677+t683+t714;
- coeff_dxx_0_p1 = t685+t722;
- coeff_dxx_p1_p1 = t683+t676+t716;
- coeff_dxx_p2_p1 = t690+t675+t717;
- coeff_dxx_m2_p2 = t676+t694+t720;
- coeff_dxx_m1_p2 = t698+t679+t718;
- coeff_dxx_0_p2 = t691+t721;
- coeff_dxx_p1_p2 = t695+t680+t718;
- coeff_dxx_p2_p2 = t686+t677+t720;
+ t677;
+ t713 = y*y;
+ t683 = RATIONAL(-1.0,49.0)*t713;
+ t708 = RATIONAL(1.0,10.0);
+ t712 = x*x;
+ t687 = t708*t712;
+ t730 = RATIONAL(71.0,2940.0)+t683+t687;
+ t696 = RATIONAL(-2.0,5.0)*t712;
+ t729 = t696+RATIONAL(226.0,735.0)+t683;
+ t728 = t696+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t713;
+ t690 = RATIONAL(2.0,49.0)*t713;
+ t727 = t690+RATIONAL(-289.0,2940.0)+t687;
+ t726 = x*y;
+ t697 = RATIONAL(1.0,5.0)*x;
+ t686 = RATIONAL(-1.0,10.0)*x;
+ t698 = RATIONAL(-1.0,5.0)*x;
+ t691 = t708*x;
+ t693 = RATIONAL(-2.0,49.0)*t713;
+ t725 = t693+RATIONAL(191.0,2940.0)+t687;
+ t699 = RATIONAL(3.0,5.0)*t712;
+ t724 = t693+t699+RATIONAL(-41.0,98.0);
+ t682 = RATIONAL(1.0,49.0)*t713;
+ t723 = t699+RATIONAL(-53.0,98.0)+t682;
+ t722 = t696+RATIONAL(166.0,735.0)+t682;
+ t695 = RATIONAL(1.0,35.0)*y;
+ t721 = t695+t730;
+ t720 = RATIONAL(-1.0,70.0)*y+t728;
+ t689 = RATIONAL(-1.0,35.0)*y;
+ t719 = t689+t730;
+ t718 = t698+t729;
+ t717 = RATIONAL(1.0,70.0)*y+t728;
+ t716 = t697+t729;
+ t678 = y*t686;
+ t715 = t678+t727;
+ t676 = y*t691;
+ t714 = t676+t727;
+ t694 = RATIONAL(2.0,35.0)*y;
+ t685 = RATIONAL(-2.0,35.0)*y;
+ t681 = y*t697;
+ t680 = y*t698;
+ t679 = RATIONAL(1.0,20.0)*t726;
+ t677 = RATIONAL(-1.0,20.0)*t726;
+ coeff_dxx_m2_m2 = t686+t685+t714;
+ coeff_dxx_m1_m2 = t680+t695+t716;
+ coeff_dxx_0_m2 = t694+t724;
+ coeff_dxx_p1_m2 = t695+t681+t718;
+ coeff_dxx_p2_m2 = t691+t685+t715;
+ coeff_dxx_m2_m1 = t679+t686+t719;
+ coeff_dxx_m1_m1 = t697+t678+t717;
+ coeff_dxx_0_m1 = t695+t723;
+ coeff_dxx_p1_m1 = t698+t676+t717;
+ coeff_dxx_p2_m1 = t691+t677+t719;
+ coeff_dxx_m2_0 = t686+t725;
+ coeff_dxx_m1_0 = t697+t722;
+ coeff_dxx_0_0 = RATIONAL(-57.0,98.0)+t690+t699;
+ coeff_dxx_p1_0 = t698+t722;
+ coeff_dxx_p2_0 = t691+t725;
+ coeff_dxx_m2_p1 = t677+t686+t721;
+ coeff_dxx_m1_p1 = t697+t676+t720;
+ coeff_dxx_0_p1 = t689+t723;
+ coeff_dxx_p1_p1 = t678+t698+t720;
+ coeff_dxx_p2_p1 = t691+t679+t721;
+ coeff_dxx_m2_p2 = t694+t686+t715;
+ coeff_dxx_m1_p2 = t689+t681+t716;
+ coeff_dxx_0_p2 = t685+t724;
+ coeff_dxx_p1_p2 = t680+t689+t718;
+ coeff_dxx_p2_p2 = t694+t691+t714;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c
index 9388906..0505d86 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c
@@ -1,119 +1,119 @@
-fp t778,
- t777,
- t796,
- t766,
- t735,
+fp t779,
+ t778,
+ t797,
t772,
- t755,
+ t774,
+ t796,
t795,
- t751,
- t748,
+ t770,
+ t739,
+ t766,
+ t752,
t794,
+ t742,
+ t759,
t793,
t792,
- t791,
t733,
+ t791,
t790,
t789,
t788,
t787,
- t765,
- t770,
+ t735,
t786,
- t734,
+ t732,
t785,
- t767,
- t771,
+ t773,
+ t777,
t784,
t783,
t782,
- t730,
t781,
t780,
- t779,
t776,
- t775,
- t773,
- t769,
- t756,
- t754,
+ t771,
+ t768,
+ t767,
+ t758,
+ t757,
+ t755,
t750,
- t749,
- t746,
- t745,
t743,
t738,
- t732,
+ t737,
+ t736,
+ t734,
t731;
- t778 = y*y;
- t777 = x*x;
- t796 = t778+t777;
- t766 = RATIONAL(-1.0,35.0);
- t735 = t766*x;
- t772 = RATIONAL(1.0,35.0);
- t755 = t772*y;
- t795 = t735+t755;
- t751 = t766*y;
- t748 = t772*x;
- t794 = t751+t748;
- t793 = t735+t751;
- t792 = t755+t748;
- t791 = x*y;
- t733 = RATIONAL(-2.0,49.0)*t791;
- t790 = t733+RATIONAL(-21.0,200.0);
- t789 = t733+RATIONAL(21.0,200.0);
- t788 = t796*RATIONAL(-1.0,20.0);
- t787 = t796*RATIONAL(1.0,20.0);
- t765 = RATIONAL(-1.0,40.0);
- t770 = RATIONAL(1.0,10.0);
- t786 = t765*t777+t770*t778+t790;
- t734 = RATIONAL(4.0,49.0)*t791;
- t785 = t734+RATIONAL(11.0,150.0)+t788;
- t767 = RATIONAL(1.0,40.0);
- t771 = RATIONAL(-1.0,10.0);
- t784 = t771*t778+t767*t777+t789;
- t783 = t771*t777+t767*t778+t789;
- t782 = t734+RATIONAL(-11.0,150.0)+t787;
- t730 = RATIONAL(1.0,49.0)*t791;
- t781 = RATIONAL(37.0,300.0)+t730+t788;
- t780 = RATIONAL(-37.0,300.0)+t730+t787;
- t779 = t765*t778+t770*t777+t790;
+ t779 = y*y;
+ t778 = x*x;
+ t797 = t778+t779;
+ t772 = RATIONAL(1.0,10.0);
+ t774 = RATIONAL(-1.0,40.0);
+ t796 = t772*t778+t774*t779;
+ t795 = t774*t778+t772*t779;
+ t770 = RATIONAL(-1.0,35.0);
+ t739 = t770*y;
+ t766 = RATIONAL(1.0,35.0);
+ t752 = t766*x;
+ t794 = t739+t752;
+ t742 = t770*x;
+ t759 = t766*y;
+ t793 = t742+t759;
+ t792 = x*y;
+ t733 = RATIONAL(-2.0,49.0)*t792;
+ t791 = t733+RATIONAL(-21.0,200.0);
+ t790 = t797*RATIONAL(-1.0,20.0);
+ t789 = t797*RATIONAL(1.0,20.0);
+ t788 = t733+RATIONAL(21.0,200.0);
+ t787 = t752+t759+t791;
+ t735 = RATIONAL(4.0,49.0)*t792;
+ t786 = t735+RATIONAL(11.0,150.0)+t790;
+ t732 = RATIONAL(1.0,49.0)*t792;
+ t785 = t732+RATIONAL(37.0,300.0)+t790;
+ t773 = RATIONAL(-1.0,10.0);
+ t777 = RATIONAL(1.0,40.0);
+ t784 = t773*t779+t777*t778+t788;
+ t783 = t735+RATIONAL(-11.0,150.0)+t789;
+ t782 = t732+RATIONAL(-37.0,300.0)+t789;
+ t781 = t773*t778+t777*t779+t788;
+ t780 = t739+t742+t791;
t776 = RATIONAL(-2.0,35.0);
- t775 = RATIONAL(-1.0,70.0);
- t773 = RATIONAL(1.0,70.0);
- t769 = RATIONAL(2.0,35.0);
- t756 = t776*y;
- t754 = t769*y;
- t750 = t776*x;
- t749 = t775*x;
- t746 = t773*y;
- t745 = t775*y;
- t743 = t769*x;
- t738 = t773*x;
- t732 = RATIONAL(-4.0,49.0)*t791;
- t731 = RATIONAL(2.0,49.0)*t791;
- coeff_dxy_m2_m2 = t756+t750+t782;
- coeff_dxy_m1_m2 = t783+t794;
- coeff_dxy_0_m2 = t743+t732;
- coeff_dxy_p1_m2 = t779+t792;
- coeff_dxy_p2_m2 = t750+t754+t785;
- coeff_dxy_m2_m1 = t784+t795;
- coeff_dxy_m1_m1 = t738+t746+t781;
- coeff_dxy_0_m1 = t748+t731;
- coeff_dxy_p1_m1 = t745+t738+t780;
- coeff_dxy_p2_m1 = t786+t793;
- coeff_dxy_m2_0 = t754+t732;
- coeff_dxy_m1_0 = t755+t731;
- coeff_dxy_0_0 = t734;
- coeff_dxy_p1_0 = t751+t731;
- coeff_dxy_p2_0 = t732+t756;
- coeff_dxy_m2_p1 = t786+t792;
- coeff_dxy_m1_p1 = t749+t746+t780;
- coeff_dxy_0_p1 = t735+t731;
- coeff_dxy_p1_p1 = t749+t745+t781;
+ t771 = RATIONAL(2.0,35.0);
+ t768 = RATIONAL(-1.0,70.0);
+ t767 = RATIONAL(1.0,70.0);
+ t758 = t771*y;
+ t757 = t768*y;
+ t755 = t776*x;
+ t750 = t771*x;
+ t743 = t767*x;
+ t738 = t767*y;
+ t737 = t776*y;
+ t736 = t768*x;
+ t734 = RATIONAL(-4.0,49.0)*t792;
+ t731 = RATIONAL(2.0,49.0)*t792;
+ coeff_dxy_m2_m2 = t755+t737+t783;
+ coeff_dxy_m1_m2 = t781+t794;
+ coeff_dxy_0_m2 = t750+t734;
+ coeff_dxy_p1_m2 = t787+t796;
+ coeff_dxy_p2_m2 = t758+t755+t786;
+ coeff_dxy_m2_m1 = t784+t793;
+ coeff_dxy_m1_m1 = t738+t743+t785;
+ coeff_dxy_0_m1 = t752+t731;
+ coeff_dxy_p1_m1 = t757+t743+t782;
+ coeff_dxy_p2_m1 = t780+t795;
+ coeff_dxy_m2_0 = t734+t758;
+ coeff_dxy_m1_0 = t731+t759;
+ coeff_dxy_0_0 = t735;
+ coeff_dxy_p1_0 = t731+t739;
+ coeff_dxy_p2_0 = t734+t737;
+ coeff_dxy_m2_p1 = t787+t795;
+ coeff_dxy_m1_p1 = t738+t736+t782;
+ coeff_dxy_0_p1 = t731+t742;
+ coeff_dxy_p1_p1 = t757+t736+t785;
coeff_dxy_p2_p1 = t784+t794;
- coeff_dxy_m2_p2 = t743+t756+t785;
- coeff_dxy_m1_p2 = t779+t793;
- coeff_dxy_0_p2 = t750+t732;
- coeff_dxy_p1_p2 = t783+t795;
- coeff_dxy_p2_p2 = t754+t743+t782;
+ coeff_dxy_m2_p2 = t750+t737+t786;
+ coeff_dxy_m1_p2 = t780+t796;
+ coeff_dxy_0_p2 = t755+t734;
+ coeff_dxy_p1_p2 = t781+t793;
+ coeff_dxy_p2_p2 = t758+t750+t783;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c
index 03c68c3..f4baa3a 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c
@@ -1,50 +1,51 @@
-fp t641,
- t611,
- t638,
- t639,
- t657,
+fp t644,
+ t609,
+ t630,
+ t643,
+ t616,
+ t675,
+ t636,
+ t621,
+ t622,
t674,
- t608,
- t656,
- t673,
t661,
- t589,
+ t601,
+ t641,
+ t673,
+ t631,
+ t600,
t672,
+ t613,
t671,
- t635,
- t640,
- t609,
+ t662,
+ t642,
+ t611,
t670,
- t574,
+ t587,
t669,
- t614,
- t598,
+ t619,
+ t590,
t668,
- t632,
- t615,
- t667,
t634,
- t619,
+ t612,
+ t667,
t666,
- t617,
- t596,
+ t610,
t665,
- t627,
- t606,
+ t632,
+ t618,
t664,
- t637,
- t616,
t663,
- t662,
t660,
- t618,
t659,
+ t574,
t658,
t586,
+ t657,
+ t623,
+ t656,
t655,
- t578,
t654,
- t620,
t653,
t652,
t651,
@@ -52,100 +53,95 @@ fp t641,
t649,
t648,
t647,
- t646,
- t645,
- t644,
- t588,
+ t592,
+ t584,
t583,
t582,
- t579,
- t576,
- t572;
- t641 = y*y;
- t611 = RATIONAL(1.0,20.0)*t641;
- t638 = x*x;
- t639 = t638*x;
- t657 = RATIONAL(1.0,60.0)*t639+x*t611;
- t674 = t657+RATIONAL(-37.0,300.0)*x;
- t608 = RATIONAL(-1.0,20.0)*t641;
- t656 = RATIONAL(-1.0,60.0)*t639+x*t608;
- t673 = t656+RATIONAL(37.0,300.0)*x;
- t661 = y*t638;
- t589 = RATIONAL(-1.0,49.0)*t661;
- t672 = t589+RATIONAL(226.0,735.0)*y;
- t671 = RATIONAL(1.0,98.0)*t661+RATIONAL(181.0,735.0)*y;
- t635 = RATIONAL(1.0,30.0);
- t640 = y*t641;
- t609 = t635*t640;
- t670 = t589+t609+RATIONAL(71.0,2940.0)*y;
- t574 = RATIONAL(2.0,49.0)*t661;
- t669 = RATIONAL(-289.0,2940.0)*y+t574+t609;
- t614 = RATIONAL(-1.0,10.0)*t641;
- t598 = RATIONAL(21.0,200.0)*x;
- t668 = x*t614+RATIONAL(1.0,120.0)*t639+t598;
- t632 = RATIONAL(1.0,35.0);
- t615 = t632*t638;
- t667 = t611+RATIONAL(-31.0,420.0)+t615;
- t634 = RATIONAL(1.0,70.0);
- t619 = t634*t638;
- t666 = t619+RATIONAL(-1.0,84.0)+t608;
- t617 = RATIONAL(1.0,10.0)*t641;
- t596 = RATIONAL(-21.0,200.0)*x;
- t665 = RATIONAL(-1.0,120.0)*t639+x*t617+t596;
- t627 = RATIONAL(-1.0,35.0);
- t606 = t627*t638;
- t664 = RATIONAL(31.0,420.0)+t606+t608;
- t637 = RATIONAL(-1.0,70.0);
- t616 = t637*t638;
- t663 = t611+t616+RATIONAL(1.0,84.0);
- t662 = x*y;
- t660 = x*t641;
- t618 = RATIONAL(-2.0,15.0)*t640;
- t659 = t617+t618;
- t658 = t618+t614;
- t586 = RATIONAL(-2.0,49.0)*t661;
- t655 = t586+t609+RATIONAL(191.0,2940.0)*y;
- t578 = RATIONAL(1.0,49.0)*t661;
- t654 = t618+RATIONAL(166.0,735.0)*y+t578;
- t620 = RATIONAL(1.0,5.0)*t640;
- t653 = RATIONAL(-53.0,98.0)*y+t620+t578;
- t652 = t620+RATIONAL(-41.0,98.0)*y+t586;
- t651 = RATIONAL(1.0,140.0)*t638+RATIONAL(-31.0,210.0)+t659+t671;
- t650 = t619+RATIONAL(11.0,105.0)+t658+t672;
- t649 = RATIONAL(-11.0,150.0)*x+t657+t669;
- t648 = RATIONAL(11.0,150.0)*x+t656+t669;
- t647 = t616+RATIONAL(-11.0,105.0)+t659+t672;
- t646 = t635*t639+RATIONAL(-1.0,40.0)*t660+t596+t670;
- t645 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t638+t658+t671;
- t644 = RATIONAL(-1.0,30.0)*t639+RATIONAL(1.0,40.0)*t660+t598+t670;
- t588 = t627*t662;
- t583 = t632*t662;
- t582 = RATIONAL(2.0,35.0)*t662;
- t579 = t634*t662;
- t576 = t637*t662;
- t572 = RATIONAL(-2.0,35.0)*t662;
- coeff_dy_m2_m2 = t572+t649+t664;
- coeff_dy_m1_m2 = t588+t644+t666;
- coeff_dy_0_m2 = RATIONAL(-17.0,420.0)+t615+t608+t655;
- coeff_dy_p1_m2 = t583+t646+t666;
- coeff_dy_p2_m2 = t582+t648+t664;
- coeff_dy_m2_m1 = t583+t647+t668;
- coeff_dy_m1_m1 = t579+t651+t673;
- coeff_dy_0_m1 = RATIONAL(-17.0,105.0)+t617+t619+t654;
- coeff_dy_p1_m1 = t576+t651+t674;
- coeff_dy_p2_m1 = t588+t647+t665;
- coeff_dy_m2_0 = t582+t652;
- coeff_dy_m1_0 = t583+t653;
- coeff_dy_0_0 = t574+RATIONAL(-57.0,98.0)*y+t620;
- coeff_dy_p1_0 = t588+t653;
- coeff_dy_p2_0 = t572+t652;
- coeff_dy_m2_p1 = t583+t650+t665;
- coeff_dy_m1_p1 = t579+t645+t674;
- coeff_dy_0_p1 = RATIONAL(17.0,105.0)+t616+t614+t654;
- coeff_dy_p1_p1 = t576+t645+t673;
- coeff_dy_p2_p1 = t588+t650+t668;
- coeff_dy_m2_p2 = t572+t648+t667;
- coeff_dy_m1_p2 = t588+t646+t663;
- coeff_dy_0_p2 = t611+RATIONAL(17.0,420.0)+t606+t655;
- coeff_dy_p1_p2 = t583+t644+t663;
- coeff_dy_p2_p2 = t582+t649+t667;
+ t578,
+ t576;
+ t644 = y*y;
+ t609 = RATIONAL(1.0,20.0)*t644;
+ t630 = RATIONAL(1.0,35.0);
+ t643 = x*x;
+ t616 = t630*t643;
+ t675 = t609+t616+RATIONAL(-31.0,420.0);
+ t636 = RATIONAL(-1.0,35.0);
+ t621 = t636*t643;
+ t622 = RATIONAL(-1.0,20.0)*t644;
+ t674 = t621+t622+RATIONAL(31.0,420.0);
+ t661 = x*t644;
+ t601 = RATIONAL(21.0,200.0)*x;
+ t641 = x*t643;
+ t673 = RATIONAL(1.0,40.0)*t661+t601+RATIONAL(-1.0,30.0)*t641;
+ t631 = RATIONAL(1.0,30.0);
+ t600 = RATIONAL(-21.0,200.0)*x;
+ t672 = t600+RATIONAL(-1.0,40.0)*t661+t631*t641;
+ t613 = RATIONAL(1.0,10.0)*t644;
+ t671 = RATIONAL(1.0,140.0)*t643+t613+RATIONAL(-31.0,210.0);
+ t662 = y*t643;
+ t642 = y*t644;
+ t611 = RATIONAL(-2.0,15.0)*t642;
+ t670 = RATIONAL(181.0,735.0)*y+t611+RATIONAL(1.0,98.0)*t662;
+ t587 = RATIONAL(-1.0,49.0)*t662;
+ t669 = RATIONAL(226.0,735.0)*y+t587+t611;
+ t619 = t631*t642;
+ t590 = RATIONAL(2.0,49.0)*t662;
+ t668 = RATIONAL(-289.0,2940.0)*y+t619+t590;
+ t634 = RATIONAL(-1.0,70.0);
+ t612 = t634*t643;
+ t667 = t613+t612+RATIONAL(-11.0,105.0);
+ t666 = RATIONAL(71.0,2940.0)*y+t587+t619;
+ t610 = RATIONAL(-1.0,10.0)*t644;
+ t665 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t643+t610;
+ t632 = RATIONAL(1.0,70.0);
+ t618 = t632*t643;
+ t664 = t610+t618+RATIONAL(11.0,105.0);
+ t663 = x*y;
+ t660 = x*t622+RATIONAL(-1.0,60.0)*t641;
+ t659 = x*t609+RATIONAL(1.0,60.0)*t641;
+ t574 = RATIONAL(1.0,49.0)*t662;
+ t658 = RATIONAL(166.0,735.0)*y+t611+t574;
+ t586 = RATIONAL(-2.0,49.0)*t662;
+ t657 = RATIONAL(191.0,2940.0)*y+t619+t586;
+ t623 = RATIONAL(1.0,5.0)*t642;
+ t656 = RATIONAL(-41.0,98.0)*y+t623+t586;
+ t655 = RATIONAL(-53.0,98.0)*y+t623+t574;
+ t654 = RATIONAL(37.0,300.0)*x+t660+t670;
+ t653 = RATIONAL(11.0,150.0)*x+t660+t668;
+ t652 = RATIONAL(-37.0,300.0)*x+t659+t670;
+ t651 = RATIONAL(-1.0,84.0)+t618+t622+t666;
+ t650 = x*t610+t601+RATIONAL(1.0,120.0)*t641+t669;
+ t649 = RATIONAL(1.0,84.0)+t609+t612+t666;
+ t648 = RATIONAL(-11.0,150.0)*x+t659+t668;
+ t647 = t600+RATIONAL(-1.0,120.0)*t641+x*t613+t669;
+ t592 = RATIONAL(-2.0,35.0)*t663;
+ t584 = t630*t663;
+ t583 = t636*t663;
+ t582 = RATIONAL(2.0,35.0)*t663;
+ t578 = t634*t663;
+ t576 = t632*t663;
+ coeff_dy_m2_m2 = t592+t648+t674;
+ coeff_dy_m1_m2 = t583+t651+t673;
+ coeff_dy_0_m2 = t616+RATIONAL(-17.0,420.0)+t622+t657;
+ coeff_dy_p1_m2 = t584+t651+t672;
+ coeff_dy_p2_m2 = t582+t653+t674;
+ coeff_dy_m2_m1 = t584+t650+t667;
+ coeff_dy_m1_m1 = t576+t654+t671;
+ coeff_dy_0_m1 = t618+t613+RATIONAL(-17.0,105.0)+t658;
+ coeff_dy_p1_m1 = t578+t652+t671;
+ coeff_dy_p2_m1 = t583+t647+t667;
+ coeff_dy_m2_0 = t582+t656;
+ coeff_dy_m1_0 = t584+t655;
+ coeff_dy_0_0 = RATIONAL(-57.0,98.0)*y+t590+t623;
+ coeff_dy_p1_0 = t583+t655;
+ coeff_dy_p2_0 = t592+t656;
+ coeff_dy_m2_p1 = t584+t647+t664;
+ coeff_dy_m1_p1 = t576+t652+t665;
+ coeff_dy_0_p1 = t612+t610+RATIONAL(17.0,105.0)+t658;
+ coeff_dy_p1_p1 = t578+t654+t665;
+ coeff_dy_p2_p1 = t583+t650+t664;
+ coeff_dy_m2_p2 = t592+t653+t675;
+ coeff_dy_m1_p2 = t583+t649+t672;
+ coeff_dy_0_p2 = t621+RATIONAL(17.0,420.0)+t609+t657;
+ coeff_dy_p1_p2 = t584+t649+t673;
+ coeff_dy_p2_p2 = t582+t648+t675;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c
index cc0ccef..9e72985 100644
--- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c
@@ -1,115 +1,115 @@
-fp t833,
- t808,
+fp t834,
+ t813,
+ t854,
t853,
- t807,
+ t830,
+ t835,
+ t814,
t852,
- t834,
- t818,
+ t816,
t851,
t850,
- t849,
- t829,
t815,
- t817,
- t806,
+ t808,
+ t818,
t820,
- t804,
+ t821,
+ t849,
t848,
+ t810,
t847,
- t812,
+ t804,
t846,
t819,
t845,
- t803,
t844,
t843,
- t798,
t842,
t841,
+ t801,
t840,
+ t798,
t839,
t838,
t837,
- t800,
t836,
- t835,
- t816,
- t814,
+ t817,
+ t812,
t811,
- t810,
t809,
+ t806,
t805,
+ t803,
t802,
- t801,
- t799,
- t797;
- t833 = x*x;
- t808 = RATIONAL(-1.0,49.0)*t833;
- t853 = t808+RATIONAL(71.0,2940.0);
- t807 = RATIONAL(2.0,49.0)*t833;
- t852 = t807+RATIONAL(-289.0,2940.0);
- t834 = y*y;
- t818 = RATIONAL(-2.0,5.0)*t834;
- t851 = t818+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t833;
- t850 = t818+t808+RATIONAL(226.0,735.0);
- t849 = x*y;
- t829 = RATIONAL(1.0,10.0);
- t815 = t829*y;
- t817 = RATIONAL(-1.0,5.0)*y;
- t806 = RATIONAL(-1.0,10.0)*y;
+ t800,
+ t799;
+ t834 = x*x;
+ t813 = RATIONAL(-1.0,49.0)*t834;
+ t854 = t813+RATIONAL(226.0,735.0);
+ t853 = RATIONAL(1.0,98.0)*t834+RATIONAL(181.0,735.0);
+ t830 = RATIONAL(1.0,10.0);
+ t835 = y*y;
+ t814 = t830*t835;
+ t852 = t813+t814+RATIONAL(71.0,2940.0);
+ t816 = RATIONAL(2.0,49.0)*t834;
+ t851 = t814+t816+RATIONAL(-289.0,2940.0);
+ t850 = x*y;
+ t815 = RATIONAL(-1.0,10.0)*y;
+ t808 = t830*y;
+ t818 = RATIONAL(-1.0,5.0)*y;
t820 = RATIONAL(1.0,5.0)*y;
- t804 = t829*t834;
- t848 = t815+t804;
- t847 = t806+t804;
- t812 = RATIONAL(1.0,49.0)*t833;
- t846 = t818+RATIONAL(166.0,735.0)+t812;
- t819 = RATIONAL(3.0,5.0)*t834;
- t845 = t819+RATIONAL(-53.0,98.0)+t812;
- t803 = RATIONAL(-2.0,49.0)*t833;
- t844 = t803+RATIONAL(191.0,2940.0)+t804;
- t843 = t803+t819+RATIONAL(-41.0,98.0);
+ t821 = RATIONAL(-2.0,5.0)*t835;
+ t849 = t818+t821;
+ t848 = t820+t821;
+ t810 = RATIONAL(-2.0,49.0)*t834;
+ t847 = t810+t814+RATIONAL(191.0,2940.0);
+ t804 = RATIONAL(1.0,49.0)*t834;
+ t846 = t804+t821+RATIONAL(166.0,735.0);
+ t819 = RATIONAL(3.0,5.0)*t835;
+ t845 = t810+t819+RATIONAL(-41.0,98.0);
+ t844 = t804+t819+RATIONAL(-53.0,98.0);
+ t843 = t849+t853;
+ t842 = t848+t853;
+ t841 = t849+t854;
+ t801 = x*t808;
+ t840 = t801+t851;
t798 = x*t815;
- t842 = t798+t851;
- t841 = t848+t853;
- t840 = t847+t853;
- t839 = t847+t852;
- t838 = t820+t850;
- t837 = t817+t850;
- t800 = x*t806;
- t836 = t800+t851;
- t835 = t848+t852;
- t816 = RATIONAL(2.0,35.0)*x;
- t814 = RATIONAL(1.0,70.0)*x;
- t811 = RATIONAL(-1.0,70.0)*x;
- t810 = RATIONAL(1.0,35.0)*x;
- t809 = RATIONAL(-1.0,35.0)*x;
- t805 = RATIONAL(-2.0,35.0)*x;
- t802 = x*t817;
- t801 = x*t820;
- t799 = RATIONAL(-1.0,20.0)*t849;
- t797 = RATIONAL(1.0,20.0)*t849;
- coeff_dyy_m2_m2 = t798+t805+t839;
- coeff_dyy_m1_m2 = t797+t809+t840;
- coeff_dyy_0_m2 = t806+t844;
- coeff_dyy_p1_m2 = t810+t799+t840;
- coeff_dyy_p2_m2 = t800+t816+t839;
- coeff_dyy_m2_m1 = t810+t802+t838;
- coeff_dyy_m1_m1 = t820+t814+t836;
+ t839 = t798+t851;
+ t838 = t808+t852;
+ t837 = t848+t854;
+ t836 = t815+t852;
+ t817 = RATIONAL(-2.0,35.0)*x;
+ t812 = RATIONAL(-1.0,35.0)*x;
+ t811 = RATIONAL(1.0,35.0)*x;
+ t809 = RATIONAL(2.0,35.0)*x;
+ t806 = RATIONAL(1.0,70.0)*x;
+ t805 = RATIONAL(-1.0,70.0)*x;
+ t803 = x*t818;
+ t802 = x*t820;
+ t800 = RATIONAL(1.0,20.0)*t850;
+ t799 = RATIONAL(-1.0,20.0)*t850;
+ coeff_dyy_m2_m2 = t815+t817+t840;
+ coeff_dyy_m1_m2 = t800+t812+t836;
+ coeff_dyy_0_m2 = t815+t847;
+ coeff_dyy_p1_m2 = t799+t811+t836;
+ coeff_dyy_p2_m2 = t815+t809+t839;
+ coeff_dyy_m2_m1 = t803+t811+t837;
+ coeff_dyy_m1_m1 = t806+t798+t842;
coeff_dyy_0_m1 = t820+t846;
- coeff_dyy_p1_m1 = t811+t820+t842;
- coeff_dyy_p2_m1 = t809+t801+t838;
- coeff_dyy_m2_0 = t816+t843;
- coeff_dyy_m1_0 = t810+t845;
- coeff_dyy_0_0 = t819+RATIONAL(-57.0,98.0)+t807;
- coeff_dyy_p1_0 = t809+t845;
- coeff_dyy_p2_0 = t805+t843;
- coeff_dyy_m2_p1 = t810+t801+t837;
- coeff_dyy_m1_p1 = t817+t814+t842;
- coeff_dyy_0_p1 = t817+t846;
- coeff_dyy_p1_p1 = t811+t817+t836;
- coeff_dyy_p2_p1 = t809+t802+t837;
- coeff_dyy_m2_p2 = t805+t800+t835;
- coeff_dyy_m1_p2 = t809+t799+t841;
- coeff_dyy_0_p2 = t815+t844;
- coeff_dyy_p1_p2 = t797+t810+t841;
- coeff_dyy_p2_p2 = t816+t798+t835;
+ coeff_dyy_p1_m1 = t805+t801+t842;
+ coeff_dyy_p2_m1 = t802+t812+t837;
+ coeff_dyy_m2_0 = t809+t845;
+ coeff_dyy_m1_0 = t811+t844;
+ coeff_dyy_0_0 = RATIONAL(-57.0,98.0)+t819+t816;
+ coeff_dyy_p1_0 = t812+t844;
+ coeff_dyy_p2_0 = t817+t845;
+ coeff_dyy_m2_p1 = t802+t811+t841;
+ coeff_dyy_m1_p1 = t806+t801+t843;
+ coeff_dyy_0_p1 = t818+t846;
+ coeff_dyy_p1_p1 = t805+t798+t843;
+ coeff_dyy_p2_p1 = t812+t803+t841;
+ coeff_dyy_m2_p2 = t817+t808+t839;
+ coeff_dyy_m1_p2 = t812+t799+t838;
+ coeff_dyy_0_p2 = t808+t847;
+ coeff_dyy_p1_p2 = t811+t800+t838;
+ coeff_dyy_p2_p2 = t808+t809+t840;