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