aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0
diff options
context:
space:
mode:
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0')
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c162
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c131
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c74
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c87
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c114
-rw-r--r--src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c71
6 files changed, 639 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c
new file mode 100644
index 0000000..cb6836d
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c
@@ -0,0 +1,162 @@
+fp t128;
+fp t116;
+fp t105;
+fp t120;
+fp t152;
+fp t127;
+fp t151;
+fp t150;
+fp t118;
+fp t149;
+fp t148;
+fp t121;
+fp t147;
+fp t146;
+fp t137;
+fp t138;
+fp t145;
+fp t139;
+fp t141;
+fp t144;
+fp t129;
+fp t130;
+fp t142;
+fp t143;
+fp t115;
+fp t104;
+fp t96;
+fp t119;
+fp t136;
+fp t110;
+fp t135;
+fp t101;
+fp t114;
+fp t134;
+fp t133;
+fp t113;
+fp t112;
+fp t111;
+fp t109;
+fp t108;
+fp t100;
+fp t98;
+fp t91;
+fp t88;
+fp t87;
+fp t86;
+fp t85;
+fp t84;
+fp t83;
+fp t82;
+ t128 = y*y;
+ t116 = RATIONAL(-3.0,40.0);
+ t105 = t116*y;
+ t120 = RATIONAL(-3.0,20.0);
+ t152 = t105+t120*t128;
+ t127 = x*x;
+ t151 = y*t127;
+ t150 = t128*x;
+ t118 = RATIONAL(-1.0,20.0);
+ t149 = t118*y;
+ t148 = RATIONAL(1.0,20.0)*y;
+ t121 = RATIONAL(3.0,20.0);
+ t147 = t121*y;
+ t146 = RATIONAL(-1.0,40.0)*y;
+ t137 = RATIONAL(1.0,8.0)*t127;
+ t138 = RATIONAL(1.0,4.0)*t127;
+ t145 = y*t137+t128*t138;
+ t139 = RATIONAL(-1.0,8.0)*t127;
+ t141 = RATIONAL(-1.0,4.0)*t127;
+ t144 = y*t139+t128*t141;
+ t129 = t127*x;
+ t130 = t128*y;
+ t142 = t129*t130;
+ t143 = RATIONAL(1.0,50.0)*t142+RATIONAL(-1.0,400.0)*t127;
+ t115 = RATIONAL(3.0,40.0);
+ t104 = t115*t128;
+ t96 = t128*t137;
+ t119 = RATIONAL(1.0,40.0);
+ t136 = t116*t128+t119*t150+t96;
+ t110 = RATIONAL(-11.0,40.0);
+ t135 = t110*t150+t104+t96;
+ t101 = t128*t139;
+ t114 = RATIONAL(11.0,40.0);
+ t134 = t114*t128+t101+x*t104;
+ t133 = RATIONAL(7.0,40.0)*t150+t101+RATIONAL(9.0,40.0)*t128;
+ t113 = RATIONAL(-1.0,300.0);
+ t112 = RATIONAL(11.0,20.0);
+ t111 = RATIONAL(-11.0,20.0);
+ t109 = RATIONAL(29.0,600.0);
+ t108 = t120*y;
+ t100 = y*t138;
+ t98 = y*t141;
+ t91 = x*t146;
+ t88 = RATIONAL(1.0,12.0)*t151;
+ t87 = RATIONAL(-1.0,12.0)*t151;
+ t86 = RATIONAL(-1.0,50.0)*t142;
+ t85 = RATIONAL(-1.0,24.0)*t151;
+ t84 = RATIONAL(1.0,24.0)*t151;
+ t83 = RATIONAL(1.0,100.0)*t142;
+ t82 = RATIONAL(-1.0,100.0)*t142;
+ coeffs_I->coeff_m1_m1 = RATIONAL(-111.0,2000.0)+RATIONAL(37.0,400.0)*t127
++t87+RATIONAL(-37.0,600.0)*t129+t149+t83+(RATIONAL(-17.0,300.0)*t127+RATIONAL(
+-2.0,125.0))*t130+(RATIONAL(481.0,6000.0)+RATIONAL(11.0,60.0)*y+RATIONAL(59.0,
+750.0)*t130)*x+t135;
+ coeffs_I->coeff_0_m1 = t88+RATIONAL(-111.0,400.0)*t127+RATIONAL(333.0,
+2000.0)+RATIONAL(-11.0,60.0)*y+RATIONAL(37.0,200.0)*t129+(RATIONAL(-481.0,
+2000.0)+t149)*x+(RATIONAL(7.0,500.0)*x+RATIONAL(-3.0,100.0)*t129+RATIONAL(13.0,
+150.0)*t127+RATIONAL(-89.0,750.0))*t130+t134;
+ coeffs_I->coeff_p1_m1 = t88+RATIONAL(-333.0,2000.0)+RATIONAL(-37.0,200.0)
+*t129+RATIONAL(111.0,400.0)*t127+t108+(RATIONAL(-7.0,60.0)*y+RATIONAL(481.0,
+2000.0))*x+(RATIONAL(3.0,100.0)*t129+t113*t127+RATIONAL(-73.0,750.0)*x+RATIONAL
+(-6.0,125.0))*t130+t133;
+ coeffs_I->coeff_p2_m1 = t82+t87+t148+RATIONAL(-37.0,400.0)*t127+RATIONAL(
+111.0,2000.0)+RATIONAL(37.0,600.0)*t129+(RATIONAL(2.0,125.0)+RATIONAL(-2.0,75.0
+)*t127)*t130+(RATIONAL(7.0,1500.0)*t130+RATIONAL(-481.0,6000.0)+RATIONAL(-1.0,
+60.0)*y)*x+t136;
+ coeffs_I->coeff_m1_0 = RATIONAL(201.0,2000.0)+RATIONAL(133.0,400.0)*t127+
+RATIONAL(-11.0,200.0)*t129+(RATIONAL(23.0,200.0)*t127+RATIONAL(1.0,150.0)*t129+
+RATIONAL(81.0,1000.0))*t130+(t114*y+RATIONAL(-957.0,2000.0)+t112*t128+RATIONAL(
+-851.0,3000.0)*t130)*x+t144+t152;
+ coeffs_I->coeff_0_0 = t86+RATIONAL(-199.0,400.0)*t127+t111*t128+t110*y+
+RATIONAL(33.0,200.0)*t129+RATIONAL(1397.0,2000.0)+(RATIONAL(-19.0,200.0)*t127+
+RATIONAL(257.0,1000.0))*t130+(RATIONAL(-129.0,2000.0)+RATIONAL(101.0,1000.0)*
+t130+t152)*x+t145;
+ coeffs_I->coeff_p1_0 = RATIONAL(-9.0,40.0)*y+RATIONAL(-33.0,200.0)*t129+
+RATIONAL(603.0,2000.0)+RATIONAL(-9.0,20.0)*t128+(RATIONAL(243.0,1000.0)+
+RATIONAL(-31.0,200.0)*t127)*t130+(RATIONAL(149.0,1000.0)*t130+RATIONAL(-7.0,
+40.0)*y+RATIONAL(-7.0,20.0)*t128+RATIONAL(1129.0,2000.0))*x+t143+t145;
+ coeffs_I->coeff_p2_0 = t115*y+RATIONAL(11.0,200.0)*t129+RATIONAL(67.0,
+400.0)*t127+RATIONAL(-43.0,2000.0)*x+t91+RATIONAL(-201.0,2000.0)+(t118*x+t121)*
+t128+(RATIONAL(-1.0,150.0)*t129+RATIONAL(101.0,3000.0)*x+RATIONAL(-81.0,1000.0)
++RATIONAL(27.0,200.0)*t127)*t130+t144;
+ coeffs_I->coeff_m1_p1 = RATIONAL(-29.0,600.0)*t129+RATIONAL(29.0,400.0)*
+t127+RATIONAL(-87.0,2000.0)+t147+t100+(t111*y+RATIONAL(377.0,6000.0))*x+(
+RATIONAL(-13.0,100.0)*t127+RATIONAL(1.0,300.0)*t129+RATIONAL(-9.0,125.0)+
+RATIONAL(203.0,750.0)*x)*t130+t135;
+ coeffs_I->coeff_0_p1 = RATIONAL(261.0,2000.0)+t98+RATIONAL(-87.0,400.0)*
+t127+t112*y+t82+RATIONAL(29.0,200.0)*t129+(RATIONAL(-71.0,250.0)+RATIONAL(7.0,
+50.0)*t127)*t130+(t147+RATIONAL(-31.0,500.0)*t130+RATIONAL(-377.0,2000.0))*x+
+t134;
+ coeffs_I->coeff_p1_p1 = RATIONAL(87.0,400.0)*t127+RATIONAL(9.0,20.0)*y+
+RATIONAL(-29.0,200.0)*t129+t98+t83+RATIONAL(-261.0,2000.0)+(RATIONAL(-27.0,
+125.0)+RATIONAL(11.0,100.0)*t127)*t130+(RATIONAL(-47.0,250.0)*t130+RATIONAL(7.0
+,20.0)*y+RATIONAL(377.0,2000.0))*x+t133;
+ coeffs_I->coeff_p2_p1 = RATIONAL(87.0,2000.0)+RATIONAL(-29.0,400.0)*t127+
+t100+t109*t129+t108+(RATIONAL(-377.0,6000.0)+t148)*x+(RATIONAL(9.0,125.0)+
+RATIONAL(-31.0,1500.0)*x+t113*t129+RATIONAL(-3.0,25.0)*t127)*t130+t136;
+ coeffs_I->coeff_m1_p2 = RATIONAL(1.0,400.0)*t127+t146+t86+RATIONAL(-1.0,
+600.0)*t129+t85+RATIONAL(-3.0,2000.0)+(RATIONAL(7.0,1000.0)+RATIONAL(43.0,600.0
+)*t127)*t130+(RATIONAL(11.0,120.0)*y+RATIONAL(-197.0,3000.0)*t130+RATIONAL(13.0
+,6000.0))*x;
+ coeffs_I->coeff_0_p2 = RATIONAL(-3.0,400.0)*t127+RATIONAL(-11.0,120.0)*y+
+t91+t84+RATIONAL(-13.0,2000.0)*x+RATIONAL(9.0,2000.0)+RATIONAL(1.0,200.0)*t129+
+(RATIONAL(-53.0,1000.0)*x+RATIONAL(3.0,50.0)*t129+RATIONAL(437.0,3000.0)+
+RATIONAL(-79.0,600.0)*t127)*t130;
+ coeffs_I->coeff_p1_p2 = RATIONAL(3.0,400.0)*t127+t105+t84+RATIONAL(-1.0,
+200.0)*t129+RATIONAL(-9.0,2000.0)+(RATIONAL(13.0,2000.0)+RATIONAL(-7.0,120.0)*y
+)*x+(RATIONAL(409.0,3000.0)*x+t109*t127+RATIONAL(21.0,1000.0)+RATIONAL(-3.0,
+50.0)*t129)*t130;
+ coeffs_I->coeff_p2_p2 = RATIONAL(3.0,2000.0)+t119*y+t85+RATIONAL(1.0,
+600.0)*t129+(RATIONAL(-7.0,1000.0)+RATIONAL(7.0,600.0)*t127)*t130+(RATIONAL(
+-53.0,3000.0)*t130+RATIONAL(-1.0,120.0)*y+RATIONAL(-13.0,6000.0))*x+t143;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c
new file mode 100644
index 0000000..a593111
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c
@@ -0,0 +1,131 @@
+fp t207;
+fp t193;
+fp t206;
+fp t205;
+fp t204;
+fp t188;
+fp t168;
+fp t203;
+fp t167;
+fp t202;
+fp t201;
+fp t191;
+fp t200;
+fp t199;
+fp t192;
+fp t194;
+fp t198;
+fp t197;
+fp t183;
+fp t196;
+fp t187;
+fp t186;
+fp t185;
+fp t184;
+fp t182;
+fp t181;
+fp t180;
+fp t179;
+fp t178;
+fp t177;
+fp t176;
+fp t175;
+fp t174;
+fp t173;
+fp t172;
+fp t171;
+fp t169;
+fp t158;
+fp t156;
+fp t155;
+fp t154;
+fp t153;
+ t207 = x*y;
+ t193 = y*y;
+ t206 = x*t193;
+ t205 = x*RATIONAL(-1.0,4.0);
+ t204 = x*RATIONAL(1.0,2.0);
+ t188 = RATIONAL(1.0,4.0);
+ t168 = t188*t206;
+ t203 = RATIONAL(-1.0,6.0)*t207+t168;
+ t167 = t193*t205;
+ t202 = RATIONAL(1.0,6.0)*t207+t167;
+ t201 = t188*t207+t193*t204;
+ t191 = RATIONAL(-1.0,2.0);
+ t200 = y*t205+t191*t206;
+ t199 = y*t204+t168;
+ t192 = x*x;
+ t194 = y*t193;
+ t198 = t192*t194;
+ t197 = t191*t207+t167;
+ t183 = RATIONAL(-1.0,200.0);
+ t196 = RATIONAL(3.0,50.0)*t198+t183*x;
+ t187 = RATIONAL(-1.0,20.0);
+ t186 = RATIONAL(3.0,200.0);
+ t185 = RATIONAL(-3.0,200.0);
+ t184 = RATIONAL(1.0,200.0);
+ t182 = RATIONAL(37.0,200.0);
+ t181 = RATIONAL(87.0,200.0);
+ t180 = RATIONAL(-37.0,200.0);
+ t179 = RATIONAL(-87.0,200.0);
+ t178 = RATIONAL(29.0,200.0);
+ t177 = RATIONAL(-29.0,200.0);
+ t176 = RATIONAL(111.0,200.0);
+ t175 = RATIONAL(-111.0,200.0);
+ t174 = RATIONAL(3.0,40.0)*t193;
+ t173 = RATIONAL(7.0,40.0)*t193;
+ t172 = RATIONAL(-1.0,40.0)*y;
+ t171 = RATIONAL(1.0,40.0)*t193;
+ t169 = RATIONAL(-11.0,40.0)*t193;
+ t158 = RATIONAL(-3.0,50.0)*t198;
+ t156 = RATIONAL(-1.0,12.0)*t207;
+ t155 = RATIONAL(1.0,12.0)*t207;
+ t154 = RATIONAL(3.0,100.0)*t198;
+ t153 = RATIONAL(-3.0,100.0)*t198;
+ coeffs_dx->coeff_m1_m1 = t180*t192+RATIONAL(481.0,6000.0)+t154+t169+
+RATIONAL(11.0,60.0)*y+RATIONAL(59.0,750.0)*t194+(RATIONAL(-17.0,150.0)*t194+
+t182)*x+t203;
+ coeffs_dx->coeff_0_m1 = t176*t192+t174+RATIONAL(-481.0,2000.0)+t175*x+
+t187*y+(RATIONAL(13.0,75.0)*x+RATIONAL(7.0,500.0)+RATIONAL(-9.0,100.0)*t192)*
+t194+t202;
+ coeffs_dx->coeff_p1_m1 = t173+t176*x+RATIONAL(-7.0,60.0)*y+t175*t192+
+RATIONAL(481.0,2000.0)+(RATIONAL(-1.0,150.0)*x+RATIONAL(-73.0,750.0)+RATIONAL(
+9.0,100.0)*t192)*t194+t202;
+ coeffs_dx->coeff_p2_m1 = RATIONAL(7.0,1500.0)*t194+t153+RATIONAL(-1.0,
+60.0)*y+t171+t182*t192+RATIONAL(-481.0,6000.0)+(RATIONAL(-4.0,75.0)*t194+t180)*
+x+t203;
+ coeffs_dx->coeff_m1_0 = RATIONAL(-33.0,200.0)*t192+RATIONAL(11.0,40.0)*y+
+RATIONAL(-957.0,2000.0)+RATIONAL(133.0,200.0)*x+RATIONAL(11.0,20.0)*t193+(
+RATIONAL(-851.0,3000.0)+RATIONAL(23.0,100.0)*x+RATIONAL(1.0,50.0)*t192)*t194+
+t200;
+ coeffs_dx->coeff_0_0 = RATIONAL(-3.0,40.0)*y+RATIONAL(-129.0,2000.0)+t158
++RATIONAL(101.0,1000.0)*t194+RATIONAL(-3.0,20.0)*t193+RATIONAL(99.0,200.0)*t192
++(RATIONAL(-19.0,100.0)*t194+RATIONAL(-199.0,200.0))*x+t201;
+ coeffs_dx->coeff_p1_0 = RATIONAL(-7.0,40.0)*y+RATIONAL(1129.0,2000.0)+
+RATIONAL(-7.0,20.0)*t193+RATIONAL(-99.0,200.0)*t192+(RATIONAL(149.0,1000.0)+
+RATIONAL(-31.0,100.0)*x)*t194+t196+t201;
+ coeffs_dx->coeff_p2_0 = t172+t187*t193+RATIONAL(33.0,200.0)*t192+RATIONAL
+(-43.0,2000.0)+RATIONAL(67.0,200.0)*x+(RATIONAL(-1.0,50.0)*t192+RATIONAL(27.0,
+100.0)*x+RATIONAL(101.0,3000.0))*t194+t200;
+ coeffs_dx->coeff_m1_p1 = RATIONAL(-11.0,20.0)*y+t169+RATIONAL(377.0,
+6000.0)+t178*x+t177*t192+(RATIONAL(-13.0,50.0)*x+RATIONAL(1.0,100.0)*t192+
+RATIONAL(203.0,750.0))*t194+t199;
+ coeffs_dx->coeff_0_p1 = t153+t181*t192+t174+RATIONAL(-31.0,500.0)*t194+
+RATIONAL(3.0,20.0)*y+RATIONAL(-377.0,2000.0)+(t179+RATIONAL(7.0,25.0)*t194)*x+
+t197;
+ coeffs_dx->coeff_p1_p1 = RATIONAL(7.0,20.0)*y+RATIONAL(-47.0,250.0)*t194+
+t179*t192+RATIONAL(377.0,2000.0)+t173+t154+(RATIONAL(11.0,50.0)*t194+t181)*x+
+t197;
+ coeffs_dx->coeff_p2_p1 = RATIONAL(1.0,20.0)*y+RATIONAL(-377.0,6000.0)+
+t178*t192+t171+t177*x+(RATIONAL(-31.0,1500.0)+RATIONAL(-6.0,25.0)*x+RATIONAL(
+-1.0,100.0)*t192)*t194+t199;
+ coeffs_dx->coeff_m1_p2 = RATIONAL(11.0,120.0)*y+t183*t192+RATIONAL(13.0,
+6000.0)+RATIONAL(-197.0,3000.0)*t194+t156+t158+(RATIONAL(43.0,300.0)*t194+t184)
+*x;
+ coeffs_dx->coeff_0_p2 = t186*t192+t172+t155+RATIONAL(-13.0,2000.0)+t185*x
++(RATIONAL(-53.0,1000.0)+RATIONAL(-79.0,300.0)*x+RATIONAL(9.0,50.0)*t192)*t194;
+ coeffs_dx->coeff_p1_p2 = t155+t185*t192+RATIONAL(-7.0,120.0)*y+t186*x+
+RATIONAL(13.0,2000.0)+(RATIONAL(409.0,3000.0)+RATIONAL(-9.0,50.0)*t192+RATIONAL
+(29.0,300.0)*x)*t194;
+ coeffs_dx->coeff_p2_p2 = RATIONAL(-1.0,120.0)*y+t184*t192+RATIONAL(-13.0,
+6000.0)+t156+(RATIONAL(-53.0,3000.0)+RATIONAL(7.0,300.0)*x)*t194+t196;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c
new file mode 100644
index 0000000..8c53891
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c
@@ -0,0 +1,74 @@
+fp t286;
+fp t287;
+fp t296;
+fp t282;
+fp t276;
+fp t284;
+fp t295;
+fp t285;
+fp t283;
+fp t294;
+fp t274;
+fp t293;
+fp t292;
+fp t291;
+fp t290;
+fp t289;
+fp t270;
+fp t269;
+fp t268;
+fp t266;
+fp t265;
+ t286 = y*y;
+ t287 = t286*y;
+ t296 = x*t287;
+ t282 = RATIONAL(1.0,4.0);
+ t276 = t282*t286;
+ t284 = RATIONAL(1.0,2.0);
+ t295 = t284*y+t276;
+ t285 = RATIONAL(-1.0,2.0);
+ t283 = RATIONAL(-1.0,4.0);
+ t294 = t283*y+t285*t286;
+ t274 = t283*t286;
+ t293 = t285*y+t274;
+ t292 = RATIONAL(3.0,25.0)*t296+RATIONAL(-1.0,200.0);
+ t291 = t284*t286+t282*y;
+ t290 = RATIONAL(-1.0,6.0)*y+t276;
+ t289 = RATIONAL(1.0,6.0)*y+t274;
+ t270 = RATIONAL(1.0,12.0)*y;
+ t269 = RATIONAL(-1.0,12.0)*y;
+ t268 = RATIONAL(-3.0,25.0)*t296;
+ t266 = RATIONAL(3.0,50.0)*t296;
+ t265 = RATIONAL(-3.0,50.0)*t296;
+ coeffs_dxx->coeff_m1_m1 = RATIONAL(-37.0,100.0)*x+RATIONAL(-17.0,150.0)*
+t287+t266+RATIONAL(37.0,200.0)+t290;
+ coeffs_dxx->coeff_0_m1 = RATIONAL(-111.0,200.0)+RATIONAL(13.0,75.0)*t287+
+(RATIONAL(-9.0,50.0)*t287+RATIONAL(111.0,100.0))*x+t289;
+ coeffs_dxx->coeff_p1_m1 = RATIONAL(-1.0,150.0)*t287+RATIONAL(111.0,200.0)
++(RATIONAL(-111.0,100.0)+RATIONAL(9.0,50.0)*t287)*x+t289;
+ coeffs_dxx->coeff_p2_m1 = RATIONAL(-37.0,200.0)+RATIONAL(-4.0,75.0)*t287+
+t265+RATIONAL(37.0,100.0)*x+t290;
+ coeffs_dxx->coeff_m1_0 = RATIONAL(23.0,100.0)*t287+RATIONAL(133.0,200.0)+
+(RATIONAL(-33.0,100.0)+RATIONAL(1.0,25.0)*t287)*x+t294;
+ coeffs_dxx->coeff_0_0 = RATIONAL(-199.0,200.0)+RATIONAL(-19.0,100.0)*t287
++RATIONAL(99.0,100.0)*x+t268+t291;
+ coeffs_dxx->coeff_p1_0 = RATIONAL(-99.0,100.0)*x+RATIONAL(-31.0,100.0)*
+t287+t291+t292;
+ coeffs_dxx->coeff_p2_0 = RATIONAL(67.0,200.0)+RATIONAL(27.0,100.0)*t287+(
+RATIONAL(-1.0,25.0)*t287+RATIONAL(33.0,100.0))*x+t294;
+ coeffs_dxx->coeff_m1_p1 = RATIONAL(-13.0,50.0)*t287+RATIONAL(29.0,200.0)+
+(RATIONAL(1.0,50.0)*t287+RATIONAL(-29.0,100.0))*x+t295;
+ coeffs_dxx->coeff_0_p1 = RATIONAL(-87.0,200.0)+t265+RATIONAL(87.0,100.0)*
+x+RATIONAL(7.0,25.0)*t287+t293;
+ coeffs_dxx->coeff_p1_p1 = RATIONAL(-87.0,100.0)*x+t266+RATIONAL(87.0,
+200.0)+RATIONAL(11.0,50.0)*t287+t293;
+ coeffs_dxx->coeff_p2_p1 = RATIONAL(-6.0,25.0)*t287+RATIONAL(-29.0,200.0)+
+(RATIONAL(-1.0,50.0)*t287+RATIONAL(29.0,100.0))*x+t295;
+ coeffs_dxx->coeff_m1_p2 = RATIONAL(43.0,300.0)*t287+RATIONAL(-1.0,100.0)*
+x+t268+t269+RATIONAL(1.0,200.0);
+ coeffs_dxx->coeff_0_p2 = RATIONAL(-79.0,300.0)*t287+RATIONAL(-3.0,200.0)+
+t270+(RATIONAL(3.0,100.0)+RATIONAL(9.0,25.0)*t287)*x;
+ coeffs_dxx->coeff_p1_p2 = t270+RATIONAL(29.0,300.0)*t287+RATIONAL(3.0,
+200.0)+(RATIONAL(-9.0,25.0)*t287+RATIONAL(-3.0,100.0))*x;
+ coeffs_dxx->coeff_p2_p2 = RATIONAL(1.0,100.0)*x+RATIONAL(7.0,300.0)*t287+
+t269+t292;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c
new file mode 100644
index 0000000..ed83ce7
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c
@@ -0,0 +1,87 @@
+fp t311;
+fp t314;
+fp t318;
+fp t301;
+fp t331;
+fp t322;
+fp t330;
+fp t315;
+fp t302;
+fp t329;
+fp t328;
+fp t316;
+fp t327;
+fp t319;
+fp t326;
+fp t323;
+fp t324;
+fp t325;
+fp t317;
+fp t310;
+fp t309;
+fp t307;
+fp t306;
+fp t300;
+fp t299;
+fp t298;
+fp t297;
+ t311 = RATIONAL(-1.0,2.0)*x;
+ t314 = RATIONAL(1.0,2.0)*x;
+ t318 = RATIONAL(7.0,20.0);
+ t301 = y*t311;
+ t331 = t318*y+t301;
+ t322 = x*y;
+ t330 = t322+RATIONAL(1.0,4.0)*x;
+ t315 = RATIONAL(-11.0,20.0);
+ t302 = y*t314;
+ t329 = t302+t315*y;
+ t328 = -t322+RATIONAL(-1.0,4.0)*x;
+ t316 = RATIONAL(1.0,20.0);
+ t327 = t302+t316*y;
+ t319 = RATIONAL(3.0,20.0);
+ t326 = t319*y+t301;
+ t323 = x*x;
+ t324 = y*y;
+ t325 = t323*t324;
+ t317 = RATIONAL(-1.0,40.0);
+ t310 = RATIONAL(-1.0,6.0)*x;
+ t309 = RATIONAL(1.0,6.0)*x;
+ t307 = RATIONAL(-1.0,12.0)*x;
+ t306 = RATIONAL(1.0,12.0)*x;
+ t300 = RATIONAL(-9.0,50.0)*t325;
+ t299 = RATIONAL(9.0,50.0)*t325;
+ t298 = RATIONAL(-9.0,100.0)*t325;
+ t297 = RATIONAL(9.0,100.0)*t325;
+ coeffs_dxy->coeff_m1_m1 = t310+RATIONAL(11.0,60.0)+t297+(RATIONAL(59.0,
+250.0)+RATIONAL(-17.0,50.0)*x)*t324+t329;
+ coeffs_dxy->coeff_0_m1 = RATIONAL(-1.0,20.0)+t309+(RATIONAL(21.0,500.0)+
+RATIONAL(13.0,25.0)*x+RATIONAL(-27.0,100.0)*t323)*t324+t326;
+ coeffs_dxy->coeff_p1_m1 = t309+RATIONAL(-7.0,60.0)+(RATIONAL(-1.0,50.0)*x
++RATIONAL(-73.0,250.0)+RATIONAL(27.0,100.0)*t323)*t324+t331;
+ coeffs_dxy->coeff_p2_m1 = t298+t310+RATIONAL(-1.0,60.0)+(RATIONAL(-4.0,
+25.0)*x+RATIONAL(7.0,500.0))*t324+t327;
+ coeffs_dxy->coeff_m1_0 = RATIONAL(11.0,10.0)*y+RATIONAL(11.0,40.0)+(
+RATIONAL(3.0,50.0)*t323+RATIONAL(-851.0,1000.0)+RATIONAL(69.0,100.0)*x)*t324+
+t328;
+ coeffs_dxy->coeff_0_0 = t300+RATIONAL(-3.0,10.0)*y+RATIONAL(-3.0,40.0)+(
+RATIONAL(303.0,1000.0)+RATIONAL(-57.0,100.0)*x)*t324+t330;
+ coeffs_dxy->coeff_p1_0 = RATIONAL(-7.0,10.0)*y+RATIONAL(-7.0,40.0)+t299+(
+RATIONAL(447.0,1000.0)+RATIONAL(-93.0,100.0)*x)*t324+t330;
+ coeffs_dxy->coeff_p2_0 = t317+RATIONAL(-1.0,10.0)*y+(RATIONAL(101.0,
+1000.0)+RATIONAL(-3.0,50.0)*t323+RATIONAL(81.0,100.0)*x)*t324+t328;
+ coeffs_dxy->coeff_m1_p1 = t315+t314+(RATIONAL(203.0,250.0)+RATIONAL(-39.0
+,50.0)*x+RATIONAL(3.0,100.0)*t323)*t324+t329;
+ coeffs_dxy->coeff_0_p1 = t319+t311+t298+(RATIONAL(-93.0,500.0)+RATIONAL(
+21.0,25.0)*x)*t324+t326;
+ coeffs_dxy->coeff_p1_p1 = t311+t297+t318+(RATIONAL(33.0,50.0)*x+RATIONAL(
+-141.0,250.0))*t324+t331;
+ coeffs_dxy->coeff_p2_p1 = t316+t314+(RATIONAL(-3.0,100.0)*t323+RATIONAL(
+-31.0,500.0)+RATIONAL(-18.0,25.0)*x)*t324+t327;
+ coeffs_dxy->coeff_m1_p2 = t307+RATIONAL(11.0,120.0)+t300+(RATIONAL(-197.0
+,1000.0)+RATIONAL(43.0,100.0)*x)*t324;
+ coeffs_dxy->coeff_0_p2 = t317+t306+(RATIONAL(-159.0,1000.0)+RATIONAL(27.0
+,50.0)*t323+RATIONAL(-79.0,100.0)*x)*t324;
+ coeffs_dxy->coeff_p1_p2 = t306+RATIONAL(-7.0,120.0)+(RATIONAL(-27.0,50.0)
+*t323+RATIONAL(29.0,100.0)*x+RATIONAL(409.0,1000.0))*t324;
+ coeffs_dxy->coeff_p2_p2 = t307+RATIONAL(-1.0,120.0)+t299+(RATIONAL(-53.0,
+1000.0)+RATIONAL(7.0,100.0)*x)*t324;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c
new file mode 100644
index 0000000..6c1ce7f
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c
@@ -0,0 +1,114 @@
+fp t264;
+fp t240;
+fp t242;
+fp t263;
+fp t247;
+fp t262;
+fp t261;
+fp t243;
+fp t226;
+fp t249;
+fp t248;
+fp t258;
+fp t232;
+fp t229;
+fp t257;
+fp t256;
+fp t234;
+fp t218;
+fp t255;
+fp t219;
+fp t254;
+fp t253;
+fp t252;
+fp t239;
+fp t251;
+fp t244;
+fp t241;
+fp t235;
+fp t228;
+fp t227;
+fp t225;
+fp t224;
+fp t223;
+fp t215;
+fp t211;
+fp t210;
+fp t208;
+ t264 = RATIONAL(7.0,20.0)*x+RATIONAL(9.0,20.0);
+ t240 = RATIONAL(1.0,20.0);
+ t242 = RATIONAL(-3.0,20.0);
+ t263 = t240*x+t242;
+ t247 = x*x;
+ t262 = y*t247;
+ t261 = RATIONAL(-11.0,20.0)*x;
+ t243 = RATIONAL(3.0,20.0);
+ t226 = t243*y;
+ t249 = t247*x;
+ t248 = y*y;
+ t258 = t249*t248;
+ t232 = RATIONAL(1.0,4.0)*t247;
+ t229 = RATIONAL(-1.0,4.0)*t247;
+ t257 = RATIONAL(-1.0,8.0)*t247+RATIONAL(-1.0,2.0)*t262;
+ t256 = RATIONAL(1.0,8.0)*t247+RATIONAL(1.0,2.0)*t262;
+ t234 = RATIONAL(11.0,20.0);
+ t218 = y*t229;
+ t255 = x*t226+t234*y+t218;
+ t219 = y*t232;
+ t254 = t226+t219+y*t261;
+ t253 = t219+t263*y;
+ t252 = t218+t264*y;
+ t239 = RATIONAL(-3.0,40.0);
+ t251 = RATIONAL(-3.0,10.0)*y+t239;
+ t244 = RATIONAL(-1.0,40.0);
+ t241 = RATIONAL(-1.0,20.0);
+ t235 = RATIONAL(-1.0,100.0);
+ t228 = RATIONAL(-1.0,24.0)*t247;
+ t227 = t244*x;
+ t225 = RATIONAL(1.0,24.0)*t247;
+ t224 = RATIONAL(-1.0,12.0)*t247;
+ t223 = RATIONAL(1.0,12.0)*t247;
+ t215 = RATIONAL(3.0,50.0)*t258;
+ t211 = RATIONAL(-3.0,50.0)*t258;
+ t210 = RATIONAL(3.0,100.0)*t258;
+ t208 = RATIONAL(-3.0,100.0)*t258;
+ coeffs_dy->coeff_m1_m1 = t241+t210+RATIONAL(11.0,60.0)*x+t224+(RATIONAL(
+-17.0,100.0)*t247+RATIONAL(-6.0,125.0)+RATIONAL(59.0,250.0)*x)*t248+t254;
+ coeffs_dy->coeff_0_m1 = t223+t241*x+RATIONAL(-11.0,60.0)+(RATIONAL(21.0,
+500.0)*x+RATIONAL(13.0,50.0)*t247+RATIONAL(-9.0,100.0)*t249+RATIONAL(-89.0,
+250.0))*t248+t255;
+ coeffs_dy->coeff_p1_m1 = t242+t223+RATIONAL(-7.0,60.0)*x+(RATIONAL(-18.0,
+125.0)+t235*t247+RATIONAL(9.0,100.0)*t249+RATIONAL(-73.0,250.0)*x)*t248+t252;
+ coeffs_dy->coeff_p2_m1 = RATIONAL(-1.0,60.0)*x+t240+t208+t224+(RATIONAL(
+-2.0,25.0)*t247+RATIONAL(6.0,125.0)+RATIONAL(7.0,500.0)*x)*t248+t253;
+ coeffs_dy->coeff_m1_0 = (RATIONAL(11.0,40.0)+RATIONAL(11.0,10.0)*y)*x+(
+RATIONAL(-851.0,1000.0)*x+RATIONAL(69.0,200.0)*t247+RATIONAL(1.0,50.0)*t249+
+RATIONAL(243.0,1000.0))*t248+t251+t257;
+ coeffs_dy->coeff_0_0 = RATIONAL(-11.0,10.0)*y+t211+RATIONAL(-11.0,40.0)+(
+RATIONAL(771.0,1000.0)+RATIONAL(-57.0,200.0)*t247)*t248+(RATIONAL(303.0,1000.0)
+*t248+t251)*x+t256;
+ coeffs_dy->coeff_p1_0 = t215+RATIONAL(-9.0,10.0)*y+RATIONAL(-9.0,40.0)+(
+RATIONAL(-93.0,200.0)*t247+RATIONAL(729.0,1000.0))*t248+(RATIONAL(-7.0,10.0)*y+
+RATIONAL(-7.0,40.0)+RATIONAL(447.0,1000.0)*t248)*x+t256;
+ coeffs_dy->coeff_p2_0 = RATIONAL(3.0,40.0)+t227+(RATIONAL(-1.0,10.0)*x+
+RATIONAL(3.0,10.0))*y+(RATIONAL(81.0,200.0)*t247+RATIONAL(101.0,1000.0)*x+
+RATIONAL(-243.0,1000.0)+RATIONAL(-1.0,50.0)*t249)*t248+t257;
+ coeffs_dy->coeff_m1_p1 = t243+t232+t261+(RATIONAL(-27.0,125.0)+RATIONAL(
+-39.0,100.0)*t247+RATIONAL(203.0,250.0)*x+RATIONAL(1.0,100.0)*t249)*t248+t254;
+ coeffs_dy->coeff_0_p1 = t229+t234+t208+t243*x+(RATIONAL(-213.0,250.0)+
+RATIONAL(21.0,50.0)*t247+RATIONAL(-93.0,500.0)*x)*t248+t255;
+ coeffs_dy->coeff_p1_p1 = t229+t210+(RATIONAL(33.0,100.0)*t247+RATIONAL(
+-141.0,250.0)*x+RATIONAL(-81.0,125.0))*t248+t252+t264;
+ coeffs_dy->coeff_p2_p1 = t232+(RATIONAL(-9.0,25.0)*t247+RATIONAL(-31.0,
+500.0)*x+t235*t249+RATIONAL(27.0,125.0))*t248+t253+t263;
+ coeffs_dy->coeff_m1_p2 = t244+t228+RATIONAL(11.0,120.0)*x+t211+(RATIONAL(
+-197.0,1000.0)*x+RATIONAL(43.0,200.0)*t247+RATIONAL(21.0,1000.0))*t248;
+ coeffs_dy->coeff_0_p2 = t225+t227+RATIONAL(-11.0,120.0)+(RATIONAL(-159.0,
+1000.0)*x+RATIONAL(9.0,50.0)*t249+RATIONAL(-79.0,200.0)*t247+RATIONAL(437.0,
+1000.0))*t248;
+ coeffs_dy->coeff_p1_p2 = t239+RATIONAL(-7.0,120.0)*x+t225+(RATIONAL(29.0,
+200.0)*t247+RATIONAL(-9.0,50.0)*t249+RATIONAL(63.0,1000.0)+RATIONAL(409.0,
+1000.0)*x)*t248;
+ coeffs_dy->coeff_p2_p2 = RATIONAL(1.0,40.0)+t215+t228+RATIONAL(-1.0,120.0
+)*x+(RATIONAL(7.0,200.0)*t247+RATIONAL(-21.0,1000.0)+RATIONAL(-53.0,1000.0)*x)*
+t248;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c
new file mode 100644
index 0000000..1b690c6
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c
@@ -0,0 +1,71 @@
+fp t350;
+fp t351;
+fp t357;
+fp t342;
+fp t356;
+fp t343;
+fp t355;
+fp t348;
+fp t354;
+fp t353;
+fp t349;
+fp t345;
+fp t341;
+fp t340;
+fp t335;
+fp t334;
+fp t333;
+fp t332;
+ t350 = x*x;
+ t351 = t350*x;
+ t357 = y*t351;
+ t342 = RATIONAL(-1.0,4.0)*t350;
+ t356 = RATIONAL(7.0,20.0)*x+t342+RATIONAL(9.0,20.0);
+ t343 = RATIONAL(1.0,4.0)*t350;
+ t355 = RATIONAL(1.0,20.0)*x+t343+RATIONAL(-3.0,20.0);
+ t348 = RATIONAL(3.0,20.0);
+ t354 = t348*x+t342+RATIONAL(11.0,20.0);
+ t353 = t343+t348+RATIONAL(-11.0,20.0)*x;
+ t349 = RATIONAL(-3.0,10.0);
+ t345 = RATIONAL(-1.0,50.0);
+ t341 = RATIONAL(-1.0,2.0)*t350;
+ t340 = RATIONAL(1.0,2.0)*t350;
+ t335 = RATIONAL(3.0,50.0)*t357;
+ t334 = RATIONAL(-3.0,50.0)*t357;
+ t333 = RATIONAL(-3.0,25.0)*t357;
+ t332 = RATIONAL(3.0,25.0)*t357;
+ coeffs_dyy->coeff_m1_m1 = t335+(RATIONAL(59.0,125.0)*x+RATIONAL(-12.0,
+125.0)+RATIONAL(-17.0,50.0)*t350)*y+t353;
+ coeffs_dyy->coeff_0_m1 = (RATIONAL(-89.0,125.0)+RATIONAL(21.0,250.0)*x+
+RATIONAL(13.0,25.0)*t350+RATIONAL(-9.0,50.0)*t351)*y+t354;
+ coeffs_dyy->coeff_p1_m1 = (RATIONAL(-73.0,125.0)*x+RATIONAL(-36.0,125.0)+
+RATIONAL(9.0,50.0)*t351+t345*t350)*y+t356;
+ coeffs_dyy->coeff_p2_m1 = t334+(RATIONAL(12.0,125.0)+RATIONAL(-4.0,25.0)*
+t350+RATIONAL(7.0,250.0)*x)*y+t355;
+ coeffs_dyy->coeff_m1_0 = RATIONAL(11.0,10.0)*x+t349+t341+(RATIONAL(243.0,
+500.0)+RATIONAL(1.0,25.0)*t351+RATIONAL(69.0,100.0)*t350+RATIONAL(-851.0,500.0)
+*x)*y;
+ coeffs_dyy->coeff_0_0 = t340+t349*x+t333+RATIONAL(-11.0,10.0)+(RATIONAL(
+-57.0,100.0)*t350+RATIONAL(771.0,500.0)+RATIONAL(303.0,500.0)*x)*y;
+ coeffs_dyy->coeff_p1_0 = t332+t340+RATIONAL(-9.0,10.0)+RATIONAL(-7.0,10.0
+)*x+(RATIONAL(729.0,500.0)+RATIONAL(-93.0,100.0)*t350+RATIONAL(447.0,500.0)*x)*
+y;
+ coeffs_dyy->coeff_p2_0 = t341+RATIONAL(3.0,10.0)+RATIONAL(-1.0,10.0)*x+(
+RATIONAL(-243.0,500.0)+RATIONAL(101.0,500.0)*x+RATIONAL(81.0,100.0)*t350+
+RATIONAL(-1.0,25.0)*t351)*y;
+ coeffs_dyy->coeff_m1_p1 = (RATIONAL(1.0,50.0)*t351+RATIONAL(-54.0,125.0)+
+RATIONAL(203.0,125.0)*x+RATIONAL(-39.0,50.0)*t350)*y+t353;
+ coeffs_dyy->coeff_0_p1 = t334+(RATIONAL(-213.0,125.0)+RATIONAL(21.0,25.0)
+*t350+RATIONAL(-93.0,250.0)*x)*y+t354;
+ coeffs_dyy->coeff_p1_p1 = t335+(RATIONAL(-141.0,125.0)*x+RATIONAL(33.0,
+50.0)*t350+RATIONAL(-162.0,125.0))*y+t356;
+ coeffs_dyy->coeff_p2_p1 = (RATIONAL(54.0,125.0)+RATIONAL(-18.0,25.0)*t350
++RATIONAL(-31.0,250.0)*x+t345*t351)*y+t355;
+ coeffs_dyy->coeff_m1_p2 = t333+(RATIONAL(21.0,500.0)+RATIONAL(-197.0,
+500.0)*x+RATIONAL(43.0,100.0)*t350)*y;
+ coeffs_dyy->coeff_0_p2 = (RATIONAL(9.0,25.0)*t351+RATIONAL(-159.0,500.0)*
+x+RATIONAL(-79.0,100.0)*t350+RATIONAL(437.0,500.0))*y;
+ coeffs_dyy->coeff_p1_p2 = (RATIONAL(29.0,100.0)*t350+RATIONAL(-9.0,25.0)*
+t351+RATIONAL(63.0,500.0)+RATIONAL(409.0,500.0)*x)*y;
+ coeffs_dyy->coeff_p2_p2 = t332+(RATIONAL(7.0,100.0)*t350+RATIONAL(-53.0,
+500.0)*x+RATIONAL(-21.0,500.0))*y;