diff options
Diffstat (limited to 'src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c')
-rw-r--r-- | src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c | 340 |
1 files changed, 176 insertions, 164 deletions
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; |