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