diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0')
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; |