diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c')
-rw-r--r-- | src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c | 394 |
1 files changed, 199 insertions, 195 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c index 49655f0..cfe391a 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c @@ -1,239 +1,243 @@ -fp t413, +fp t417, t418, - t466, - t403, - t417, - t358, + t468, t414, + t411, + t361, + t416, t373, - t438, - t444, - t465, - t405, - t364, + t435, + t446, + t467, + t410, t367, - t439, + t368, + t437, + t466, + t465, + t441, t464, - t436, t463, - t410, - t359, - t442, - t409, - t318, - t392, - t462, - t411, - t360, - t441, + t395, + t445, t400, - t298, - t461, + t403, t366, - t443, - t395, + t462, + t397, + t404, + t362, t401, + t461, + t389, + t409, + t316, t460, - t440, - t296, - t399, + t402, + t312, + t384, t459, - t372, + t412, + t443, + t407, + t309, t458, - t333, + t295, + t399, + t442, t457, - t396, - t412, - t327, + t406, + t405, + t319, t456, - t407, - t331, - t391, + t444, + t370, t455, - t402, - t313, - t397, + t315, t454, t325, t453, - t324, + t358, + t398, + t294, t452, - t311, t451, + t413, + t364, t450, - t387, - t398, + t357, + t380, + t381, t449, + t332, t415, - t371, - t416, - t357, - t383, - t382, - t448, - t317, t356, + t448, t447, - t362, - t446, - t445, - t437, - t435, - t393, - t404, - t369, - t376, + t440, + t360, + t439, + t438, + t436, + t390, + t408, + t363, + t392, t434, - t365, - t384, t386, + t374, t433, + t372, t432, t431, + t317, t430, t429, + t326, t428, t427, - t406, + t385, t426, t425, - t408, + t388, t424, t423, t394, - t389, - t388, - t385; - t413 = x*x; + t393, + t391, + t387; + t417 = x*x; t418 = y*y; - t466 = t418+t413; - t403 = RATIONAL(-1.0,60.0); - t417 = y*t418; - t358 = t403*t417; - t414 = t413*x; - t373 = t403*t414; - t438 = y*t373+x*t358; - t444 = x*y; - t465 = t438+RATIONAL(37.0,300.0)*t444; - t405 = RATIONAL(1.0,60.0); - t364 = t405*t414; - t367 = t405*t417; - t439 = y*t364+x*t367; - t464 = t439+RATIONAL(-37.0,300.0)*t444; - t436 = t413*t418; - t463 = RATIONAL(-99.0,1225.0)+RATIONAL(-1.0,98.0)*t436; - t410 = RATIONAL(-1.0,30.0); - t359 = t410*t417; - t442 = y*t413; - t409 = RATIONAL(1.0,70.0); - t318 = t409*t442; - t392 = RATIONAL(11.0,105.0); - t462 = t359+t318+t392*y; - t411 = RATIONAL(1.0,30.0); - t360 = t411*t417; - t441 = y*t414; - t400 = RATIONAL(-1.0,120.0); - t298 = RATIONAL(-21.0,200.0)*t444; - t461 = x*t360+t298+t400*t441; - t366 = t410*t414; - t443 = x*t418; - t395 = RATIONAL(31.0,210.0); + t468 = t418+t417; + t414 = t417*x; + t411 = RATIONAL(1.0,60.0); + t361 = t411*t414; + t416 = y*t418; + t373 = t411*t416; + t435 = x*t373+y*t361; + t446 = x*y; + t467 = RATIONAL(-11.0,150.0)*t446+t435; + t410 = RATIONAL(-1.0,60.0); + t367 = t410*t416; + t368 = t410*t414; + t437 = x*t367+y*t368; + t466 = RATIONAL(37.0,300.0)*t446+t437; + t465 = t437+RATIONAL(11.0,150.0)*t446; + t441 = t417*t418; + t464 = RATIONAL(-1.0,98.0)*t441+RATIONAL(-99.0,1225.0); + t463 = RATIONAL(-37.0,300.0)*t446+t435; + t395 = RATIONAL(-31.0,210.0); + t445 = x*t418; + t400 = RATIONAL(1.0,140.0); + t403 = RATIONAL(1.0,30.0); + t366 = t403*t414; + t462 = t400*t445+t395*x+t366; + t397 = RATIONAL(31.0,210.0); + t404 = RATIONAL(-1.0,30.0); + t362 = t404*t414; t401 = RATIONAL(-1.0,140.0); - t460 = t395*x+t366+t401*t443; - t440 = x*t417; - t296 = RATIONAL(21.0,200.0)*t444; + t461 = t362+t401*t445+t397*x; + t389 = RATIONAL(31.0,420.0); + t409 = RATIONAL(-1.0,35.0); + t316 = t409*t445; + t460 = t389*x+t316+t368; + t402 = RATIONAL(1.0,35.0); + t312 = t402*t445; + t384 = RATIONAL(-31.0,420.0); + t459 = t312+t384*x+t361; + t412 = RATIONAL(1.0,84.0); + t443 = y*t417; + t407 = RATIONAL(-1.0,70.0); + t309 = t407*t443; + t458 = t309+t412*y+t373; + t295 = RATIONAL(21.0,200.0)*t446; t399 = RATIONAL(1.0,120.0); - t459 = t296+y*t366+t399*t440; - t372 = t411*t414; - t458 = t400*t440+y*t372+t298; - t333 = t409*t443; - t457 = t366+t333+t392*x; - t396 = RATIONAL(-11.0,105.0); - t412 = RATIONAL(-1.0,70.0); - t327 = t412*t443; - t456 = t396*x+t372+t327; - t407 = RATIONAL(1.0,35.0); - t331 = t407*t442; - t391 = RATIONAL(-31.0,420.0); - t455 = t331+t391*y+t367; - t402 = RATIONAL(-1.0,35.0); - t313 = t402*t442; - t397 = RATIONAL(31.0,420.0); - t454 = t358+t397*y+t313; - t325 = t402*t443; - t453 = t397*x+t325+t373; - t324 = t407*t443; - t452 = t364+t324+t391*x; - t311 = t412*t442; - t451 = t360+t311+t396*y; - t450 = x*t359+t296+t399*t441; - t387 = RATIONAL(-31.0,210.0); - t398 = RATIONAL(1.0,140.0); - t449 = t372+t387*x+t398*t443; - t415 = t413*t413; - t371 = t410*t415; - t416 = t418*t418; - t357 = t399*t416; - t383 = RATIONAL(71.0,5880.0); - t382 = RATIONAL(113.0,735.0); - t448 = t371+t382*t413+t383*t418+t357+t463; - t317 = RATIONAL(1.0,49.0)*t436; + t442 = x*t416; + t457 = t295+t399*t442+y*t362; + t406 = RATIONAL(-1.0,84.0); + t405 = RATIONAL(1.0,70.0); + t319 = t405*t443; + t456 = t406*y+t319+t367; + t444 = t414*y; + t370 = t404*t416; + t455 = t295+t399*t444+x*t370; + t315 = t405*t445; + t454 = t406*x+t315+t368; + t325 = t407*t445; + t453 = t361+t412*x+t325; + t358 = t403*t416; + t398 = RATIONAL(-1.0,120.0); + t294 = RATIONAL(-21.0,200.0)*t446; + t452 = t398*t444+t294+x*t358; + t451 = y*t366+t398*t442+t294; + t413 = t417*t417; + t364 = t404*t413; + t450 = RATIONAL(-24.0,1225.0)+t364+RATIONAL(1.0,196.0)*t441+t468*RATIONAL +(181.0,1470.0); + t357 = t399*t413; + t380 = RATIONAL(71.0,5880.0); + t381 = RATIONAL(113.0,735.0); + t449 = t357+t380*t417+t381*t418+t464; + t332 = RATIONAL(1.0,49.0)*t441; + t415 = t418*t418; t356 = t399*t415; - t447 = t317+RATIONAL(51.0,1225.0)+t356+t357+t466*RATIONAL(-289.0,5880.0); - t362 = t410*t416; - t446 = t382*t418+t356+t383*t413+t362+t463; - t445 = t362+RATIONAL(1.0,196.0)*t436+t371+RATIONAL(-24.0,1225.0)+t466* -RATIONAL(181.0,1470.0); - t437 = RATIONAL(-1.0,49.0)*t436+RATIONAL(96.0,1225.0); - t435 = RATIONAL(1.0,98.0)*t436+RATIONAL(246.0,1225.0); - t393 = RATIONAL(-41.0,196.0); - t404 = RATIONAL(1.0,20.0); - t369 = t404*t415; - t376 = RATIONAL(191.0,5880.0); - t434 = t376*t418+t369+t357+t393*t413+t437; - t365 = t404*t416; - t384 = RATIONAL(-53.0,196.0); - t386 = RATIONAL(83.0,735.0); - t433 = t365+t371+t384*t418+t386*t413+t435; - t432 = t365+t376*t413+t356+t393*t418+t437; - t431 = t369+t384*t413+t362+t386*t418+t435; - t430 = RATIONAL(-11.0,150.0)*t444+t439+t447; - t429 = t359+t395*y+t401*t442+t445; - t428 = t360+t387*y+t398*t442+t445; - t427 = RATIONAL(11.0,150.0)*t444+t438+t447; - t406 = RATIONAL(1.0,84.0); - t426 = t367+t311+t406*y+t448; - t425 = t364+t406*x+t327+t446; - t408 = RATIONAL(-1.0,84.0); - t424 = t358+t318+t408*y+t448; - t423 = t333+t408*x+t373+t446; - t394 = RATIONAL(-17.0,420.0); - t389 = RATIONAL(17.0,105.0); - t388 = RATIONAL(17.0,420.0); - t385 = RATIONAL(-17.0,105.0); - coeff_I_m2_m2 = t430+t453+t454; - coeff_I_m1_m2 = t424+t456+t459; - coeff_I_0_m2 = t331+t394*y+t358+t434; - coeff_I_p1_m2 = t424+t457+t458; - coeff_I_p2_m2 = t427+t452+t454; - coeff_I_m2_m1 = t423+t450+t451; - coeff_I_m1_m1 = t428+t449+t465; - coeff_I_0_m1 = t360+t385*y+t318+t431; - coeff_I_p1_m1 = t428+t460+t464; - coeff_I_p2_m1 = t425+t451+t461; - coeff_I_m2_0 = t373+t394*x+t324+t432; - coeff_I_m1_0 = t385*x+t333+t372+t433; - coeff_I_0_0 = t317+t365+t369+RATIONAL(541.0,1225.0)+t466*RATIONAL(-57.0, + t448 = RATIONAL(51.0,1225.0)+t357+t332+t356+t468*RATIONAL(-289.0,5880.0); + t447 = t380*t418+t364+t381*t417+t356+t464; + t440 = RATIONAL(1.0,98.0)*t441+RATIONAL(246.0,1225.0); + t360 = t404*t415; + t439 = t360+t370; + t438 = t358+t360; + t436 = RATIONAL(-1.0,49.0)*t441+RATIONAL(96.0,1225.0); + t390 = RATIONAL(83.0,735.0); + t408 = RATIONAL(1.0,20.0); + t363 = t408*t413; + t392 = RATIONAL(-53.0,196.0); + t434 = t360+t390*t418+t363+t392*t417+t440; + t386 = RATIONAL(-41.0,196.0); + t374 = RATIONAL(191.0,5880.0); + t433 = t386*t417+t374*t418+t363+t356+t436; + t372 = t408*t415; + t432 = t357+t374*t417+t386*t418+t372+t436; + t431 = t364+t392*t418+t390*t417+t372+t440; + t317 = t402*t443; + t430 = t384*y+t317+t373+t448; + t429 = t395*y+t400*t443+t438+t450; + t326 = t409*t443; + t428 = t389*y+t367+t326+t448; + t427 = t397*y+t401*t443+t439+t450; + t385 = RATIONAL(11.0,105.0); + t426 = t385*x+t315+t362+t447; + t425 = t385*y+t319+t439+t449; + t388 = RATIONAL(-11.0,105.0); + t424 = t388*y+t309+t438+t449; + t423 = t388*x+t366+t325+t447; + t394 = RATIONAL(17.0,420.0); + t393 = RATIONAL(17.0,105.0); + t391 = RATIONAL(-17.0,105.0); + t387 = RATIONAL(-17.0,420.0); + coeff_I_m2_m2 = t428+t460+t467; + coeff_I_m1_m2 = t423+t456+t457; + coeff_I_0_m2 = t317+t367+t387*y+t433; + coeff_I_p1_m2 = t426+t451+t456; + coeff_I_p2_m2 = t428+t459+t465; + coeff_I_m2_m1 = t424+t454+t455; + coeff_I_m1_m1 = t429+t462+t466; + coeff_I_0_m1 = t319+t391*y+t358+t434; + coeff_I_p1_m1 = t429+t461+t463; + coeff_I_p2_m1 = t424+t452+t453; + coeff_I_m2_0 = t368+t312+t387*x+t432; + coeff_I_m1_0 = t366+t315+t391*x+t431; + coeff_I_0_0 = RATIONAL(541.0,1225.0)+t363+t372+t332+t468*RATIONAL(-57.0, 196.0); - coeff_I_p1_0 = t389*x+t327+t366+t433; - coeff_I_p2_0 = t325+t364+t388*x+t432; - coeff_I_m2_p1 = t423+t461+t462; - coeff_I_m1_p1 = t429+t449+t464; - coeff_I_0_p1 = t359+t389*y+t311+t431; - coeff_I_p1_p1 = t429+t460+t465; - coeff_I_p2_p1 = t425+t450+t462; - coeff_I_m2_p2 = t427+t453+t455; - coeff_I_m1_p2 = t426+t456+t458; - coeff_I_0_p2 = t313+t388*y+t367+t434; - coeff_I_p1_p2 = t426+t457+t459; - coeff_I_p2_p2 = t430+t452+t455; + coeff_I_p1_0 = t325+t362+t393*x+t431; + coeff_I_p2_0 = t394*x+t316+t361+t432; + coeff_I_m2_p1 = t425+t452+t454; + coeff_I_m1_p1 = t427+t462+t463; + coeff_I_0_p1 = t309+t370+t393*y+t434; + coeff_I_p1_p1 = t427+t461+t466; + coeff_I_p2_p1 = t425+t453+t455; + coeff_I_m2_p2 = t430+t460+t465; + coeff_I_m1_p2 = t423+t451+t458; + coeff_I_0_p2 = t394*y+t326+t373+t433; + coeff_I_p1_p2 = t426+t457+t458; + coeff_I_p2_p2 = t430+t459+t467; |