diff options
Diffstat (limited to 'src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c')
-rw-r--r-- | src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c | 245 |
1 files changed, 245 insertions, 0 deletions
diff --git a/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c b/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..d629ea0 --- /dev/null +++ b/src/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; |