diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs')
12 files changed, 2337 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-I.compute.c new file mode 100644 index 0000000..18caffa --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-I.compute.c @@ -0,0 +1,130 @@ +fp t38; +fp t41; +fp t65; +fp t40; +fp t39; +fp t63; +fp t48; +fp t49; +fp t36; +fp t68; +fp t29; +fp t67; +fp t66; +fp t32; +fp t27; +fp t53; +fp t7; +fp t51; +fp t64; +fp t58; +fp t62; +fp t35; +fp t33; +fp t55; +fp t8; +fp t61; +fp t37; +fp t57; +fp t60; +fp t59; +fp t56; +fp t54; +fp t52; +fp t31; +fp t50; +fp t26; +fp t28; +fp t47; +fp t45; +fp t34; +fp t30; +fp t25; +fp t24; +fp t23; +fp t22; +fp t21; +fp t20; +fp t19; +fp t17; +fp t15; +fp t14; +fp t13; +fp t12; +fp t11; +fp t10; +fp t4; + t38 = y*y; + t41 = y*t38; + t65 = -t41-y; + t40 = x*x; + t39 = x*t40; + t63 = t39+x; + t48 = t63*t38; + t49 = t65*t40; + t36 = RATIONAL(1.0,4.0); + t68 = y*t36; + t29 = RATIONAL(-1.0,4.0); + t67 = t29*x; + t66 = t29*y; + t32 = RATIONAL(5.0,4.0); + t27 = RATIONAL(3.0,4.0); + t53 = t27*t41; + t7 = t39*t53; + t51 = t40*t38; + t64 = t32*t51+t7; + t58 = t39*t41; + t62 = RATIONAL(-9.0,4.0)*t58+RATIONAL(-5.0,1.0)*t51; + t35 = RATIONAL(-5.0,2.0); + t33 = RATIONAL(-3.0,4.0); + t55 = t33*t41; + t8 = t39*t55; + t61 = t8+t35*t51; + t37 = RATIONAL(-1.0,2.0); + t57 = t29*t39; + t60 = t41*t57+t37*t51; + t59 = t40*t41; + t56 = t33*t39; + t54 = t27*t39; + t52 = t39*t38; + t31 = RATIONAL(2.0,1.0); + t50 = t31*t51+x*t66+t7; + t26 = RATIONAL(-5.0,4.0); + t28 = RATIONAL(1.0,2.0); + t47 = t26*t38+t28; + t45 = t28+t26*t40; + t34 = RATIONAL(3.0,2.0); + t30 = RATIONAL(-3.0,2.0); + t25 = RATIONAL(3.0,1.0); + t24 = RATIONAL(-3.0,1.0); + t23 = RATIONAL(15.0,4.0); + t22 = RATIONAL(-15.0,4.0); + t21 = x*t53; + t20 = x*t68; + t19 = y*t54; + t17 = x*t55; + t15 = RATIONAL(9.0,4.0)*t58; + t14 = y*t56; + t13 = t41*t67; + t12 = t36*x*t41; + t11 = t39*t68; + t10 = t36*t58; + t4 = y*t57; + coeffs_I->coeff_m1_m1 = t20+t12+t11+t10+t51+(-t49+t48)*t37; + coeffs_I->coeff_0_m1 = t38+t14+t34*t52+t61-t65*(t32*t40+t37); + coeffs_I->coeff_p1_m1 = t13+t19+t49+(t30*t39+t28*x)*t38+t50; + coeffs_I->coeff_p2_m1 = t28*t52+t4-t36*t49+t60; + coeffs_I->coeff_m1_0 = t17+t40+t34*t59+t61+t63*(t32*t38+t37); + coeffs_I->coeff_0_0 = t15+RATIONAL(1.0,1.0)+(t35+t22*t41)*t40+(t41+t39)* +t34+(t35+t22*t39+RATIONAL(25.0,4.0)*t40)*t38; + coeffs_I->coeff_p1_0 = t21+(t31+t25*t41)*t40+(t30+t23*t38)*t39+t47*x+t62; + coeffs_I->coeff_p2_0 = (t37+t55)*t40+t47*t39+t64; + coeffs_I->coeff_m1_p1 = t4+t21+(t30*t41+t28*y)*t40-t48+t50; + coeffs_I->coeff_0_p1 = t19+(t23*t40+t30)*t41+(t25*t39+t31)*t38+t45*y+t62; + coeffs_I->coeff_p1_p1 = t20+t17+t14+t15+(t24*t41+y)*t40+(RATIONAL(4.0,1.0 +)*t40+x+t24*t39)*t38; + coeffs_I->coeff_p2_p1 = t8+t11+t52+(t53+t66-t38)*t40; + coeffs_I->coeff_m1_p2 = t28*t59+t13+t36*t48+t60; + coeffs_I->coeff_0_p2 = t45*t41+(t56+t37)*t38+t64; + coeffs_I->coeff_p1_p2 = t8+t59+t12+(t54+t67-t40)*t38; + coeffs_I->coeff_p2_p2 = t10+t36*t51+(t59+t52)*t29; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dx.compute.c new file mode 100644 index 0000000..b7bcbbe --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dx.compute.c @@ -0,0 +1,125 @@ +fp t110; +fp t127; +fp t112; +fp t126; +fp t125; +fp t90; +fp t96; +fp t111; +fp t116; +fp t73; +fp t124; +fp t123; +fp t98; +fp t122; +fp t92; +fp t99; +fp t76; +fp t121; +fp t120; +fp t119; +fp t105; +fp t118; +fp t100; +fp t117; +fp t95; +fp t102; +fp t115; +fp t114; +fp t109; +fp t108; +fp t107; +fp t106; +fp t104; +fp t103; +fp t101; +fp t97; +fp t94; +fp t93; +fp t91; +fp t89; +fp t88; +fp t86; +fp t85; +fp t84; +fp t81; +fp t80; +fp t79; +fp t77; +fp t75; +fp t74; +fp t71; + t110 = y*y; + t127 = x*t110; + t112 = t110*y; + t126 = x*t112; + t125 = t112+y; + t90 = RATIONAL(-2.0,1.0); + t96 = RATIONAL(-9.0,4.0); + t111 = x*x; + t116 = t111*t112; + t73 = t96*t116; + t124 = t90*t127+t73; + t123 = RATIONAL(-10.0,1.0)*t127+RATIONAL(-27.0,4.0)*t116; + t98 = RATIONAL(-5.0,1.0); + t122 = t73+t98*t127; + t92 = RATIONAL(5.0,2.0); + t99 = RATIONAL(9.0,4.0); + t76 = t99*t116; + t121 = t92*t127+t76; + t120 = t99*t111; + t119 = t111*t110; + t105 = RATIONAL(-3.0,4.0); + t118 = t105*t111; + t100 = RATIONAL(3.0,4.0); + t117 = t100*t111; + t95 = RATIONAL(4.0,1.0); + t102 = RATIONAL(-1.0,4.0); + t115 = t102*y+t95*t127+t76; + t114 = x*t125; + t109 = RATIONAL(-1.0,2.0); + t108 = RATIONAL(1.0,4.0); + t107 = RATIONAL(-5.0,2.0); + t106 = RATIONAL(3.0,2.0); + t104 = RATIONAL(2.0,1.0); + t103 = RATIONAL(-3.0,2.0); + t101 = RATIONAL(1.0,2.0); + t97 = RATIONAL(3.0,1.0); + t94 = RATIONAL(-3.0,1.0); + t93 = RATIONAL(9.0,2.0); + t91 = RATIONAL(-9.0,2.0); + t89 = t108*y; + t88 = t108*t112; + t86 = t102*t112; + t85 = t105*t112; + t84 = t100*t112; + t81 = y*t117; + t80 = t100*t116; + t79 = t96*t111*y; + t77 = y*t120; + t75 = y*t118; + t74 = t105*t116; + t71 = RATIONAL(27.0,4.0)*t116; + coeffs_dx->coeff_m1_m1 = t89+t81+t88+t80+(t103*t111+t109)*t110+(t104*t110 +-t125)*x; + coeffs_dx->coeff_0_m1 = t79+t93*t119+t92*t114+t122; + coeffs_dx->coeff_p1_m1 = t86+t77+(t101+t91*t111)*t110+t90*t114+t115; + coeffs_dx->coeff_p2_m1 = t106*t119+t75+t74+(-t110+t125*t101)*x; + coeffs_dx->coeff_m1_0 = t109+RATIONAL(5.0,4.0)*t110+t85+(t103+RATIONAL( +15.0,4.0)*t110)*t111+(t97*t112+t104)*x+t122; + coeffs_dx->coeff_0_0 = t71+(t93+RATIONAL(-45.0,4.0)*t110)*t111+(RATIONAL( +25.0,2.0)*t110+RATIONAL(-15.0,2.0)*t112+t98)*x; + coeffs_dx->coeff_p1_0 = RATIONAL(-5.0,4.0)*t110+t101+t84+(RATIONAL(45.0, +4.0)*t110+t91)*t111+(t95+RATIONAL(6.0,1.0)*t112)*x+t123; + coeffs_dx->coeff_p2_0 = -x+t103*t126+(t106+RATIONAL(-15.0,4.0)*t110)*t111 ++t121; + coeffs_dx->coeff_m1_p1 = t84-t110+t94*t119+t75+(t94*t112+y)*x+t115; + coeffs_dx->coeff_0_p1 = t77+RATIONAL(9.0,1.0)*t119+(RATIONAL(15.0,2.0)* +t112+t107*y)*x+t123; + coeffs_dx->coeff_p1_p1 = t85+t110+t89+RATIONAL(-9.0,1.0)*t119+t71+t79+( +RATIONAL(8.0,1.0)*t110+RATIONAL(-6.0,1.0)*t112+t104*y)*x; + coeffs_dx->coeff_p2_p1 = t97*t119+t81+(t106*t112+t109*y)*x+t124; + coeffs_dx->coeff_m1_p2 = t126+t74+t86+(t117-x+t108)*t110; + coeffs_dx->coeff_0_p2 = t96*t119+t107*t126+t121; + coeffs_dx->coeff_p1_p2 = t88+t104*t126+(t102+t120)*t110+t124; + coeffs_dx->coeff_p2_p2 = t80+t110*t118+(t101*t110+t109*t112)*x; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dxx.compute.c new file mode 100644 index 0000000..656e57f --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dxx.compute.c @@ -0,0 +1,106 @@ +fp t238; +fp t216; +fp t237; +fp t223; +fp t224; +fp t236; +fp t235; +fp t234; +fp t208; +fp t209; +fp t192; +fp t233; +fp t232; +fp t206; +fp t207; +fp t195; +fp t231; +fp t212; +fp t230; +fp t214; +fp t229; +fp t218; +fp t228; +fp t227; +fp t226; +fp t222; +fp t221; +fp t220; +fp t219; +fp t217; +fp t215; +fp t213; +fp t211; +fp t210; +fp t205; +fp t204; +fp t203; +fp t197; +fp t196; +fp t194; +fp t193; +fp t190; +fp t189; + t238 = x*y; + t216 = RATIONAL(-3.0,2.0); + t237 = t216*x; + t223 = y*y; + t224 = t223*y; + t236 = y+t224; + t235 = x*t224; + t234 = x*t223; + t208 = RATIONAL(5.0,2.0); + t209 = RATIONAL(9.0,2.0); + t192 = t209*t235; + t233 = t208*t223+t192; + t232 = RATIONAL(-10.0,1.0)*t223+RATIONAL(-27.0,2.0)*t235; + t206 = RATIONAL(-2.0,1.0); + t207 = RATIONAL(-9.0,2.0); + t195 = t207*t235; + t231 = t206*t223+t195; + t212 = RATIONAL(4.0,1.0); + t230 = t212*t223+t192; + t214 = RATIONAL(-5.0,1.0); + t229 = t214*t223+t195; + t218 = RATIONAL(3.0,2.0); + t228 = t218*t224; + t227 = t216*t224; + t226 = -t223+x*t227; + t222 = RATIONAL(-9.0,1.0); + t221 = RATIONAL(-6.0,1.0); + t220 = RATIONAL(-1.0,2.0); + t219 = RATIONAL(-5.0,2.0); + t217 = RATIONAL(2.0,1.0); + t215 = RATIONAL(1.0,2.0); + t213 = RATIONAL(3.0,1.0); + t211 = RATIONAL(-3.0,1.0); + t210 = RATIONAL(9.0,1.0); + t205 = RATIONAL(6.0,1.0); + t204 = RATIONAL(15.0,2.0); + t203 = RATIONAL(-15.0,2.0); + t197 = x*t228; + t196 = t218*t238; + t194 = t207*t238; + t193 = t209*t238; + t190 = y*t237; + t189 = RATIONAL(27.0,2.0)*t235; + coeffs_dxx->coeff_m1_m1 = t196+t197+(t217+t211*x)*t223-t236; + coeffs_dxx->coeff_0_m1 = t210*t234+t194+t236*t208+t229; + coeffs_dxx->coeff_p1_m1 = t222*t234+t193+t236*t206+t230; + coeffs_dxx->coeff_p2_m1 = t190+t213*t234+t236*t215+t226; + coeffs_dxx->coeff_m1_0 = t213*t224+t217+(t204*t223+t211)*x+t229; + coeffs_dxx->coeff_0_0 = t214+t203*t224+RATIONAL(25.0,2.0)*t223+t189+(t210 ++RATIONAL(-45.0,2.0)*t223)*x; + coeffs_dxx->coeff_p1_0 = t212+t205*t224+(t222+RATIONAL(45.0,2.0)*t223)*x+ +t232; + coeffs_dxx->coeff_p2_0 = RATIONAL(-1.0,1.0)+t227+(t213+t203*t223)*x+t233; + coeffs_dxx->coeff_m1_p1 = y+t221*t234+t211*t224+t190+t230; + coeffs_dxx->coeff_0_p1 = RATIONAL(18.0,1.0)*t234+t219*y+t193+t204*t224+ +t232; + coeffs_dxx->coeff_p1_p1 = t189+t194+t217*y+t221*t224+(RATIONAL(-18.0,1.0) +*x+RATIONAL(8.0,1.0))*t223; + coeffs_dxx->coeff_p2_p1 = t196+t220*y+t228+t205*t234+t231; + coeffs_dxx->coeff_m1_p2 = t218*t234+t224+t226; + coeffs_dxx->coeff_0_p2 = t207*t234+t219*t224+t233; + coeffs_dxx->coeff_p1_p2 = t217*t224+t209*t234+t231; + coeffs_dxx->coeff_p2_p2 = t197+t220*t224+(t215+t237)*t223; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dxy.compute.c new file mode 100644 index 0000000..3ccbd01 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dxy.compute.c @@ -0,0 +1,121 @@ +fp t298; +fp t288; +fp t297; +fp t287; +fp t296; +fp t291; +fp t295; +fp t239; +fp t294; +fp t266; +fp t271; +fp t293; +fp t240; +fp t292; +fp t290; +fp t289; +fp t286; +fp t285; +fp t284; +fp t283; +fp t282; +fp t281; +fp t280; +fp t279; +fp t278; +fp t277; +fp t276; +fp t275; +fp t274; +fp t273; +fp t272; +fp t270; +fp t269; +fp t268; +fp t267; +fp t264; +fp t263; +fp t262; +fp t261; +fp t260; +fp t259; +fp t258; +fp t257; +fp t256; +fp t255; +fp t254; +fp t253; +fp t252; +fp t251; +fp t246; +fp t244; + t298 = x*y; + t288 = x*x; + t297 = t288*y; + t287 = y*y; + t296 = x*t287; + t291 = t288*t287; + t295 = RATIONAL(-81.0,4.0)*t291+RATIONAL(-20.0,1.0)*t298; + t239 = RATIONAL(27.0,4.0)*t291; + t294 = t239+RATIONAL(5.0,1.0)*t298; + t266 = RATIONAL(-9.0,4.0); + t271 = RATIONAL(-2.0,1.0); + t293 = t271*t298+t266*t291; + t240 = RATIONAL(-27.0,4.0)*t291; + t292 = t240+RATIONAL(-4.0,1.0)*t298; + t290 = RATIONAL(-10.0,1.0)*t298+t240; + t289 = RATIONAL(-1.0,4.0)+t239+RATIONAL(8.0,1.0)*t298; + t286 = RATIONAL(-9.0,1.0); + t285 = RATIONAL(-6.0,1.0); + t284 = RATIONAL(-3.0,4.0); + t283 = RATIONAL(-1.0,2.0); + t282 = RATIONAL(-5.0,2.0); + t281 = RATIONAL(3.0,2.0); + t280 = RATIONAL(2.0,1.0); + t279 = RATIONAL(-3.0,2.0); + t278 = RATIONAL(1.0,2.0); + t277 = RATIONAL(3.0,1.0); + t276 = RATIONAL(-3.0,1.0); + t275 = RATIONAL(9.0,1.0); + t274 = RATIONAL(9.0,2.0); + t273 = RATIONAL(5.0,2.0); + t272 = RATIONAL(-9.0,2.0); + t270 = RATIONAL(3.0,4.0); + t269 = RATIONAL(9.0,4.0); + t268 = RATIONAL(1.0,4.0); + t267 = RATIONAL(6.0,1.0); + t264 = RATIONAL(45.0,2.0); + t263 = RATIONAL(-45.0,2.0); + t262 = RATIONAL(18.0,1.0); + t261 = RATIONAL(15.0,2.0); + t260 = RATIONAL(-18.0,1.0); + t259 = RATIONAL(-15.0,2.0); + t258 = t284*t287; + t257 = t284*t288; + t256 = t266*t287; + t255 = t269*t287; + t254 = t270*t288; + t253 = t270*t287; + t252 = t266*t288; + t251 = t269*t288; + t246 = t269*t291; + t244 = RATIONAL(81.0,4.0)*t291; + coeffs_dxy->coeff_m1_m1 = t246+t268+t276*t297+t254+t253-y+(t276*t287-1.0+ +RATIONAL(4.0,1.0)*y)*x; + coeffs_dxy->coeff_0_m1 = t275*t297+t252+(t273+t261*t287)*x+t290; + coeffs_dxy->coeff_p1_m1 = t251+t286*t297+y+t258+(t271+t285*t287)*x+t289; + coeffs_dxy->coeff_p2_m1 = t257+t277*t297+(t278+t281*t287)*x+t293; + coeffs_dxy->coeff_m1_0 = t275*t296+t256+(t273+t261*t288)*y+t290; + coeffs_dxy->coeff_0_0 = t244+t263*t297+(RATIONAL(25.0,1.0)*y+t263*t287)*x +; + coeffs_dxy->coeff_p1_0 = t262*t296+t255+(t264*t288+t282)*y+t295; + coeffs_dxy->coeff_p2_0 = t259*t297+t272*t296+t294; + coeffs_dxy->coeff_m1_p1 = x+t255+t286*t296+t257+(t271+t285*t288)*y+t289; + coeffs_dxy->coeff_0_p1 = t251+t262*t297+(t264*t287+t282)*x+t295; + coeffs_dxy->coeff_p1_p1 = t256+t244+t268+t252+(t260*t288+t280)*y+(t260* +t287+RATIONAL(16.0,1.0)*y+t280)*x; + coeffs_dxy->coeff_p2_p1 = t267*t297+t254+(t283+t274*t287)*x+t292; + coeffs_dxy->coeff_m1_p2 = t277*t296+t258+(t281*t288+t278)*y+t293; + coeffs_dxy->coeff_0_p2 = t272*t297+t259*t296+t294; + coeffs_dxy->coeff_p1_p2 = t267*t296+t253+(t283+t274*t288)*y+t292; + coeffs_dxy->coeff_p2_p2 = t246+t279*t297+(y+t279*t287)*x; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dy.compute.c new file mode 100644 index 0000000..cc35204 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dy.compute.c @@ -0,0 +1,123 @@ +fp t171; +fp t187; +fp t170; +fp t186; +fp t185; +fp t169; +fp t184; +fp t157; +fp t177; +fp t141; +fp t183; +fp t179; +fp t182; +fp t151; +fp t175; +fp t136; +fp t181; +fp t149; +fp t180; +fp t159; +fp t178; +fp t143; +fp t164; +fp t176; +fp t174; +fp t154; +fp t161; +fp t173; +fp t168; +fp t167; +fp t166; +fp t165; +fp t163; +fp t162; +fp t160; +fp t156; +fp t153; +fp t152; +fp t150; +fp t148; +fp t147; +fp t146; +fp t144; +fp t142; +fp t139; +fp t137; +fp t135; +fp t134; +fp t131; +fp t130; + t171 = x*x; + t187 = y*t171; + t170 = x*t171; + t186 = t170*y; + t185 = t170+x; + t169 = y*y; + t184 = t171*t169; + t157 = RATIONAL(-5.0,1.0); + t177 = RATIONAL(-9.0,4.0)*t169; + t141 = t170*t177; + t183 = t157*t187+t141; + t179 = t169*t170; + t182 = RATIONAL(-10.0,1.0)*t187+RATIONAL(-27.0,4.0)*t179; + t151 = RATIONAL(5.0,2.0); + t175 = RATIONAL(9.0,4.0)*t169; + t136 = t170*t175; + t181 = t151*t187+t136; + t149 = RATIONAL(-2.0,1.0); + t180 = t149*t187+t141; + t159 = RATIONAL(3.0,4.0); + t178 = t159*t169; + t143 = t159*t170; + t164 = RATIONAL(-3.0,4.0); + t176 = t164*t169; + t174 = y*t185; + t154 = RATIONAL(4.0,1.0); + t161 = RATIONAL(-1.0,4.0); + t173 = t136+t154*t187+t161*x; + t168 = RATIONAL(-1.0,2.0); + t167 = RATIONAL(1.0,4.0); + t166 = RATIONAL(-5.0,2.0); + t165 = RATIONAL(3.0,2.0); + t163 = RATIONAL(2.0,1.0); + t162 = RATIONAL(-3.0,2.0); + t160 = RATIONAL(1.0,2.0); + t156 = RATIONAL(3.0,1.0); + t153 = RATIONAL(-3.0,1.0); + t152 = RATIONAL(9.0,2.0); + t150 = RATIONAL(-9.0,2.0); + t148 = t167*x; + t147 = t167*t170; + t146 = t164*t170; + t144 = t161*t170; + t142 = t169*t143; + t139 = x*t176; + t137 = t170*t176; + t135 = x*t177; + t134 = x*t175; + t131 = x*t178; + t130 = RATIONAL(27.0,4.0)*t179; + coeffs_dy->coeff_m1_m1 = t148+t147+t142+t131+(t168+t162*t169)*t171+(t163* +t171-t185)*y; + coeffs_dy->coeff_0_m1 = t146+RATIONAL(5.0,4.0)*t171+t168+(RATIONAL(15.0, +4.0)*t171+t162)*t169+(t156*t170+t163)*y+t183; + coeffs_dy->coeff_p1_m1 = -t171+t153*t184+t143+t139+(x+t153*t170)*y+t173; + coeffs_dy->coeff_p2_m1 = t137+t144+t186+(t167-y+t178)*t171; + coeffs_dy->coeff_m1_0 = t135+t152*t184+t151*t174+t183; + coeffs_dy->coeff_0_0 = t130+(t152+RATIONAL(-45.0,4.0)*t171)*t169+( +RATIONAL(-15.0,2.0)*t170+RATIONAL(25.0,2.0)*t171+t157)*y; + coeffs_dy->coeff_p1_0 = RATIONAL(9.0,1.0)*t184+t134+(RATIONAL(15.0,2.0)* +t170+t166*x)*y+t182; + coeffs_dy->coeff_p2_0 = t166*t186+t171*t177+t181; + coeffs_dy->coeff_m1_p1 = t144+t134+(t160+t150*t169)*t171+t149*t174+t173; + coeffs_dy->coeff_0_p1 = t143+t160+RATIONAL(-5.0,4.0)*t171+(t150+RATIONAL( +45.0,4.0)*t171)*t169+(t154+RATIONAL(6.0,1.0)*t170)*y+t182; + coeffs_dy->coeff_p1_p1 = t135+RATIONAL(-9.0,1.0)*t184+t146+t171+t130+t148 ++(t163*x+RATIONAL(8.0,1.0)*t171+RATIONAL(-6.0,1.0)*t170)*y; + coeffs_dy->coeff_p2_p1 = t163*t186+t147+(t175+t161)*t171+t180; + coeffs_dy->coeff_m1_p2 = t137+t165*t184+t139+(-t171+t185*t160)*y; + coeffs_dy->coeff_0_p2 = -y+t162*t186+(RATIONAL(-15.0,4.0)*t171+t165)*t169 ++t181; + coeffs_dy->coeff_p1_p2 = t156*t184+t131+(t168*x+t165*t170)*y+t180; + coeffs_dy->coeff_p2_p2 = t142+t171*t176+(t168*t170+t160*t171)*y; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dyy.compute.c new file mode 100644 index 0000000..dadfd2d --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-dyy.compute.c @@ -0,0 +1,108 @@ +fp t335; +fp t334; +fp t351; +fp t350; +fp t318; +fp t349; +fp t320; +fp t348; +fp t347; +fp t346; +fp t327; +fp t345; +fp t306; +fp t317; +fp t344; +fp t329; +fp t343; +fp t342; +fp t323; +fp t341; +fp t340; +fp t303; +fp t324; +fp t339; +fp t338; +fp t319; +fp t337; +fp t333; +fp t332; +fp t331; +fp t330; +fp t328; +fp t326; +fp t325; +fp t322; +fp t316; +fp t315; +fp t314; +fp t308; +fp t307; +fp t305; +fp t304; +fp t301; +fp t299; + t335 = x*x; + t334 = t335*x; + t351 = t334*y; + t350 = t334+x; + t318 = RATIONAL(-9.0,2.0); + t349 = y*t318; + t320 = RATIONAL(9.0,2.0); + t348 = y*t320; + t347 = t335*y; + t346 = RATIONAL(9.0,1.0)*y; + t327 = RATIONAL(-3.0,2.0); + t345 = t327*y; + t306 = t334*t349; + t317 = RATIONAL(-2.0,1.0); + t344 = t306+t317*t335; + t329 = RATIONAL(3.0,2.0); + t343 = t329*t334; + t342 = t327*t334; + t323 = RATIONAL(-5.0,1.0); + t341 = t306+t323*t335; + t340 = y*t342-t335; + t303 = t334*t348; + t324 = RATIONAL(4.0,1.0); + t339 = t303+t324*t335; + t338 = RATIONAL(-10.0,1.0)*t335+RATIONAL(-27.0,2.0)*t351; + t319 = RATIONAL(5.0,2.0); + t337 = t319*t335+t303; + t333 = RATIONAL(-9.0,1.0); + t332 = RATIONAL(-6.0,1.0); + t331 = RATIONAL(-1.0,2.0); + t330 = RATIONAL(-5.0,2.0); + t328 = RATIONAL(2.0,1.0); + t326 = RATIONAL(1.0,2.0); + t325 = RATIONAL(3.0,1.0); + t322 = RATIONAL(-3.0,1.0); + t316 = RATIONAL(6.0,1.0); + t315 = RATIONAL(15.0,2.0); + t314 = RATIONAL(-15.0,2.0); + t308 = x*t345; + t307 = x*t348; + t305 = y*t343; + t304 = t329*x*y; + t301 = x*t349; + t299 = RATIONAL(27.0,2.0)*t351; + coeffs_dyy->coeff_m1_m1 = t305+t304+(t322*y+t328)*t335-t350; + coeffs_dyy->coeff_0_m1 = t328+t325*t334+(t322+t315*t335)*y+t341; + coeffs_dyy->coeff_p1_m1 = t322*t334+t308+x+t332*t347+t339; + coeffs_dyy->coeff_p2_m1 = t329*t347+t334+t340; + coeffs_dyy->coeff_m1_0 = t335*t346+t301+t350*t319+t341; + coeffs_dyy->coeff_0_0 = t299+t323+t314*t334+t346+(RATIONAL(-45.0,2.0)*y+ +RATIONAL(25.0,2.0))*t335; + coeffs_dyy->coeff_p1_0 = t330*x+t307+t315*t334+RATIONAL(18.0,1.0)*t347+ +t338; + coeffs_dyy->coeff_p2_0 = t318*t347+t330*t334+t337; + coeffs_dyy->coeff_m1_p1 = t333*t347+t307+t350*t317+t339; + coeffs_dyy->coeff_0_p1 = t324+t316*t334+(RATIONAL(45.0,2.0)*t335+t333)*y+ +t338; + coeffs_dyy->coeff_p1_p1 = t328*x+t301+t332*t334+t299+(RATIONAL(8.0,1.0)+ +RATIONAL(-18.0,1.0)*y)*t335; + coeffs_dyy->coeff_p2_p1 = t320*t347+t328*t334+t344; + coeffs_dyy->coeff_m1_p2 = t308+t325*t347+t350*t326+t340; + coeffs_dyy->coeff_0_p2 = RATIONAL(-1.0,1.0)+t342+(t314*t335+t325)*y+t337; + coeffs_dyy->coeff_p1_p2 = t316*t347+t343+t331*x+t304+t344; + coeffs_dyy->coeff_p2_p2 = t305+t331*t334+(t326+t345)*t335; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-I.compute.c new file mode 100644 index 0000000..a5b6360 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-I.compute.c @@ -0,0 +1,337 @@ +fp t469; +fp t472; +fp t517; +fp t471; +fp t470; +fp t516; +fp t515; +fp t514; +fp t409; +fp t513; +fp t416; +fp t512; +fp t511; +fp t435; +fp t510; +fp t509; +fp t423; +fp t508; +fp t507; +fp t506; +fp t413; +fp t487; +fp t368; +fp t489; +fp t445; +fp t505; +fp t410; +fp t504; +fp t499; +fp t503; +fp t502; +fp t449; +fp t501; +fp t453; +fp t491; +fp t401; +fp t441; +fp t500; +fp t442; +fp t490; +fp t408; +fp t498; +fp t497; +fp t421; +fp t366; +fp t496; +fp t495; +fp t403; +fp t494; +fp t404; +fp t493; +fp t424; +fp t492; +fp t434; +fp t488; +fp t450; +fp t371; +fp t486; +fp t485; +fp t484; +fp t389; +fp t483; +fp t444; +fp t482; +fp t481; +fp t480; +fp t479; +fp t426; +fp t428; +fp t478; +fp t475; +fp t468; +fp t467; +fp t466; +fp t465; +fp t464; +fp t463; +fp t462; +fp t461; +fp t460; +fp t459; +fp t458; +fp t457; +fp t456; +fp t454; +fp t452; +fp t451; +fp t448; +fp t447; +fp t446; +fp t443; +fp t440; +fp t439; +fp t438; +fp t437; +fp t436; +fp t433; +fp t432; +fp t431; +fp t430; +fp t429; +fp t427; +fp t425; +fp t422; +fp t420; +fp t419; +fp t418; +fp t417; +fp t415; +fp t414; +fp t412; +fp t411; +fp t407; +fp t406; +fp t405; +fp t402; +fp t400; +fp t396; +fp t395; +fp t394; +fp t393; +fp t391; +fp t386; +fp t383; +fp t382; +fp t380; +fp t379; +fp t374; +fp t373; +fp t372; +fp t369; +fp t364; +fp t362; +fp t361; +fp t359; +fp t358; +fp t357; +fp t356; +fp t355; +fp t354; +fp t352; + t469 = y*y; + t472 = t469*y; + t517 = t472+y; + t471 = x*x; + t470 = x*t471; + t516 = t470+x; + t515 = x*y; + t514 = x*t472; + t409 = RATIONAL(1.0,144.0); + t513 = x*t409; + t416 = RATIONAL(-1.0,144.0); + t512 = t416*x; + t511 = y*t470; + t435 = RATIONAL(1.0,18.0); + t510 = t435*x; + t509 = t435*y; + t423 = RATIONAL(-1.0,18.0); + t508 = t423*x; + t507 = t470*t469; + t506 = t471*t472; + t413 = RATIONAL(-7.0,144.0); + t487 = t413*t470; + t368 = t472*t487; + t489 = t469*t471; + t445 = RATIONAL(-1.0,24.0); + t505 = t445*t489+t368; + t410 = RATIONAL(7.0,144.0); + t504 = t410*t472; + t499 = t470*t472; + t503 = RATIONAL(-35.0,9.0)*t489+RATIONAL(-16.0,9.0)*t499; + t502 = t410*t470; + t449 = RATIONAL(1.0,9.0); + t501 = t449*t472; + t453 = RATIONAL(-1.0,9.0); + t491 = t453*t472; + t401 = t470*t491; + t441 = RATIONAL(-7.0,36.0); + t500 = t441*t489+t401; + t442 = RATIONAL(7.0,18.0); + t490 = t449*t470; + t408 = t472*t490; + t498 = t442*t489+t408; + t497 = t413*t472; + t421 = RATIONAL(5.0,48.0); + t366 = t410*t499; + t496 = t421*t489+t366; + t495 = t453*t470; + t403 = RATIONAL(7.0,9.0)*t499; + t494 = RATIONAL(7.0,6.0)*t489+t403; + t404 = RATIONAL(-7.0,9.0)*t499; + t493 = t404+RATIONAL(-35.0,12.0)*t489; + t424 = RATIONAL(-1.0,72.0); + t492 = t424*t489+t416*t499; + t434 = RATIONAL(5.0,36.0); + t488 = t408+t434*t489; + t450 = RATIONAL(-5.0,6.0); + t371 = y*t508; + t486 = t450*t489+t404+t371; + t485 = RATIONAL(-4.0,9.0)*t515+t403+RATIONAL(25.0,12.0)*t489; + t484 = t368+RATIONAL(-5.0,24.0)*t489+t371; + t389 = x*t509; + t483 = RATIONAL(-5.0,18.0)*t489+t389+t401; + t444 = RATIONAL(1.0,12.0); + t482 = y*t512+t444*t489+t366; + t481 = t389+RATIONAL(-49.0,144.0)*t499+RATIONAL(-5.0,8.0)*t489; + t480 = t516*t469; + t479 = t517*t471; + t426 = RATIONAL(7.0,36.0); + t428 = RATIONAL(-1.0,12.0); + t478 = t428+t426*t471; + t475 = t426*t469+t428; + t468 = RATIONAL(5.0,3.0); + t467 = RATIONAL(-8.0,9.0); + t466 = RATIONAL(5.0,4.0); + t465 = RATIONAL(-2.0,3.0); + t464 = RATIONAL(8.0,9.0); + t463 = RATIONAL(-5.0,3.0); + t462 = RATIONAL(5.0,6.0); + t461 = RATIONAL(1.0,3.0); + t460 = RATIONAL(-1.0,6.0); + t459 = RATIONAL(-1.0,2.0); + t458 = RATIONAL(-2.0,9.0); + t457 = RATIONAL(-1.0,3.0); + t456 = RATIONAL(2.0,9.0); + t454 = RATIONAL(-7.0,3.0); + t452 = RATIONAL(-4.0,3.0); + t451 = RATIONAL(2.0,3.0); + t448 = RATIONAL(1.0,72.0); + t447 = RATIONAL(14.0,9.0); + t446 = RATIONAL(-7.0,12.0); + t443 = RATIONAL(-10.0,9.0); + t440 = RATIONAL(1.0,24.0); + t439 = RATIONAL(-7.0,24.0); + t438 = RATIONAL(-7.0,18.0); + t437 = RATIONAL(-5.0,48.0); + t436 = RATIONAL(7.0,24.0); + t433 = RATIONAL(-5.0,36.0); + t432 = RATIONAL(-28.0,9.0); + t431 = RATIONAL(28.0,9.0); + t430 = RATIONAL(20.0,9.0); + t429 = RATIONAL(-14.0,9.0); + t427 = RATIONAL(7.0,12.0); + t425 = RATIONAL(-7.0,72.0); + t422 = RATIONAL(10.0,9.0); + t420 = RATIONAL(7.0,72.0); + t419 = RATIONAL(-20.0,9.0); + t418 = RATIONAL(-35.0,48.0); + t417 = RATIONAL(49.0,36.0); + t415 = RATIONAL(-35.0,36.0); + t414 = RATIONAL(35.0,48.0); + t412 = RATIONAL(35.0,36.0); + t411 = RATIONAL(-49.0,36.0); + t407 = RATIONAL(4.0,9.0)*t515; + t406 = y*t490; + t405 = t467*t511; + t402 = y*t495; + t400 = t464*t511; + t396 = x*t501; + t395 = t467*t514; + t394 = t464*t514; + t393 = x*t491; + t391 = t442*t511; + t386 = t472*t510; + t383 = t438*t514; + t382 = t438*t511; + t380 = t470*t509; + t379 = RATIONAL(16.0,9.0)*t499; + t374 = t423*t511; + t373 = t472*t508; + t372 = t442*t514; + t369 = y*t487; + t364 = y*t502; + t362 = t409*t511; + t361 = t472*t512; + t359 = x*t504; + t358 = t416*t511; + t357 = t472*t513; + t356 = x*t497; + t355 = t409*t499; + t354 = y*t513; + t352 = RATIONAL(49.0,144.0)*t499; + coeffs_I->coeff_m2_m2 = t354+t357+t362+t355+RATIONAL(1.0,36.0)*t489+(t479 ++t480)*t424; + coeffs_I->coeff_m1_m2 = t369+t373+t421*t479+(t420*t470+t449*x)*t469+t484; + coeffs_I->coeff_0_m2 = t406+(t458*t470+t460)*t469+t498+t517*(t444+t441* +t471); + coeffs_I->coeff_p1_m2 = t386+t402+t434*t479+(t456*t470+t453*x)*t469+t483; + coeffs_I->coeff_p2_m2 = t361+t364+t445*t479+(t425*t470+t448*x)*t469+t482; + coeffs_I->coeff_p3_m2 = t358+t448*t507+t409*t479+t492; + coeffs_I->coeff_m2_m1 = t356+t374+(t420*t472+t449*y)*t471+t421*t480+t484; + coeffs_I->coeff_m1_m1 = t372+t407+t352+t391+(t450*y+t418*t472)*t471+(t450 +*x+RATIONAL(25.0,16.0)*t471+t418*t470)*t469; + coeffs_I->coeff_0_m1 = t405+(t417*t471+t446)*t472+(t468*t470+t466)*t469+( +t447*t471+t465)*y+t493; + coeffs_I->coeff_p1_m1 = t400+t383+(t415*t472+t443*y)*t471+(t462*x+t463* +t470)*t469+t485; + coeffs_I->coeff_p2_m1 = t359+t382+(t436*t472+t461*y)*t471+(t414*t470+t437 +*x)*t469+t481; + coeffs_I->coeff_p3_m1 = t437*t507+t380+(t423*y+t497)*t471+t496; + coeffs_I->coeff_m2_0 = t396+(t460+t458*t472)*t471+t498+t516*(t444+t441* +t469); + coeffs_I->coeff_m1_0 = t395+(t468*t472+t466)*t471+(t417*t469+t446)*t470+( +t465+t447*t469)*x+t493; + coeffs_I->coeff_0_0 = t379+RATIONAL(1.0,1.0)+(t470+t472)*RATIONAL(4.0,3.0 +)+(t454+t432*t472)*t471+(RATIONAL(49.0,9.0)*t471+t454+t432*t470)*t469; + coeffs_I->coeff_p1_0 = t394+(t430*t472+t468)*t471+(t452+t431*t469)*t470+( +t429*t469+t451)*x+t503; + coeffs_I->coeff_p2_0 = t393+(t465*t472+t459)*t471+(t411*t469+t427)*t470+ +t475*x+t494; + coeffs_I->coeff_p3_0 = (t444+t501)*t471+t475*t470+t500; + coeffs_I->coeff_m2_p1 = t393+t380+(t456*t472+t453*y)*t471+t434*t480+t483; + coeffs_I->coeff_m1_p1 = t382+t394+(t462*y+t463*t472)*t471+(t443*x+t415* +t470)*t469+t485; + coeffs_I->coeff_0_p1 = t400+(t452+t431*t471)*t472+(t468+t430*t470)*t469+( +t451+t429*t471)*y+t503; + coeffs_I->coeff_p1_p1 = t379+t405+t407+t395+(t422*y+t419*t472)*t471+(t422 +*x+t419*t470+RATIONAL(25.0,9.0)*t471)*t469; + coeffs_I->coeff_p2_p1 = t391+t396+(t457*y+t451*t472)*t471+(t412*t470+t433 +*x)*t469+t486; + coeffs_I->coeff_p3_p1 = t433*t507+t374+(t509+t491)*t471+t488; + coeffs_I->coeff_m2_p2 = t359+t358+(t425*t472+t448*y)*t471+t445*t480+t482; + coeffs_I->coeff_m1_p2 = t383+t364+(t414*t472+t437*y)*t471+(t461*x+t436* +t470)*t469+t481; + coeffs_I->coeff_0_p2 = t402+(t411*t471+t427)*t472+(t459+t465*t470)*t469+ +t478*y+t494; + coeffs_I->coeff_p1_p2 = t372+t406+(t412*t472+t433*y)*t471+(t457*x+t451* +t470)*t469+t486; + coeffs_I->coeff_p2_p2 = t369+t354+t352+t356+(t439*t472+t440*y)*t471+(t440 +*x+t439*t470+RATIONAL(1.0,4.0)*t471)*t469; + coeffs_I->coeff_p3_p2 = t440*t507+t362+(t416*y+t504)*t471+t505; + coeffs_I->coeff_m2_p3 = t361+t448*t506+t409*t480+t492; + coeffs_I->coeff_m1_p3 = t386+t437*t506+(t508+t487)*t469+t496; + coeffs_I->coeff_0_p3 = t478*t472+(t490+t444)*t469+t500; + coeffs_I->coeff_p1_p3 = t373+t433*t506+(t510+t495)*t469+t488; + coeffs_I->coeff_p2_p3 = t440*t506+t357+(t512+t502)*t469+t505; + coeffs_I->coeff_p3_p3 = t355+t409*t489+(t507+t506)*t416; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dx.compute.c new file mode 100644 index 0000000..8d83f9f --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dx.compute.c @@ -0,0 +1,276 @@ +fp t616; +fp t618; +fp t648; +fp t617; +fp t647; +fp t646; +fp t645; +fp t644; +fp t635; +fp t608; +fp t548; +fp t643; +fp t551; +fp t642; +fp t597; +fp t567; +fp t610; +fp t628; +fp t540; +fp t641; +fp t577; +fp t640; +fp t607; +fp t631; +fp t539; +fp t592; +fp t639; +fp t627; +fp t580; +fp t638; +fp t579; +fp t637; +fp t636; +fp t634; +fp t633; +fp t535; +fp t578; +fp t632; +fp t630; +fp t588; +fp t520; +fp t629; +fp t626; +fp t583; +fp t564; +fp t612; +fp t625; +fp t624; +fp t591; +fp t563; +fp t623; +fp t611; +fp t622; +fp t621; +fp t614; +fp t576; +fp t620; +fp t615; +fp t613; +fp t609; +fp t606; +fp t605; +fp t604; +fp t603; +fp t602; +fp t601; +fp t600; +fp t599; +fp t598; +fp t595; +fp t594; +fp t593; +fp t589; +fp t587; +fp t586; +fp t585; +fp t584; +fp t575; +fp t574; +fp t573; +fp t572; +fp t571; +fp t570; +fp t568; +fp t566; +fp t565; +fp t562; +fp t560; +fp t559; +fp t558; +fp t557; +fp t556; +fp t553; +fp t552; +fp t550; +fp t547; +fp t544; +fp t543; +fp t542; +fp t536; +fp t534; +fp t532; +fp t524; +fp t522; +fp t521; +fp t519; + t616 = y*y; + t618 = t616*y; + t648 = x*t618; + t617 = x*x; + t647 = t617*y; + t646 = y+t618; + t645 = x*t616; + t644 = t617*t616; + t635 = t617*t618; + t608 = RATIONAL(7.0,3.0); + t548 = t608*t635; + t643 = t548+t608*t645; + t551 = RATIONAL(-7.0,3.0)*t635; + t642 = RATIONAL(-35.0,6.0)*t645+t551; + t597 = RATIONAL(7.0,18.0); + t567 = t597*t618; + t610 = RATIONAL(1.0,3.0); + t628 = t617*t610; + t540 = t618*t628; + t641 = t540+RATIONAL(7.0,9.0)*t645; + t577 = RATIONAL(5.0,18.0); + t640 = t577*t645+t540; + t607 = RATIONAL(-1.0,3.0); + t631 = t607*t617; + t539 = t618*t631; + t592 = RATIONAL(-7.0,18.0); + t639 = t592*t645+t539; + t627 = t617*RATIONAL(-1.0,48.0); + t580 = RATIONAL(-1.0,36.0); + t638 = t618*t627+t580*t645; + t579 = RATIONAL(1.0,48.0); + t637 = t579*t617; + t636 = RATIONAL(-70.0,9.0)*t645+RATIONAL(-16.0,3.0)*t635; + t634 = t616*RATIONAL(-5.0,12.0); + t633 = RATIONAL(7.0,48.0)*t617; + t535 = t618*t633; + t578 = RATIONAL(5.0,24.0); + t632 = t578*t645+t535; + t630 = RATIONAL(-7.0,48.0)*t617; + t588 = RATIONAL(-1.0,12.0); + t520 = t618*t630; + t629 = t520+t588*t645; + t626 = x*t646; + t583 = RATIONAL(-1.0,18.0); + t564 = t583*y; + t612 = RATIONAL(-5.0,3.0); + t625 = t551+t564+t612*t645; + t624 = t564+t520+x*t634; + t591 = RATIONAL(1.0,18.0); + t563 = t591*y; + t623 = t563+RATIONAL(-49.0,48.0)*t635+RATIONAL(-5.0,4.0)*t645; + t611 = RATIONAL(-4.0,9.0); + t622 = RATIONAL(25.0,6.0)*t645+t548+t611*y; + t621 = t563+t539+RATIONAL(-5.0,9.0)*t645; + t614 = RATIONAL(1.0,6.0); + t576 = RATIONAL(-1.0,144.0); + t620 = t535+t576*y+t614*t645; + t615 = RATIONAL(4.0,9.0); + t613 = RATIONAL(-2.0,3.0); + t609 = RATIONAL(-2.0,9.0); + t606 = RATIONAL(2.0,9.0); + t605 = RATIONAL(-1.0,9.0); + t604 = RATIONAL(2.0,3.0); + t603 = RATIONAL(1.0,9.0); + t602 = RATIONAL(1.0,72.0); + t601 = RATIONAL(-5.0,24.0); + t600 = RATIONAL(-7.0,12.0); + t599 = RATIONAL(-1.0,24.0); + t598 = RATIONAL(1.0,12.0); + t595 = RATIONAL(-7.0,36.0); + t594 = RATIONAL(1.0,24.0); + t593 = RATIONAL(10.0,3.0); + t589 = RATIONAL(-5.0,18.0); + t587 = RATIONAL(7.0,12.0); + t586 = RATIONAL(7.0,36.0); + t585 = RATIONAL(-1.0,72.0); + t584 = RATIONAL(1.0,36.0); + t575 = RATIONAL(1.0,144.0); + t574 = t603*t618; + t573 = t615*y; + t572 = RATIONAL(-8.0,9.0)*t618; + t571 = RATIONAL(8.0,9.0)*t618; + t570 = t605*t618; + t568 = t591*t618; + t566 = t592*t618; + t565 = t583*t618; + t562 = t576*t618; + t560 = RATIONAL(-7.0,144.0)*t618; + t559 = RATIONAL(7.0,144.0)*t618; + t558 = t575*y; + t557 = t575*t618; + t556 = y*t631; + t553 = RATIONAL(7.0,6.0)*t647; + t552 = RATIONAL(-1.0,6.0)*t647; + t550 = RATIONAL(-8.0,3.0)*t647; + t547 = RATIONAL(8.0,3.0)*t647; + t544 = RATIONAL(-7.0,6.0)*t647; + t543 = t614*t647; + t542 = y*t628; + t536 = y*t633; + t534 = y*t627; + t532 = RATIONAL(16.0,3.0)*t635; + t524 = t579*t635; + t522 = y*t630; + t521 = y*t637; + t519 = RATIONAL(49.0,48.0)*t635; + coeffs_dx->coeff_m2_m2 = t558+t521+t557+t524+(t599*t617+t585)*t616+(t591* +t616+t646*t580)*x; + coeffs_dx->coeff_m1_m2 = t522+t565+(RATIONAL(7.0,24.0)*t617+t603)*t616+ +t578*t626+t624; + coeffs_dx->coeff_0_m2 = t613*t644+t542+t592*t626+t641; + coeffs_dx->coeff_p1_m2 = t568+t556+(t605+t604*t617)*t616+t577*t626+t621; + coeffs_dx->coeff_p2_m2 = t562+t536+(t602+RATIONAL(-7.0,24.0)*t617)*t616+ +t588*t626+t620; + coeffs_dx->coeff_p3_m2 = t534+t594*t644+t602*t626+t638; + coeffs_dx->coeff_m2_m1 = t560+t552+(RATIONAL(5.0,48.0)+RATIONAL(5.0,16.0) +*t617)*t616+(t606*y+t586*t618)*x+t624; + coeffs_dx->coeff_m1_m1 = t567+t553+t573+t519+(RATIONAL(-35.0,16.0)*t617+ +RATIONAL(-5.0,6.0))*t616+(RATIONAL(-35.0,24.0)*t618+t612*y+RATIONAL(25.0,8.0)* +t616)*x; + coeffs_dx->coeff_0_m1 = RATIONAL(5.0,1.0)*t644+t550+(RATIONAL(49.0,18.0)* +t618+RATIONAL(28.0,9.0)*y)*x+t642; + coeffs_dx->coeff_p1_m1 = t547+t566+(RATIONAL(-5.0,1.0)*t617+RATIONAL(5.0, +6.0))*t616+(RATIONAL(-20.0,9.0)*y+RATIONAL(-35.0,18.0)*t618)*x+t622; + coeffs_dx->coeff_p2_m1 = t544+t559+(RATIONAL(35.0,16.0)*t617+RATIONAL( +-5.0,48.0))*t616+(t587*t618+t604*y)*x+t623; + coeffs_dx->coeff_p3_m1 = t543+RATIONAL(-5.0,16.0)*t644+(t605*y+RATIONAL( +-7.0,72.0)*t618)*x+t632; + coeffs_dx->coeff_m2_0 = t595*t616+t574+t598+(t600*t616+RATIONAL(1.0,4.0)) +*t617+(t611*t618+t607)*x+t641; + coeffs_dx->coeff_m1_0 = t613+t572+RATIONAL(14.0,9.0)*t616+(RATIONAL(-7.0, +4.0)+RATIONAL(49.0,12.0)*t616)*t617+(RATIONAL(5.0,2.0)+t593*t618)*x+t642; + coeffs_dx->coeff_0_0 = t532+(RATIONAL(-28.0,3.0)*t616+RATIONAL(4.0,1.0))* +t617+(RATIONAL(-14.0,3.0)+RATIONAL(-56.0,9.0)*t618+RATIONAL(98.0,9.0)*t616)*x; + coeffs_dx->coeff_p1_0 = t571+t604+RATIONAL(-14.0,9.0)*t616+(RATIONAL(-4.0 +,1.0)+RATIONAL(28.0,3.0)*t616)*t617+(RATIONAL(40.0,9.0)*t618+t593)*x+t636; + coeffs_dx->coeff_p2_0 = t570-x+t588+RATIONAL(-4.0,3.0)*t648+t586*t616+( +RATIONAL(-49.0,12.0)*t616+RATIONAL(7.0,4.0))*t617+t643; + coeffs_dx->coeff_p3_0 = (t587*t616+RATIONAL(-1.0,4.0))*t617+(t606*t618+ +t614)*x+t639; + coeffs_dx->coeff_m2_p1 = t543+t570+(RATIONAL(5.0,12.0)*t617+RATIONAL(5.0, +36.0))*t616+(t615*t618+t609*y)*x+t621; + coeffs_dx->coeff_m1_p1 = t571+t544+(RATIONAL(-35.0,12.0)*t617+RATIONAL( +-10.0,9.0))*t616+(RATIONAL(5.0,3.0)*y+RATIONAL(-10.0,3.0)*t618)*x+t622; + coeffs_dx->coeff_0_p1 = t547+RATIONAL(20.0,3.0)*t644+(RATIONAL(56.0,9.0)* +t618+RATIONAL(-28.0,9.0)*y)*x+t636; + coeffs_dx->coeff_p1_p1 = t550+t572+t573+t532+(RATIONAL(-20.0,3.0)*t617+ +RATIONAL(10.0,9.0))*t616+(RATIONAL(20.0,9.0)*y+RATIONAL(50.0,9.0)*t616+RATIONAL +(-40.0,9.0)*t618)*x; + coeffs_dx->coeff_p2_p1 = t574+t553+(RATIONAL(35.0,12.0)*t617+RATIONAL( +-5.0,36.0))*t616+(t613*y+RATIONAL(4.0,3.0)*t618)*x+t625; + coeffs_dx->coeff_p3_p1 = t552+t617*t634+(t603*y+t609*t618)*x+t640; + coeffs_dx->coeff_m2_p2 = t559+t534+(t599+RATIONAL(-1.0,8.0)*t617)*t616+( +t584*y+t595*t618)*x+t620; + coeffs_dx->coeff_m1_p2 = t536+t566+(t610+RATIONAL(7.0,8.0)*t617)*t616+( +RATIONAL(35.0,24.0)*t618+t601*y)*x+t623; + coeffs_dx->coeff_0_p2 = t556+RATIONAL(-2.0,1.0)*t644+(RATIONAL(-49.0,18.0 +)*t618+t597*y)*x+t643; + coeffs_dx->coeff_p1_p2 = t542+t567+(t607+RATIONAL(2.0,1.0)*t617)*t616+( +t589*y+RATIONAL(35.0,18.0)*t618)*x+t625; + coeffs_dx->coeff_p2_p2 = t522+t558+t560+t519+(t594+RATIONAL(-7.0,8.0)* +t617)*t616+(t600*t618+RATIONAL(1.0,2.0)*t616+t598*y)*x; + coeffs_dx->coeff_p3_p2 = t521+RATIONAL(1.0,8.0)*t644+(t585*y+RATIONAL(7.0 +,72.0)*t618)*x+t629; + coeffs_dx->coeff_m2_p3 = t562+t584*t648+(t575+t637)*t616+t638; + coeffs_dx->coeff_m1_p3 = t601*t648+t568+(t583+t630)*t616+t632; + coeffs_dx->coeff_0_p3 = t616*t628+x*t567+t639; + coeffs_dx->coeff_p1_p3 = t565+t589*t648+(t631+t591)*t616+t640; + coeffs_dx->coeff_p2_p3 = t557+t598*t648+(t633+t576)*t616+t629; + coeffs_dx->coeff_p3_p3 = t524+t616*t627+(t585*t618+t602*t616)*x; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dxx.compute.c new file mode 100644 index 0000000..7650439 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dxx.compute.c @@ -0,0 +1,215 @@ +fp t875; +fp t840; +fp t874; +fp t851; +fp t852; +fp t873; +fp t836; +fp t872; +fp t849; +fp t871; +fp t870; +fp t869; +fp t848; +fp t830; +fp t795; +fp t868; +fp t829; +fp t785; +fp t825; +fp t867; +fp t810; +fp t866; +fp t819; +fp t796; +fp t865; +fp t828; +fp t799; +fp t864; +fp t788; +fp t863; +fp t862; +fp t861; +fp t860; +fp t821; +fp t859; +fp t858; +fp t850; +fp t827; +fp t789; +fp t857; +fp t845; +fp t856; +fp t855; +fp t820; +fp t854; +fp t847; +fp t846; +fp t844; +fp t843; +fp t842; +fp t841; +fp t839; +fp t838; +fp t837; +fp t835; +fp t834; +fp t832; +fp t831; +fp t826; +fp t824; +fp t823; +fp t822; +fp t803; +fp t802; +fp t801; +fp t800; +fp t798; +fp t797; +fp t794; +fp t793; +fp t791; +fp t787; +fp t786; +fp t784; +fp t783; +fp t782; +fp t780; + t875 = x*y; + t840 = RATIONAL(2.0,3.0); + t874 = t840*y; + t851 = y*y; + t852 = t851*y; + t873 = y+t852; + t836 = RATIONAL(-1.0,24.0); + t872 = t836*x; + t849 = RATIONAL(-2.0,3.0); + t871 = t849*y; + t870 = x*t852; + t869 = x*t851; + t848 = RATIONAL(-5.0,3.0); + t830 = RATIONAL(-14.0,3.0); + t795 = t830*t870; + t868 = t848*t851+t795; + t829 = RATIONAL(-7.0,24.0); + t785 = t829*t870; + t825 = RATIONAL(-1.0,12.0); + t867 = t785+t825*t851; + t810 = RATIONAL(-35.0,6.0)*t851; + t866 = t810+t795; + t819 = RATIONAL(5.0,18.0); + t796 = t840*t870; + t865 = t819*t851+t796; + t828 = RATIONAL(-7.0,18.0); + t799 = t849*t870; + t864 = t828*t851+t799; + t788 = RATIONAL(14.0,3.0)*t870; + t863 = t788+RATIONAL(25.0,6.0)*t851; + t862 = RATIONAL(-5.0,12.0)*t851+t785; + t861 = RATIONAL(-32.0,3.0)*t870+RATIONAL(-70.0,9.0)*t851; + t860 = RATIONAL(-5.0,9.0)*t851+t799; + t821 = RATIONAL(-1.0,36.0); + t859 = t836*t870+t821*t851; + t858 = RATIONAL(7.0,9.0)*t851+t796; + t850 = RATIONAL(1.0,6.0); + t827 = RATIONAL(7.0,24.0); + t789 = t827*t870; + t857 = t850*t851+t789; + t845 = RATIONAL(7.0,3.0); + t856 = t788+t845*t851; + t855 = RATIONAL(-49.0,24.0)*t870+RATIONAL(-5.0,4.0)*t851; + t820 = RATIONAL(5.0,24.0); + t854 = t789+t820*t851; + t847 = RATIONAL(-2.0,9.0); + t846 = RATIONAL(1.0,2.0); + t844 = RATIONAL(-1.0,3.0); + t843 = RATIONAL(2.0,9.0); + t842 = RATIONAL(4.0,3.0); + t841 = RATIONAL(-4.0,3.0); + t839 = RATIONAL(1.0,72.0); + t838 = RATIONAL(-5.0,24.0); + t837 = RATIONAL(-7.0,12.0); + t835 = RATIONAL(1.0,12.0); + t834 = RATIONAL(7.0,18.0); + t832 = RATIONAL(1.0,24.0); + t831 = RATIONAL(10.0,3.0); + t826 = RATIONAL(-5.0,18.0); + t824 = RATIONAL(7.0,12.0); + t823 = RATIONAL(-1.0,72.0); + t822 = RATIONAL(1.0,36.0); + t803 = RATIONAL(-7.0,3.0)*t875; + t802 = RATIONAL(1.0,3.0)*t875; + t801 = t845*t875; + t800 = t844*t875; + t798 = x*t871; + t797 = x*t874; + t794 = RATIONAL(32.0,3.0)*t870; + t793 = RATIONAL(-16.0,3.0)*t875; + t791 = y*t872; + t787 = t827*t875; + t786 = RATIONAL(16.0,3.0)*t875; + t784 = t829*t875; + t783 = t832*t870; + t782 = t832*t875; + t780 = RATIONAL(49.0,24.0)*t870; + coeffs_dxx->coeff_m2_m2 = t782+t783+t873*t821+(t825*x+RATIONAL(1.0,18.0)) +*t851; + coeffs_dxx->coeff_m1_m2 = t784+t824*t869+t873*t820+t862; + coeffs_dxx->coeff_0_m2 = t797+t841*t869+t873*t828+t858; + coeffs_dxx->coeff_p1_m2 = t798+t842*t869+t873*t819+t860; + coeffs_dxx->coeff_p2_m2 = t837*t869+t787+t873*t825+t857; + coeffs_dxx->coeff_p3_m2 = t835*t869+t791+t873*t839+t859; + coeffs_dxx->coeff_m2_m1 = t843*y+RATIONAL(7.0,36.0)*t852+RATIONAL(5.0,8.0 +)*t869+t800+t862; + coeffs_dxx->coeff_m1_m1 = RATIONAL(-35.0,24.0)*t852+t780+t801+t848*y+( +RATIONAL(25.0,8.0)+RATIONAL(-35.0,8.0)*x)*t851; + coeffs_dxx->coeff_0_m1 = RATIONAL(28.0,9.0)*y+t793+RATIONAL(49.0,18.0)* +t852+RATIONAL(10.0,1.0)*t869+t866; + coeffs_dxx->coeff_p1_m1 = t786+RATIONAL(-10.0,1.0)*t869+RATIONAL(-20.0, +9.0)*y+RATIONAL(-35.0,18.0)*t852+t863; + coeffs_dxx->coeff_p2_m1 = RATIONAL(35.0,8.0)*t869+t874+t803+t824*t852+ +t855; + coeffs_dxx->coeff_p3_m1 = RATIONAL(-7.0,72.0)*t852+RATIONAL(-1.0,9.0)*y+ +t802+RATIONAL(-5.0,8.0)*t869+t854; + coeffs_dxx->coeff_m2_0 = RATIONAL(-4.0,9.0)*t852+t844+(t846+RATIONAL(-7.0 +,6.0)*t851)*x+t858; + coeffs_dxx->coeff_m1_0 = t831*t852+RATIONAL(5.0,2.0)+(RATIONAL(-7.0,2.0)+ +RATIONAL(49.0,6.0)*t851)*x+t866; + coeffs_dxx->coeff_0_0 = RATIONAL(-56.0,9.0)*t852+t830+RATIONAL(98.0,9.0)* +t851+t794+(RATIONAL(8.0,1.0)+RATIONAL(-56.0,3.0)*t851)*x; + coeffs_dxx->coeff_p1_0 = t831+RATIONAL(40.0,9.0)*t852+(RATIONAL(56.0,3.0) +*t851+RATIONAL(-8.0,1.0))*x+t861; + coeffs_dxx->coeff_p2_0 = t841*t852+RATIONAL(-1.0,1.0)+(RATIONAL(7.0,2.0)+ +RATIONAL(-49.0,6.0)*t851)*x+t856; + coeffs_dxx->coeff_p3_0 = t843*t852+t850+(RATIONAL(-1.0,2.0)+RATIONAL(7.0, +6.0)*t851)*x+t864; + coeffs_dxx->coeff_m2_p1 = RATIONAL(5.0,6.0)*t869+RATIONAL(4.0,9.0)*t852+ +t802+t847*y+t860; + coeffs_dxx->coeff_m1_p1 = RATIONAL(5.0,3.0)*y+RATIONAL(-10.0,3.0)*t852+ +t803+x*t810+t863; + coeffs_dxx->coeff_0_p1 = RATIONAL(56.0,9.0)*t852+RATIONAL(-28.0,9.0)*y+ +RATIONAL(40.0,3.0)*t869+t786+t861; + coeffs_dxx->coeff_p1_p1 = t793+t794+RATIONAL(20.0,9.0)*y+RATIONAL(-40.0, +9.0)*t852+(RATIONAL(-40.0,3.0)*x+RATIONAL(50.0,9.0))*t851; + coeffs_dxx->coeff_p2_p1 = t801+RATIONAL(35.0,6.0)*t869+t842*t852+t871+ +t868; + coeffs_dxx->coeff_p3_p1 = t800+t847*t852+RATIONAL(1.0,9.0)*y+RATIONAL( +-5.0,6.0)*t869+t865; + coeffs_dxx->coeff_m2_p2 = RATIONAL(-7.0,36.0)*t852+t791+RATIONAL(-1.0,4.0 +)*t869+t822*y+t857; + coeffs_dxx->coeff_m1_p2 = t838*y+RATIONAL(35.0,24.0)*t852+RATIONAL(7.0, +4.0)*t869+t787+t855; + coeffs_dxx->coeff_0_p2 = RATIONAL(-49.0,18.0)*t852+t834*y+t798+RATIONAL( +-4.0,1.0)*t869+t856; + coeffs_dxx->coeff_p1_p2 = t797+RATIONAL(4.0,1.0)*t869+t826*y+RATIONAL( +35.0,18.0)*t852+t868; + coeffs_dxx->coeff_p2_p2 = t837*t852+t835*y+t784+t780+(t846+RATIONAL(-7.0, +4.0)*x)*t851; + coeffs_dxx->coeff_p3_p2 = RATIONAL(7.0,72.0)*t852+t823*y+t782+RATIONAL( +1.0,4.0)*t869+t867; + coeffs_dxx->coeff_m2_p3 = t832*t869+t822*t852+t859; + coeffs_dxx->coeff_m1_p3 = t838*t852+t829*t869+t854; + coeffs_dxx->coeff_0_p3 = t834*t852+t840*t869+t864; + coeffs_dxx->coeff_p1_p3 = t849*t869+t826*t852+t865; + coeffs_dxx->coeff_p2_p3 = t827*t869+t835*t852+t867; + coeffs_dxx->coeff_p3_p3 = t823*t852+t783+(t839+t872)*t851; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dxy.compute.c new file mode 100644 index 0000000..7cc7307 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dxy.compute.c @@ -0,0 +1,299 @@ +fp t1018; +fp t999; +fp t1017; +fp t998; +fp t1016; +fp t972; +fp t1015; +fp t936; +fp t997; +fp t1014; +fp t973; +fp t921; +fp t971; +fp t915; +fp t1013; +fp t916; +fp t919; +fp t896; +fp t1012; +fp t894; +fp t1011; +fp t1010; +fp t1009; +fp t1008; +fp t886; +fp t1007; +fp t890; +fp t970; +fp t1006; +fp t988; +fp t1005; +fp t1004; +fp t1003; +fp t952; +fp t1002; +fp t1001; +fp t1000; +fp t996; +fp t995; +fp t994; +fp t993; +fp t992; +fp t991; +fp t989; +fp t987; +fp t986; +fp t985; +fp t984; +fp t983; +fp t982; +fp t981; +fp t980; +fp t979; +fp t978; +fp t977; +fp t976; +fp t975; +fp t974; +fp t969; +fp t968; +fp t967; +fp t966; +fp t965; +fp t964; +fp t963; +fp t962; +fp t961; +fp t960; +fp t959; +fp t958; +fp t957; +fp t956; +fp t955; +fp t954; +fp t953; +fp t951; +fp t950; +fp t949; +fp t948; +fp t947; +fp t946; +fp t945; +fp t944; +fp t943; +fp t942; +fp t941; +fp t940; +fp t939; +fp t938; +fp t937; +fp t935; +fp t934; +fp t933; +fp t932; +fp t931; +fp t930; +fp t929; +fp t928; +fp t927; +fp t926; +fp t924; +fp t923; +fp t922; +fp t920; +fp t918; +fp t917; +fp t914; +fp t913; +fp t912; +fp t911; +fp t910; +fp t909; +fp t908; +fp t907; +fp t906; +fp t905; +fp t904; +fp t903; +fp t902; +fp t901; +fp t889; +fp t883; +fp t878; + t1018 = x*y; + t999 = x*x; + t1017 = y*t999; + t998 = y*y; + t1016 = t998*x; + t972 = RATIONAL(-5.0,6.0); + t1015 = x*t972; + t936 = RATIONAL(-1.0,18.0); + t997 = t999*t998; + t1014 = t936*t1018+RATIONAL(-1.0,16.0)*t997; + t973 = RATIONAL(7.0,6.0); + t921 = t973*t998; + t971 = RATIONAL(-7.0,6.0); + t915 = t971*t998; + t1013 = RATIONAL(-7.0,9.0)*t1018-t997; + t916 = t971*t999; + t919 = t973*t999; + t896 = RATIONAL(7.0,1.0)*t997; + t1012 = t896+RATIONAL(14.0,3.0)*t1018; + t894 = RATIONAL(-7.0,1.0)*t997; + t1011 = RATIONAL(-35.0,3.0)*t1018+t894; + t1010 = RATIONAL(5.0,9.0)*t1018+t997; + t1009 = RATIONAL(-16.0,1.0)*t997+RATIONAL(-140.0,9.0)*t1018; + t1008 = RATIONAL(14.0,9.0)*t1018+t997; + t886 = RATIONAL(7.0,16.0)*t997; + t1007 = RATIONAL(5.0,12.0)*t1018+t886; + t890 = RATIONAL(-7.0,16.0)*t997; + t970 = RATIONAL(-1.0,6.0); + t1006 = t970*t1018+t890; + t988 = RATIONAL(1.0,3.0); + t1005 = RATIONAL(-1.0,144.0)+t988*t1018+t886; + t1004 = y*t1015+t936+t890; + t1003 = t894+t936+RATIONAL(-10.0,3.0)*t1018; + t952 = RATIONAL(1.0,18.0); + t1002 = RATIONAL(-5.0,2.0)*t1018+t952+RATIONAL(-49.0,16.0)*t997; + t1001 = t952-t997+RATIONAL(-10.0,9.0)*t1018; + t1000 = t896+RATIONAL(25.0,3.0)*t1018+RATIONAL(-4.0,9.0); + t996 = RATIONAL(4.0,9.0); + t995 = RATIONAL(5.0,3.0); + t994 = RATIONAL(1.0,6.0); + t993 = RATIONAL(-2.0,3.0); + t992 = RATIONAL(-5.0,3.0); + t991 = RATIONAL(5.0,6.0); + t989 = RATIONAL(-7.0,4.0); + t987 = RATIONAL(4.0,1.0); + t986 = RATIONAL(5.0,8.0); + t985 = RATIONAL(-2.0,9.0); + t984 = RATIONAL(-4.0,1.0); + t983 = RATIONAL(7.0,4.0); + t982 = RATIONAL(1.0,4.0); + t981 = RATIONAL(-1.0,4.0); + t980 = RATIONAL(-5.0,8.0); + t979 = RATIONAL(-1.0,3.0); + t978 = RATIONAL(2.0,9.0); + t977 = RATIONAL(4.0,3.0); + t976 = RATIONAL(-1.0,9.0); + t975 = RATIONAL(-4.0,3.0); + t974 = RATIONAL(2.0,3.0); + t969 = RATIONAL(1.0,9.0); + t968 = RATIONAL(-8.0,3.0); + t967 = RATIONAL(8.0,3.0); + t966 = RATIONAL(1.0,72.0); + t965 = RATIONAL(-5.0,24.0); + t964 = RATIONAL(-7.0,12.0); + t963 = RATIONAL(-1.0,24.0); + t962 = RATIONAL(1.0,12.0); + t961 = RATIONAL(7.0,18.0); + t960 = RATIONAL(-35.0,6.0); + t959 = RATIONAL(1.0,24.0); + t958 = RATIONAL(-7.0,24.0); + t957 = RATIONAL(-7.0,18.0); + t956 = RATIONAL(7.0,24.0); + t955 = RATIONAL(-35.0,8.0); + t954 = RATIONAL(10.0,1.0); + t953 = RATIONAL(35.0,8.0); + t951 = RATIONAL(-10.0,1.0); + t950 = RATIONAL(-56.0,3.0); + t949 = RATIONAL(-28.0,9.0); + t948 = RATIONAL(49.0,6.0); + t947 = RATIONAL(28.0,9.0); + t946 = RATIONAL(56.0,3.0); + t945 = RATIONAL(-5.0,18.0); + t944 = RATIONAL(20.0,9.0); + t943 = RATIONAL(-49.0,6.0); + t942 = RATIONAL(-1.0,12.0); + t941 = RATIONAL(7.0,12.0); + t940 = RATIONAL(1.0,48.0); + t939 = RATIONAL(40.0,3.0); + t938 = RATIONAL(-40.0,3.0); + t937 = RATIONAL(-7.0,48.0); + t935 = RATIONAL(35.0,6.0); + t934 = RATIONAL(-1.0,48.0); + t933 = RATIONAL(7.0,48.0); + t932 = RATIONAL(-1.0,72.0); + t931 = RATIONAL(1.0,36.0); + t930 = RATIONAL(-20.0,9.0); + t929 = RATIONAL(-1.0,36.0); + t928 = RATIONAL(5.0,24.0); + t927 = RATIONAL(5.0,18.0); + t926 = RATIONAL(1.0,144.0); + t924 = t979*t999; + t923 = t970*t998; + t922 = t988*t999; + t920 = t970*t999; + t918 = t994*t998; + t917 = t968*t999; + t914 = t967*t999; + t913 = t968*t998; + t912 = t988*t998; + t911 = t994*t999; + t910 = t967*t998; + t909 = t979*t998; + t908 = t940*t998; + t907 = t940*t999; + t906 = t937*t999; + t905 = t934*t998; + t904 = t933*t999; + t903 = t934*t999; + t902 = t937*t998; + t901 = t933*t998; + t889 = RATIONAL(1.0,16.0)*t997; + t883 = RATIONAL(16.0,1.0)*t997; + t878 = RATIONAL(49.0,16.0)*t997; + coeffs_dxy->coeff_m2_m2 = t907+t908+t889+t926+(t942*t999+t929)*y+(t942* +t998+t969*y+t929)*x; + coeffs_dxy->coeff_m1_m2 = t906+t923+(t978+t941*t999)*y+(t928+t986*t998)*x ++t1004; + coeffs_dxy->coeff_0_m2 = t922+t975*t1017+(t957+t915)*x+t1008; + coeffs_dxy->coeff_p1_m2 = t918+t924+(t985+t977*t999)*y+(t991*t998+t927)*x ++t1001; + coeffs_dxy->coeff_p2_m2 = t905+t904+(t964*t999+t931)*y+(t981*t998+t942)*x ++t1005; + coeffs_dxy->coeff_p3_m2 = t903+t962*t1017+(t966+t959*t998)*x+t1014; + coeffs_dxy->coeff_m2_m1 = t902+t920+(t986*t999+t928)*y+(t978+t941*t998)*x ++t1004; + coeffs_dxy->coeff_m1_m1 = t996+t921+t919+t878+(t955*t999+t992)*y+(t992+ +RATIONAL(25.0,4.0)*y+t955*t998)*x; + coeffs_dxy->coeff_0_m1 = t954*t1017+t917+(t948*t998+t947)*x+t1011; + coeffs_dxy->coeff_p1_m1 = t915+t914+(t951*t999+t995)*y+(t930+t960*t998)*x ++t1000; + coeffs_dxy->coeff_p2_m1 = t916+t901+(t953*t999+t965)*y+(t974+t983*t998)*x ++t1002; + coeffs_dxy->coeff_p3_m1 = t980*t1017+t911+(t958*t998+t976)*x+t1007; + coeffs_dxy->coeff_m2_0 = t912+t975*t1016+(t957+t916)*y+t1008; + coeffs_dxy->coeff_m1_0 = t954*t1016+t913+(t948*t999+t947)*y+t1011; + coeffs_dxy->coeff_0_0 = t883+t950*t1017+(t950*t998+RATIONAL(196.0,9.0)*y) +*x; + coeffs_dxy->coeff_p1_0 = t939*t1016+t910+(t946*t999+t949)*y+t1009; + coeffs_dxy->coeff_p2_0 = t909+t984*t1016+(t961+t943*t999)*y+t1012; + coeffs_dxy->coeff_p3_0 = t974*t1016+y*t919+t1013; + coeffs_dxy->coeff_m2_p1 = t911+t909+(t991*t999+t927)*y+(t977*t998+t985)*x ++t1001; + coeffs_dxy->coeff_m1_p1 = t916+t910+(t960*t999+t930)*y+(t995+t951*t998)*x ++t1000; + coeffs_dxy->coeff_0_p1 = t939*t1017+t914+(t949+t946*t998)*x+t1009; + coeffs_dxy->coeff_p1_p1 = t883+t917+t996+t913+(t944+t938*t999)*y+(t938* +t998+RATIONAL(100.0,9.0)*y+t944)*x; + coeffs_dxy->coeff_p2_p1 = t912+t919+(t945+t935*t999)*y+(t987*t998+t993)*x ++t1003; + coeffs_dxy->coeff_p3_p1 = t920+t972*t1017+(t969+t993*t998)*x+t1010; + coeffs_dxy->coeff_m2_p2 = t903+t901+(t981*t999+t942)*y+(t964*t998+t931)*x ++t1005; + coeffs_dxy->coeff_m1_p2 = t904+t915+(t974+t983*t999)*y+(t965+t953*t998)*x ++t1002; + coeffs_dxy->coeff_0_p2 = t984*t1017+t924+(t943*t998+t961)*x+t1012; + coeffs_dxy->coeff_p1_p2 = t922+t921+(t993+t987*t999)*y+(t935*t998+t945)*x ++t1003; + coeffs_dxy->coeff_p2_p2 = t878+t902+t906+t926+(t962+t989*t999)*y+(t989* +t998+t962+y)*x; + coeffs_dxy->coeff_p3_p2 = t907+t982*t1017+(t956*t998+t932)*x+t1006; + coeffs_dxy->coeff_m2_p3 = t905+t962*t1016+(t959*t999+t966)*y+t1014; + coeffs_dxy->coeff_m1_p3 = t918+t980*t1016+(t976+t958*t999)*y+t1007; + coeffs_dxy->coeff_0_p3 = t974*t1017+x*t921+t1013; + coeffs_dxy->coeff_p1_p3 = t923+t998*t1015+(t969+t993*t999)*y+t1010; + coeffs_dxy->coeff_p2_p3 = t908+t982*t1016+(t956*t999+t932)*y+t1006; + coeffs_dxy->coeff_p3_p3 = t963*t1017+t889+(t963*t998+t931*y)*x; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dy.compute.c new file mode 100644 index 0000000..b2d27c0 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dy.compute.c @@ -0,0 +1,278 @@ +fp t749; +fp t748; +fp t779; +fp t778; +fp t777; +fp t747; +fp t776; +fp t775; +fp t741; +fp t759; +fp t684; +fp t708; +fp t774; +fp t727; +fp t773; +fp t738; +fp t762; +fp t681; +fp t723; +fp t772; +fp t713; +fp t771; +fp t770; +fp t769; +fp t664; +fp t709; +fp t768; +fp t767; +fp t728; +fp t696; +fp t739; +fp t676; +fp t766; +fp t764; +fp t711; +fp t765; +fp t678; +fp t763; +fp t719; +fp t666; +fp t761; +fp t760; +fp t758; +fp t757; +fp t722; +fp t698; +fp t756; +fp t742; +fp t755; +fp t745; +fp t707; +fp t754; +fp t714; +fp t699; +fp t753; +fp t743; +fp t752; +fp t751; +fp t746; +fp t744; +fp t740; +fp t737; +fp t736; +fp t735; +fp t734; +fp t733; +fp t732; +fp t731; +fp t730; +fp t729; +fp t726; +fp t725; +fp t724; +fp t720; +fp t718; +fp t717; +fp t716; +fp t715; +fp t706; +fp t705; +fp t704; +fp t703; +fp t702; +fp t701; +fp t697; +fp t695; +fp t694; +fp t693; +fp t692; +fp t691; +fp t690; +fp t689; +fp t686; +fp t683; +fp t679; +fp t677; +fp t674; +fp t673; +fp t672; +fp t671; +fp t663; +fp t662; +fp t660; +fp t659; +fp t656; +fp t654; +fp t650; + t749 = x*x; + t748 = t749*x; + t779 = x+t748; + t778 = y*t749; + t777 = y*t748; + t747 = y*y; + t776 = x*t747; + t775 = t747*t749; + t741 = RATIONAL(1.0,3.0); + t759 = t741*t747; + t684 = t748*t759; + t708 = RATIONAL(5.0,18.0); + t774 = t684+t708*t778; + t727 = RATIONAL(7.0,48.0); + t773 = t727*t747; + t738 = RATIONAL(-1.0,3.0); + t762 = t738*t747; + t681 = t748*t762; + t723 = RATIONAL(-7.0,18.0); + t772 = t723*t778+t681; + t713 = RATIONAL(-7.0,48.0); + t771 = t713*t747; + t770 = t748*t747; + t769 = t749*RATIONAL(-5.0,12.0); + t664 = t727*t770; + t709 = RATIONAL(5.0,24.0); + t768 = t709*t778+t664; + t767 = RATIONAL(-16.0,3.0)*t770+RATIONAL(-70.0,9.0)*t778; + t728 = RATIONAL(7.0,18.0); + t696 = t728*t748; + t739 = RATIONAL(7.0,3.0); + t676 = t739*t770; + t766 = t676+t739*t778; + t764 = t747*RATIONAL(-1.0,48.0); + t711 = RATIONAL(-1.0,36.0); + t765 = t748*t764+t711*t778; + t678 = RATIONAL(-7.0,3.0)*t770; + t763 = RATIONAL(-35.0,6.0)*t778+t678; + t719 = RATIONAL(-1.0,12.0); + t666 = t713*t770; + t761 = t719*t778+t666; + t760 = RATIONAL(1.0,48.0)*t747; + t758 = t684+RATIONAL(7.0,9.0)*t778; + t757 = t779*y; + t722 = RATIONAL(1.0,18.0); + t698 = t722*x; + t756 = RATIONAL(-49.0,48.0)*t770+RATIONAL(-5.0,4.0)*t778+t698; + t742 = RATIONAL(-4.0,9.0); + t755 = t676+RATIONAL(25.0,6.0)*t778+t742*x; + t745 = RATIONAL(1.0,6.0); + t707 = RATIONAL(-1.0,144.0); + t754 = t707*x+t745*t778+t664; + t714 = RATIONAL(-1.0,18.0); + t699 = t714*x; + t753 = t699+y*t769+t666; + t743 = RATIONAL(-5.0,3.0); + t752 = t743*t778+t678+t699; + t751 = RATIONAL(-5.0,9.0)*t778+t681+t698; + t746 = RATIONAL(4.0,9.0); + t744 = RATIONAL(-2.0,3.0); + t740 = RATIONAL(-2.0,9.0); + t737 = RATIONAL(2.0,9.0); + t736 = RATIONAL(-1.0,9.0); + t735 = RATIONAL(2.0,3.0); + t734 = RATIONAL(1.0,9.0); + t733 = RATIONAL(1.0,72.0); + t732 = RATIONAL(-5.0,24.0); + t731 = RATIONAL(-7.0,12.0); + t730 = RATIONAL(-1.0,24.0); + t729 = RATIONAL(1.0,12.0); + t726 = RATIONAL(-7.0,36.0); + t725 = RATIONAL(1.0,24.0); + t724 = RATIONAL(10.0,3.0); + t720 = RATIONAL(-5.0,18.0); + t718 = RATIONAL(7.0,12.0); + t717 = RATIONAL(7.0,36.0); + t716 = RATIONAL(-1.0,72.0); + t715 = RATIONAL(1.0,36.0); + t706 = RATIONAL(1.0,144.0); + t705 = t736*t748; + t704 = t734*t748; + t703 = t746*x; + t702 = RATIONAL(-8.0,9.0)*t748; + t701 = RATIONAL(8.0,9.0)*t748; + t697 = t714*t748; + t695 = t723*t748; + t694 = t722*t748; + t693 = t706*t748; + t692 = t706*x; + t691 = RATIONAL(-7.0,144.0)*t748; + t690 = RATIONAL(7.0,144.0)*t748; + t689 = t707*t748; + t686 = t745*t776; + t683 = RATIONAL(-1.0,6.0)*t776; + t679 = RATIONAL(7.0,6.0)*t776; + t677 = RATIONAL(-7.0,6.0)*t776; + t674 = x*t759; + t673 = RATIONAL(-8.0,3.0)*t776; + t672 = RATIONAL(8.0,3.0)*t776; + t671 = x*t762; + t663 = x*t764; + t662 = t748*t760; + t660 = x*t760; + t659 = x*t771; + t656 = x*t773; + t654 = RATIONAL(16.0,3.0)*t770; + t650 = RATIONAL(49.0,48.0)*t770; + coeffs_dy->coeff_m2_m2 = t693+t692+t660+t662+(t716+t730*t747)*t749+(t722* +t749+t779*t711)*y; + coeffs_dy->coeff_m1_m2 = t691+t683+(RATIONAL(5.0,48.0)+RATIONAL(5.0,16.0) +*t747)*t749+(t737*x+t717*t748)*y+t753; + coeffs_dy->coeff_0_m2 = t729+t726*t749+t704+(t731*t749+RATIONAL(1.0,4.0)) +*t747+(t742*t748+t738)*y+t758; + coeffs_dy->coeff_p1_m2 = t705+t686+(RATIONAL(5.0,12.0)*t747+RATIONAL(5.0, +36.0))*t749+(t740*x+t746*t748)*y+t751; + coeffs_dy->coeff_p2_m2 = t663+t690+(t730+RATIONAL(-1.0,8.0)*t747)*t749+( +t715*x+t726*t748)*y+t754; + coeffs_dy->coeff_p3_m2 = t715*t777+t689+(t706+t760)*t749+t765; + coeffs_dy->coeff_m2_m1 = t659+t697+(RATIONAL(7.0,24.0)*t747+t734)*t749+ +t709*t757+t753; + coeffs_dy->coeff_m1_m1 = t696+t650+t703+t679+(RATIONAL(-5.0,6.0)+RATIONAL +(-35.0,16.0)*t747)*t749+(RATIONAL(25.0,8.0)*t749+t743*x+RATIONAL(-35.0,24.0)* +t748)*y; + coeffs_dy->coeff_0_m1 = t702+RATIONAL(14.0,9.0)*t749+t744+(RATIONAL(-7.0, +4.0)+RATIONAL(49.0,12.0)*t749)*t747+(RATIONAL(5.0,2.0)+t724*t748)*y+t763; + coeffs_dy->coeff_p1_m1 = t677+t701+(RATIONAL(-35.0,12.0)*t747+RATIONAL( +-10.0,9.0))*t749+(RATIONAL(-10.0,3.0)*t748+RATIONAL(5.0,3.0)*x)*y+t755; + coeffs_dy->coeff_p2_m1 = t695+t656+(RATIONAL(7.0,8.0)*t747+t741)*t749+( +t732*x+RATIONAL(35.0,24.0)*t748)*y+t756; + coeffs_dy->coeff_p3_m1 = t732*t777+t694+(t714+t771)*t749+t768; + coeffs_dy->coeff_m2_0 = t744*t775+t674+t723*t757+t758; + coeffs_dy->coeff_m1_0 = RATIONAL(5.0,1.0)*t775+t673+(RATIONAL(28.0,9.0)*x ++RATIONAL(49.0,18.0)*t748)*y+t763; + coeffs_dy->coeff_0_0 = t654+(RATIONAL(4.0,1.0)+RATIONAL(-28.0,3.0)*t749)* +t747+(RATIONAL(-14.0,3.0)+RATIONAL(-56.0,9.0)*t748+RATIONAL(98.0,9.0)*t749)*y; + coeffs_dy->coeff_p1_0 = RATIONAL(20.0,3.0)*t775+t672+(RATIONAL(56.0,9.0)* +t748+RATIONAL(-28.0,9.0)*x)*y+t767; + coeffs_dy->coeff_p2_0 = t671+RATIONAL(-2.0,1.0)*t775+(RATIONAL(-49.0,18.0 +)*t748+t728*x)*y+t766; + coeffs_dy->coeff_p3_0 = y*t696+t749*t759+t772; + coeffs_dy->coeff_m2_p1 = t694+t671+(t735*t747+t736)*t749+t708*t757+t751; + coeffs_dy->coeff_m1_p1 = t672+t695+(RATIONAL(5.0,6.0)+RATIONAL(-5.0,1.0)* +t747)*t749+(RATIONAL(-35.0,18.0)*t748+RATIONAL(-20.0,9.0)*x)*y+t755; + coeffs_dy->coeff_0_p1 = RATIONAL(-14.0,9.0)*t749+t701+t735+(RATIONAL(-4.0 +,1.0)+RATIONAL(28.0,3.0)*t749)*t747+(RATIONAL(40.0,9.0)*t748+t724)*y+t767; + coeffs_dy->coeff_p1_p1 = t703+t702+t673+t654+(RATIONAL(-20.0,3.0)*t747+ +RATIONAL(10.0,9.0))*t749+(RATIONAL(-40.0,9.0)*t748+RATIONAL(20.0,9.0)*x+ +RATIONAL(50.0,9.0)*t749)*y; + coeffs_dy->coeff_p2_p1 = t696+t674+(t738+RATIONAL(2.0,1.0)*t747)*t749+( +RATIONAL(35.0,18.0)*t748+t720*x)*y+t752; + coeffs_dy->coeff_p3_p1 = t697+t720*t777+(t762+t722)*t749+t774; + coeffs_dy->coeff_m2_p2 = t689+t656+(t733+RATIONAL(-7.0,24.0)*t747)*t749+ +t719*t757+t754; + coeffs_dy->coeff_m1_p2 = t677+t690+(RATIONAL(-5.0,48.0)+RATIONAL(35.0, +16.0)*t747)*t749+(t718*t748+t735*x)*y+t756; + coeffs_dy->coeff_0_p2 = t719+RATIONAL(-4.0,3.0)*t777+t717*t749+t705-y+( +RATIONAL(7.0,4.0)+RATIONAL(-49.0,12.0)*t749)*t747+t766; + coeffs_dy->coeff_p1_p2 = t679+t704+(RATIONAL(35.0,12.0)*t747+RATIONAL( +-5.0,36.0))*t749+(t744*x+RATIONAL(4.0,3.0)*t748)*y+t752; + coeffs_dy->coeff_p2_p2 = t659+t650+t692+t691+(t725+RATIONAL(-7.0,8.0)* +t747)*t749+(t729*x+RATIONAL(1.0,2.0)*t749+t731*t748)*y; + coeffs_dy->coeff_p3_p2 = t693+t729*t777+(t773+t707)*t749+t761; + coeffs_dy->coeff_m2_p3 = t663+t725*t775+t733*t757+t765; + coeffs_dy->coeff_m1_p3 = t686+RATIONAL(-5.0,16.0)*t775+(RATIONAL(-7.0, +72.0)*t748+t736*x)*y+t768; + coeffs_dy->coeff_0_p3 = (t718*t749+RATIONAL(-1.0,4.0))*t747+(t737*t748+ +t745)*y+t772; + coeffs_dy->coeff_p1_p3 = t683+t747*t769+(t740*t748+t734*x)*y+t774; + coeffs_dy->coeff_p2_p3 = t660+RATIONAL(1.0,8.0)*t775+(RATIONAL(7.0,72.0)* +t748+t716*x)*y+t761; + coeffs_dy->coeff_p3_p3 = t662+t749*t764+(t716*t748+t733*t749)*y; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dyy.compute.c new file mode 100644 index 0000000..b75626c --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order3/coeffs-dyy.compute.c @@ -0,0 +1,219 @@ +fp t1116; +fp t1091; +fp t1115; +fp t1071; +fp t1114; +fp t1070; +fp t1113; +fp t1080; +fp t1112; +fp t1088; +fp t1111; +fp t1075; +fp t1110; +fp t1090; +fp t1109; +fp t1108; +fp t1032; +fp t1107; +fp t1050; +fp t1064; +fp t1028; +fp t1106; +fp t1066; +fp t1105; +fp t1104; +fp t1040; +fp t1058; +fp t1103; +fp t1059; +fp t1068; +fp t1034; +fp t1102; +fp t1101; +fp t1087; +fp t1100; +fp t1027; +fp t1099; +fp t1089; +fp t1098; +fp t1097; +fp t1042; +fp t1096; +fp t1060; +fp t1095; +fp t1082; +fp t1094; +fp t1069; +fp t1093; +fp t1086; +fp t1085; +fp t1084; +fp t1083; +fp t1081; +fp t1079; +fp t1078; +fp t1077; +fp t1076; +fp t1074; +fp t1073; +fp t1067; +fp t1065; +fp t1063; +fp t1062; +fp t1061; +fp t1041; +fp t1039; +fp t1038; +fp t1037; +fp t1036; +fp t1035; +fp t1033; +fp t1031; +fp t1030; +fp t1029; +fp t1025; +fp t1024; +fp t1023; +fp t1022; +fp t1020; + t1116 = x*y; + t1091 = x*x; + t1115 = y*t1091; + t1071 = RATIONAL(1.0,24.0); + t1114 = y*t1071; + t1070 = RATIONAL(-7.0,24.0); + t1113 = y*t1070; + t1080 = RATIONAL(2.0,3.0); + t1112 = y*t1080; + t1088 = RATIONAL(-2.0,3.0); + t1111 = y*t1088; + t1075 = RATIONAL(-1.0,24.0); + t1110 = t1075*y; + t1090 = t1091*x; + t1109 = y*t1090; + t1108 = x+t1090; + t1032 = t1070*t1109; + t1107 = RATIONAL(-5.0,12.0)*t1091+t1032; + t1050 = RATIONAL(-35.0,6.0)*t1091; + t1064 = RATIONAL(-14.0,3.0); + t1028 = t1064*t1109; + t1106 = t1050+t1028; + t1066 = RATIONAL(-1.0,12.0); + t1105 = t1066*t1091+t1032; + t1104 = RATIONAL(-70.0,9.0)*t1091+RATIONAL(-32.0,3.0)*t1109; + t1040 = t1080*t1109; + t1058 = RATIONAL(5.0,18.0); + t1103 = t1040+t1058*t1091; + t1059 = RATIONAL(5.0,24.0); + t1068 = RATIONAL(7.0,24.0); + t1034 = t1068*t1109; + t1102 = t1059*t1091+t1034; + t1101 = RATIONAL(-5.0,4.0)*t1091+RATIONAL(-49.0,24.0)*t1109; + t1087 = RATIONAL(-5.0,3.0); + t1100 = t1087*t1091+t1028; + t1027 = RATIONAL(14.0,3.0)*t1109; + t1099 = RATIONAL(25.0,6.0)*t1091+t1027; + t1089 = RATIONAL(1.0,6.0); + t1098 = t1089*t1091+t1034; + t1097 = t1040+RATIONAL(7.0,9.0)*t1091; + t1042 = t1088*t1109; + t1096 = t1042+RATIONAL(-5.0,9.0)*t1091; + t1060 = RATIONAL(-1.0,36.0); + t1095 = t1060*t1091+t1075*t1109; + t1082 = RATIONAL(7.0,3.0); + t1094 = t1027+t1082*t1091; + t1069 = RATIONAL(-7.0,18.0); + t1093 = t1042+t1069*t1091; + t1086 = RATIONAL(-2.0,9.0); + t1085 = RATIONAL(-1.0,3.0); + t1084 = RATIONAL(2.0,9.0); + t1083 = RATIONAL(4.0,3.0); + t1081 = RATIONAL(-4.0,3.0); + t1079 = RATIONAL(1.0,2.0); + t1078 = RATIONAL(1.0,72.0); + t1077 = RATIONAL(-5.0,24.0); + t1076 = RATIONAL(-7.0,12.0); + t1074 = RATIONAL(1.0,12.0); + t1073 = RATIONAL(7.0,18.0); + t1067 = RATIONAL(-5.0,18.0); + t1065 = RATIONAL(7.0,12.0); + t1063 = RATIONAL(10.0,3.0); + t1062 = RATIONAL(-1.0,72.0); + t1061 = RATIONAL(1.0,36.0); + t1041 = RATIONAL(1.0,3.0)*t1116; + t1039 = t1085*t1116; + t1038 = t1082*t1116; + t1037 = RATIONAL(-7.0,3.0)*t1116; + t1036 = x*t1112; + t1035 = x*t1111; + t1033 = x*t1110; + t1031 = t1071*t1109; + t1030 = x*t1114; + t1029 = x*t1113; + t1025 = t1068*t1116; + t1024 = RATIONAL(-16.0,3.0)*t1116; + t1023 = RATIONAL(32.0,3.0)*t1109; + t1022 = RATIONAL(16.0,3.0)*t1116; + t1020 = RATIONAL(49.0,24.0)*t1109; + coeffs_dyy->coeff_m2_m2 = t1030+t1031+(t1066*y+RATIONAL(1.0,18.0))*t1091+ +t1108*t1060; + coeffs_dyy->coeff_m1_m2 = RATIONAL(5.0,8.0)*t1115+RATIONAL(7.0,36.0)* +t1090+t1039+t1084*x+t1107; + coeffs_dyy->coeff_0_m2 = t1085+RATIONAL(-4.0,9.0)*t1090+(t1079+RATIONAL( +-7.0,6.0)*t1091)*y+t1097; + coeffs_dyy->coeff_p1_m2 = RATIONAL(5.0,6.0)*t1115+t1041+RATIONAL(4.0,9.0) +*t1090+t1086*x+t1096; + coeffs_dyy->coeff_p2_m2 = t1033+RATIONAL(-1.0,4.0)*t1115+t1061*x+RATIONAL +(-7.0,36.0)*t1090+t1098; + coeffs_dyy->coeff_p3_m2 = t1091*t1114+t1061*t1090+t1095; + coeffs_dyy->coeff_m2_m1 = t1065*t1115+t1029+t1108*t1059+t1107; + coeffs_dyy->coeff_m1_m1 = t1020+t1038+t1087*x+RATIONAL(-35.0,24.0)*t1090+ +(RATIONAL(-35.0,8.0)*y+RATIONAL(25.0,8.0))*t1091; + coeffs_dyy->coeff_0_m1 = RATIONAL(5.0,2.0)+t1063*t1090+(RATIONAL(49.0,6.0 +)*t1091+RATIONAL(-7.0,2.0))*y+t1106; + coeffs_dyy->coeff_p1_m1 = RATIONAL(5.0,3.0)*x+t1037+RATIONAL(-10.0,3.0)* +t1090+y*t1050+t1099; + coeffs_dyy->coeff_p2_m1 = RATIONAL(7.0,4.0)*t1115+t1077*x+RATIONAL(35.0, +24.0)*t1090+t1025+t1101; + coeffs_dyy->coeff_p3_m1 = t1091*t1113+t1077*t1090+t1102; + coeffs_dyy->coeff_m2_0 = t1036+t1081*t1115+t1108*t1069+t1097; + coeffs_dyy->coeff_m1_0 = RATIONAL(10.0,1.0)*t1115+RATIONAL(28.0,9.0)*x+ +t1024+RATIONAL(49.0,18.0)*t1090+t1106; + coeffs_dyy->coeff_0_0 = t1064+RATIONAL(-56.0,9.0)*t1090+RATIONAL(98.0,9.0 +)*t1091+t1023+(RATIONAL(-56.0,3.0)*t1091+RATIONAL(8.0,1.0))*y; + coeffs_dyy->coeff_p1_0 = t1022+RATIONAL(-28.0,9.0)*x+RATIONAL(56.0,9.0)* +t1090+RATIONAL(40.0,3.0)*t1115+t1104; + coeffs_dyy->coeff_p2_0 = t1035+RATIONAL(-49.0,18.0)*t1090+t1073*x+ +RATIONAL(-4.0,1.0)*t1115+t1094; + coeffs_dyy->coeff_p3_0 = t1073*t1090+t1091*t1112+t1093; + coeffs_dyy->coeff_m2_p1 = t1083*t1115+t1035+t1108*t1058+t1096; + coeffs_dyy->coeff_m1_p1 = RATIONAL(-35.0,18.0)*t1090+RATIONAL(-20.0,9.0)* +x+t1022+RATIONAL(-10.0,1.0)*t1115+t1099; + coeffs_dyy->coeff_0_p1 = t1063+RATIONAL(40.0,9.0)*t1090+(RATIONAL(56.0, +3.0)*t1091+RATIONAL(-8.0,1.0))*y+t1104; + coeffs_dyy->coeff_p1_p1 = t1023+t1024+RATIONAL(20.0,9.0)*x+RATIONAL(-40.0 +,9.0)*t1090+(RATIONAL(50.0,9.0)+RATIONAL(-40.0,3.0)*y)*t1091; + coeffs_dyy->coeff_p2_p1 = t1036+RATIONAL(35.0,18.0)*t1090+RATIONAL(4.0, +1.0)*t1115+t1067*x+t1100; + coeffs_dyy->coeff_p3_p1 = t1091*t1111+t1067*t1090+t1103; + coeffs_dyy->coeff_m2_p2 = t1076*t1115+t1025+t1108*t1066+t1098; + coeffs_dyy->coeff_m1_p2 = RATIONAL(35.0,8.0)*t1115+t1037+t1080*x+t1065* +t1090+t1101; + coeffs_dyy->coeff_0_p2 = t1081*t1090+RATIONAL(-1.0,1.0)+(RATIONAL(-49.0, +6.0)*t1091+RATIONAL(7.0,2.0))*y+t1094; + coeffs_dyy->coeff_p1_p2 = RATIONAL(35.0,6.0)*t1115+t1083*t1090+t1088*x+ +t1038+t1100; + coeffs_dyy->coeff_p2_p2 = t1020+t1029+t1076*t1090+t1074*x+(RATIONAL(-7.0, +4.0)*y+t1079)*t1091; + coeffs_dyy->coeff_p3_p2 = t1074*t1090+t1068*t1115+t1105; + coeffs_dyy->coeff_m2_p3 = t1074*t1115+t1033+t1108*t1078+t1095; + coeffs_dyy->coeff_m1_p3 = RATIONAL(-7.0,72.0)*t1090+t1041+RATIONAL(-1.0, +9.0)*x+RATIONAL(-5.0,8.0)*t1115+t1102; + coeffs_dyy->coeff_0_p3 = t1084*t1090+t1089+(RATIONAL(7.0,6.0)*t1091+ +RATIONAL(-1.0,2.0))*y+t1093; + coeffs_dyy->coeff_p1_p3 = RATIONAL(1.0,9.0)*x+t1086*t1090+RATIONAL(-5.0, +6.0)*t1115+t1039+t1103; + coeffs_dyy->coeff_p2_p3 = t1030+RATIONAL(1.0,4.0)*t1115+t1062*x+RATIONAL( +7.0,72.0)*t1090+t1105; + coeffs_dyy->coeff_p3_p3 = t1062*t1090+t1031+(t1110+t1078)*t1091; |