diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0')
10 files changed, 1672 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..d629ea0 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c @@ -0,0 +1,245 @@ +fp t82; +fp t84; +fp t114; +fp t83; +fp t104; +fp t26; +fp t106; +fp t34; +fp t154; +fp t85; +fp t129; +fp t105; +fp t25; +fp t29; +fp t153; +fp t81; +fp t130; +fp t107; +fp t36; +fp t152; +fp t108; +fp t38; +fp t46; +fp t151; +fp t109; +fp t45; +fp t150; +fp t27; +fp t131; +fp t110; +fp t35; +fp t149; +fp t42; +fp t48; +fp t148; +fp t147; +fp t77; +fp t111; +fp t51; +fp t78; +fp t128; +fp t146; +fp t123; +fp t60; +fp t127; +fp t145; +fp t115; +fp t64; +fp t31; +fp t119; +fp t144; +fp t52; +fp t143; +fp t124; +fp t142; +fp t70; +fp t33; +fp t141; +fp t63; +fp t44; +fp t140; +fp t139; +fp t103; +fp t41; +fp t68; +fp t138; +fp t47; +fp t137; +fp t136; +fp t50; +fp t135; +fp t126; +fp t65; +fp t134; +fp t133; +fp t24; +fp t132; +fp t125; +fp t121; +fp t80; +fp t79; +fp t76; +fp t32; +fp t122; +fp t120; +fp t118; +fp t117; +fp t116; +fp t113; +fp t112; +fp t102; +fp t101; +fp t100; +fp t99; +fp t98; +fp t59; +fp t62; +fp t97; +fp t96; +fp t95; +fp t94; +fp t93; +fp t92; +fp t91; +fp t90; +fp t89; +fp t88; +fp t87; +fp t86; + t82 = RATIONAL(1.0,8.0); + t84 = z*z; + t114 = t82*t84; + t83 = x*x; + t104 = t83*t114; + t26 = y*t104; + t106 = x*t114; + t34 = y*t106; + t154 = t26+t34; + t85 = y*y; + t129 = z*t85; + t105 = t82*t129; + t25 = t83*t105; + t29 = x*t105; + t153 = t25+t29; + t81 = RATIONAL(-1.0,8.0); + t130 = x*t81; + t107 = t84*t130; + t36 = y*t107; + t152 = t26+t36; + t108 = t81*t129; + t38 = x*t108; + t46 = t83*t108; + t151 = t38+t46; + t109 = y*t81*t83; + t45 = t84*t109; + t150 = t34+t45; + t27 = z*t109; + t131 = y*z; + t110 = t82*t131; + t35 = x*t110; + t149 = t27+t35; + t42 = t130*t131; + t48 = t83*t110; + t148 = t42+t48; + t147 = t36+t45; + t77 = RATIONAL(1.0,4.0); + t111 = t77*t83; + t51 = z*t111; + t78 = RATIONAL(-1.0,4.0); + t128 = z*t78; + t146 = (t128+t51)*y; + t123 = t78*t84; + t60 = x*t123; + t127 = y*t77; + t145 = x*t127+y*t60; + t115 = t78*t85; + t64 = x*t115; + t31 = t84*t64; + t119 = t77*t85; + t144 = x*t119+t31; + t52 = z*t119; + t143 = (t128+t52)*x; + t124 = t77*t84; + t142 = x*t124+t31; + t70 = t84*t119; + t33 = x*t70; + t141 = t60+t33; + t63 = t85*t111; + t44 = z*t63; + t140 = z*t115+t44; + t139 = (t64+t77*x)*z; + t103 = t83*t115; + t41 = t84*t103; + t68 = t84*t111; + t138 = t41+t68; + t47 = t85*t104; + t137 = t85*t106+t47; + t136 = t85*t107+t47; + t50 = y*t124; + t135 = (t78*y+t50)*x; + t126 = t78*t83; + t65 = y*t126; + t134 = (t127+t65)*z; + t133 = t64+t33; + t24 = z*t103; + t132 = t24+t52; + t125 = t85+t84; + t121 = t83*t84; + t80 = RATIONAL(-1.0,2.0); + t79 = RATIONAL(1.0,2.0); + t76 = t79*t85; + t32 = t76*t121; + t122 = t80*t121+t32; + t120 = t85*t84; + t118 = y*t68+t41; + t117 = t84*t65+t41; + t116 = t84+t83; + t113 = t83+t85; + t112 = t80*t85; + t102 = 1.0+t120; + t101 = t121+1.0; + t100 = t83*t85+1.0; + t99 = t29+t46+t136; + t98 = z*t126+t44+t138; + t59 = t84*t112; + t62 = t83*t112; + t97 = t59+t32+t76+t62; + t96 = t79*t83+t62+t122; + t95 = t25+t38+t136; + t94 = y*t123+t70+t118; + t93 = t50+t70+t117; + t92 = t59+t79*t84+t122; + t91 = t63+t65+t118; + t90 = t63+y*t111+t117; + t89 = t35+t48+t137; + t88 = t27+t42+t137; + t87 = t51+t24+t138; + t86 = -t120+t125; + coeffs_I->coeff_m1_m1_m1 = t99+t148+t150; + coeffs_I->coeff_0_m1_m1 = t94+t134+t140; + coeffs_I->coeff_p1_m1_m1 = t89+t147+t151; + coeffs_I->coeff_m1_0_m1 = t98+t139+t141; + coeffs_I->coeff_0_0_m1 = (t100*t80+t113*t79)*z+t92; + coeffs_I->coeff_p1_0_m1 = t98+t142+t143; + coeffs_I->coeff_m1_p1_m1 = t99+t149+t152; + coeffs_I->coeff_0_p1_m1 = t93+t140+t146; + coeffs_I->coeff_p1_p1_m1 = t88+t151+t154; + coeffs_I->coeff_m1_m1_0 = t91+t133+t145; + coeffs_I->coeff_0_m1_0 = (t101*t80+t116*t79)*y+t97; + coeffs_I->coeff_p1_m1_0 = t91+t135+t144; + coeffs_I->coeff_m1_0_0 = (t102*t80+t125*t79)*x+t96; + coeffs_I->coeff_0_0_0 = RATIONAL(1.0,1.0)+(-1.0+t86)*t83-t86; + coeffs_I->coeff_p1_0_0 = (t102*t79+t125*t80)*x+t96; + coeffs_I->coeff_m1_p1_0 = t90+t133+t135; + coeffs_I->coeff_0_p1_0 = (t101*t79+t116*t80)*y+t97; + coeffs_I->coeff_p1_p1_0 = t90+t144+t145; + coeffs_I->coeff_m1_m1_p1 = t95+t149+t150; + coeffs_I->coeff_0_m1_p1 = t94+t132+t146; + coeffs_I->coeff_p1_m1_p1 = t88+t147+t153; + coeffs_I->coeff_m1_0_p1 = t87+t141+t143; + coeffs_I->coeff_0_0_p1 = (t100*t79+t113*t80)*z+t92; + coeffs_I->coeff_p1_0_p1 = t87+t139+t142; + coeffs_I->coeff_m1_p1_p1 = t95+t148+t152; + coeffs_I->coeff_0_p1_p1 = t93+t132+t134; + coeffs_I->coeff_p1_p1_p1 = t89+t153+t154; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..52dc6d6 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c @@ -0,0 +1,186 @@ +fp t202; +fp t206; +fp t207; +fp t223; +fp t205; +fp t160; +fp t201; +fp t231; +fp t174; +fp t242; +fp t218; +fp t179; +fp t241; +fp t198; +fp t194; +fp t178; +fp t197; +fp t240; +fp t228; +fp t180; +fp t156; +fp t239; +fp t165; +fp t238; +fp t196; +fp t175; +fp t237; +fp t172; +fp t161; +fp t236; +fp t193; +fp t188; +fp t235; +fp t166; +fp t182; +fp t234; +fp t189; +fp t177; +fp t233; +fp t232; +fp t230; +fp t229; +fp t227; +fp t226; +fp t225; +fp t224; +fp t222; +fp t203; +fp t221; +fp t220; +fp t204; +fp t219; +fp t171; +fp t187; +fp t217; +fp t216; +fp t167; +fp t168; +fp t158; +fp t162; +fp t215; +fp t214; +fp t213; +fp t212; +fp t211; +fp t181; +fp t170; +fp t159; +fp t210; +fp t163; +fp t184; +fp t186; +fp t209; +fp t173; +fp t185; +fp t155; +fp t208; +fp t183; +fp t176; + t202 = RATIONAL(-1.0,2.0); + t206 = z*z; + t207 = y*y; + t223 = t207*t206; + t205 = x*t223; + t160 = t202*t205; + t201 = RATIONAL(1.0,2.0); + t231 = x*t201; + t174 = t206*t231; + t242 = t160+t174; + t218 = t201*t207; + t179 = x*t218; + t241 = t179+t160; + t198 = RATIONAL(-1.0,4.0); + t194 = t198*t206; + t178 = y*t194; + t197 = RATIONAL(1.0,4.0); + t240 = t178+t197*y; + t228 = x*t202; + t180 = y*t228; + t156 = y*t174; + t239 = t180+t156; + t165 = t206*t180; + t238 = y*t231+t165; + t196 = t197*t206; + t175 = y*t196; + t237 = t198*y+t175; + t172 = z*t228; + t161 = z*t179; + t236 = t172+t161; + t193 = t198*t207; + t188 = t206*t193; + t235 = t188+t196; + t166 = t207*t172; + t182 = z*t231; + t234 = t166+t182; + t189 = t197*t207; + t177 = t206*t189; + t233 = t177+t194; + t232 = y*z; + t230 = t207*z; + t229 = t202-x; + t227 = y*t206; + t226 = t201-x; + t225 = x+t205; + t224 = x*t232; + t222 = z*t180+t160; + t203 = RATIONAL(-1.0,8.0); + t221 = t203*t207; + t220 = t160+y*t182; + t204 = RATIONAL(1.0,8.0); + t219 = t204*t206; + t171 = z*t193; + t187 = z*t189; + t217 = -y+t227; + t216 = -t230+z; + t167 = t198*t224; + t168 = x*t175; + t158 = x*t171; + t162 = x*t177; + t215 = t167+t168+t158+t162; + t214 = t171+t197*z+t242; + t213 = t187+t198*z+t242; + t212 = t188+t189+t241; + t211 = t177+t193+t241; + t181 = y*t219; + t170 = t204*t230; + t159 = t197*t224; + t210 = t181+t170+t159+t162; + t163 = x*t178; + t184 = t204*t232; + t186 = z*t221; + t209 = t163+t184+t186+t162; + t173 = t203*t232; + t185 = t203*t227; + t155 = x*t187; + t208 = t173+t185+t155+t162; + t183 = t207*t219; + t176 = t206*t221; + coeffs_dx->coeff_m1_m1_m1 = t163+t176+t158+t173+t210; + coeffs_dx->coeff_0_m1_m1 = t156+t161+t222; + coeffs_dx->coeff_p1_m1_m1 = t158+t159+t185+t183+t209; + coeffs_dx->coeff_m1_0_m1 = t214+t233+t236; + coeffs_dx->coeff_0_0_m1 = t205+(-t206+t216)*x; + coeffs_dx->coeff_p1_0_m1 = t213+t235+t236; + coeffs_dx->coeff_m1_p1_m1 = t184+t176+t185+t170+t215; + coeffs_dx->coeff_0_p1_m1 = t165+t161+t220; + coeffs_dx->coeff_p1_p1_m1 = t173+t183+t181+t186+t215; + coeffs_dx->coeff_m1_m1_0 = t211+t239+t240; + coeffs_dx->coeff_0_m1_0 = t205+(-t207-t217)*x; + coeffs_dx->coeff_p1_m1_0 = t212+t237+t239; + coeffs_dx->coeff_m1_0_0 = t202+t226*t207+(t202*t207+t226)*t206+t225; + coeffs_dx->coeff_0_0_0 = ((1.0+t223)*RATIONAL(-2.0,1.0)+(t207+t206)* +RATIONAL(2.0,1.0))*x; + coeffs_dx->coeff_p1_0_0 = t201+t229*t207+(t218+t229)*t206+t225; + coeffs_dx->coeff_m1_p1_0 = t211+t237+t238; + coeffs_dx->coeff_0_p1_0 = t205+(-t207+t217)*x; + coeffs_dx->coeff_p1_p1_0 = t212+t238+t240; + coeffs_dx->coeff_m1_m1_p1 = t167+t181+t155+t176+t209; + coeffs_dx->coeff_0_m1_p1 = t156+t166+t220; + coeffs_dx->coeff_p1_m1_p1 = t163+t183+t167+t170+t208; + coeffs_dx->coeff_m1_0_p1 = t213+t233+t234; + coeffs_dx->coeff_0_0_p1 = t205+(-t206-t216)*x; + coeffs_dx->coeff_p1_0_p1 = t214+t234+t235; + coeffs_dx->coeff_m1_p1_p1 = t168+t176+t159+t186+t208; + coeffs_dx->coeff_0_p1_p1 = t165+t166+t222; + coeffs_dx->coeff_p1_p1_p1 = t168+t155+t183+t184+t210; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..e8ab5c4 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,120 @@ +fp t465; +fp t485; +fp t464; +fp t484; +fp t467; +fp t483; +fp t468; +fp t482; +fp t460; +fp t481; +fp t461; +fp t480; +fp t466; +fp t479; +fp t478; +fp t477; +fp t473; +fp t448; +fp t476; +fp t459; +fp t458; +fp t452; +fp t453; +fp t475; +fp t474; +fp t444; +fp t472; +fp t471; +fp t470; +fp t469; +fp t455; +fp t454; +fp t451; +fp t450; +fp t447; +fp t446; +fp t445; +fp t443; +fp t442; +fp t441; +fp t440; +fp t439; +fp t438; +fp t437; +fp t436; +fp t435; + t465 = RATIONAL(-1.0,2.0); + t485 = t465*y; + t464 = RATIONAL(1.0,2.0); + t484 = t464*y; + t467 = z*z; + t483 = y*t467; + t468 = y*y; + t482 = t468*z; + t460 = RATIONAL(1.0,4.0); + t481 = y*t460; + t461 = RATIONAL(-1.0,4.0); + t480 = z*t461; + t466 = t467*t468; + t479 = -t468+t466; + t478 = t466-t467; + t477 = t468+t467; + t473 = t460*t468; + t448 = t467*t473; + t476 = z*t473+t448; + t459 = t464*t467; + t458 = t464*t468; + t452 = z*t458; + t453 = t465*t466; + t475 = t452+t453; + t474 = t448+t468*t480; + t444 = t465*t483; + t472 = t444+t453; + t471 = t453+t465*t482; + t470 = t482-z; + t469 = t483-y; + t455 = t461*t483; + t454 = y*t459; + t451 = z*t485; + t450 = y*t480; + t447 = t467*t481; + t446 = z*t484; + t445 = z*t481; + t443 = t466+RATIONAL(1.0,1.0)-t477; + t442 = t464*z+t459+t471; + t441 = t458+t484+t472; + t440 = t465*z+t459+t475; + t439 = t485+t454+t453+t458; + t438 = t447+t445+t476; + t437 = t455+t450+t476; + t436 = t455+t445+t474; + t435 = t450+t447+t474; + coeffs_dxx->coeff_m1_m1_m1 = t436; + coeffs_dxx->coeff_0_m1_m1 = t451+t454+t475; + coeffs_dxx->coeff_p1_m1_m1 = t436; + coeffs_dxx->coeff_m1_0_m1 = t440; + coeffs_dxx->coeff_0_0_m1 = -t470+t478; + coeffs_dxx->coeff_p1_0_m1 = t440; + coeffs_dxx->coeff_m1_p1_m1 = t435; + coeffs_dxx->coeff_0_p1_m1 = t452+t446+t472; + coeffs_dxx->coeff_p1_p1_m1 = t435; + coeffs_dxx->coeff_m1_m1_0 = t439; + coeffs_dxx->coeff_0_m1_0 = -t469+t479; + coeffs_dxx->coeff_p1_m1_0 = t439; + coeffs_dxx->coeff_m1_0_0 = t443; + coeffs_dxx->coeff_0_0_0 = (t466+1.0)*RATIONAL(-2.0,1.0)+t477*RATIONAL(2.0 +,1.0); + coeffs_dxx->coeff_p1_0_0 = t443; + coeffs_dxx->coeff_m1_p1_0 = t441; + coeffs_dxx->coeff_0_p1_0 = t469+t479; + coeffs_dxx->coeff_p1_p1_0 = t441; + coeffs_dxx->coeff_m1_m1_p1 = t437; + coeffs_dxx->coeff_0_m1_p1 = t446+t454+t471; + coeffs_dxx->coeff_p1_m1_p1 = t437; + coeffs_dxx->coeff_m1_0_p1 = t442; + coeffs_dxx->coeff_0_0_p1 = t470+t478; + coeffs_dxx->coeff_p1_0_p1 = t442; + coeffs_dxx->coeff_m1_p1_p1 = t438; + coeffs_dxx->coeff_0_p1_p1 = t444+t451+t471; + coeffs_dxx->coeff_p1_p1_p1 = t438; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxy.compute.c new file mode 100644 index 0000000..d5b8bd9 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxy.compute.c @@ -0,0 +1,170 @@ +fp t521; +fp t527; +fp t541; +fp t496; +fp t522; +fp t513; +fp t560; +fp t518; +fp t515; +fp t517; +fp t559; +fp t542; +fp t498; +fp t511; +fp t558; +fp t524; +fp t516; +fp t523; +fp t510; +fp t557; +fp t545; +fp t500; +fp t504; +fp t556; +fp t512; +fp t499; +fp t555; +fp t507; +fp t508; +fp t554; +fp t492; +fp t506; +fp t553; +fp t546; +fp t503; +fp t552; +fp t493; +fp t505; +fp t551; +fp t495; +fp t489; +fp t550; +fp t509; +fp t491; +fp t549; +fp t548; +fp t547; +fp t526; +fp t514; +fp t544; +fp t525; +fp t519; +fp t540; +fp t486; +fp t494; +fp t543; +fp t539; +fp t538; +fp t537; +fp t536; +fp t535; +fp t534; +fp t533; +fp t532; +fp t531; +fp t530; +fp t529; +fp t490; +fp t528; +fp t501; +fp t497; + t521 = RATIONAL(1.0,2.0); + t527 = z*z; + t541 = t521*t527; + t496 = y*t541; + t522 = RATIONAL(-1.0,2.0); + t513 = t522*y; + t560 = t496+t513; + t518 = RATIONAL(-1.0,4.0); + t515 = t518*t527; + t517 = RATIONAL(1.0,4.0); + t559 = t515+t517; + t542 = t522*t527; + t498 = y*t542; + t511 = t521*y; + t558 = t498+t511; + t524 = RATIONAL(-1.0,8.0); + t516 = t524*z; + t523 = RATIONAL(1.0,8.0); + t510 = t523*t527; + t557 = t516+t510; + t545 = z*t518; + t500 = y*t545; + t504 = y*t515; + t556 = t500+t504; + t512 = t517*t527; + t499 = y*t512; + t555 = t499+t500; + t507 = t523*z; + t508 = t524*t527; + t554 = t507+t508; + t492 = x*t545; + t506 = x*t515; + t553 = t492+t506; + t546 = z*t517; + t503 = y*t546; + t552 = t499+t503; + t493 = x*t512; + t505 = x*t546; + t551 = t493+t505; + t495 = z*t513; + t489 = x*t496; + t550 = x*t495+t489; + t509 = t521*x; + t491 = z*t509; + t549 = y*t491+t489; + t548 = t504+t503; + t547 = t518+t512; + t526 = y*x; + t514 = t522*x; + t544 = y*t527; + t525 = z*t526; + t519 = RATIONAL(2.0,1.0); + t540 = x*t544; + t486 = t519*t540; + t494 = RATIONAL(-2.0,1.0)*t526; + t543 = t494+t486; + t539 = (1.0-t527)*x; + t538 = y-t544; + t537 = -t525-t540; + t536 = t525-t540; + t535 = x*t541-t540; + t534 = t495-t540; + t533 = t516+t508+t549; + t532 = t505+t506+t550; + t531 = t507+t510+t549; + t530 = t492+t493+t550; + t529 = t514+t526+t535; + t490 = x*t542; + t528 = t490+t526-t540+t509; + t501 = z*t514; + t497 = z*t511; + coeffs_dxy->coeff_m1_m1_m1 = t532+t548+t557; + coeffs_dxy->coeff_0_m1_m1 = t501+t525+t535; + coeffs_dxy->coeff_p1_m1_m1 = t532+t554+t555; + coeffs_dxy->coeff_m1_0_m1 = t525+t496+t534; + coeffs_dxy->coeff_0_0_m1 = t486+z*t494; + coeffs_dxy->coeff_p1_0_m1 = t497+t498+t536; + coeffs_dxy->coeff_m1_p1_m1 = t530+t548+t554; + coeffs_dxy->coeff_0_p1_m1 = t491+t490+t536; + coeffs_dxy->coeff_p1_p1_m1 = t530+t555+t557; + coeffs_dxy->coeff_m1_m1_0 = t529+t559+t560; + coeffs_dxy->coeff_0_m1_0 = t539+t543; + coeffs_dxy->coeff_p1_m1_0 = t529+t547+t558; + coeffs_dxy->coeff_m1_0_0 = t538+t543; + coeffs_dxy->coeff_0_0_0 = (RATIONAL(-4.0,1.0)*t527+RATIONAL(4.0,1.0))* +t526; + coeffs_dxy->coeff_p1_0_0 = -t538+t543; + coeffs_dxy->coeff_m1_p1_0 = t528+t547+t560; + coeffs_dxy->coeff_0_p1_0 = -t539+t543; + coeffs_dxy->coeff_p1_p1_0 = t528+t558+t559; + coeffs_dxy->coeff_m1_m1_p1 = t531+t553+t556; + coeffs_dxy->coeff_0_m1_p1 = t491-t525+t535; + coeffs_dxy->coeff_p1_m1_p1 = t533+t552+t553; + coeffs_dxy->coeff_m1_0_p1 = t497+t496+t537; + coeffs_dxy->coeff_0_0_p1 = t486+t519*t525; + coeffs_dxy->coeff_p1_0_p1 = -t525+t498+t534; + coeffs_dxy->coeff_m1_p1_p1 = t533+t551+t556; + coeffs_dxy->coeff_0_p1_p1 = t490+t501+t537; + coeffs_dxy->coeff_p1_p1_p1 = t531+t551+t552; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxz.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxz.compute.c new file mode 100644 index 0000000..66fe899 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxz.compute.c @@ -0,0 +1,152 @@ +fp t596; +fp t602; +fp t597; +fp t618; +fp t565; +fp t629; +fp t593; +fp t587; +fp t592; +fp t628; +fp t617; +fp t576; +fp t627; +fp t583; +fp t626; +fp t624; +fp t616; +fp t601; +fp t625; +fp t623; +fp t622; +fp t621; +fp t620; +fp t571; +fp t615; +fp t562; +fp t619; +fp t581; +fp t567; +fp t614; +fp t613; +fp t600; +fp t612; +fp t611; +fp t610; +fp t609; +fp t598; +fp t589; +fp t590; +fp t564; +fp t561; +fp t608; +fp t599; +fp t591; +fp t573; +fp t569; +fp t607; +fp t588; +fp t577; +fp t566; +fp t606; +fp t578; +fp t580; +fp t563; +fp t605; +fp t572; +fp t604; +fp t570; +fp t603; +fp t579; +fp t575; +fp t574; + t596 = RATIONAL(1.0,2.0); + t602 = y*y; + t597 = RATIONAL(-1.0,2.0); + t618 = t602*t597; + t565 = x*t618; + t629 = t565+t596*x; + t593 = RATIONAL(-1.0,4.0); + t587 = t593*t602; + t592 = RATIONAL(1.0,4.0); + t628 = t587+t592; + t617 = t602*t596; + t576 = x*t617; + t627 = t576+t597*x; + t583 = t592*t602; + t626 = t593+t583; + t624 = t602*x; + t616 = z*t624; + t601 = x*z; + t625 = -t616+t601; + t623 = y*t593; + t622 = y*t592; + t621 = y*t597; + t620 = y*t596; + t571 = RATIONAL(-2.0,1.0)*t601; + t615 = RATIONAL(2.0,1.0)*t601; + t562 = t602*t615; + t619 = t571+t562; + t581 = z*t620; + t567 = z*t621; + t614 = (1.0-t602)*z; + t613 = x-t624; + t600 = y*t601; + t612 = -t600-t616; + t611 = t600-t616; + t610 = -t616+x*t620; + t609 = -t616+t567; + t598 = RATIONAL(1.0,8.0); + t589 = t598*y; + t590 = t598*t602; + t564 = z*t576; + t561 = x*t581; + t608 = t589+t590+t564+t561; + t599 = RATIONAL(-1.0,8.0); + t591 = t599*y; + t573 = z*t583; + t569 = x*t623; + t607 = t591+t573+t569+t564; + t588 = t599*t602; + t577 = z*t587; + t566 = x*t583; + t606 = t588+t577+t564+t566; + t578 = x*t622; + t580 = x*t587; + t563 = x*t567; + t605 = t564+t578+t580+t563; + t572 = z*t617; + t604 = t572+t597*z+t625; + t570 = z*t618; + t603 = t570+t596*z+t625; + t579 = z*t622; + t575 = x*t621; + t574 = z*t623; + coeffs_dxz->coeff_m1_m1_m1 = t579+t591+t590+t577+t605; + coeffs_dxz->coeff_0_m1_m1 = t575+t576+t611; + coeffs_dxz->coeff_p1_m1_m1 = t589+t574+t588+t573+t605; + coeffs_dxz->coeff_m1_0_m1 = t604+t627+t628; + coeffs_dxz->coeff_0_0_m1 = t613+t619; + coeffs_dxz->coeff_p1_0_m1 = t603+t626+t627; + coeffs_dxz->coeff_m1_p1_m1 = t580+t574+t569+t577+t608; + coeffs_dxz->coeff_0_p1_m1 = -t600+t576+t610; + coeffs_dxz->coeff_p1_p1_m1 = t588+t580+t561+t579+t607; + coeffs_dxz->coeff_m1_m1_0 = t572+t600+t609; + coeffs_dxz->coeff_0_m1_0 = t562+y*t571; + coeffs_dxz->coeff_p1_m1_0 = t581+t570+t611; + coeffs_dxz->coeff_m1_0_0 = t614+t619; + coeffs_dxz->coeff_0_0_0 = (RATIONAL(-4.0,1.0)*t602+RATIONAL(4.0,1.0))* +t601; + coeffs_dxz->coeff_p1_0_0 = -t614+t619; + coeffs_dxz->coeff_m1_p1_0 = t581+t572+t612; + coeffs_dxz->coeff_0_p1_0 = t562+y*t615; + coeffs_dxz->coeff_p1_p1_0 = t570-t600+t609; + coeffs_dxz->coeff_m1_m1_p1 = t579+t569+t563+t589+t606; + coeffs_dxz->coeff_0_m1_p1 = t565+t600+t610; + coeffs_dxz->coeff_p1_m1_p1 = t566+t574+t563+t590+t607; + coeffs_dxz->coeff_m1_0_p1 = t604+t626+t629; + coeffs_dxz->coeff_0_0_p1 = -t613+t619; + coeffs_dxz->coeff_p1_0_p1 = t603+t628+t629; + coeffs_dxz->coeff_m1_p1_p1 = t574+t591+t561+t578+t606; + coeffs_dxz->coeff_0_p1_p1 = t575+t565+t612; + coeffs_dxz->coeff_p1_p1_p1 = t579+t566+t578+t573+t608; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dy.compute.c new file mode 100644 index 0000000..9bc1495 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dy.compute.c @@ -0,0 +1,203 @@ +fp t285; +fp t295; +fp t286; +fp t282; +fp t262; +fp t339; +fp t294; +fp t284; +fp t307; +fp t265; +fp t338; +fp t289; +fp t313; +fp t263; +fp t244; +fp t290; +fp t319; +fp t273; +fp t337; +fp t272; +fp t254; +fp t336; +fp t270; +fp t278; +fp t335; +fp t280; +fp t257; +fp t334; +fp t333; +fp t332; +fp t291; +fp t306; +fp t246; +fp t331; +fp t292; +fp t316; +fp t268; +fp t320; +fp t311; +fp t251; +fp t330; +fp t259; +fp t256; +fp t329; +fp t258; +fp t243; +fp t328; +fp t252; +fp t317; +fp t269; +fp t327; +fp t326; +fp t314; +fp t260; +fp t264; +fp t255; +fp t325; +fp t324; +fp t293; +fp t253; +fp t310; +fp t274; +fp t323; +fp t271; +fp t245; +fp t266; +fp t322; +fp t321; +fp t318; +fp t315; +fp t312; +fp t261; +fp t309; +fp t308; +fp t305; +fp t304; +fp t303; +fp t302; +fp t301; +fp t300; +fp t299; +fp t298; +fp t297; +fp t296; +fp t288; +fp t287; +fp t249; +fp t248; + t285 = RATIONAL(1.0,4.0); + t295 = z*z; + t286 = RATIONAL(-1.0,4.0); + t282 = t286*t295; + t262 = x*t282; + t339 = t262+t285*x; + t294 = x*x; + t284 = t285*t294; + t307 = t294*t295; + t265 = t286*t307; + t338 = t284+t265; + t289 = RATIONAL(1.0,2.0); + t313 = y*t289; + t263 = t294*t313; + t244 = z*t263; + t290 = RATIONAL(-1.0,2.0); + t319 = y*t290; + t273 = z*t319; + t337 = t244+t273; + t272 = z*t313; + t254 = t294*t273; + t336 = t272+t254; + t270 = t285*t307; + t278 = t286*t294; + t335 = t270+t278; + t280 = t285*t295; + t257 = x*t280; + t334 = t286*x+t257; + t333 = t265+t280; + t332 = t282+t270; + t291 = RATIONAL(-1.0,8.0); + t306 = t294*t291; + t246 = y*t270; + t331 = t246+t295*t306; + t292 = RATIONAL(1.0,8.0); + t316 = x*t292; + t268 = z*t316; + t320 = x*z; + t311 = y*t320; + t251 = t286*t311; + t330 = t268+t251; + t259 = t295*t316; + t256 = y*t262; + t329 = t259+t256; + t258 = t291*t320; + t243 = t285*t311; + t328 = t258+t243; + t252 = y*t257; + t317 = x*t295; + t269 = t291*t317; + t327 = t252+t269; + t326 = t243+t268; + t314 = t294*z; + t260 = t292*t314; + t264 = z*t284; + t255 = y*t264; + t325 = t260+t255; + t324 = t258+t251; + t293 = y*t307; + t253 = t290*t293; + t310 = t289*t295; + t274 = y*t310; + t323 = t253+t274; + t271 = z*t278; + t245 = y*t271; + t266 = z*t306; + t322 = t245+t266; + t321 = t246+t292*t307; + t318 = -y+t290; + t315 = t293+y; + t312 = t289-y; + t261 = x*t319; + t309 = t295*t261+t253; + t308 = x*t274+t253; + t305 = -x+t317; + t304 = -z+t314; + t303 = t252+t259+t321; + t302 = t269+t256+t321; + t301 = t260+t245+t331; + t300 = t285*z+t271+t323; + t299 = t266+t255+t331; + t298 = t263+x*t313+t309; + t297 = t264+t286*z+t323; + t296 = t261+t263+t308; + t288 = RATIONAL(-2.0,1.0); + t287 = RATIONAL(2.0,1.0); + t249 = z*t261; + t248 = x*t272; + coeffs_dy->coeff_m1_m1_m1 = t301+t328+t329; + coeffs_dy->coeff_0_m1_m1 = t300+t332+t337; + coeffs_dy->coeff_p1_m1_m1 = t301+t327+t330; + coeffs_dy->coeff_m1_0_m1 = t249+t244+t308; + coeffs_dy->coeff_0_0_m1 = t293+(-t295-t304)*y; + coeffs_dy->coeff_p1_0_m1 = t248+t244+t309; + coeffs_dy->coeff_m1_p1_m1 = t302+t322+t326; + coeffs_dy->coeff_0_p1_m1 = t297+t333+t337; + coeffs_dy->coeff_p1_p1_m1 = t303+t322+t324; + coeffs_dy->coeff_m1_m1_0 = t296+t335+t339; + coeffs_dy->coeff_0_m1_0 = t290+t312*t295+(t290*t295+t312)*t294+t315; + coeffs_dy->coeff_p1_m1_0 = t298+t334+t335; + coeffs_dy->coeff_m1_0_0 = t293+(-t294-t305)*y; + coeffs_dy->coeff_0_0_0 = (t287*t295+t288+(t288*t295+t287)*t294)*y; + coeffs_dy->coeff_p1_0_0 = t293+(-t294+t305)*y; + coeffs_dy->coeff_m1_p1_0 = t296+t334+t338; + coeffs_dy->coeff_0_p1_0 = t289+t318*t295+(t310+t318)*t294+t315; + coeffs_dy->coeff_p1_p1_0 = t298+t338+t339; + coeffs_dy->coeff_m1_m1_p1 = t299+t329+t330; + coeffs_dy->coeff_0_m1_p1 = t297+t332+t336; + coeffs_dy->coeff_p1_m1_p1 = t299+t327+t328; + coeffs_dy->coeff_m1_0_p1 = t254+t248+t308; + coeffs_dy->coeff_0_0_p1 = t293+(-t295+t304)*y; + coeffs_dy->coeff_p1_0_p1 = t254+t249+t309; + coeffs_dy->coeff_m1_p1_p1 = t302+t324+t325; + coeffs_dy->coeff_0_p1_p1 = t300+t333+t336; + coeffs_dy->coeff_p1_p1_p1 = t303+t325+t326; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyy.compute.c new file mode 100644 index 0000000..5328b28 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyy.compute.c @@ -0,0 +1,116 @@ +fp t678; +fp t659; +fp t677; +fp t660; +fp t676; +fp t663; +fp t675; +fp t662; +fp t674; +fp t661; +fp t673; +fp t672; +fp t639; +fp t647; +fp t671; +fp t653; +fp t648; +fp t670; +fp t654; +fp t669; +fp t655; +fp t652; +fp t668; +fp t667; +fp t656; +fp t666; +fp t665; +fp t664; +fp t650; +fp t649; +fp t646; +fp t645; +fp t643; +fp t642; +fp t641; +fp t638; +fp t637; +fp t636; +fp t635; +fp t634; +fp t633; +fp t632; +fp t631; +fp t630; + t678 = x*z; + t659 = RATIONAL(1.0,2.0); + t677 = t659*x; + t660 = RATIONAL(-1.0,2.0); + t676 = t660*z; + t663 = z*z; + t675 = t663*x; + t662 = x*x; + t674 = z*t662; + t661 = t663*t662; + t673 = t661-t662; + t672 = t661-t663; + t639 = t660*t674; + t647 = t660*t661; + t671 = t639+t647; + t653 = t659*t662; + t648 = z*t653; + t670 = t647+t648; + t654 = t659*t663; + t669 = x*t654+t647; + t655 = RATIONAL(1.0,4.0); + t652 = t655*t661; + t668 = t655*t674+t652; + t667 = t663+t662; + t656 = RATIONAL(-1.0,4.0); + t666 = t652+t656*t674; + t665 = x-t675; + t664 = -z+t674; + t650 = t655*t678; + t649 = t656*t675; + t646 = t660*t675; + t645 = z*t677; + t643 = x*t676; + t642 = t656*t678; + t641 = t655*t675; + t638 = t661+RATIONAL(1.0,1.0)-t667; + t637 = t659*z+t654+t671; + t636 = t647+t653+t646+t677; + t635 = t653+t660*x+t669; + t634 = t654+t676+t670; + t633 = t649+t642+t668; + t632 = t650+t641+t668; + t631 = t641+t642+t666; + t630 = t649+t650+t666; + coeffs_dyy->coeff_m1_m1_m1 = t630; + coeffs_dyy->coeff_0_m1_m1 = t634; + coeffs_dyy->coeff_p1_m1_m1 = t631; + coeffs_dyy->coeff_m1_0_m1 = t643+t648+t669; + coeffs_dyy->coeff_0_0_m1 = -t664+t672; + coeffs_dyy->coeff_p1_0_m1 = t645+t646+t670; + coeffs_dyy->coeff_m1_p1_m1 = t630; + coeffs_dyy->coeff_0_p1_m1 = t634; + coeffs_dyy->coeff_p1_p1_m1 = t631; + coeffs_dyy->coeff_m1_m1_0 = t635; + coeffs_dyy->coeff_0_m1_0 = t638; + coeffs_dyy->coeff_p1_m1_0 = t636; + coeffs_dyy->coeff_m1_0_0 = t665+t673; + coeffs_dyy->coeff_0_0_0 = (1.0+t661)*RATIONAL(-2.0,1.0)+t667*RATIONAL(2.0 +,1.0); + coeffs_dyy->coeff_p1_0_0 = -t665+t673; + coeffs_dyy->coeff_m1_p1_0 = t635; + coeffs_dyy->coeff_0_p1_0 = t638; + coeffs_dyy->coeff_p1_p1_0 = t636; + coeffs_dyy->coeff_m1_m1_p1 = t633; + coeffs_dyy->coeff_0_m1_p1 = t637; + coeffs_dyy->coeff_p1_m1_p1 = t632; + coeffs_dyy->coeff_m1_0_p1 = t639+t645+t669; + coeffs_dyy->coeff_0_0_p1 = t664+t672; + coeffs_dyy->coeff_p1_0_p1 = t646+t643+t671; + coeffs_dyy->coeff_m1_p1_p1 = t633; + coeffs_dyy->coeff_0_p1_p1 = t637; + coeffs_dyy->coeff_p1_p1_p1 = t632; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyz.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyz.compute.c new file mode 100644 index 0000000..a17d581 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyz.compute.c @@ -0,0 +1,164 @@ +fp t716; +fp t720; +fp t714; +fp t737; +fp t694; +fp t682; +fp t756; +fp t711; +fp t740; +fp t693; +fp t710; +fp t742; +fp t696; +fp t755; +fp t715; +fp t754; +fp t705; +fp t753; +fp t691; +fp t752; +fp t717; +fp t738; +fp t692; +fp t681; +fp t708; +fp t751; +fp t699; +fp t750; +fp t704; +fp t739; +fp t684; +fp t680; +fp t749; +fp t697; +fp t701; +fp t688; +fp t748; +fp t695; +fp t747; +fp t746; +fp t698; +fp t686; +fp t745; +fp t735; +fp t687; +fp t744; +fp t743; +fp t719; +fp t741; +fp t733; +fp t679; +fp t690; +fp t736; +fp t734; +fp t732; +fp t731; +fp t730; +fp t729; +fp t728; +fp t727; +fp t726; +fp t725; +fp t724; +fp t723; +fp t722; +fp t721; +fp t718; +fp t689; + t716 = RATIONAL(1.0,8.0); + t720 = x*x; + t714 = RATIONAL(1.0,2.0); + t737 = t714*t720; + t694 = y*t737; + t682 = z*t694; + t756 = t682+t716*t720; + t711 = RATIONAL(-1.0,4.0); + t740 = x*t711; + t693 = y*t740; + t710 = RATIONAL(1.0,4.0); + t742 = x*t710; + t696 = z*t742; + t755 = t693+t696; + t715 = RATIONAL(-1.0,8.0); + t754 = t715*t720+t682; + t705 = t711*t720; + t753 = t705+t710; + t691 = z*t740; + t752 = t693+t691; + t717 = RATIONAL(-1.0,2.0); + t738 = x*t717; + t692 = y*t738; + t681 = z*t692; + t708 = t715*x; + t751 = t681+t708; + t699 = y*t742; + t750 = t696+t699; + t704 = t716*x; + t739 = x*t714; + t684 = z*t739; + t680 = y*t684; + t749 = t704+t680; + t697 = y*t705; + t701 = t710*t720; + t688 = z*t701; + t748 = t697+t688; + t695 = z*t737; + t747 = t717*z+t695; + t746 = t711+t701; + t698 = z*t705; + t686 = y*t701; + t745 = t698+t686; + t735 = t717*t720; + t687 = z*t735; + t744 = t714*z+t687; + t743 = t699+t691; + t719 = y*z; + t741 = y*t720; + t733 = RATIONAL(2.0,1.0)*t719; + t679 = t720*t733; + t690 = RATIONAL(-2.0,1.0)*t719; + t736 = t679+t690; + t734 = z*t741; + t732 = (t720-1.0)*z; + t731 = -y+t741; + t730 = -t734+z*t738; + t729 = -t734+t684; + t728 = t694-t734; + t727 = y*t735-t734; + t726 = t704+t681+t754; + t725 = t686+t688+t756; + t724 = t697+t698+t756; + t723 = t680+t708+t754; + t722 = t719+t714*y+t727; + t721 = t719+t717*y+t728; + t718 = x*t719; + t689 = y*t739; + coeffs_dyz->coeff_m1_m1_m1 = t724+t750+t751; + coeffs_dyz->coeff_0_m1_m1 = t721+t747+t753; + coeffs_dyz->coeff_p1_m1_m1 = t724+t749+t752; + coeffs_dyz->coeff_m1_0_m1 = t692+t718+t728; + coeffs_dyz->coeff_0_0_m1 = -t731+t736; + coeffs_dyz->coeff_p1_0_m1 = t689-t718+t728; + coeffs_dyz->coeff_m1_p1_m1 = t726+t743+t748; + coeffs_dyz->coeff_0_p1_m1 = t721+t744+t746; + coeffs_dyz->coeff_p1_p1_m1 = t723+t748+t755; + coeffs_dyz->coeff_m1_m1_0 = t695+t718+t730; + coeffs_dyz->coeff_0_m1_0 = -t732+t736; + coeffs_dyz->coeff_p1_m1_0 = t695-t718+t729; + coeffs_dyz->coeff_m1_0_0 = x*t690+t679; + coeffs_dyz->coeff_0_0_0 = (RATIONAL(4.0,1.0)+RATIONAL(-4.0,1.0)*t720)* +t719; + coeffs_dyz->coeff_p1_0_0 = x*t733+t679; + coeffs_dyz->coeff_m1_p1_0 = t687+t718+t729; + coeffs_dyz->coeff_0_p1_0 = t732+t736; + coeffs_dyz->coeff_p1_p1_0 = t687-t718+t730; + coeffs_dyz->coeff_m1_m1_p1 = t726+t745+t755; + coeffs_dyz->coeff_0_m1_p1 = t722+t746+t747; + coeffs_dyz->coeff_p1_m1_p1 = t723+t743+t745; + coeffs_dyz->coeff_m1_0_p1 = t689+t718+t727; + coeffs_dyz->coeff_0_0_p1 = t731+t736; + coeffs_dyz->coeff_p1_0_p1 = t692-t718+t727; + coeffs_dyz->coeff_m1_p1_p1 = t725+t751+t752; + coeffs_dyz->coeff_0_p1_p1 = t722+t744+t753; + coeffs_dyz->coeff_p1_p1_p1 = t725+t749+t750; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dz.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dz.compute.c new file mode 100644 index 0000000..19e7ca2 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dz.compute.c @@ -0,0 +1,196 @@ +fp t383; +fp t382; +fp t392; +fp t375; +fp t359; +fp t434; +fp t386; +fp t414; +fp t369; +fp t387; +fp t391; +fp t403; +fp t390; +fp t348; +fp t433; +fp t365; +fp t353; +fp t388; +fp t432; +fp t389; +fp t431; +fp t408; +fp t373; +fp t347; +fp t430; +fp t409; +fp t358; +fp t379; +fp t371; +fp t351; +fp t429; +fp t368; +fp t374; +fp t356; +fp t345; +fp t428; +fp t363; +fp t378; +fp t364; +fp t352; +fp t427; +fp t426; +fp t404; +fp t360; +fp t349; +fp t410; +fp t370; +fp t425; +fp t344; +fp t424; +fp t354; +fp t423; +fp t422; +fp t421; +fp t420; +fp t415; +fp t367; +fp t407; +fp t350; +fp t419; +fp t418; +fp t417; +fp t366; +fp t341; +fp t416; +fp t413; +fp t412; +fp t411; +fp t406; +fp t362; +fp t405; +fp t361; +fp t402; +fp t401; +fp t400; +fp t399; +fp t398; +fp t397; +fp t396; +fp t395; +fp t394; +fp t393; +fp t343; +fp t342; + t383 = RATIONAL(-1.0,4.0); + t382 = RATIONAL(1.0,4.0); + t392 = y*y; + t375 = t382*t392; + t359 = x*t375; + t434 = t383*x+t359; + t386 = RATIONAL(1.0,2.0); + t414 = z*t386; + t369 = t392*t414; + t387 = RATIONAL(-1.0,2.0); + t391 = x*x; + t403 = t392*t391; + t390 = z*t403; + t348 = t387*t390; + t433 = t369+t348; + t365 = t382*t403; + t353 = z*t365; + t388 = RATIONAL(-1.0,8.0); + t432 = t353+t388*t403; + t389 = RATIONAL(1.0,8.0); + t431 = t353+t389*t403; + t408 = x*t392; + t373 = t389*t408; + t347 = z*t359; + t430 = t373+t347; + t409 = t391*y; + t358 = t388*t409; + t379 = t383*t391; + t371 = y*t379; + t351 = z*t371; + t429 = t358+t351; + t368 = t389*t409; + t374 = t382*t391; + t356 = y*t374; + t345 = z*t356; + t428 = t368+t345; + t363 = t388*t408; + t378 = t383*t392; + t364 = x*t378; + t352 = z*t364; + t427 = t363+t352; + t426 = t378+t365; + t404 = t386*t391; + t360 = z*t404; + t349 = y*t360; + t410 = z*t387; + t370 = y*t410; + t425 = t349+t370; + t344 = t391*t370; + t424 = t344+y*t414; + t354 = t383*t403; + t423 = t354+t375; + t422 = t354+t374; + t421 = t382*x+t364; + t420 = t379+t365; + t415 = y*x; + t367 = t388*t415; + t407 = z*t415; + t350 = t382*t407; + t419 = t367+t350; + t418 = t345+t358; + t417 = t368+t351; + t366 = t389*t415; + t341 = t383*t407; + t416 = t366+t341; + t413 = z+t390; + t412 = t386-z; + t411 = -z+t387; + t406 = t348+x*t369; + t362 = x*t410; + t405 = t348+t392*t362; + t361 = x*t414; + t402 = t409-y; + t401 = -x+t408; + t400 = t373+t352+t432; + t399 = t360+t362+t406; + t398 = t371+t382*y+t433; + t397 = t363+t347+t432; + t396 = t360+t361+t405; + t395 = t366+t350+t431; + t394 = t356+t383*y+t433; + t393 = t367+t341+t431; + t343 = x*t370; + t342 = y*t361; + coeffs_dz->coeff_m1_m1_m1 = t400+t417+t419; + coeffs_dz->coeff_0_m1_m1 = t398+t425+t426; + coeffs_dz->coeff_p1_m1_m1 = t397+t416+t417; + coeffs_dz->coeff_m1_0_m1 = t399+t420+t421; + coeffs_dz->coeff_0_0_m1 = t387+t412*t391+(t387*t391+t412)*t392+t413; + coeffs_dz->coeff_p1_0_m1 = t396+t420+t434; + coeffs_dz->coeff_m1_p1_m1 = t400+t416+t418; + coeffs_dz->coeff_0_p1_m1 = t394+t424+t426; + coeffs_dz->coeff_p1_p1_m1 = t397+t418+t419; + coeffs_dz->coeff_m1_m1_0 = t343+t349+t406; + coeffs_dz->coeff_0_m1_0 = t390+(-t392-t402)*z; + coeffs_dz->coeff_p1_m1_0 = t342+t349+t405; + coeffs_dz->coeff_m1_0_0 = t390+(-t391-t401)*z; + coeffs_dz->coeff_0_0_0 = ((1.0+t403)*RATIONAL(-2.0,1.0)+(t392+t391)* +RATIONAL(2.0,1.0))*z; + coeffs_dz->coeff_p1_0_0 = t390+(-t391+t401)*z; + coeffs_dz->coeff_m1_p1_0 = t344+t342+t406; + coeffs_dz->coeff_0_p1_0 = t390+(-t392+t402)*z; + coeffs_dz->coeff_p1_p1_0 = t344+t343+t405; + coeffs_dz->coeff_m1_m1_p1 = t395+t427+t429; + coeffs_dz->coeff_0_m1_p1 = t394+t423+t425; + coeffs_dz->coeff_p1_m1_p1 = t393+t429+t430; + coeffs_dz->coeff_m1_0_p1 = t399+t422+t434; + coeffs_dz->coeff_0_0_p1 = t386+t411*t391+(t404+t411)*t392+t413; + coeffs_dz->coeff_p1_0_p1 = t396+t421+t422; + coeffs_dz->coeff_m1_p1_p1 = t393+t427+t428; + coeffs_dz->coeff_0_p1_p1 = t398+t423+t424; + coeffs_dz->coeff_p1_p1_p1 = t395+t428+t430; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dzz.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dzz.compute.c new file mode 100644 index 0000000..c4ac417 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dzz.compute.c @@ -0,0 +1,120 @@ +fp t786; +fp t808; +fp t787; +fp t807; +fp t790; +fp t806; +fp t783; +fp t805; +fp t782; +fp t804; +fp t789; +fp t788; +fp t803; +fp t793; +fp t776; +fp t802; +fp t801; +fp t800; +fp t794; +fp t771; +fp t799; +fp t798; +fp t797; +fp t781; +fp t796; +fp t780; +fp t795; +fp t792; +fp t791; +fp t778; +fp t775; +fp t774; +fp t772; +fp t770; +fp t769; +fp t768; +fp t767; +fp t765; +fp t764; +fp t763; +fp t762; +fp t761; +fp t760; +fp t759; +fp t758; +fp t757; + t786 = RATIONAL(1.0,2.0); + t808 = t786*y; + t787 = RATIONAL(-1.0,2.0); + t807 = t787*x; + t790 = y*y; + t806 = t790*x; + t783 = RATIONAL(-1.0,4.0); + t805 = y*t783; + t782 = RATIONAL(1.0,4.0); + t804 = y*t782; + t789 = x*x; + t788 = t790*t789; + t803 = -t790+t788; + t793 = t787*t789; + t776 = t790*t793; + t802 = t776+y*t807; + t801 = -t789+t788; + t800 = t789+t790; + t794 = t782*t790; + t771 = t789*t794; + t799 = t783*t806+t771; + t798 = t771+x*t794; + t797 = t776+x*t808; + t781 = t786*t790; + t796 = t776+t781; + t780 = t786*t789; + t795 = t776+t780; + t792 = x-t806; + t791 = (1.0-t789)*y; + t778 = x*t805; + t775 = y*t780; + t774 = x*t804; + t772 = t789*t805; + t770 = x*t781; + t769 = t789*t804; + t768 = t787*t806; + t767 = y*t793; + t765 = RATIONAL(1.0,1.0)+t788-t800; + t764 = t775+t787*y+t796; + t763 = t770+t807+t795; + t762 = t786*x+t768+t795; + t761 = t808+t767+t796; + t760 = t772+t774+t799; + t759 = t772+t778+t798; + t758 = t769+t778+t799; + t757 = t769+t774+t798; + coeffs_dzz->coeff_m1_m1_m1 = t760; + coeffs_dzz->coeff_0_m1_m1 = t764; + coeffs_dzz->coeff_p1_m1_m1 = t759; + coeffs_dzz->coeff_m1_0_m1 = t763; + coeffs_dzz->coeff_0_0_m1 = t765; + coeffs_dzz->coeff_p1_0_m1 = t762; + coeffs_dzz->coeff_m1_p1_m1 = t758; + coeffs_dzz->coeff_0_p1_m1 = t761; + coeffs_dzz->coeff_p1_p1_m1 = t757; + coeffs_dzz->coeff_m1_m1_0 = t775+t770+t802; + coeffs_dzz->coeff_0_m1_0 = t791+t803; + coeffs_dzz->coeff_p1_m1_0 = t775+t768+t797; + coeffs_dzz->coeff_m1_0_0 = t792+t801; + coeffs_dzz->coeff_0_0_0 = (1.0+t788)*RATIONAL(-2.0,1.0)+t800*RATIONAL(2.0 +,1.0); + coeffs_dzz->coeff_p1_0_0 = -t792+t801; + coeffs_dzz->coeff_m1_p1_0 = t770+t767+t797; + coeffs_dzz->coeff_0_p1_0 = -t791+t803; + coeffs_dzz->coeff_p1_p1_0 = t768+t767+t802; + coeffs_dzz->coeff_m1_m1_p1 = t760; + coeffs_dzz->coeff_0_m1_p1 = t764; + coeffs_dzz->coeff_p1_m1_p1 = t759; + coeffs_dzz->coeff_m1_0_p1 = t763; + coeffs_dzz->coeff_0_0_p1 = t765; + coeffs_dzz->coeff_p1_0_p1 = t762; + coeffs_dzz->coeff_m1_p1_p1 = t758; + coeffs_dzz->coeff_0_p1_p1 = t761; + coeffs_dzz->coeff_p1_p1_p1 = t757; |