diff options
author | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-08-18 15:12:49 +0000 |
---|---|---|
committer | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-08-18 15:12:49 +0000 |
commit | 15c98b3abee2525b24ca89932256ff2cda63c9bc (patch) | |
tree | fba7c114f0b5e7ffdad3e60254fd0e8dd1e42f97 /src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0 | |
parent | 7be0a94bec4b4e89c433e5380fce2cf84e5e0681 (diff) |
This commit reorganizes the Maple scripts which generate the interpolation
coefficients, and the coefficient files themselves, to properly support
multiple interpolation operators, and adds some skeleton support for
Hermite interpolation. I will add full support for Hermite in a future
checkin.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@81 df1f8a13-aa1d-4dd4-9681-27ded5b42416
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0')
10 files changed, 3619 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-I.compute.c new file mode 100644 index 0000000..8736aff --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-I.compute.c @@ -0,0 +1,906 @@ +fp t484, + t501, + t340, + t244, + t502, + t505, + t493, + t359, + t730, + t429, + t617, + t405, + t729, + t486, + t597, + t272, + t582, + t265, + t606, + t506, + t503, + t491, + t371, + t728, + t477, + t339, + t243, + t586, + t727, + t435, + t441, + t635, + t726, + t488, + t343, + t251, + t504, + t500, + t373, + t725, + t362, + t724, + t242, + t367, + t723, + t331, + t223, + t722, + t469, + t631, + t210, + t159, + t721, + t645, + t610, + t157, + t720, + t332, + t269, + t498, + t357, + t719, + t718, + t368, + t268, + t717, + t716, + t346, + t232, + t358, + t715, + t467, + t391, + t714, + t495, + t379, + t376, + t713, + t341, + t227, + t712, + t464, + t614, + t197, + t154, + t711, + t334, + t221, + t710, + t158, + t709, + t708, + t349, + t707, + t648, + t392, + t481, + t706, + t705, + t245, + t590, + t704, + t470, + t703, + t702, + t701, + t700, + t699, + t698, + t697, + t696, + t499, + t608, + t262, + t695, + t153, + t694, + t228, + t693, + t692, + t691, + t344, + t254, + t690, + t419, + t629, + t400, + t689, + t363, + t641, + t688, + t401, + t420, + t687, + t686, + t475, + t612, + t388, + t685, + t476, + t489, + t450, + t627, + t169, + t684, + t461, + t459, + t327, + t398, + t683, + t471, + t487, + t347, + t395, + t682, + t428, + t465, + t490, + t681, + t348, + t640, + t680, + t679, + t647, + t678, + t619, + t328, + t677, + t611, + t676, + t496, + t354, + t433, + t675, + t438, + t439, + t497, + t353, + t674, + t415, + t448, + t673, + t636, + t672, + t455, + t451, + t638, + t492, + t671, + t166, + t670, + t453, + t458, + t669, + t668, + t479, + t478, + t336, + t407, + t667, + t445, + t436, + t183, + t494, + t375, + t666, + t337, + t644, + t665, + t383, + t642, + t184, + t664, + t364, + t624, + t663, + t662, + t372, + t483, + t412, + t661, + t660, + t618, + t338, + t659, + t325, + t658, + t360, + t190, + t657, + t656, + t381, + t637, + t655, + t378, + t654, + t653, + t634, + t652, + t351, + t651, + t621, + t650, + t649, + t646, + t425, + t643, + t446, + t639, + t430, + t633, + t413, + t632, + t485, + t630, + t473, + t628, + t434, + t626, + t456, + t625, + t431, + t623, + t474, + t622, + t472, + t620, + t424, + t616, + t403, + t615, + t426, + t613, + t266, + t609, + t386, + t607, + t605, + t604, + t247, + t603, + t226, + t602, + t601, + t257, + t249, + t600, + t377, + t382, + t598, + t596, + t595, + t259, + t594, + t380, + t593, + t256, + t592, + t591, + t589, + t588, + t587, + t585, + t584, + t583, + t581, + t263, + t580, + t579, + t208, + t191, + t578, + t577, + t576, + t575, + t574, + t573, + t156, + t572, + t571, + t570, + t569, + t568, + t567, + t566, + t155, + t565, + t564, + t563, + t562, + t561, + t560, + t559, + t558, + t557, + t556, + t152, + t555, + t554, + t553, + t552, + t551, + t550, + t549, + t548, + t547, + t546, + t545, + t544, + t543, + t542, + t541, + t540, + t402, + t539, + t538, + t537, + t536, + t192, + t535, + t534, + t345, + t533, + t532, + t531, + t530, + t529, + t528, + t527, + t437, + t526, + t525, + t449, + t524, + t523, + t399, + t522, + t521, + t447, + t520, + t389, + t519, + t411, + t518, + t517, + t516, + t515, + t514, + t513, + t512, + t511, + t510, + t480, + t466, + t463, + t462, + t460, + t457, + t452, + t443, + t442, + t417, + t416, + t414, + t410, + t409, + t408, + t404, + t396, + t394, + t393, + t390, + t385, + t370, + t178, + t177, + t161; + t484 = RATIONAL(1.0,160.0); + t501 = z*z; + t340 = t484*t501; + t244 = y*t340; + t502 = y*y; + t505 = t502*y; + t493 = RATIONAL(-1.0,32.0); + t359 = t493*t505; + t730 = t244+t359; + t429 = RATIONAL(-31.0,1000.0); + t617 = z*t429; + t405 = RATIONAL(141.0,2000.0); + t729 = x*t617+t405*y; + t486 = RATIONAL(-3.0,160.0); + t597 = t486*t502; + t272 = x*t597; + t582 = t486*t501; + t265 = x*t582; + t606 = t272+t265; + t506 = x*x; + t503 = t506*x; + t491 = RATIONAL(1.0,96.0); + t371 = t491*t503; + t728 = t606+t371; + t477 = RATIONAL(3.0,160.0); + t339 = t477*t502; + t243 = x*t339; + t586 = t243+t265; + t727 = t586+t371; + t435 = RATIONAL(-97.0,3000.0); + t441 = RATIONAL(-37.0,1000.0); + t635 = z*t441; + t726 = x*t635+t435*y; + t488 = RATIONAL(-1.0,160.0); + t343 = t488*t502; + t251 = z*t343; + t504 = t501*z; + t500 = RATIONAL(1.0,32.0); + t373 = t500*t504; + t725 = t251+t373; + t362 = t493*t504; + t724 = t362+t251; + t242 = y*t582; + t367 = t491*t505; + t723 = t242+t367; + t331 = t484*t502; + t223 = z*t331; + t722 = t362+t223; + t469 = RATIONAL(-9.0,1000.0); + t631 = x*t469; + t210 = z*t631; + t159 = y*t210; + t721 = RATIONAL(71.0,2000.0)+t159; + t645 = x*y; + t610 = z*t645; + t157 = RATIONAL(3.0,1000.0)*t610; + t720 = RATIONAL(143.0,2000.0)+t157; + t332 = t477*t501; + t269 = y*t332; + t498 = RATIONAL(-1.0,96.0); + t357 = t498*t505; + t719 = t269+t357; + t718 = t441*t645+t435*z; + t368 = t498*t504; + t268 = z*t597; + t717 = t368+t268; + t716 = RATIONAL(-59.0,2000.0)+t159; + t346 = t488*t501; + t232 = y*t346; + t358 = t500*t505; + t715 = t232+t358; + t467 = RATIONAL(11.0,500.0); + t391 = RATIONAL(-167.0,12000.0); + t714 = t467*t645+t391*z; + t495 = RATIONAL(3.0,80.0); + t379 = t495*t501; + t376 = t495*t502; + t713 = t379+t376; + t341 = t488*t506; + t227 = y*t341; + t712 = t227+t730; + t464 = RATIONAL(-3.0,1000.0); + t614 = y*t464; + t197 = z*t614; + t154 = x*t197; + t711 = RATIONAL(49.0,4000.0)+t154; + t334 = t484*t506; + t221 = y*t334; + t710 = t221+t359; + t158 = RATIONAL(9.0,1000.0)*t610; + t709 = RATIONAL(-97.0,4000.0)+t158; + t708 = t242+t357; + t349 = t491*t504; + t707 = t349+t268; + t648 = y*z; + t392 = RATIONAL(281.0,12000.0); + t481 = RATIONAL(-9.0,250.0); + t706 = t392*x+t481*t648; + t705 = t223+t373; + t245 = x*t332; + t590 = t243+t245; + t704 = t371+t590; + t470 = RATIONAL(-41.0,750.0); + t703 = t470*y+t210; + t702 = RATIONAL(-67.0,4000.0)+t158; + t701 = RATIONAL(73.0,4000.0)+t158; + t700 = y*t631+t470*z; + t699 = t429*t645+t405*z; + t698 = t244+t358; + t697 = RATIONAL(-23.0,500.0)+t157; + t696 = RATIONAL(49.0,1000.0)+t157; + t499 = RATIONAL(1.0,20.0); + t608 = t506*t486; + t262 = z*t608; + t695 = t499*t506+t262; + t153 = RATIONAL(27.0,1000.0)*t610; + t694 = t153+RATIONAL(-59.0,1000.0); + t228 = z*t339; + t693 = t368+t228; + t692 = t232+t359; + t691 = RATIONAL(7.0,250.0)+t159; + t344 = t477*t506; + t254 = z*t344; + t690 = t244+t254; + t419 = RATIONAL(-13.0,2000.0); + t629 = x*t419; + t400 = RATIONAL(323.0,4000.0); + t689 = t332+y*t629+t400*z; + t363 = t500*t501; + t641 = x*t481; + t688 = t363+y*t641+t392*z; + t401 = RATIONAL(53.0,12000.0); + t420 = RATIONAL(49.0,2000.0); + t687 = t420*t645+t332+t401*z; + t686 = t344+t400*x+t419*t648; + t475 = RATIONAL(-7.0,250.0); + t612 = y*t475; + t388 = RATIONAL(-137.0,12000.0); + t685 = t388*z+x*t612+t340; + t476 = RATIONAL(-4.0,375.0); + t489 = RATIONAL(-1.0,80.0); + t450 = RATIONAL(31.0,2000.0); + t627 = z*t450; + t169 = y*t627; + t684 = t476*x+t489*t506+t169; + t461 = RATIONAL(-17.0,500.0); + t459 = RATIONAL(-11.0,160.0); + t327 = t459*t506; + t398 = RATIONAL(-107.0,4000.0); + t683 = t398*x+t461*t648+t327; + t471 = RATIONAL(-1.0,2000.0); + t487 = RATIONAL(7.0,160.0); + t347 = t487*t501; + t395 = RATIONAL(-607.0,12000.0); + t682 = t347+t471*t645+t395*z; + t428 = RATIONAL(-49.0,2000.0); + t465 = RATIONAL(-7.0,2000.0); + t490 = RATIONAL(-3.0,40.0); + t681 = t490*t501+t465*t645+t428*z; + t348 = t487*t502; + t640 = z*t471; + t680 = t348+x*t640+t395*y; + t679 = z*t629+t400*y+t339; + t647 = x*z; + t678 = t388*y+t475*t647+t331; + t619 = x*t461; + t328 = t459*t501; + t677 = y*t619+t398*z+t328; + t611 = z*t465; + t676 = t428*y+x*t611+t490*t502; + t496 = RATIONAL(-1.0,16.0); + t354 = t496*t501; + t433 = RATIONAL(-87.0,2000.0); + t675 = x*t614+t433*z+t354; + t438 = RATIONAL(13.0,1000.0); + t439 = RATIONAL(-29.0,6000.0); + t497 = RATIONAL(1.0,40.0); + t353 = t497*t502; + t674 = t439*y+t438*t647+t353; + t415 = RATIONAL(183.0,4000.0); + t448 = RATIONAL(37.0,2000.0); + t673 = t348+t448*t647+t415*y; + t636 = z*t420; + t672 = t344+t401*x+y*t636; + t455 = RATIONAL(-79.0,2000.0); + t451 = RATIONAL(43.0,2000.0); + t638 = x*t451; + t492 = RATIONAL(-1.0,20.0); + t671 = t492*t502+t455*y+z*t638; + t166 = x*t627; + t670 = t166+t489*t502+t476*y; + t453 = RATIONAL(19.0,1000.0); + t458 = RATIONAL(53.0,1000.0); + t669 = t379+t458*z+t453*t645; + t668 = y*t638+t455*z+t492*t501; + t479 = RATIONAL(2.0,125.0); + t478 = RATIONAL(-9.0,160.0); + t336 = t478*t502; + t407 = RATIONAL(-137.0,4000.0); + t667 = t407*y+t479*t647+t336; + t445 = RATIONAL(-17.0,3000.0); + t436 = RATIONAL(-19.0,2000.0); + t183 = t436*t648; + t494 = RATIONAL(1.0,80.0); + t375 = t494*t506; + t666 = t445*x+t183+t375; + t337 = t487*t506; + t644 = y*t448; + t665 = t415*x+t337+z*t644; + t383 = t494*t501; + t642 = x*t436; + t184 = y*t642; + t664 = t445*z+t383+t184; + t364 = t495*t506; + t624 = z*t453; + t663 = t364+y*t624+t458*x; + t662 = t455*x+t492*t506+t451*t648; + t372 = t500*t502; + t483 = RATIONAL(-1.0,500.0); + t412 = RATIONAL(299.0,4000.0); + t661 = t483*t647+t372+t412*y; + t660 = x*t624+t376+t458*y; + t618 = y*t479; + t338 = t478*t501; + t659 = t407*z+x*t618+t338; + t325 = t459*t502; + t658 = z*t619+t398*y+t325; + t360 = t494*t502; + t190 = z*t642; + t657 = t360+t445*y+t190; + t656 = x*t636+t401*y+t339; + t381 = t497*t501; + t637 = y*t438; + t655 = t381+x*t637+t439*z; + t378 = t496*t502; + t654 = t464*t647+t378+t433*y; + t653 = t334+z*t612+t388*x; + t634 = y*t483; + t652 = t363+x*t634+t412*z; + t351 = t497*t506; + t651 = t439*x+t351+z*t637; + t621 = z*t467; + t650 = y*t621+t341+t371+t391*x; + t649 = t381+t351+RATIONAL(-61.0,4000.0)+t353; + t646 = z+y; + t425 = RATIONAL(-39.0,2000.0); + t643 = x*t425; + t446 = RATIONAL(27.0,2000.0); + t639 = y*t446; + t430 = RATIONAL(-23.0,2000.0); + t633 = x*t430; + t413 = RATIONAL(111.0,2000.0); + t632 = z*t413; + t485 = RATIONAL(-9.0,500.0); + t630 = x*t485; + t473 = RATIONAL(-3.0,2000.0); + t628 = x*t473; + t434 = RATIONAL(-21.0,2000.0); + t626 = y*t434; + t456 = RATIONAL(-27.0,1000.0); + t625 = z*t456; + t431 = RATIONAL(93.0,2000.0); + t623 = y*t431; + t474 = RATIONAL(-1.0,250.0); + t622 = y*t474; + t472 = RATIONAL(-1.0,1000.0); + t620 = x*t472; + t424 = RATIONAL(33.0,2000.0); + t616 = z*t424; + t403 = RATIONAL(129.0,2000.0); + t615 = x*t403; + t426 = RATIONAL(-17.0,2000.0); + t613 = z*t426; + t266 = y*t608; + t609 = t266+t254; + t386 = t500*t503; + t607 = t496*t506+t386; + t605 = t243+t268; + t604 = t272+t245; + t247 = z*t334; + t603 = t221+t247; + t226 = z*t341; + t602 = t226+t227; + t601 = t272+t228; + t257 = x*t346; + t249 = x*t331; + t600 = t257+t249; + t377 = t493*t503; + t382 = t500*t506; + t598 = t377+t382; + t596 = t221+t226; + t595 = t226+t223; + t259 = y*t344; + t594 = t269+t259; + t380 = t498*t503; + t593 = t380+t382; + t256 = x*t340; + t592 = t256+t249; + t591 = t228+t254; + t589 = t259+t262; + t588 = t242+t266; + t587 = t269+t266; + t585 = t232+t265; + t584 = t247+t251; + t583 = t247+t227; + t581 = t372+t357; + t263 = x*t343; + t580 = t256+t263; + t579 = t257+t263; + t208 = z*t620; + t191 = x*t625; + t578 = t259+t725; + t577 = t266+t725; + t576 = t247+t723; + t575 = t259+t722; + t574 = t226+t719; + t573 = t227+t693; + t156 = RATIONAL(1.0,1000.0)*t610; + t572 = t377+t156+t579; + t571 = t380+t349+t265+t254; + t570 = t368+t262+t245+t371; + t569 = t269+t367+t717; + t568 = t349+t228+t708; + t567 = t707+t715; + t566 = t380+t154+t590; + t155 = y*t208; + t565 = t155+t386+t592; + t564 = t377+t367+t249+t266; + t563 = t698+t705; + t562 = t247+t376+t719; + t561 = t698+t724; + t560 = t221+t379+t693; + t559 = t692+t705; + t558 = t272+t371+t710; + t557 = t349+t346+t262+t263; + t556 = t708+t717; + t152 = y*t191; + t555 = t152+t380+t606; + t554 = t692+t724; + t553 = t554+t602; + t552 = t561+t596; + t551 = t569+t609; + t550 = t563+t603; + t549 = t380+y*t635+t364+t435*x+t709; + t548 = t386+t600+t712; + t547 = t559+t583; + t546 = t559+t596; + t545 = t571+t605; + t544 = t377+t221+t580+t715; + t543 = t377+t373+t584+t600; + t542 = t570+t601; + t541 = t568+t589; + t540 = t259+t380+t604+t723; + t402 = RATIONAL(469.0,4000.0); + t539 = t402+t428*x+y*t611+t490*t506+t154+t386; + t538 = t349+t367+t591+t594; + t537 = t563+t602; + t536 = t347+t415*z+t362+x*t644+t697; + t192 = t450*t645; + t535 = t192+t349+t489*t501+t476*z+t702; + t534 = t362+t386+t580+t595; + t345 = t478*t506; + t533 = z*t618+t345+t386+t407*x+t720; + t532 = t569+t589; + t531 = t561+t583; + t530 = t367+t391*y+t343+x*t621+t691; + t529 = t357+t587+t727; + t528 = t259+t254+t556; + t527 = t380+t337+y*t640+t395*x+t721; + t437 = RATIONAL(79.0,3000.0); + t526 = t437*x+t472*t648+t711+t728; + t525 = t262+t437*z+y*t620+t707+t711; + t449 = RATIONAL(-59.0,4000.0); + t524 = t470*x+t152+t449+t469*t648+t704; + t523 = t437*y+t208+t367+t588+t711; + t399 = RATIONAL(191.0,4000.0); + t522 = t399*x+t485*t648+t592+t598+t716; + t521 = t363+y*t630+t247+t399*z+t716+t722; + t447 = RATIONAL(-33.0,2000.0); + t520 = y*t625+t447*x+t579+t607+t701; + t389 = RATIONAL(-691.0,12000.0); + t519 = t389*x+z*t622+t590+t593+t696; + t411 = RATIONAL(397.0,4000.0); + t518 = t411+t197+t433*x+t156+t592+t607; + t517 = t378+t447*y+t191+t227+t701+t715; + t516 = t377+y*t617+t154+t405*x+t649; + t515 = t226+t456*t645+t447*z+t354+t701+t725; + t514 = t392*y+z*t641+t581+t588+t694; + t513 = t389*y+t474*t647+t581+t594+t696; + t512 = t363+t389*z+t368+x*t622+t591+t696; + t511 = t244+t399*y+t372+z*t630+t710+t716; + t510 = t155+t420+z*t634+t412*x+t579+t598; + t480 = RATIONAL(-7.0,160.0); + t466 = RATIONAL(21.0,250.0); + t463 = RATIONAL(-17.0,375.0); + t462 = RATIONAL(13.0,500.0); + t460 = RATIONAL(-13.0,160.0); + t457 = RATIONAL(29.0,2000.0); + t452 = RATIONAL(-47.0,1500.0); + t443 = RATIONAL(-57.0,2000.0); + t442 = RATIONAL(-39.0,4000.0); + t417 = RATIONAL(147.0,2000.0); + t416 = RATIONAL(133.0,2000.0); + t414 = RATIONAL(307.0,6000.0); + t410 = RATIONAL(157.0,2000.0); + t409 = RATIONAL(-129.0,4000.0); + t408 = RATIONAL(-191.0,4000.0); + t404 = RATIONAL(-437.0,6000.0); + t396 = RATIONAL(-539.0,12000.0); + t394 = RATIONAL(319.0,12000.0); + t393 = RATIONAL(-449.0,12000.0); + t390 = RATIONAL(349.0,12000.0); + t385 = t499*t502; + t370 = t499*t501; + t178 = z*t626; + t177 = x*t626; + t161 = t434*t647; + coeff_I_m1_m1_m1 = RATIONAL(-269.0,4000.0)+t266+t452*z+t385+t370+(t417*z+ +t452)*y+(t452+t646*t417)*x+t555+t556+t695; + coeff_I_0_m1_m1 = t528+t442*x+t403*t648+t460*t506+t386+t579+t680+t682+ +t721; + coeff_I_p1_m1_m1 = t377+t528+t410*x+y*t632+t375+t592+t709+t713+t718+t726; + coeff_I_p2_m1_m1 = t344+t393*x+t514+z*t623+t570+t605+t688; + coeff_I_m1_0_m1 = z*t615+t527+t460*t502+t442*y+t243+t262+t358+t573+t585+ +t682; + coeff_I_0_0_m1 = t426*t645+t404*z+t539+t254+t560+t600+t676+t715; + coeff_I_p1_0_m1 = t512+t544+t658+t686; + coeff_I_p2_0_m1 = t517+t542+t655+t666; + coeff_I_m1_p1_m1 = t360+t549+t410*y+t262+x*t632+t560+t586+t718+t730; + coeff_I_0_p1_m1 = t512+t548+t679+t683; + coeff_I_p1_p1_m1 = t463*z+t516+t424*t645+t359+t573+t580+t690+t729; + coeff_I_p2_p1_m1 = t511+t542+t653+t687; + coeff_I_m1_p2_m1 = t393*y+t532+t431*t647+t339+t593+t606+t688+t694+t706; + coeff_I_0_p2_m1 = t520+t551+t655+t657; + coeff_I_p1_p2_m1 = t522+t551+t678+t687; + coeff_I_p2_p2_m1 = t524+t414*z+y*t628+t383+t532+t703; + coeff_I_m1_m1_0 = t527+y*t615+t442*z+t460*t501+t574+t577+t604+t680; + coeff_I_0_m1_0 = x*t613+t404*y+t539+t562+t578+t580+t681; + coeff_I_p1_m1_0 = t513+t543+t677+t686; + coeff_I_p2_m1_0 = t515+t529+t666+t674; + coeff_I_m1_0_0 = t402+y*t613+t404*x+t364+t537+t566+t676+t681; + coeff_I_0_0_0 = t550+RATIONAL(177.0,1000.0)+t325+t190+(x+t646)*RATIONAL( +-161.0,4000.0)+t183+t327+t328+t184+t565; + coeff_I_p1_0_0 = t550+t411+t178+t351+t466*x+t572+t654+t675; + coeff_I_p2_0_0 = t334+t430*t648+t390*x+t537+t659+t667+t720+t728; + coeff_I_m1_p1_0 = t519+t546+t677+t679; + coeff_I_0_p1_0 = t518+t353+t466*y+t161+t547+t675; + coeff_I_p1_p1_0 = t408*z+t192+t510+t338+t547+t661; + coeff_I_p2_p1_0 = t526+t546+t660+t668; + coeff_I_m1_p2_0 = t515+t540+t651+t657; + coeff_I_0_p2_0 = t331+t533+z*t633+t390*y+t576+t577+t580+t659; + coeff_I_p1_p2_0 = t523+t543+t663+t668; + coeff_I_p2_p2_0 = t242+t530+t177+t409*z+t226+t480*t501+t578+t586+t650; + coeff_I_m1_m1_p1 = t549+t410*z+t413*t645+t383+t266+t562+t604+t722+t726; + coeff_I_0_m1_p1 = t513+t534+t683+t689; + coeff_I_p1_m1_p1 = t463*y+t516+x*t616+t574+t575+t600+t699; + coeff_I_p2_m1_p1 = t521+t529+t653+t656; + coeff_I_m1_0_p1 = t519+t531+t658+t689; + coeff_I_0_0_p1 = t381+t177+t518+t466*z+t552+t654; + coeff_I_p1_0_p1 = t408*y+t166+t336+t510+t552+t652; + coeff_I_p2_0_p1 = t526+t531+t669+t671; + coeff_I_m1_p1_p1 = t554+t463*x+y*t616+t566+t603+t649+t699+t729; + coeff_I_0_p1_p1 = t345+t408*x+t553+t420+t169+t565+t652+t661; + coeff_I_p1_p1_p1 = RATIONAL(-193.0,4000.0)+t364+t416*z+(t416+t457*z)*y+( +t416+t646*t457)*x+t553+t572+t713; + coeff_I_p2_p1_p1 = t536+z*t639+t394*x+t341+t558+t584+t585+t673; + coeff_I_m1_p2_p1 = t521+t540+t672+t678; + coeff_I_0_p2_p1 = t523+t534+t662+t669; + coeff_I_p1_p2_p1 = t242+t257+t536+t446*t647+t394*y+t343+t564+t595+t665; + coeff_I_p2_p2_p1 = t370+y*t643+t462*z+t575+t576+t670+t684+t702+t727; + coeff_I_m1_m1_p2 = t332+t382+t393*z+t514+x*t623+t571+t601+t706; + coeff_I_0_m1_p2 = t520+t541+t664+t674; + coeff_I_p1_m1_p2 = t522+t541+t656+t685; + coeff_I_p2_m1_p2 = t360+z*t628+t524+t414*y+t568+t609+t700; + coeff_I_m1_0_p2 = t517+t545+t651+t664; + coeff_I_0_0_p2 = t221+t340+t262+t533+y*t633+t390*z+t567+t600+t667; + coeff_I_p1_0_p2 = t525+t544+t663+t671; + coeff_I_p2_0_p2 = t409*y+t480*t502+t346+t567+t227+t254+t161+t604+t650+ +t691+t714; + coeff_I_m1_p1_p2 = t511+t545+t672+t685; + coeff_I_0_p1_p2 = t525+t548+t660+t662; + coeff_I_p1_p1_p2 = t377+t557+t256+x*t639+t394*z+t268+t665+t673+t697+t712; + coeff_I_p2_p1_p2 = t462*y+t385+t245+t268+z*t643+t535+t558+t684+t690; + coeff_I_m1_p2_p2 = t538+t414*x+t449+t473*t648+t375+t555+t700+t703; + coeff_I_0_p2_p2 = t257+t530+t409*x+t480*t506+t178+t386+t228+t557+t587+ +t714; + coeff_I_p1_p2_p2 = t269+t256+t462*x+t228+t425*t648+t535+t564+t670+t695; + coeff_I_p2_p2_p2 = RATIONAL(77.0,2000.0)+t396*z+t153+(t396+t443*z)*y+( +t502+t501+t506)*t486+(t396+t646*t443)*x+t538+t704; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dx.compute.c new file mode 100644 index 0000000..7374de1 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dx.compute.c @@ -0,0 +1,507 @@ +fp t893, + t888, + t801, + t894, + t804, + t960, + t895, + t808, + t887, + t795, + t995, + t885, + t796, + t798, + t994, + t807, + t993, + t816, + t992, + t991, + t822, + t990, + t938, + t809, + t989, + t881, + t787, + t793, + t936, + t988, + t806, + t987, + t986, + t820, + t985, + t984, + t817, + t983, + t819, + t982, + t866, + t963, + t736, + t981, + t815, + t980, + t979, + t791, + t954, + t978, + t800, + t977, + t813, + t976, + t975, + t974, + t973, + t951, + t972, + t971, + t970, + t969, + t746, + t968, + t873, + t777, + t735, + t886, + t967, + t966, + t965, + t964, + t962, + t870, + t772, + t754, + t961, + t876, + t778, + t740, + t959, + t957, + t744, + t739, + t955, + t953, + t952, + t741, + t743, + t950, + t949, + t814, + t947, + t742, + t946, + t945, + t944, + t745, + t943, + t942, + t941, + t935, + t738, + t934, + t933, + t737, + t878, + t932, + t931, + t930, + t929, + t928, + t927, + t926, + t925, + t734, + t924, + t923, + t922, + t921, + t920, + t919, + t918, + t733, + t917, + t916, + t732, + t915, + t914, + t818, + t913, + t912, + t911, + t731, + t910, + t909, + t908, + t907, + t906, + t905, + t904, + t903, + t902, + t901, + t900, + t823, + t899, + t898, + t897, + t896, + t884, + t883, + t882, + t879, + t877, + t875, + t872, + t871, + t869, + t864, + t863, + t862, + t861, + t859, + t856, + t854, + t852, + t850, + t848, + t847, + t846, + t844, + t843, + t842, + t841, + t840, + t837, + t833, + t802, + t797, + t794, + t792, + t790, + t789, + t788, + t786, + t785, + t784, + t783, + t782, + t781, + t780, + t779, + t776, + t775, + t774, + t773, + t770, + t769, + t768, + t766, + t765, + t764, + t763, + t762, + t761, + t760, + t759, + t758, + t757, + t756, + t755, + t753, + t752, + t751, + t750, + t749, + t748, + t747; + t893 = z*z; + t888 = RATIONAL(-1.0,160.0); + t801 = t888*t893; + t894 = y*y; + t804 = t888*t894; + t960 = t801+t804; + t895 = x*x; + t808 = RATIONAL(1.0,32.0)*t895; + t887 = RATIONAL(-3.0,160.0); + t795 = t887*t893; + t995 = t808+t795; + t885 = RATIONAL(1.0,160.0); + t796 = t885*t894; + t798 = t885*t893; + t994 = t796+t798; + t807 = RATIONAL(-1.0,32.0)*t895; + t993 = t795+t807; + t816 = RATIONAL(1.0,20.0)*x; + t992 = t816+RATIONAL(141.0,2000.0); + t991 = RATIONAL(-79.0,2000.0)+RATIONAL(-1.0,10.0)*x; + t822 = RATIONAL(7.0,80.0)*x; + t990 = t822+RATIONAL(183.0,4000.0); + t938 = t801+t796; + t809 = RATIONAL(-3.0,32.0)*t895; + t989 = t938+t809; + t881 = RATIONAL(3.0,160.0); + t787 = t881*t893; + t793 = t887*t894; + t936 = t787+t793; + t988 = t936+t807; + t806 = RATIONAL(3.0,80.0)*x; + t987 = RATIONAL(323.0,4000.0)+t806; + t986 = t816+RATIONAL(-29.0,6000.0); + t820 = RATIONAL(-1.0,80.0)*x; + t985 = t820+RATIONAL(-167.0,12000.0); + t984 = t793+t995; + t817 = RATIONAL(3.0,40.0)*x; + t983 = RATIONAL(53.0,1000.0)+t817; + t819 = RATIONAL(3.0,32.0)*t895; + t982 = t819+t994; + t866 = RATIONAL(-27.0,1000.0); + t963 = y*z; + t736 = RATIONAL(9.0,1000.0)*t963; + t981 = t736+t866*z; + t815 = RATIONAL(1.0,80.0)*x; + t980 = t815+RATIONAL(-137.0,12000.0); + t979 = t819+t960; + t791 = t881*t894; + t954 = t787+t791; + t978 = t954+t807; + t800 = RATIONAL(-11.0,80.0)*x; + t977 = RATIONAL(-107.0,4000.0)+t800; + t813 = RATIONAL(1.0,40.0)*x; + t976 = t813+RATIONAL(-17.0,3000.0); + t975 = t817+RATIONAL(-97.0,3000.0); + t974 = t793+t993; + t973 = t822+RATIONAL(-607.0,12000.0); + t951 = t809+t798; + t972 = t951+t804; + t971 = t808+t954; + t970 = t791+t995; + t969 = t808+t936; + t746 = y*t815; + t968 = t938+t746; + t873 = RATIONAL(-9.0,1000.0); + t777 = t873*z; + t735 = y*t777; + t886 = RATIONAL(-9.0,500.0); + t967 = t735+t886*z; + t966 = RATIONAL(-1.0,40.0)*x+RATIONAL(-4.0,375.0); + t965 = RATIONAL(53.0,12000.0)+t806; + t964 = t791+t993; + t962 = z+y; + t870 = RATIONAL(-3.0,1000.0); + t772 = t870*z; + t754 = t866*y; + t961 = RATIONAL(-3.0,80.0)*x; + t876 = RATIONAL(-1.0,1000.0); + t778 = t876*z; + t740 = z*t820; + t959 = t735+t740; + t957 = t740+t746; + t744 = y*t806; + t739 = z*t815; + t955 = t744+t739; + t953 = t739+t746; + t952 = t804+t798; + t741 = y*t961; + t743 = z*t806; + t950 = t741+t743; + t949 = t739+t741; + t814 = RATIONAL(1.0,16.0)*x; + t947 = t809+t814; + t742 = z*t961; + t946 = t744+t742; + t945 = t741+t742; + t944 = t736+t739; + t745 = y*t820; + t943 = t739+t745; + t942 = t745+t740; + t941 = t744+t743; + t935 = RATIONAL(1.0,10.0)*x+t945; + t738 = y*t772; + t934 = t738+t742+t876*y; + t933 = t817+t942; + t737 = RATIONAL(3.0,1000.0)*t963; + t878 = RATIONAL(-1.0,250.0); + t932 = t737+t878*y+t743; + t931 = t816+t953; + t930 = t735+t739+t886*y; + t929 = t754+t736+t740; + t928 = t813+t941; + t927 = t742+t969; + t926 = t736+t796+t951; + t925 = t819+t745+t938; + t734 = y*t778; + t924 = t734+t982; + t923 = t738+t989; + t922 = t819+t740+t952; + t921 = t744+t970; + t920 = t746+t972; + t919 = t735+t979; + t918 = t737+t989; + t733 = RATIONAL(1.0,1000.0)*t963; + t917 = t733+t809+t960; + t916 = t741+t970; + t732 = RATIONAL(27.0,1000.0)*t963; + t915 = t732+t971; + t914 = t745+t972; + t818 = RATIONAL(-9.0,80.0)*x; + t913 = t819+t737+RATIONAL(-137.0,4000.0)+t818; + t912 = t737+t984; + t911 = t741+t988; + t731 = z*t754; + t910 = t731+t974; + t909 = t738+t978; + t908 = t745+t964; + t907 = t746+t964; + t906 = t743+t969; + t905 = t744+t988; + t904 = t819+RATIONAL(-3.0,20.0)*x+t738+RATIONAL(-49.0,2000.0); + t903 = t738+RATIONAL(79.0,3000.0)+t984; + t902 = t731+RATIONAL(-41.0,750.0)+t971; + t901 = t732+RATIONAL(281.0,12000.0)+t814+t974; + t900 = RATIONAL(-691.0,12000.0)+t737+t814+t978; + t823 = RATIONAL(-1.0,8.0)*x; + t899 = t823+t736+RATIONAL(-33.0,2000.0)+t979; + t898 = t733+t823+RATIONAL(-87.0,2000.0)+t982; + t897 = t735+RATIONAL(191.0,4000.0)+t947+t994; + t896 = t734+RATIONAL(299.0,4000.0)+t947+t960; + t884 = RATIONAL(-1.0,500.0); + t883 = RATIONAL(-9.0,250.0); + t882 = RATIONAL(2.0,125.0); + t879 = RATIONAL(-7.0,250.0); + t877 = RATIONAL(-3.0,2000.0); + t875 = RATIONAL(-1.0,2000.0); + t872 = RATIONAL(11.0,500.0); + t871 = RATIONAL(-7.0,2000.0); + t869 = RATIONAL(-17.0,500.0); + t864 = RATIONAL(19.0,1000.0); + t863 = RATIONAL(43.0,2000.0); + t862 = RATIONAL(31.0,2000.0); + t861 = RATIONAL(37.0,2000.0); + t859 = RATIONAL(27.0,2000.0); + t856 = RATIONAL(-37.0,1000.0); + t854 = RATIONAL(13.0,1000.0); + t852 = RATIONAL(-19.0,2000.0); + t850 = RATIONAL(-21.0,2000.0); + t848 = RATIONAL(93.0,2000.0); + t847 = RATIONAL(-23.0,2000.0); + t846 = RATIONAL(-31.0,1000.0); + t844 = RATIONAL(-17.0,2000.0); + t843 = RATIONAL(-39.0,2000.0); + t842 = RATIONAL(33.0,2000.0); + t841 = RATIONAL(49.0,2000.0); + t840 = RATIONAL(-13.0,2000.0); + t837 = RATIONAL(111.0,2000.0); + t833 = RATIONAL(129.0,2000.0); + t802 = t884*y; + t797 = t879*z; + t794 = t884*z; + t792 = t883*z; + t790 = t878*z; + t789 = t879*y; + t788 = t883*y; + t786 = t882*y; + t785 = t882*z; + t784 = t871*z; + t783 = t872*z; + t782 = t875*z; + t781 = t869*z; + t780 = t872*y; + t779 = t875*y; + t776 = t871*y; + t775 = t869*y; + t774 = t873*y; + t773 = t870*y; + t770 = t841*z; + t769 = t862*z; + t768 = t854*y; + t766 = t864*y; + t765 = t863*z; + t764 = t840*z; + t763 = t846*z; + t762 = t841*y; + t761 = t856*z; + t760 = t861*y; + t759 = t861*z; + t758 = t856*y; + t757 = t852*z; + t756 = t850*z; + t755 = t840*y; + t753 = t862*y; + t752 = t854*z; + t751 = t863*y; + t750 = t864*z; + t749 = t846*y; + t748 = t852*y; + t747 = t850*y; + coeff_dx_m1_m1_m1 = RATIONAL(-47.0,1500.0)+t962*RATIONAL(147.0,2000.0)+ +t910+t935; + coeff_dx_0_m1_m1 = t782+RATIONAL(-39.0,4000.0)+t779+RATIONAL(-13.0,80.0)* +x+t919+t941; + coeff_dx_p1_m1_m1 = RATIONAL(157.0,2000.0)+t761+t758+t926+t928; + coeff_dx_p2_m1_m1 = t792+RATIONAL(-449.0,12000.0)+t788+t806+t915+t945; + coeff_dx_m1_0_m1 = t833*z+t742+t779+t735+t908+t973; + coeff_dx_0_0_m1 = t784+t743+t844*y+t904+t968; + coeff_dx_p1_0_m1 = t781+t920+t932+t987; + coeff_dx_p2_0_m1 = t768+t745+t927+t976+t981; + coeff_dx_m1_p1_m1 = t837*z+t736+t758+t742+t907+t975; + coeff_dx_0_p1_m1 = t764+t925+t932+t977; + coeff_dx_p1_p1_m1 = t743+t738+t842*y+t763+t914+t992; + coeff_dx_p2_p1_m1 = t762+t746+t927+t967+t980; + coeff_dx_m1_p2_m1 = t788+t848*z+t901+t946; + coeff_dx_0_p2_m1 = t768+t757+t899+t950; + coeff_dx_p1_p2_m1 = t797+t762+t897+t950; + coeff_dx_p2_p2_m1 = t777+t877*y+t902+t946; + coeff_dx_m1_m1_0 = t833*y+t782+t911+t959+t973; + coeff_dx_0_m1_0 = t844*z+t776+t904+t952+t955; + coeff_dx_p1_m1_0 = t790+t775+t918+t955+t987; + coeff_dx_p2_m1_0 = t752+t916+t929+t976; + coeff_dx_m1_0_0 = RATIONAL(-437.0,6000.0)+t784+t776+t909+t933; + coeff_dx_0_0_0 = t800+t757+t748+RATIONAL(-161.0,4000.0)+t924+t953; + coeff_dx_p1_0_0 = RATIONAL(21.0,250.0)+t773+t772+t917+t931; + coeff_dx_p2_0_0 = t785+t815+t786+RATIONAL(349.0,12000.0)+t912+t942; + coeff_dx_m1_p1_0 = t775+t764+t900+t957; + coeff_dx_0_p1_0 = t773+t756+t898+t943; + coeff_dx_p1_p1_0 = t794+t753+t896+t943; + coeff_dx_p2_p1_0 = t750+t751+t903+t957; + coeff_dx_m1_p2_0 = t757+t905+t929+t986; + coeff_dx_0_p2_0 = t847*z+t786+t913+t949+t952; + coeff_dx_p1_p2_0 = t778+t751+t923+t949+t983; + coeff_dx_p2_p2_0 = t783+t747+t921+t959+t985; + coeff_dx_m1_m1_p1 = t761+t837*y+t911+t944+t975; + coeff_dx_0_m1_p1 = t737+t790+t755+t744+t922+t977; + coeff_dx_p1_m1_p1 = t842*z+t744+t749+t740+t923+t992; + coeff_dx_p2_m1_p1 = t770+t916+t930+t980; + coeff_dx_m1_0_p1 = t755+t781+t900+t943; + coeff_dx_0_0_p1 = t772+t747+t898+t957; + coeff_dx_p1_0_p1 = t802+t769+t896+t957; + coeff_dx_p2_0_p1 = t765+t766+t903+t943; + coeff_dx_m1_p1_p1 = t749+RATIONAL(-17.0,375.0)+t763+t909+t931; + coeff_dx_0_p1_p1 = RATIONAL(-191.0,4000.0)+t818+t794+t802+t924+t942; + coeff_dx_p1_p1_p1 = RATIONAL(133.0,2000.0)+t962*RATIONAL(29.0,2000.0)+ +t917+t933; + coeff_dx_p2_p1_p1 = t820+t759+RATIONAL(319.0,12000.0)+t760+t912+t953; + coeff_dx_m1_p2_p1 = t797+t905+t930+t965; + coeff_dx_0_p2_p1 = t741+t778+t738+t766+t922+t991; + coeff_dx_p1_p2_p1 = t740+t859*z+t741+t760+t918+t990; + coeff_dx_p2_p2_p1 = t769+t843*y+t921+t944+t966; + coeff_dx_m1_m1_p2 = t848*y+t792+t901+t950; + coeff_dx_0_m1_p2 = t752+t748+t899+t946; + coeff_dx_p1_m1_p2 = t770+t789+t897+t946; + coeff_dx_p2_m1_p2 = t774+t877*z+t902+t950; + coeff_dx_m1_0_p2 = t743+t748+t908+t981+t986; + coeff_dx_0_0_p2 = t847*y+t785+t742+t913+t968; + coeff_dx_p1_0_p2 = t765+t920+t934+t983; + coeff_dx_p2_0_p2 = t756+t745+t780+t735+t906+t985; + coeff_dx_m1_p1_p2 = t743+t789+t907+t965+t967; + coeff_dx_0_p1_p2 = t750+t925+t934+t991; + coeff_dx_p1_p1_p2 = t737+t759+t859*y+t742+t914+t990; + coeff_dx_p2_p1_p2 = t753+t746+t736+t843*z+t906+t966; + coeff_dx_m1_p2_p2 = t777+RATIONAL(307.0,6000.0)+t774+t910+t928; + coeff_dx_0_p2_p2 = t783+RATIONAL(-129.0,4000.0)+t780+RATIONAL(-7.0,80.0)* +x+t919+t945; + coeff_dx_p1_p2_p2 = RATIONAL(13.0,500.0)+t753+t769+t926+t935; + coeff_dx_p2_p2_p2 = RATIONAL(-539.0,12000.0)+t961+t962*RATIONAL(-57.0, +2000.0)+t915+t941; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dxx.compute.c new file mode 100644 index 0000000..7a88287 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dxx.compute.c @@ -0,0 +1,182 @@ +fp t1538, + t1531, + t1539, + t1534, + t1584, + t1536, + t1545, + t1530, + t1583, + t1582, + t1527, + t1533, + t1581, + t1528, + t1537, + t1580, + t1579, + t1578, + t1552, + t1577, + t1544, + t1535, + t1576, + t1525, + t1575, + t1574, + t1573, + t1572, + t1543, + t1526, + t1571, + t1570, + t1569, + t1568, + t1567, + t1548, + t1566, + t1565, + t1564, + t1563, + t1546, + t1562, + t1561, + t1560, + t1547, + t1559, + t1558, + t1557, + t1551, + t1556, + t1555, + t1554, + t1553, + t1550, + t1549, + t1542, + t1540, + t1532, + t1529; + t1538 = RATIONAL(-3.0,80.0); + t1531 = t1538*y; + t1539 = RATIONAL(-1.0,80.0); + t1534 = t1539*z; + t1584 = t1531+t1534; + t1536 = t1538*z; + t1545 = RATIONAL(3.0,80.0); + t1530 = t1545*y; + t1583 = t1536+t1530; + t1582 = t1536+t1531; + t1527 = t1539*y; + t1533 = t1545*z; + t1581 = t1527+t1533; + t1528 = RATIONAL(3.0,16.0)*x; + t1537 = RATIONAL(-11.0,80.0); + t1580 = t1528+t1537; + t1579 = t1528+t1534; + t1578 = t1530+t1534; + t1552 = RATIONAL(-1.0,8.0); + t1577 = t1528+t1552; + t1544 = RATIONAL(1.0,80.0); + t1535 = t1544*z; + t1576 = t1528+t1535; + t1525 = t1544*y; + t1575 = t1528+t1525; + t1574 = t1531+t1533; + t1573 = t1530+t1533; + t1572 = t1536+t1527; + t1543 = RATIONAL(1.0,16.0); + t1526 = t1543*x; + t1571 = t1527+t1526; + t1570 = t1526+t1535; + t1569 = t1525+t1526; + t1568 = t1543+t1525+t1534; + t1567 = t1527+t1543+t1535; + t1548 = RATIONAL(1.0,20.0); + t1566 = t1548+t1525+t1535; + t1565 = t1533+t1545+t1525; + t1564 = t1543+t1583; + t1563 = t1530+t1545+t1535; + t1546 = RATIONAL(1.0,40.0); + t1562 = t1546+t1573; + t1561 = t1543+t1574; + t1560 = t1548+t1581; + t1547 = RATIONAL(3.0,40.0); + t1559 = t1536+t1525+t1547; + t1558 = t1527+t1547+t1534; + t1557 = t1531+t1547+t1535; + t1551 = RATIONAL(7.0,80.0); + t1556 = t1551+t1572; + t1555 = t1551+t1584; + t1554 = RATIONAL(1.0,10.0)+t1582; + t1553 = t1548+t1578; + t1550 = RATIONAL(-3.0,20.0); + t1549 = RATIONAL(-9.0,80.0); + t1542 = RATIONAL(-1.0,10.0); + t1540 = RATIONAL(-1.0,40.0); + t1532 = RATIONAL(-1.0,16.0)*x; + t1529 = RATIONAL(-3.0,16.0)*x; + coeff_dxx_m1_m1_m1 = t1532+t1554; + coeff_dxx_0_m1_m1 = RATIONAL(-13.0,80.0)+t1528+t1573; + coeff_dxx_p1_m1_m1 = t1529+t1562; + coeff_dxx_p2_m1_m1 = t1545+t1526+t1582; + coeff_dxx_m1_0_m1 = t1532+t1556; + coeff_dxx_0_0_m1 = t1533+t1550+t1575; + coeff_dxx_p1_0_m1 = t1529+t1565; + coeff_dxx_p2_0_m1 = t1536+t1546+t1571; + coeff_dxx_m1_p1_m1 = t1532+t1559; + coeff_dxx_0_p1_m1 = t1580+t1581; + coeff_dxx_p1_p1_m1 = t1529+t1560; + coeff_dxx_p2_p1_m1 = t1536+t1544+t1569; + coeff_dxx_m1_p2_m1 = t1532+t1564; + coeff_dxx_0_p2_m1 = t1574+t1577; + coeff_dxx_p1_p2_m1 = t1529+t1561; + coeff_dxx_p2_p2_m1 = t1526+t1583; + coeff_dxx_m1_m1_0 = t1532+t1555; + coeff_dxx_0_m1_0 = t1530+t1550+t1576; + coeff_dxx_p1_m1_0 = t1529+t1563; + coeff_dxx_p2_m1_0 = t1546+t1526+t1584; + coeff_dxx_m1_0_0 = t1532+t1558; + coeff_dxx_0_0_0 = t1535+t1537+t1575; + coeff_dxx_p1_0_0 = t1529+t1566; + coeff_dxx_p2_0_0 = t1534+t1544+t1571; + coeff_dxx_m1_p1_0 = t1532+t1568; + coeff_dxx_0_p1_0 = t1552+t1527+t1576; + coeff_dxx_p1_p1_0 = t1529+t1567; + coeff_dxx_p2_p1_0 = t1534+t1569; + coeff_dxx_m1_p2_0 = t1532+t1553; + coeff_dxx_0_p2_0 = t1531+t1549+t1576; + coeff_dxx_p1_p2_0 = t1529+t1557; + coeff_dxx_p2_p2_0 = t1539+t1526+t1578; + coeff_dxx_m1_m1_p1 = t1532+t1557; + coeff_dxx_0_m1_p1 = t1578+t1580; + coeff_dxx_p1_m1_p1 = t1529+t1553; + coeff_dxx_p2_m1_p1 = t1544+t1531+t1570; + coeff_dxx_m1_0_p1 = t1532+t1567; + coeff_dxx_0_0_p1 = t1552+t1534+t1575; + coeff_dxx_p1_0_p1 = t1529+t1568; + coeff_dxx_p2_0_p1 = t1527+t1570; + coeff_dxx_m1_p1_p1 = t1532+t1566; + coeff_dxx_0_p1_p1 = t1527+t1549+t1579; + coeff_dxx_p1_p1_p1 = t1529+t1558; + coeff_dxx_p2_p1_p1 = t1539+t1535+t1569; + coeff_dxx_m1_p2_p1 = t1532+t1563; + coeff_dxx_0_p2_p1 = t1542+t1531+t1579; + coeff_dxx_p1_p2_p1 = t1529+t1555; + coeff_dxx_p2_p2_p1 = t1530+t1540+t1570; + coeff_dxx_m1_m1_p2 = t1532+t1561; + coeff_dxx_0_m1_p2 = t1577+t1583; + coeff_dxx_p1_m1_p2 = t1529+t1564; + coeff_dxx_p2_m1_p2 = t1526+t1574; + coeff_dxx_m1_0_p2 = t1532+t1560; + coeff_dxx_0_0_p2 = t1536+t1549+t1575; + coeff_dxx_p1_0_p2 = t1529+t1559; + coeff_dxx_p2_0_p2 = t1533+t1539+t1571; + coeff_dxx_m1_p1_p2 = t1532+t1565; + coeff_dxx_0_p1_p2 = t1542+t1528+t1572; + coeff_dxx_p1_p1_p2 = t1529+t1556; + coeff_dxx_p2_p1_p2 = t1533+t1540+t1569; + coeff_dxx_m1_p2_p2 = t1532+t1562; + coeff_dxx_0_p2_p2 = RATIONAL(-7.0,80.0)+t1528+t1582; + coeff_dxx_p1_p2_p2 = t1529+t1554; + coeff_dxx_p2_p2_p2 = t1538+t1526+t1573; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dxy.compute.c new file mode 100644 index 0000000..f470308 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dxy.compute.c @@ -0,0 +1,204 @@ +fp t1614, + t1588, + t1662, + t1626, + t1596, + t1628, + t1597, + t1661, + t1590, + t1660, + t1629, + t1594, + t1599, + t1659, + t1627, + t1598, + t1658, + t1657, + t1595, + t1656, + t1593, + t1655, + t1654, + t1653, + t1652, + t1651, + t1619, + t1591, + t1650, + t1649, + t1602, + t1585, + t1648, + t1586, + t1647, + t1600, + t1646, + t1645, + t1592, + t1644, + t1643, + t1642, + t1641, + t1640, + t1639, + t1638, + t1637, + t1636, + t1635, + t1634, + t1587, + t1633, + t1632, + t1615, + t1589, + t1631, + t1630, + t1625, + t1618, + t1616, + t1612, + t1611, + t1610, + t1609, + t1607, + t1606, + t1605, + t1604, + t1603, + t1601; + t1614 = RATIONAL(-9.0,1000.0); + t1588 = t1614*z; + t1662 = t1588+RATIONAL(-7.0,250.0); + t1626 = RATIONAL(-3.0,80.0); + t1596 = t1626*y; + t1628 = RATIONAL(1.0,80.0); + t1597 = t1628*x; + t1661 = t1596+t1597; + t1590 = RATIONAL(3.0,1000.0)*z; + t1660 = RATIONAL(-13.0,2000.0)+t1590; + t1629 = RATIONAL(3.0,80.0); + t1594 = t1629*x; + t1599 = t1629*y; + t1659 = t1594+t1599; + t1627 = RATIONAL(-1.0,80.0); + t1598 = t1627*x; + t1658 = t1598+t1599; + t1657 = RATIONAL(-9.0,500.0)+t1588; + t1595 = t1626*x; + t1656 = t1588+t1595; + t1593 = t1627*y; + t1655 = t1593+t1595; + t1654 = RATIONAL(2.0,125.0)+t1590; + t1653 = RATIONAL(-1.0,2000.0)+t1588; + t1652 = t1593+t1598; + t1651 = t1596+t1598; + t1619 = RATIONAL(-3.0,1000.0); + t1591 = t1619*z; + t1650 = t1591+t1598; + t1649 = t1595+t1599; + t1602 = RATIONAL(-27.0,1000.0); + t1585 = t1602*z; + t1648 = t1585+t1596; + t1586 = RATIONAL(27.0,1000.0)*z; + t1647 = t1586+t1596; + t1600 = t1628*y; + t1646 = t1594+t1600; + t1645 = t1593+t1597; + t1592 = RATIONAL(9.0,1000.0)*z; + t1644 = t1592+t1594; + t1643 = t1586+t1599; + t1642 = t1593+t1594; + t1641 = RATIONAL(11.0,500.0)+t1588; + t1640 = t1598+t1600; + t1639 = t1591+t1600; + t1638 = t1597+t1599; + t1637 = t1592+t1599; + t1636 = t1597+t1600; + t1635 = t1590+t1597; + t1634 = t1591+t1594; + t1587 = RATIONAL(1.0,1000.0)*z; + t1633 = t1587+t1598; + t1632 = t1590+RATIONAL(-1.0,250.0); + t1615 = RATIONAL(-1.0,1000.0); + t1589 = t1615*z; + t1631 = t1589+RATIONAL(-1.0,500.0); + t1630 = t1592+t1595; + t1625 = RATIONAL(-9.0,250.0); + t1618 = RATIONAL(-7.0,2000.0); + t1616 = RATIONAL(-17.0,500.0); + t1612 = RATIONAL(49.0,2000.0); + t1611 = RATIONAL(31.0,2000.0); + t1610 = RATIONAL(43.0,2000.0); + t1609 = RATIONAL(-21.0,2000.0); + t1607 = RATIONAL(-31.0,1000.0); + t1606 = RATIONAL(19.0,1000.0); + t1605 = RATIONAL(37.0,2000.0); + t1604 = RATIONAL(13.0,1000.0); + t1603 = RATIONAL(-19.0,2000.0); + t1601 = RATIONAL(-37.0,1000.0); + coeff_dxy_m1_m1_m1 = RATIONAL(147.0,2000.0)+t1595+t1648; + coeff_dxy_0_m1_m1 = t1642+t1653; + coeff_dxy_p1_m1_m1 = t1601+t1600+t1644; + coeff_dxy_p2_m1_m1 = t1595+t1625+t1643; + coeff_dxy_m1_0_m1 = t1653+t1658; + coeff_dxy_0_0_m1 = RATIONAL(-17.0,2000.0)+t1591+t1636; + coeff_dxy_p1_0_m1 = t1632+t1645; + coeff_dxy_p2_0_m1 = t1592+t1604+t1651; + coeff_dxy_m1_p1_m1 = t1601+t1597+t1637; + coeff_dxy_0_p1_m1 = t1632+t1640; + coeff_dxy_p1_p1_m1 = t1593+RATIONAL(33.0,2000.0)+t1650; + coeff_dxy_p2_p1_m1 = t1612+t1588+t1661; + coeff_dxy_m1_p2_m1 = t1594+t1625+t1647; + coeff_dxy_0_p2_m1 = t1593+t1604+t1630; + coeff_dxy_p1_p2_m1 = t1600+t1612+t1656; + coeff_dxy_p2_p2_m1 = RATIONAL(-3.0,2000.0)+t1585+t1659; + coeff_dxy_m1_m1_0 = RATIONAL(129.0,2000.0)+t1596+t1656; + coeff_dxy_0_m1_0 = t1618+t1593+t1634; + coeff_dxy_p1_m1_0 = t1590+t1616+t1646; + coeff_dxy_p2_m1_0 = t1599+t1602+t1630; + coeff_dxy_m1_0_0 = t1599+t1618+t1650; + coeff_dxy_0_0_0 = t1603+t1589+t1636; + coeff_dxy_p1_0_0 = t1587+t1619+t1645; + coeff_dxy_p2_0_0 = t1651+t1654; + coeff_dxy_m1_p1_0 = t1616+t1599+t1635; + coeff_dxy_0_p1_0 = t1600+t1619+t1633; + coeff_dxy_p1_p1_0 = t1611+t1589+t1652; + coeff_dxy_p2_p1_0 = t1610+t1591+t1661; + coeff_dxy_m1_p2_0 = t1596+t1602+t1644; + coeff_dxy_0_p2_0 = t1654+t1655; + coeff_dxy_p1_p2_0 = t1610+t1595+t1639; + coeff_dxy_p2_p2_0 = t1588+t1609+t1659; + coeff_dxy_m1_m1_p1 = t1596+RATIONAL(111.0,2000.0)+t1630; + coeff_dxy_0_m1_p1 = t1642+t1660; + coeff_dxy_p1_m1_p1 = t1607+t1600+t1634; + coeff_dxy_p2_m1_p1 = t1649+t1657; + coeff_dxy_m1_0_p1 = t1658+t1660; + coeff_dxy_0_0_p1 = t1609+t1587+t1636; + coeff_dxy_p1_0_p1 = t1631+t1645; + coeff_dxy_p2_0_p1 = t1596+t1606+t1650; + coeff_dxy_m1_p1_p1 = t1591+t1607+t1638; + coeff_dxy_0_p1_p1 = t1631+t1640; + coeff_dxy_p1_p1_p1 = t1593+RATIONAL(29.0,2000.0)+t1633; + coeff_dxy_p2_p1_p1 = t1596+t1605+t1635; + coeff_dxy_m1_p2_p1 = t1596+t1594+t1657; + coeff_dxy_0_p2_p1 = t1606+t1591+t1655; + coeff_dxy_p1_p2_p1 = t1600+t1595+t1590+t1605; + coeff_dxy_p2_p2_p1 = t1594+RATIONAL(-39.0,2000.0)+t1637; + coeff_dxy_m1_m1_p2 = t1595+RATIONAL(93.0,2000.0)+t1647; + coeff_dxy_0_m1_p2 = t1603+t1592+t1642; + coeff_dxy_p1_m1_p2 = t1646+t1662; + coeff_dxy_p2_m1_p2 = t1614+t1585+t1649; + coeff_dxy_m1_0_p2 = t1603+t1598+t1637; + coeff_dxy_0_0_p2 = RATIONAL(-23.0,2000.0)+t1600+t1635; + coeff_dxy_p1_0_p2 = t1615+t1591+t1645; + coeff_dxy_p2_0_p2 = t1641+t1651; + coeff_dxy_m1_p1_p2 = t1638+t1662; + coeff_dxy_0_p1_p2 = t1598+t1615+t1639; + coeff_dxy_p1_p1_p2 = t1590+RATIONAL(27.0,2000.0)+t1652; + coeff_dxy_p2_p1_p2 = t1611+t1592+t1661; + coeff_dxy_m1_p2_p2 = t1594+t1614+t1648; + coeff_dxy_0_p2_p2 = t1641+t1655; + coeff_dxy_p1_p2_p2 = t1600+t1611+t1630; + coeff_dxy_p2_p2_p2 = t1594+RATIONAL(-57.0,2000.0)+t1643; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dxz.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dxz.compute.c new file mode 100644 index 0000000..0da29f3 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dxz.compute.c @@ -0,0 +1,212 @@ +fp t1706, + t1672, + t1675, + t1748, + t1697, + t1667, + t1747, + t1707, + t1674, + t1676, + t1746, + t1704, + t1678, + t1745, + t1668, + t1744, + t1677, + t1743, + t1705, + t1671, + t1742, + t1693, + t1741, + t1673, + t1740, + t1739, + t1666, + t1738, + t1737, + t1669, + t1736, + t1665, + t1735, + t1692, + t1670, + t1734, + t1733, + t1732, + t1731, + t1730, + t1729, + t1728, + t1727, + t1663, + t1726, + t1725, + t1680, + t1664, + t1724, + t1723, + t1722, + t1721, + t1720, + t1719, + t1718, + t1717, + t1716, + t1715, + t1681, + t1714, + t1713, + t1712, + t1711, + t1710, + t1689, + t1709, + t1708, + t1702, + t1701, + t1700, + t1699, + t1698, + t1690, + t1687; + t1706 = RATIONAL(1.0,80.0); + t1672 = t1706*x; + t1675 = t1706*z; + t1748 = t1672+t1675; + t1697 = RATIONAL(-3.0,1000.0); + t1667 = t1697*y; + t1747 = t1667+RATIONAL(43.0,2000.0); + t1707 = RATIONAL(3.0,80.0); + t1674 = t1707*x; + t1676 = t1707*z; + t1746 = t1674+t1676; + t1704 = RATIONAL(-3.0,80.0); + t1678 = t1704*x; + t1745 = t1675+t1678; + t1668 = RATIONAL(3.0,1000.0)*y; + t1744 = t1668+RATIONAL(-17.0,500.0); + t1677 = t1704*z; + t1743 = t1674+t1677; + t1705 = RATIONAL(-1.0,80.0); + t1671 = t1705*x; + t1742 = t1671+t1675; + t1693 = RATIONAL(-1.0,1000.0); + t1741 = t1667+t1693; + t1673 = t1705*z; + t1740 = t1672+t1673; + t1739 = t1671+t1673; + t1666 = t1693*y; + t1738 = t1666+t1673; + t1737 = t1675+t1674; + t1669 = RATIONAL(9.0,1000.0)*y; + t1736 = t1669+RATIONAL(-37.0,1000.0); + t1665 = RATIONAL(1.0,1000.0)*y; + t1735 = t1697+t1665; + t1692 = RATIONAL(-9.0,1000.0); + t1670 = t1692*y; + t1734 = t1670+t1674; + t1733 = t1676+t1678; + t1732 = t1667+RATIONAL(-7.0,2000.0); + t1731 = t1670+RATIONAL(-1.0,2000.0); + t1730 = t1672+t1676; + t1729 = t1667+RATIONAL(-31.0,1000.0); + t1728 = t1670+RATIONAL(11.0,500.0); + t1727 = t1677+t1678; + t1663 = RATIONAL(27.0,1000.0)*y; + t1726 = t1663+RATIONAL(-9.0,250.0); + t1725 = t1668+RATIONAL(-13.0,2000.0); + t1680 = RATIONAL(-27.0,1000.0); + t1664 = t1680*y; + t1724 = t1664+t1677; + t1723 = t1667+RATIONAL(19.0,1000.0); + t1722 = t1670+t1678; + t1721 = t1669+RATIONAL(13.0,1000.0); + t1720 = t1673+t1678; + t1719 = t1671+t1677; + t1718 = t1668+t1673; + t1717 = t1668+t1671; + t1716 = t1671+t1676; + t1715 = t1672+t1677; + t1681 = RATIONAL(-19.0,2000.0); + t1714 = t1669+t1681; + t1713 = t1668+RATIONAL(37.0,2000.0); + t1712 = t1670+t1672; + t1711 = t1666+t1675; + t1710 = t1673+t1674; + t1689 = RATIONAL(31.0,2000.0); + t1709 = t1669+t1689; + t1708 = t1669+t1680; + t1702 = RATIONAL(2.0,125.0); + t1701 = RATIONAL(-9.0,500.0); + t1700 = RATIONAL(-1.0,250.0); + t1699 = RATIONAL(-1.0,500.0); + t1698 = RATIONAL(-7.0,250.0); + t1690 = RATIONAL(49.0,2000.0); + t1687 = RATIONAL(-21.0,2000.0); + coeff_dxz_m1_m1_m1 = t1678+RATIONAL(147.0,2000.0)+t1724; + coeff_dxz_0_m1_m1 = t1710+t1731; + coeff_dxz_p1_m1_m1 = t1736+t1737; + coeff_dxz_p2_m1_m1 = t1726+t1733; + coeff_dxz_m1_0_m1 = RATIONAL(129.0,2000.0)+t1677+t1722; + coeff_dxz_0_0_m1 = t1710+t1732; + coeff_dxz_p1_0_m1 = t1737+t1744; + coeff_dxz_p2_0_m1 = t1708+t1733; + coeff_dxz_m1_p1_m1 = RATIONAL(111.0,2000.0)+t1669+t1727; + coeff_dxz_0_p1_m1 = t1710+t1725; + coeff_dxz_p1_p1_m1 = t1729+t1737; + coeff_dxz_p2_p1_m1 = t1676+t1701+t1722; + coeff_dxz_m1_p2_m1 = RATIONAL(93.0,2000.0)+t1663+t1727; + coeff_dxz_0_p2_m1 = t1710+t1714; + coeff_dxz_p1_p2_m1 = t1675+t1698+t1734; + coeff_dxz_p2_p2_m1 = t1664+t1692+t1733; + coeff_dxz_m1_m1_0 = t1716+t1731; + coeff_dxz_0_m1_0 = t1667+RATIONAL(-17.0,2000.0)+t1748; + coeff_dxz_p1_m1_0 = t1672+t1700+t1718; + coeff_dxz_p2_m1_0 = t1719+t1721; + coeff_dxz_m1_0_0 = t1716+t1732; + coeff_dxz_0_0_0 = t1672+t1681+t1711; + coeff_dxz_p1_0_0 = t1735+t1740; + coeff_dxz_p2_0_0 = t1677+t1702+t1717; + coeff_dxz_m1_p1_0 = t1716+t1725; + coeff_dxz_0_p1_0 = t1665+t1687+t1748; + coeff_dxz_p1_p1_0 = t1699+t1672+t1738; + coeff_dxz_p2_p1_0 = t1719+t1723; + coeff_dxz_m1_p2_0 = t1714+t1716; + coeff_dxz_0_p2_0 = RATIONAL(-23.0,2000.0)+t1668+t1748; + coeff_dxz_p1_p2_0 = t1740+t1741; + coeff_dxz_p2_p2_0 = t1719+t1728; + coeff_dxz_m1_m1_p1 = t1730+t1736; + coeff_dxz_0_m1_p1 = t1675+t1700+t1717; + coeff_dxz_p1_m1_p1 = RATIONAL(33.0,2000.0)+t1667+t1739; + coeff_dxz_p2_m1_p1 = t1690+t1677+t1712; + coeff_dxz_m1_0_p1 = t1730+t1744; + coeff_dxz_0_0_p1 = t1735+t1742; + coeff_dxz_p1_0_p1 = t1671+t1689+t1738; + coeff_dxz_p2_0_p1 = t1715+t1747; + coeff_dxz_m1_p1_p1 = t1729+t1730; + coeff_dxz_0_p1_p1 = t1699+t1671+t1711; + coeff_dxz_p1_p1_p1 = t1665+RATIONAL(29.0,2000.0)+t1739; + coeff_dxz_p2_p1_p1 = t1713+t1715; + coeff_dxz_m1_p2_p1 = t1698+t1676+t1712; + coeff_dxz_0_p2_p1 = t1741+t1742; + coeff_dxz_p1_p2_p1 = RATIONAL(27.0,2000.0)+t1673+t1717; + coeff_dxz_p2_p2_p1 = t1709+t1715; + coeff_dxz_m1_m1_p2 = t1726+t1743; + coeff_dxz_0_m1_p2 = t1720+t1721; + coeff_dxz_p1_m1_p2 = t1690+t1675+t1722; + coeff_dxz_p2_m1_p2 = RATIONAL(-3.0,2000.0)+t1664+t1746; + coeff_dxz_m1_0_p2 = t1708+t1743; + coeff_dxz_0_0_p2 = t1702+t1678+t1718; + coeff_dxz_p1_0_p2 = t1745+t1747; + coeff_dxz_p2_0_p2 = t1687+t1676+t1734; + coeff_dxz_m1_p1_p2 = t1677+t1701+t1734; + coeff_dxz_0_p1_p2 = t1720+t1723; + coeff_dxz_p1_p1_p2 = t1713+t1745; + coeff_dxz_p2_p1_p2 = t1669+RATIONAL(-39.0,2000.0)+t1746; + coeff_dxz_m1_p2_p2 = t1692+t1674+t1724; + coeff_dxz_0_p2_p2 = t1720+t1728; + coeff_dxz_p1_p2_p2 = t1709+t1745; + coeff_dxz_p2_p2_p2 = RATIONAL(-57.0,2000.0)+t1663+t1746; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dy.compute.c new file mode 100644 index 0000000..ab2c1ef --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dy.compute.c @@ -0,0 +1,519 @@ +fp t1153, + t1158, + t1054, + t1160, + t1055, + t1212, + t1146, + t1060, + t1159, + t1078, + t1264, + t1076, + t1263, + t1150, + t1069, + t1087, + t1206, + t1053, + t1262, + t1152, + t1062, + t1261, + t1260, + t1061, + t1201, + t1259, + t1258, + t1057, + t1199, + t1257, + t1075, + t1256, + t1255, + t1254, + t1253, + t1252, + t1251, + t1083, + t1250, + t1086, + t1249, + t1135, + t1044, + t1001, + t1248, + t1082, + t1247, + t1080, + t1246, + t1245, + t1244, + t1077, + t1243, + t1081, + t1242, + t1238, + t1003, + t1241, + t1240, + t1239, + t1237, + t1131, + t1034, + t1138, + t1049, + t1141, + t1043, + t1236, + t1007, + t998, + t1235, + t999, + t1234, + t1006, + t1233, + t1004, + t1232, + t1231, + t1230, + t1009, + t1229, + t1073, + t1228, + t1008, + t1227, + t1225, + t1223, + t1222, + t1010, + t1221, + t1220, + t1219, + t1217, + t1216, + t1215, + t1011, + t1214, + t1211, + t1005, + t1210, + t1208, + t1115, + t1207, + t1204, + t1203, + t1202, + t1198, + t1151, + t1197, + t1196, + t1195, + t1194, + t1193, + t1192, + t1191, + t997, + t1190, + t996, + t1189, + t1188, + t1187, + t1186, + t1002, + t1185, + t1184, + t1085, + t1183, + t1182, + t1181, + t1180, + t1179, + t1178, + t1177, + t1176, + t1175, + t1174, + t1173, + t1000, + t1172, + t1171, + t1051, + t1170, + t1169, + t1168, + t1167, + t1166, + t1074, + t1165, + t1164, + t1163, + t1162, + t1161, + t1149, + t1148, + t1147, + t1144, + t1143, + t1142, + t1140, + t1137, + t1136, + t1134, + t1129, + t1128, + t1127, + t1126, + t1124, + t1121, + t1119, + t1117, + t1113, + t1112, + t1111, + t1109, + t1108, + t1107, + t1106, + t1105, + t1102, + t1098, + t1070, + t1068, + t1067, + t1065, + t1064, + t1063, + t1059, + t1058, + t1056, + t1052, + t1050, + t1048, + t1047, + t1046, + t1045, + t1042, + t1041, + t1040, + t1039, + t1038, + t1037, + t1036, + t1035, + t1033, + t1032, + t1030, + t1029, + t1028, + t1027, + t1026, + t1025, + t1024, + t1023, + t1022, + t1021, + t1020, + t1019, + t1018, + t1017, + t1016, + t1015, + t1014, + t1013, + t1012; + t1153 = RATIONAL(-1.0,160.0); + t1158 = z*z; + t1054 = t1153*t1158; + t1160 = x*x; + t1055 = t1153*t1160; + t1212 = t1054+t1055; + t1146 = RATIONAL(3.0,160.0); + t1060 = t1146*t1160; + t1159 = y*y; + t1078 = RATIONAL(-1.0,32.0)*t1159; + t1264 = t1060+t1078; + t1076 = RATIONAL(1.0,32.0)*t1159; + t1263 = t1060+t1076; + t1150 = RATIONAL(1.0,160.0); + t1069 = t1150*t1158; + t1087 = RATIONAL(-3.0,32.0)*t1159; + t1206 = t1069+t1087; + t1053 = t1150*t1160; + t1262 = t1206+t1053; + t1152 = RATIONAL(-3.0,160.0); + t1062 = t1152*t1158; + t1261 = t1062+t1264; + t1260 = RATIONAL(-4.0,375.0)+RATIONAL(-1.0,40.0)*y; + t1061 = t1152*t1160; + t1201 = t1061+t1062; + t1259 = t1076+t1201; + t1258 = t1206+t1055; + t1057 = t1146*t1158; + t1199 = t1061+t1057; + t1257 = t1199+t1078; + t1075 = RATIONAL(3.0,40.0)*y; + t1256 = RATIONAL(-97.0,3000.0)+t1075; + t1255 = t1078+t1201; + t1254 = t1057+t1264; + t1253 = t1087+t1212; + t1252 = t1053+t1087+t1054; + t1251 = t1062+t1263; + t1083 = RATIONAL(3.0,80.0)*y; + t1250 = RATIONAL(53.0,12000.0)+t1083; + t1086 = RATIONAL(1.0,40.0)*y; + t1249 = RATIONAL(-17.0,3000.0)+t1086; + t1135 = RATIONAL(-3.0,1000.0); + t1044 = t1135*z; + t1001 = x*t1044; + t1248 = RATIONAL(53.0,1000.0)+t1075+t1001; + t1082 = RATIONAL(1.0,20.0)*y; + t1247 = t1082+RATIONAL(-29.0,6000.0); + t1080 = RATIONAL(-1.0,80.0)*y; + t1246 = t1080+RATIONAL(-167.0,12000.0); + t1245 = t1057+t1263; + t1244 = t1199+t1076; + t1077 = RATIONAL(7.0,80.0)*y; + t1243 = t1077+RATIONAL(-607.0,12000.0); + t1081 = RATIONAL(1.0,80.0)*y; + t1242 = RATIONAL(-137.0,12000.0)+t1081; + t1238 = x*z; + t1003 = RATIONAL(3.0,1000.0)*t1238; + t1241 = t1003+RATIONAL(183.0,4000.0)+t1077; + t1240 = RATIONAL(141.0,2000.0)+t1001+t1082; + t1239 = RATIONAL(323.0,4000.0)+t1003+t1083; + t1237 = x+z; + t1131 = RATIONAL(-27.0,1000.0); + t1034 = t1131*z; + t1138 = RATIONAL(-9.0,1000.0); + t1049 = t1138*x; + t1141 = RATIONAL(-1.0,1000.0); + t1043 = t1141*z; + t1236 = RATIONAL(-3.0,80.0)*y; + t1007 = x*t1083; + t998 = z*t1049; + t1235 = t1007+t998; + t999 = RATIONAL(9.0,1000.0)*t1238; + t1234 = t1007+t999; + t1006 = x*t1081; + t1233 = t1006+t999; + t1004 = x*t1236; + t1232 = t1004+t999; + t1231 = t1004+t998; + t1230 = t1069+t1007; + t1009 = z*t1081; + t1229 = t1009+t1055; + t1073 = RATIONAL(3.0,32.0)*t1159; + t1228 = t1073+t1001; + t1008 = z*t1083; + t1227 = t1008+t1007; + t1225 = t1008+t1053; + t1223 = t1009+t1006; + t1222 = t1006+t1054; + t1010 = z*t1236; + t1221 = t1010+t1053; + t1220 = t1010+t1004; + t1219 = t1003+t1073; + t1217 = t1069+t1004; + t1216 = t1069+t1053; + t1215 = t1073+RATIONAL(-1.0,8.0)*y; + t1011 = z*t1080; + t1214 = t1011+t1006; + t1211 = t1010+t1007; + t1005 = x*t1080; + t1210 = t1009+t1005; + t1208 = t1011+t1055; + t1115 = RATIONAL(-21.0,2000.0); + t1207 = t1115*z+t1005; + t1204 = t1054+t1005; + t1203 = t1008+t1004; + t1202 = t1011+t1005; + t1198 = t1034+t1005+t999; + t1151 = RATIONAL(-9.0,500.0); + t1197 = t1006+t1151*z+t998; + t1196 = t1086+t1227; + t1195 = RATIONAL(1.0,10.0)*y+t1220; + t1194 = t1075+t1202; + t1193 = t1082+t1223; + t1192 = t999+t1262; + t1191 = t1073+t998+t1212; + t997 = RATIONAL(27.0,1000.0)*t1238; + t1190 = t997+t1245; + t996 = x*t1034; + t1189 = t996+t1255; + t1188 = RATIONAL(-3.0,20.0)*y+RATIONAL(-49.0,2000.0)+t1228; + t1187 = t1010+t1261; + t1186 = t1009+t1252; + t1002 = x*t1043; + t1185 = t1002+t1073+t1216; + t1184 = t1008+t1261; + t1085 = RATIONAL(-9.0,80.0)*y; + t1183 = RATIONAL(-137.0,4000.0)+t1085+t1219; + t1182 = t1011+t1252; + t1181 = t1001+t1254; + t1180 = t1008+t1244; + t1179 = t1011+t1251; + t1178 = t1009+t1251; + t1177 = t1003+t1259; + t1176 = t1010+t1244; + t1175 = t1011+t1257; + t1174 = t1006+t1258; + t1173 = RATIONAL(-79.0,2000.0)+RATIONAL(-1.0,10.0)*y+t1228; + t1000 = RATIONAL(1.0,1000.0)*t1238; + t1172 = t1000+t1253; + t1171 = t1009+t1257; + t1051 = RATIONAL(-11.0,80.0)*y; + t1170 = RATIONAL(-107.0,4000.0)+t1051+t1219; + t1169 = t1005+t1258; + t1168 = RATIONAL(-41.0,750.0)+t996+t1245; + t1167 = t1001+RATIONAL(79.0,3000.0)+t1259; + t1166 = RATIONAL(-33.0,2000.0)+t999+t1212+t1215; + t1074 = RATIONAL(1.0,16.0)*y; + t1165 = t1074+t997+RATIONAL(281.0,12000.0)+t1255; + t1164 = t1074+RATIONAL(191.0,4000.0)+t998+t1262; + t1163 = t1002+RATIONAL(299.0,4000.0)+t1074+t1253; + t1162 = RATIONAL(-87.0,2000.0)+t1000+t1215+t1216; + t1161 = t1003+t1074+RATIONAL(-691.0,12000.0)+t1254; + t1149 = RATIONAL(-1.0,500.0); + t1148 = RATIONAL(-9.0,250.0); + t1147 = RATIONAL(2.0,125.0); + t1144 = RATIONAL(-7.0,250.0); + t1143 = RATIONAL(-1.0,250.0); + t1142 = RATIONAL(-3.0,2000.0); + t1140 = RATIONAL(-1.0,2000.0); + t1137 = RATIONAL(11.0,500.0); + t1136 = RATIONAL(-7.0,2000.0); + t1134 = RATIONAL(-17.0,500.0); + t1129 = RATIONAL(19.0,1000.0); + t1128 = RATIONAL(43.0,2000.0); + t1127 = RATIONAL(31.0,2000.0); + t1126 = RATIONAL(37.0,2000.0); + t1124 = RATIONAL(27.0,2000.0); + t1121 = RATIONAL(-37.0,1000.0); + t1119 = RATIONAL(13.0,1000.0); + t1117 = RATIONAL(-19.0,2000.0); + t1113 = RATIONAL(93.0,2000.0); + t1112 = RATIONAL(-23.0,2000.0); + t1111 = RATIONAL(-31.0,1000.0); + t1109 = RATIONAL(-17.0,2000.0); + t1108 = RATIONAL(-39.0,2000.0); + t1107 = RATIONAL(33.0,2000.0); + t1106 = RATIONAL(49.0,2000.0); + t1105 = RATIONAL(-13.0,2000.0); + t1102 = RATIONAL(111.0,2000.0); + t1098 = RATIONAL(129.0,2000.0); + t1070 = t1144*x; + t1068 = t1147*x; + t1067 = t1143*z; + t1065 = t1149*z; + t1064 = t1147*z; + t1063 = t1144*z; + t1059 = t1151*x; + t1058 = t1148*x; + t1056 = t1148*z; + t1052 = t1149*x; + t1050 = t1143*x; + t1048 = t1135*x; + t1047 = t1141*x; + t1046 = t1137*x; + t1045 = t1134*z; + t1042 = t1136*x; + t1041 = t1137*z; + t1040 = t1140*x; + t1039 = t1134*x; + t1038 = t1138*z; + t1037 = t1140*z; + t1036 = t1136*z; + t1035 = t1126*x; + t1033 = t1121*z; + t1032 = t1117*x; + t1030 = t1127*x; + t1029 = t1127*z; + t1028 = t1106*z; + t1027 = t1128*x; + t1026 = t1119*z; + t1025 = t1129*z; + t1024 = t1111*z; + t1023 = t1115*x; + t1022 = t1105*x; + t1021 = t1106*x; + t1020 = t1126*z; + t1019 = t1121*x; + t1018 = t1111*x; + t1017 = t1105*z; + t1016 = t1131*x; + t1015 = t1117*z; + t1014 = t1128*z; + t1013 = t1119*x; + t1012 = t1129*x; + coeff_dy_m1_m1_m1 = RATIONAL(-47.0,1500.0)+t1237*RATIONAL(147.0,2000.0)+ +t1189+t1195; + coeff_dy_0_m1_m1 = t1040+t1098*z+t1005+t998+t1187+t1243; + coeff_dy_p1_m1_m1 = t1019+t1102*z+t1187+t1233+t1256; + coeff_dy_p2_m1_m1 = t1113*z+t1058+t1165+t1211; + coeff_dy_m1_0_m1 = t1040+RATIONAL(-39.0,4000.0)+RATIONAL(-13.0,80.0)*y+ +t1037+t1191+t1227; + coeff_dy_0_0_m1 = t1109*x+t1036+t1188+t1222+t1225; + coeff_dy_p1_0_m1 = t1017+t1050+t1170+t1204+t1225; + coeff_dy_p2_0_m1 = t1015+t1013+t1166+t1203; + coeff_dy_m1_p1_m1 = t1019+t1033+RATIONAL(157.0,2000.0)+t1192+t1196; + coeff_dy_0_p1_m1 = t1045+t1008+t1050+t1174+t1239; + coeff_dy_p1_p1_m1 = t1107*x+t1008+t1024+t1169+t1240; + coeff_dy_p2_p1_m1 = t1063+t1021+t1164+t1203; + coeff_dy_m1_p2_m1 = RATIONAL(-449.0,12000.0)+t1056+t1058+t1083+t1190+ +t1220; + coeff_dy_0_p2_m1 = t1013+t1176+t1198+t1249; + coeff_dy_p1_p2_m1 = t1021+t1176+t1197+t1242; + coeff_dy_p2_p2_m1 = t1038+t1142*x+t1168+t1211; + coeff_dy_m1_m1_0 = t1037+t1098*x+t1175+t1231+t1243; + coeff_dy_0_m1_0 = RATIONAL(-437.0,6000.0)+t1036+t1042+t1181+t1194; + coeff_dy_p1_m1_0 = t1039+t1017+t1161+t1214; + coeff_dy_p2_m1_0 = t1015+t1016+t1175+t1234+t1247; + coeff_dy_m1_0_0 = t1109*z+t1042+t1188+t1229+t1230; + coeff_dy_0_0_0 = t1051+t1032+RATIONAL(-161.0,4000.0)+t1015+t1185+t1223; + coeff_dy_p1_0_0 = t1048+t1009+t1162+t1207; + coeff_dy_p2_0_0 = t1068+t1112*z+t1183+t1217+t1229; + coeff_dy_m1_p1_0 = t1039+t1007+t1067+t1186+t1239; + coeff_dy_0_p1_0 = RATIONAL(21.0,250.0)+t1048+t1044+t1172+t1193; + coeff_dy_p1_p1_0 = t1065+t1030+t1163+t1210; + coeff_dy_p2_p1_0 = t1043+t1004+t1027+t1186+t1248; + coeff_dy_m1_p2_0 = t1016+t1026+t1179+t1232+t1249; + coeff_dy_0_p2_0 = t1068+t1064+t1081+RATIONAL(349.0,12000.0)+t1177+t1202; + coeff_dy_p1_p2_0 = t1027+t1025+t1167+t1214; + coeff_dy_p2_p2_0 = t1041+t1023+t1179+t1235+t1246; + coeff_dy_m1_m1_p1 = t1102*x+t1033+t1171+t1232+t1256; + coeff_dy_0_m1_p1 = t1022+t1045+t1161+t1210; + coeff_dy_p1_m1_p1 = t1018+RATIONAL(-17.0,375.0)+t1024+t1181+t1193; + coeff_dy_p2_m1_p1 = t1059+t1063+t1171+t1235+t1250; + coeff_dy_m1_0_p1 = t1022+t1067+t1170+t1208+t1230; + coeff_dy_0_0_p1 = t1044+t1023+t1162+t1214; + coeff_dy_p1_0_p1 = RATIONAL(-191.0,4000.0)+t1052+t1085+t1065+t1185+t1202; + coeff_dy_p2_0_p1 = t1043+t1012+t1173+t1208+t1217; + coeff_dy_m1_p1_p1 = t1007+t1107*z+t1018+t1182+t1240; + coeff_dy_0_p1_p1 = t1052+t1029+t1163+t1214; + coeff_dy_p1_p1_p1 = RATIONAL(133.0,2000.0)+t1237*RATIONAL(29.0,2000.0)+ +t1172+t1194; + coeff_dy_p2_p1_p1 = t1035+t1124*z+t1004+t1182+t1241; + coeff_dy_m1_p2_p1 = t1059+t1028+t1178+t1231+t1242; + coeff_dy_0_p2_p1 = t1012+t1014+t1167+t1210; + coeff_dy_p1_p2_p1 = RATIONAL(319.0,12000.0)+t1035+t1020+t1080+t1177+t1223 +; + coeff_dy_p2_p2_p1 = t1108*x+t1029+t1178+t1234+t1260; + coeff_dy_m1_m1_p2 = t1113*x+t1056+t1165+t1203; + coeff_dy_0_m1_p2 = t1032+t1184+t1198+t1247; + coeff_dy_p1_m1_p2 = t1070+t1184+t1197+t1250; + coeff_dy_p2_m1_p2 = t1038+t1049+RATIONAL(307.0,6000.0)+t1189+t1196; + coeff_dy_m1_0_p2 = t1026+t1032+t1166+t1211; + coeff_dy_0_0_p2 = t1064+t1112*x+t1183+t1221+t1222; + coeff_dy_p1_0_p2 = t1047+t1025+t1173+t1204+t1221; + coeff_dy_p2_0_p2 = t1046+RATIONAL(-7.0,80.0)*y+t1041+RATIONAL(-129.0, +4000.0)+t1191+t1220; + coeff_dy_m1_p1_p2 = t1070+t1028+t1164+t1211; + coeff_dy_0_p1_p2 = t1047+t1010+t1014+t1174+t1248; + coeff_dy_p1_p1_p2 = t1010+t1124*x+t1020+t1169+t1241; + coeff_dy_p2_p1_p2 = RATIONAL(13.0,500.0)+t1030+t1029+t1192+t1195; + coeff_dy_m1_p2_p2 = t1049+t1142*z+t1168+t1203; + coeff_dy_0_p2_p2 = t1046+t998+t1180+t1207+t1246; + coeff_dy_p1_p2_p2 = t1030+t1108*z+t1180+t1233+t1260; + coeff_dy_p2_p2_p2 = RATIONAL(-539.0,12000.0)+t1236+t1237*RATIONAL(-57.0, +2000.0)+t1190+t1227; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dyy.compute.c new file mode 100644 index 0000000..3b9d0c2 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dyy.compute.c @@ -0,0 +1,186 @@ +fp t1775, + t1751, + t1772, + t1756, + t1812, + t1774, + t1755, + t1758, + t1811, + t1765, + t1749, + t1810, + t1773, + t1753, + t1809, + t1750, + t1808, + t1807, + t1806, + t1805, + t1804, + t1766, + t1754, + t1757, + t1803, + t1802, + t1801, + t1800, + t1799, + t1752, + t1798, + t1797, + t1796, + t1795, + t1794, + t1793, + t1792, + t1791, + t1768, + t1790, + t1789, + t1767, + t1788, + t1787, + t1786, + t1785, + t1784, + t1783, + t1782, + t1781, + t1780, + t1770, + t1779, + t1769, + t1778, + t1777, + t1776, + t1763, + t1762, + t1760, + t1759; + t1775 = RATIONAL(3.0,80.0); + t1751 = t1775*z; + t1772 = RATIONAL(-3.0,80.0); + t1756 = t1772*x; + t1812 = t1751+t1756; + t1774 = RATIONAL(1.0,80.0); + t1755 = t1774*x; + t1758 = t1772*z; + t1811 = t1755+t1758; + t1765 = RATIONAL(-9.0,80.0); + t1749 = RATIONAL(3.0,16.0)*y; + t1810 = t1765+t1749; + t1773 = RATIONAL(-1.0,80.0); + t1753 = t1773*z; + t1809 = t1753+t1756; + t1750 = t1774*z; + t1808 = t1750+t1756; + t1807 = t1756+t1758; + t1806 = RATIONAL(-11.0,80.0)+t1749; + t1805 = t1750+t1755; + t1804 = RATIONAL(-3.0,20.0)+t1749; + t1766 = RATIONAL(1.0,16.0); + t1754 = t1766*y; + t1757 = t1773*x; + t1803 = t1754+t1757; + t1802 = t1754+t1758; + t1801 = t1749+t1757; + t1800 = t1751+t1757; + t1799 = t1753+t1749; + t1752 = t1775*x; + t1798 = t1752+t1749; + t1797 = t1753+t1752; + t1796 = t1754+t1755; + t1795 = t1751+t1755; + t1794 = t1751+t1754; + t1793 = t1752+t1750; + t1792 = t1766+t1812; + t1791 = t1752+t1766+t1758; + t1768 = RATIONAL(1.0,20.0); + t1790 = t1768+t1805; + t1789 = t1768+t1797; + t1767 = RATIONAL(3.0,40.0); + t1788 = t1767+t1811; + t1787 = t1750+t1766+t1757; + t1786 = t1775+t1795; + t1785 = t1767+t1808; + t1784 = t1775+t1793; + t1783 = t1753+t1767+t1757; + t1782 = t1768+t1800; + t1781 = t1753+t1755+t1766; + t1780 = RATIONAL(1.0,10.0)+t1807; + t1770 = RATIONAL(1.0,40.0); + t1779 = t1752+t1751+t1770; + t1769 = RATIONAL(7.0,80.0); + t1778 = t1769+t1757+t1758; + t1777 = t1769+t1809; + t1776 = RATIONAL(-1.0,8.0); + t1763 = RATIONAL(-1.0,10.0); + t1762 = RATIONAL(-1.0,40.0); + t1760 = RATIONAL(-1.0,16.0)*y; + t1759 = RATIONAL(-3.0,16.0)*y; + coeff_dyy_m1_m1_m1 = t1760+t1780; + coeff_dyy_0_m1_m1 = t1760+t1778; + coeff_dyy_p1_m1_m1 = t1760+t1788; + coeff_dyy_p2_m1_m1 = t1760+t1791; + coeff_dyy_m1_0_m1 = RATIONAL(-13.0,80.0)+t1751+t1798; + coeff_dyy_0_0_m1 = t1795+t1804; + coeff_dyy_p1_0_m1 = t1800+t1806; + coeff_dyy_p2_0_m1 = t1749+t1776+t1812; + coeff_dyy_m1_p1_m1 = t1759+t1779; + coeff_dyy_0_p1_m1 = t1759+t1786; + coeff_dyy_p1_p1_m1 = t1759+t1782; + coeff_dyy_p2_p1_m1 = t1759+t1792; + coeff_dyy_m1_p2_m1 = t1756+t1775+t1802; + coeff_dyy_0_p2_m1 = t1757+t1770+t1802; + coeff_dyy_p1_p2_m1 = t1758+t1774+t1796; + coeff_dyy_p2_p2_m1 = t1752+t1802; + coeff_dyy_m1_m1_0 = t1760+t1777; + coeff_dyy_0_m1_0 = t1760+t1783; + coeff_dyy_p1_m1_0 = t1760+t1781; + coeff_dyy_p2_m1_0 = t1760+t1789; + coeff_dyy_m1_0_0 = t1793+t1804; + coeff_dyy_0_0_0 = t1805+t1806; + coeff_dyy_p1_0_0 = t1750+t1776+t1801; + coeff_dyy_p2_0_0 = t1808+t1810; + coeff_dyy_m1_p1_0 = t1759+t1784; + coeff_dyy_0_p1_0 = t1759+t1790; + coeff_dyy_p1_p1_0 = t1759+t1787; + coeff_dyy_p2_p1_0 = t1759+t1785; + coeff_dyy_m1_p2_0 = t1770+t1754+t1809; + coeff_dyy_0_p2_0 = t1753+t1774+t1803; + coeff_dyy_p1_p2_0 = t1753+t1796; + coeff_dyy_p2_p2_0 = t1773+t1754+t1797; + coeff_dyy_m1_m1_p1 = t1760+t1785; + coeff_dyy_0_m1_p1 = t1760+t1787; + coeff_dyy_p1_m1_p1 = t1760+t1790; + coeff_dyy_p2_m1_p1 = t1760+t1784; + coeff_dyy_m1_0_p1 = t1797+t1806; + coeff_dyy_0_0_p1 = t1776+t1755+t1799; + coeff_dyy_p1_0_p1 = t1757+t1765+t1799; + coeff_dyy_p2_0_p1 = t1763+t1756+t1799; + coeff_dyy_m1_p1_p1 = t1759+t1789; + coeff_dyy_0_p1_p1 = t1759+t1781; + coeff_dyy_p1_p1_p1 = t1759+t1783; + coeff_dyy_p2_p1_p1 = t1759+t1777; + coeff_dyy_m1_p2_p1 = t1774+t1754+t1808; + coeff_dyy_0_p2_p1 = t1750+t1803; + coeff_dyy_p1_p2_p1 = t1750+t1773+t1796; + coeff_dyy_p2_p2_p1 = t1754+t1762+t1793; + coeff_dyy_m1_m1_p2 = t1760+t1792; + coeff_dyy_0_m1_p2 = t1760+t1782; + coeff_dyy_p1_m1_p2 = t1760+t1786; + coeff_dyy_p2_m1_p2 = t1760+t1779; + coeff_dyy_m1_0_p2 = t1776+t1758+t1798; + coeff_dyy_0_0_p2 = t1810+t1811; + coeff_dyy_p1_0_p2 = t1758+t1763+t1801; + coeff_dyy_p2_0_p2 = RATIONAL(-7.0,80.0)+t1749+t1807; + coeff_dyy_m1_p1_p2 = t1759+t1791; + coeff_dyy_0_p1_p2 = t1759+t1788; + coeff_dyy_p1_p1_p2 = t1759+t1778; + coeff_dyy_p2_p1_p2 = t1759+t1780; + coeff_dyy_m1_p2_p2 = t1756+t1794; + coeff_dyy_0_p2_p2 = t1757+t1773+t1794; + coeff_dyy_p1_p2_p2 = t1762+t1755+t1794; + coeff_dyy_p2_p2_p2 = t1772+t1752+t1794; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dyz.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dyz.compute.c new file mode 100644 index 0000000..089bfc4 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dyz.compute.c @@ -0,0 +1,208 @@ +fp t1847, + t1819, + t1891, + t1857, + t1823, + t1854, + t1827, + t1890, + t1845, + t1817, + t1889, + t1828, + t1888, + t1826, + t1887, + t1836, + t1814, + t1886, + t1856, + t1821, + t1855, + t1825, + t1885, + t1822, + t1824, + t1884, + t1842, + t1883, + t1816, + t1882, + t1881, + t1813, + t1880, + t1815, + t1879, + t1878, + t1877, + t1876, + t1875, + t1874, + t1818, + t1873, + t1872, + t1871, + t1820, + t1870, + t1869, + t1868, + t1867, + t1866, + t1865, + t1864, + t1863, + t1862, + t1861, + t1860, + t1859, + t1858, + t1853, + t1850, + t1849, + t1846, + t1844, + t1843, + t1841, + t1840, + t1837, + t1835, + t1833, + t1832, + t1831, + t1829; + t1847 = RATIONAL(-9.0,1000.0); + t1819 = t1847*x; + t1891 = t1819+RATIONAL(49.0,2000.0); + t1857 = RATIONAL(3.0,80.0); + t1823 = t1857*z; + t1854 = RATIONAL(-3.0,80.0); + t1827 = t1854*y; + t1890 = t1823+t1827; + t1845 = RATIONAL(-3.0,1000.0); + t1817 = t1845*x; + t1889 = RATIONAL(43.0,2000.0)+t1817; + t1828 = t1854*z; + t1888 = t1828+t1827; + t1826 = t1857*y; + t1887 = t1826+t1828; + t1836 = RATIONAL(-27.0,1000.0); + t1814 = t1836*x; + t1886 = t1814+t1827; + t1856 = RATIONAL(1.0,80.0); + t1821 = t1856*y; + t1855 = RATIONAL(-1.0,80.0); + t1825 = t1855*z; + t1885 = t1821+t1825; + t1822 = t1855*y; + t1824 = t1856*z; + t1884 = t1822+t1824; + t1842 = RATIONAL(-1.0,1000.0); + t1883 = t1842+t1817; + t1816 = t1842*x; + t1882 = RATIONAL(-1.0,500.0)+t1816; + t1881 = t1823+t1826; + t1813 = RATIONAL(27.0,1000.0)*x; + t1880 = t1813+t1828; + t1815 = RATIONAL(1.0,1000.0)*x; + t1879 = t1845+t1815; + t1878 = t1819+RATIONAL(-9.0,500.0); + t1877 = t1819+RATIONAL(-7.0,250.0); + t1876 = t1823+t1813; + t1875 = t1814+t1826; + t1874 = RATIONAL(19.0,1000.0)+t1817; + t1818 = RATIONAL(3.0,1000.0)*x; + t1873 = t1818+t1821; + t1872 = t1819+t1822; + t1871 = t1825+t1827; + t1820 = RATIONAL(9.0,1000.0)*x; + t1870 = t1820+t1822; + t1869 = t1824+t1827; + t1868 = t1821+t1828; + t1867 = t1820+t1836; + t1866 = t1818+t1828; + t1865 = t1818+RATIONAL(-13.0,2000.0); + t1864 = t1826+t1825; + t1863 = t1822+t1825; + t1862 = t1820+t1821; + t1861 = t1824+t1826; + t1860 = t1820+t1824; + t1859 = t1821+t1824; + t1858 = t1823+t1817; + t1853 = RATIONAL(-9.0,250.0); + t1850 = RATIONAL(-1.0,250.0); + t1849 = RATIONAL(2.0,125.0); + t1846 = RATIONAL(-1.0,2000.0); + t1844 = RATIONAL(-7.0,2000.0); + t1843 = RATIONAL(-17.0,500.0); + t1841 = RATIONAL(11.0,500.0); + t1840 = RATIONAL(31.0,2000.0); + t1837 = RATIONAL(37.0,2000.0); + t1835 = RATIONAL(-19.0,2000.0); + t1833 = RATIONAL(-37.0,1000.0); + t1832 = RATIONAL(-31.0,1000.0); + t1831 = RATIONAL(-21.0,2000.0); + t1829 = RATIONAL(13.0,1000.0); + coeff_dyz_m1_m1_m1 = RATIONAL(147.0,2000.0)+t1828+t1886; + coeff_dyz_0_m1_m1 = t1819+RATIONAL(129.0,2000.0)+t1888; + coeff_dyz_p1_m1_m1 = t1820+RATIONAL(111.0,2000.0)+t1888; + coeff_dyz_p2_m1_m1 = t1827+RATIONAL(93.0,2000.0)+t1880; + coeff_dyz_m1_0_m1 = t1819+t1846+t1864; + coeff_dyz_0_0_m1 = t1817+t1844+t1864; + coeff_dyz_p1_0_m1 = t1864+t1865; + coeff_dyz_p2_0_m1 = t1835+t1820+t1864; + coeff_dyz_m1_p1_m1 = t1826+t1833+t1860; + coeff_dyz_0_p1_m1 = t1818+t1843+t1861; + coeff_dyz_p1_p1_m1 = t1832+t1817+t1861; + coeff_dyz_p2_p1_m1 = t1861+t1877; + coeff_dyz_m1_p2_m1 = t1853+t1827+t1876; + coeff_dyz_0_p2_m1 = t1867+t1890; + coeff_dyz_p1_p2_m1 = t1878+t1890; + coeff_dyz_p2_p2_m1 = t1847+t1823+t1886; + coeff_dyz_m1_m1_0 = t1823+t1846+t1872; + coeff_dyz_0_m1_0 = t1844+t1822+t1858; + coeff_dyz_p1_m1_0 = t1822+t1823+t1865; + coeff_dyz_p2_m1_0 = t1835+t1823+t1870; + coeff_dyz_m1_0_0 = t1817+RATIONAL(-17.0,2000.0)+t1859; + coeff_dyz_0_0_0 = t1816+t1835+t1859; + coeff_dyz_p1_0_0 = t1831+t1815+t1859; + coeff_dyz_p2_0_0 = t1818+RATIONAL(-23.0,2000.0)+t1859; + coeff_dyz_m1_p1_0 = t1825+t1850+t1873; + coeff_dyz_0_p1_0 = t1879+t1885; + coeff_dyz_p1_p1_0 = t1882+t1885; + coeff_dyz_p2_p1_0 = t1883+t1885; + coeff_dyz_m1_p2_0 = t1829+t1828+t1870; + coeff_dyz_0_p2_0 = t1822+t1849+t1866; + coeff_dyz_p1_p2_0 = t1828+t1822+t1874; + coeff_dyz_p2_p2_0 = t1841+t1828+t1872; + coeff_dyz_m1_m1_p1 = t1833+t1823+t1862; + coeff_dyz_0_m1_p1 = t1843+t1823+t1873; + coeff_dyz_p1_m1_p1 = t1832+t1821+t1858; + coeff_dyz_p2_m1_p1 = t1821+t1823+t1877; + coeff_dyz_m1_0_p1 = t1818+t1850+t1884; + coeff_dyz_0_0_p1 = t1879+t1884; + coeff_dyz_p1_0_p1 = t1882+t1884; + coeff_dyz_p2_0_p1 = t1883+t1884; + coeff_dyz_m1_p1_p1 = t1817+RATIONAL(33.0,2000.0)+t1863; + coeff_dyz_0_p1_p1 = t1840+t1816+t1863; + coeff_dyz_p1_p1_p1 = t1815+RATIONAL(29.0,2000.0)+t1863; + coeff_dyz_p2_p1_p1 = t1818+RATIONAL(27.0,2000.0)+t1863; + coeff_dyz_m1_p2_p1 = t1868+t1891; + coeff_dyz_0_p2_p1 = t1868+t1889; + coeff_dyz_p1_p2_p1 = t1821+t1837+t1866; + coeff_dyz_p2_p2_p1 = t1828+t1840+t1862; + coeff_dyz_m1_m1_p2 = t1853+t1826+t1880; + coeff_dyz_0_m1_p2 = t1867+t1887; + coeff_dyz_p1_m1_p2 = t1878+t1887; + coeff_dyz_p2_m1_p2 = t1828+t1847+t1875; + coeff_dyz_m1_0_p2 = t1820+t1829+t1871; + coeff_dyz_0_0_p2 = t1818+t1849+t1871; + coeff_dyz_p1_0_p2 = t1871+t1874; + coeff_dyz_p2_0_p2 = t1819+t1841+t1871; + coeff_dyz_m1_p1_p2 = t1869+t1891; + coeff_dyz_0_p1_p2 = t1869+t1889; + coeff_dyz_p1_p1_p2 = t1818+t1837+t1869; + coeff_dyz_p2_p1_p2 = t1840+t1827+t1860; + coeff_dyz_m1_p2_p2 = RATIONAL(-3.0,2000.0)+t1823+t1875; + coeff_dyz_0_p2_p2 = t1831+t1819+t1881; + coeff_dyz_p1_p2_p2 = t1820+RATIONAL(-39.0,2000.0)+t1881; + coeff_dyz_p2_p2_p2 = t1826+RATIONAL(-57.0,2000.0)+t1876; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dz.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dz.compute.c new file mode 100644 index 0000000..e78181a --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dz.compute.c @@ -0,0 +1,507 @@ +fp t1415, + t1428, + t1329, + t1421, + t1429, + t1333, + t1497, + t1422, + t1327, + t1419, + t1325, + t1524, + t1332, + t1427, + t1343, + t1523, + t1345, + t1328, + t1479, + t1522, + t1342, + t1273, + t1400, + t1292, + t1521, + t1331, + t1472, + t1520, + t1519, + t1518, + t1344, + t1517, + t1516, + t1348, + t1515, + t1354, + t1466, + t1514, + t1469, + t1324, + t1513, + t1512, + t1511, + t1420, + t1346, + t1274, + t1510, + t1352, + t1509, + t1508, + t1340, + t1280, + t1507, + t1335, + t1506, + t1498, + t1275, + t1505, + t1349, + t1500, + t1270, + t1504, + t1407, + t1313, + t1269, + t1503, + t1502, + t1501, + t1499, + t1410, + t1316, + t1404, + t1305, + t1268, + t1496, + t1278, + t1495, + t1276, + t1494, + t1277, + t1493, + t1492, + t1267, + t1491, + t1490, + t1489, + t1488, + t1279, + t1487, + t1347, + t1350, + t1484, + t1483, + t1482, + t1481, + t1480, + t1478, + t1477, + t1473, + t1470, + t1468, + t1467, + t1465, + t1464, + t1463, + t1462, + t1461, + t1460, + t1265, + t1459, + t1458, + t1457, + t1456, + t1271, + t1455, + t1454, + t1453, + t1452, + t1272, + t1451, + t1450, + t1449, + t1448, + t1266, + t1447, + t1446, + t1445, + t1444, + t1443, + t1351, + t1442, + t1441, + t1440, + t1439, + t1438, + t1437, + t1436, + t1357, + t1435, + t1434, + t1433, + t1432, + t1431, + t1430, + t1418, + t1417, + t1416, + t1413, + t1412, + t1411, + t1409, + t1406, + t1405, + t1403, + t1398, + t1397, + t1396, + t1395, + t1393, + t1390, + t1388, + t1386, + t1384, + t1382, + t1381, + t1380, + t1378, + t1377, + t1376, + t1375, + t1374, + t1371, + t1367, + t1339, + t1338, + t1337, + t1336, + t1334, + t1330, + t1326, + t1323, + t1322, + t1320, + t1319, + t1318, + t1317, + t1315, + t1314, + t1312, + t1311, + t1310, + t1309, + t1308, + t1307, + t1306, + t1304, + t1303, + t1302, + t1301, + t1300, + t1299, + t1298, + t1297, + t1296, + t1295, + t1294, + t1293, + t1291, + t1290, + t1289, + t1288, + t1287, + t1286, + t1285, + t1284, + t1283, + t1282, + t1281; + t1415 = RATIONAL(3.0,160.0); + t1428 = y*y; + t1329 = t1415*t1428; + t1421 = RATIONAL(-3.0,160.0); + t1429 = x*x; + t1333 = t1421*t1429; + t1497 = t1329+t1333; + t1422 = RATIONAL(-1.0,160.0); + t1327 = t1422*t1428; + t1419 = RATIONAL(1.0,160.0); + t1325 = t1419*t1429; + t1524 = t1327+t1325; + t1332 = t1421*t1428; + t1427 = z*z; + t1343 = RATIONAL(1.0,32.0)*t1427; + t1523 = t1332+t1343; + t1345 = RATIONAL(3.0,32.0)*t1427; + t1328 = t1422*t1429; + t1479 = t1327+t1328; + t1522 = t1345+t1479; + t1342 = RATIONAL(-1.0,80.0)*z; + t1273 = y*t1342; + t1400 = RATIONAL(-27.0,1000.0); + t1292 = t1400*x; + t1521 = t1273+t1292; + t1331 = t1415*t1429; + t1472 = t1329+t1331; + t1520 = t1343+t1472; + t1519 = t1343+t1497; + t1518 = t1331+t1523; + t1344 = RATIONAL(1.0,20.0)*z; + t1517 = t1344+RATIONAL(-29.0,6000.0); + t1516 = RATIONAL(-49.0,2000.0)+RATIONAL(-3.0,20.0)*z; + t1348 = RATIONAL(3.0,40.0)*z; + t1515 = RATIONAL(-97.0,3000.0)+t1348; + t1354 = RATIONAL(-1.0,32.0)*t1427; + t1466 = t1354+t1333; + t1514 = t1332+t1466; + t1469 = t1345+t1325; + t1324 = t1419*t1428; + t1513 = t1469+t1324; + t1512 = t1348+RATIONAL(53.0,1000.0); + t1511 = t1472+t1354; + t1420 = RATIONAL(-9.0,500.0); + t1346 = RATIONAL(1.0,80.0)*z; + t1274 = y*t1346; + t1510 = t1420*x+t1274; + t1352 = RATIONAL(7.0,80.0)*z; + t1509 = t1352+RATIONAL(183.0,4000.0); + t1508 = t1333+t1523; + t1340 = RATIONAL(3.0,80.0)*z; + t1280 = x*t1340; + t1507 = t1280+t1524; + t1335 = RATIONAL(-11.0,80.0)*z; + t1506 = t1335+RATIONAL(-107.0,4000.0); + t1498 = RATIONAL(-3.0,80.0)*z; + t1275 = y*t1498; + t1505 = t1275+t1524; + t1349 = RATIONAL(1.0,40.0)*z; + t1500 = x*y; + t1270 = RATIONAL(9.0,1000.0)*t1500; + t1504 = RATIONAL(-17.0,3000.0)+t1349+t1270; + t1407 = RATIONAL(-9.0,1000.0); + t1313 = t1407*y; + t1269 = x*t1313; + t1503 = t1346+RATIONAL(-137.0,12000.0)+t1269; + t1502 = RATIONAL(-4.0,375.0)+RATIONAL(-1.0,40.0)*z+t1270; + t1501 = t1342+t1269+RATIONAL(-167.0,12000.0); + t1499 = x+y; + t1410 = RATIONAL(-1.0,1000.0); + t1316 = t1410*x; + t1404 = RATIONAL(-3.0,1000.0); + t1305 = t1404*y; + t1268 = x*t1305; + t1496 = t1268+t1274; + t1278 = x*t1346; + t1495 = t1278+t1274; + t1276 = x*t1498; + t1494 = t1276+t1274; + t1277 = x*t1342; + t1493 = t1277+t1274; + t1492 = t1328+t1324; + t1267 = RATIONAL(3.0,1000.0)*t1500; + t1491 = t1277+t1267; + t1490 = t1276+t1273; + t1489 = t1273+t1267; + t1488 = t1324+t1325; + t1279 = y*t1340; + t1487 = t1280+t1279; + t1347 = RATIONAL(1.0,16.0)*z; + t1350 = RATIONAL(-3.0,32.0)*t1427; + t1484 = t1347+t1350; + t1483 = t1277+t1279; + t1482 = t1332+t1331; + t1481 = t1278+t1268; + t1480 = t1275+t1276; + t1478 = t1275+t1280; + t1477 = t1277+t1273; + t1473 = t1276+t1279; + t1470 = t1278+t1279; + t1468 = t1354+t1269; + t1467 = t1273+t1278; + t1465 = RATIONAL(1.0,10.0)*z+t1480; + t1464 = t1344+t1495; + t1463 = t1348+t1477; + t1462 = t1349+t1487; + t1461 = t1329+t1270+t1466; + t1460 = t1267+t1508; + t1265 = y*t1292; + t1459 = t1265+t1514; + t1458 = t1268+t1511; + t1457 = t1275+t1350+t1492; + t1456 = t1280+t1345+t1492; + t1271 = RATIONAL(1.0,1000.0)*t1500; + t1455 = t1350+t1271+t1479; + t1454 = t1276+t1350+t1524; + t1453 = t1269+t1522; + t1452 = t1350+t1270+t1488; + t1272 = y*t1316; + t1451 = t1272+t1513; + t1450 = t1276+t1518; + t1449 = t1354+t1270+t1482; + t1448 = t1277+t1519; + t1266 = RATIONAL(27.0,1000.0)*t1500; + t1447 = t1266+t1520; + t1446 = t1344+t1350+t1268+RATIONAL(141.0,2000.0); + t1445 = t1278+t1519; + t1444 = t1340+t1350+t1267+RATIONAL(323.0,4000.0); + t1443 = t1340+RATIONAL(53.0,12000.0)+t1468; + t1351 = RATIONAL(-9.0,80.0)*z; + t1442 = t1345+t1351+RATIONAL(-137.0,4000.0)+t1267; + t1441 = t1352+RATIONAL(-607.0,12000.0)+t1468; + t1440 = t1280+t1518; + t1439 = t1327+t1279+t1469; + t1438 = RATIONAL(-79.0,2000.0)+t1345+RATIONAL(-1.0,10.0)*z+t1268; + t1437 = RATIONAL(79.0,3000.0)+t1268+t1508; + t1436 = RATIONAL(-41.0,750.0)+t1265+t1520; + t1357 = RATIONAL(-1.0,8.0)*z; + t1435 = RATIONAL(-87.0,2000.0)+t1357+t1271+t1513; + t1434 = t1347+t1267+RATIONAL(-691.0,12000.0)+t1511; + t1433 = t1272+RATIONAL(299.0,4000.0)+t1479+t1484; + t1432 = t1357+RATIONAL(-33.0,2000.0)+t1270+t1522; + t1431 = t1266+t1347+RATIONAL(281.0,12000.0)+t1514; + t1430 = t1269+RATIONAL(191.0,4000.0)+t1484+t1488; + t1418 = RATIONAL(-1.0,500.0); + t1417 = RATIONAL(-9.0,250.0); + t1416 = RATIONAL(2.0,125.0); + t1413 = RATIONAL(-7.0,250.0); + t1412 = RATIONAL(-1.0,250.0); + t1411 = RATIONAL(-3.0,2000.0); + t1409 = RATIONAL(-1.0,2000.0); + t1406 = RATIONAL(11.0,500.0); + t1405 = RATIONAL(-7.0,2000.0); + t1403 = RATIONAL(-17.0,500.0); + t1398 = RATIONAL(19.0,1000.0); + t1397 = RATIONAL(43.0,2000.0); + t1396 = RATIONAL(31.0,2000.0); + t1395 = RATIONAL(37.0,2000.0); + t1393 = RATIONAL(27.0,2000.0); + t1390 = RATIONAL(-37.0,1000.0); + t1388 = RATIONAL(13.0,1000.0); + t1386 = RATIONAL(-19.0,2000.0); + t1384 = RATIONAL(-21.0,2000.0); + t1382 = RATIONAL(93.0,2000.0); + t1381 = RATIONAL(-23.0,2000.0); + t1380 = RATIONAL(-31.0,1000.0); + t1378 = RATIONAL(-17.0,2000.0); + t1377 = RATIONAL(-39.0,2000.0); + t1376 = RATIONAL(33.0,2000.0); + t1375 = RATIONAL(49.0,2000.0); + t1374 = RATIONAL(-13.0,2000.0); + t1371 = RATIONAL(111.0,2000.0); + t1367 = RATIONAL(129.0,2000.0); + t1339 = t1412*x; + t1338 = t1420*y; + t1337 = t1418*y; + t1336 = t1418*x; + t1334 = t1416*y; + t1330 = t1417*x; + t1326 = t1416*x; + t1323 = t1412*y; + t1322 = t1413*y; + t1320 = t1417*y; + t1319 = t1413*x; + t1318 = t1410*y; + t1317 = t1404*x; + t1315 = t1406*x; + t1314 = t1406*y; + t1312 = t1409*x; + t1311 = t1407*x; + t1310 = t1405*y; + t1309 = t1409*y; + t1308 = t1405*x; + t1307 = t1403*x; + t1306 = t1403*y; + t1304 = t1398*x; + t1303 = t1388*x; + t1302 = t1388*y; + t1301 = t1398*y; + t1300 = t1384*y; + t1299 = t1395*y; + t1298 = t1390*x; + t1297 = t1375*x; + t1296 = t1374*y; + t1295 = t1396*x; + t1294 = t1397*x; + t1293 = t1386*y; + t1291 = t1390*y; + t1290 = t1374*x; + t1289 = t1380*x; + t1288 = t1396*y; + t1287 = t1380*y; + t1286 = t1384*x; + t1285 = t1395*x; + t1284 = t1375*y; + t1283 = t1400*y; + t1282 = t1386*x; + t1281 = t1397*y; + coeff_dz_m1_m1_m1 = RATIONAL(-47.0,1500.0)+t1499*RATIONAL(147.0,2000.0)+ +t1459+t1465; + coeff_dz_0_m1_m1 = t1275+t1367*y+t1312+t1277+t1441+t1482; + coeff_dz_p1_m1_m1 = t1298+t1278+t1371*y+t1275+t1449+t1515; + coeff_dz_p2_m1_m1 = t1330+t1382*y+t1431+t1478; + coeff_dz_m1_0_m1 = t1309+t1367*x+t1441+t1490+t1497; + coeff_dz_0_0_m1 = t1310+t1308+RATIONAL(-437.0,6000.0)+t1458+t1463; + coeff_dz_p1_0_m1 = t1307+t1296+t1434+t1467; + coeff_dz_p2_0_m1 = t1280+t1293+t1461+t1517+t1521; + coeff_dz_m1_p1_m1 = t1371*x+t1291+t1461+t1494+t1515; + coeff_dz_0_p1_m1 = t1290+t1306+t1434+t1493; + coeff_dz_p1_p1_m1 = RATIONAL(-17.0,375.0)+t1289+t1287+t1458+t1464; + coeff_dz_p2_p1_m1 = t1322+t1280+t1443+t1497+t1510; + coeff_dz_m1_p2_m1 = t1320+t1382*x+t1431+t1473; + coeff_dz_0_p2_m1 = t1283+t1282+t1449+t1483+t1517; + coeff_dz_p1_p2_m1 = t1319+t1338+t1443+t1470+t1482; + coeff_dz_p2_p2_m1 = t1313+RATIONAL(307.0,6000.0)+t1311+t1459+t1462; + coeff_dz_m1_m1_0 = t1309+t1312+RATIONAL(-13.0,80.0)*z+RATIONAL(-39.0, +4000.0)+t1453+t1487; + coeff_dz_0_m1_0 = t1378*x+t1310+t1439+t1481+t1516; + coeff_dz_p1_m1_0 = t1339+t1296+t1439+t1491+t1506; + coeff_dz_p2_m1_0 = t1293+t1303+t1432+t1473; + coeff_dz_m1_0_0 = t1308+t1378*y+t1456+t1496+t1516; + coeff_dz_0_0_0 = t1335+RATIONAL(-161.0,4000.0)+t1282+t1293+t1451+t1495; + coeff_dz_p1_0_0 = t1317+t1300+t1435+t1493; + coeff_dz_p2_0_0 = t1381*y+t1326+t1442+t1492+t1494; + coeff_dz_m1_p1_0 = t1290+t1323+t1456+t1489+t1506; + coeff_dz_0_p1_0 = t1305+t1286+t1435+t1467; + coeff_dz_p1_p1_0 = RATIONAL(-191.0,4000.0)+t1336+t1351+t1337+t1451+t1477; + coeff_dz_p2_p1_0 = t1318+t1304+t1438+t1490+t1492; + coeff_dz_m1_p2_0 = t1282+t1302+t1432+t1478; + coeff_dz_0_p2_0 = t1334+t1381*x+t1278+t1442+t1505; + coeff_dz_p1_p2_0 = t1301+t1316+t1277+t1438+t1505; + coeff_dz_p2_p2_0 = RATIONAL(-129.0,4000.0)+t1314+RATIONAL(-7.0,80.0)*z+ +t1315+t1453+t1480; + coeff_dz_m1_m1_p1 = t1291+RATIONAL(157.0,2000.0)+t1298+t1452+t1462; + coeff_dz_0_m1_p1 = t1306+t1339+t1444+t1470+t1492; + coeff_dz_p1_m1_p1 = t1376*x+t1287+t1446+t1483+t1492; + coeff_dz_p2_m1_p1 = t1297+t1322+t1430+t1473; + coeff_dz_m1_0_p1 = t1274+t1307+t1323+t1444+t1507; + coeff_dz_0_0_p1 = RATIONAL(21.0,250.0)+t1305+t1317+t1455+t1464; + coeff_dz_p1_0_p1 = t1295+t1337+t1433+t1493; + coeff_dz_p2_0_p1 = t1318+t1294+t1454+t1496+t1512; + coeff_dz_m1_p1_p1 = t1273+t1376*y+t1289+t1446+t1507; + coeff_dz_0_p1_p1 = t1336+t1288+t1433+t1467; + coeff_dz_p1_p1_p1 = RATIONAL(133.0,2000.0)+t1499*RATIONAL(29.0,2000.0)+ +t1455+t1463; + coeff_dz_p2_p1_p1 = t1285+t1393*y+t1454+t1489+t1509; + coeff_dz_m1_p2_p1 = t1284+t1319+t1430+t1478; + coeff_dz_0_p2_p1 = t1316+t1281+t1457+t1481+t1512; + coeff_dz_p1_p2_p1 = t1393*x+t1299+t1457+t1491+t1509; + coeff_dz_p2_p2_p1 = RATIONAL(13.0,500.0)+t1288+t1295+t1452+t1465; + coeff_dz_m1_m1_p2 = t1330+t1340+t1320+RATIONAL(-449.0,12000.0)+t1447+ +t1480; + coeff_dz_0_m1_p2 = t1275+t1303+t1283+t1448+t1504; + coeff_dz_p1_m1_p2 = t1338+t1297+t1275+t1445+t1503; + coeff_dz_p2_m1_p2 = t1313+t1411*x+t1436+t1478; + coeff_dz_m1_0_p2 = t1302+t1450+t1504+t1521; + coeff_dz_0_0_p2 = t1326+t1334+RATIONAL(349.0,12000.0)+t1346+t1460+t1477; + coeff_dz_p1_0_p2 = t1294+t1301+t1437+t1467; + coeff_dz_p2_0_p2 = t1273+t1286+t1314+t1440+t1501; + coeff_dz_m1_p1_p2 = t1284+t1450+t1503+t1510; + coeff_dz_0_p1_p2 = t1304+t1281+t1437+t1493; + coeff_dz_p1_p1_p2 = t1285+t1299+t1342+RATIONAL(319.0,12000.0)+t1460+t1495 +; + coeff_dz_p2_p1_p2 = t1274+t1377*x+t1288+t1440+t1502; + coeff_dz_m1_p2_p2 = t1311+t1411*y+t1436+t1473; + coeff_dz_0_p2_p2 = t1279+t1315+t1300+t1448+t1501; + coeff_dz_p1_p2_p2 = t1377*y+t1295+t1279+t1445+t1502; + coeff_dz_p2_p2_p2 = t1498+RATIONAL(-539.0,12000.0)+t1499*RATIONAL(-57.0, +2000.0)+t1447+t1487; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dzz.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dzz.compute.c new file mode 100644 index 0000000..60580ca --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeff-dzz.compute.c @@ -0,0 +1,188 @@ +fp t1914, + t1899, + t1901, + t1956, + t1916, + t1893, + t1915, + t1896, + t1955, + t1892, + t1954, + t1895, + t1953, + t1910, + t1903, + t1952, + t1951, + t1919, + t1950, + t1917, + t1894, + t1949, + t1897, + t1948, + t1947, + t1946, + t1898, + t1945, + t1944, + t1943, + t1942, + t1941, + t1940, + t1939, + t1938, + t1937, + t1936, + t1909, + t1935, + t1934, + t1933, + t1932, + t1931, + t1912, + t1930, + t1929, + t1907, + t1928, + t1927, + t1911, + t1926, + t1925, + t1924, + t1923, + t1922, + t1921, + t1920, + t1918, + t1908, + t1904, + t1902, + t1900; + t1914 = RATIONAL(-3.0,80.0); + t1899 = t1914*y; + t1901 = t1914*x; + t1956 = t1899+t1901; + t1916 = RATIONAL(1.0,80.0); + t1893 = t1916*y; + t1915 = RATIONAL(-1.0,80.0); + t1896 = t1915*x; + t1955 = t1893+t1896; + t1892 = RATIONAL(3.0,16.0)*z; + t1954 = t1892+RATIONAL(-1.0,10.0); + t1895 = t1915*y; + t1953 = t1895+t1901; + t1910 = RATIONAL(1.0,16.0); + t1903 = t1910*z; + t1952 = t1903+t1916; + t1951 = t1896+t1899; + t1919 = RATIONAL(-1.0,8.0); + t1950 = t1892+t1919; + t1917 = RATIONAL(3.0,80.0); + t1894 = t1917*y; + t1949 = t1896+t1894; + t1897 = t1916*x; + t1948 = t1893+t1897; + t1947 = t1895+t1897; + t1946 = t1901+t1894; + t1898 = t1917*x; + t1945 = t1898+t1903; + t1944 = t1897+t1903; + t1943 = t1892+t1898; + t1942 = t1892+t1894; + t1941 = t1893+t1901; + t1940 = t1892+RATIONAL(-9.0,80.0); + t1939 = t1898+t1894; + t1938 = t1895+t1896; + t1937 = t1897+t1899; + t1936 = t1899+t1903; + t1909 = RATIONAL(1.0,20.0); + t1935 = t1909+t1948; + t1934 = t1895+t1898+t1909; + t1933 = t1910+t1955; + t1932 = t1897+t1894+t1917; + t1931 = t1893+t1898+t1917; + t1912 = RATIONAL(3.0,40.0); + t1930 = t1912+t1941; + t1929 = t1910+t1947; + t1907 = RATIONAL(1.0,40.0); + t1928 = t1907+t1939; + t1927 = t1898+t1899+t1910; + t1911 = RATIONAL(7.0,80.0); + t1926 = t1911+t1951; + t1925 = t1911+t1953; + t1924 = t1909+t1949; + t1923 = t1910+t1946; + t1922 = RATIONAL(1.0,10.0)+t1956; + t1921 = t1912+t1938; + t1920 = t1912+t1937; + t1918 = RATIONAL(-1.0,40.0); + t1908 = RATIONAL(-3.0,20.0); + t1904 = RATIONAL(-11.0,80.0); + t1902 = RATIONAL(-3.0,16.0)*z; + t1900 = RATIONAL(-1.0,16.0)*z; + coeff_dzz_m1_m1_m1 = t1900+t1922; + coeff_dzz_0_m1_m1 = t1900+t1926; + coeff_dzz_p1_m1_m1 = t1900+t1920; + coeff_dzz_p2_m1_m1 = t1900+t1927; + coeff_dzz_m1_0_m1 = t1900+t1925; + coeff_dzz_0_0_m1 = t1900+t1921; + coeff_dzz_p1_0_m1 = t1900+t1929; + coeff_dzz_p2_0_m1 = t1900+t1934; + coeff_dzz_m1_p1_m1 = t1900+t1930; + coeff_dzz_0_p1_m1 = t1900+t1933; + coeff_dzz_p1_p1_m1 = t1900+t1935; + coeff_dzz_p2_p1_m1 = t1900+t1931; + coeff_dzz_m1_p2_m1 = t1900+t1923; + coeff_dzz_0_p2_m1 = t1900+t1924; + coeff_dzz_p1_p2_m1 = t1900+t1932; + coeff_dzz_p2_p2_m1 = t1900+t1928; + coeff_dzz_m1_m1_0 = t1892+RATIONAL(-13.0,80.0)+t1939; + coeff_dzz_0_m1_0 = t1908+t1897+t1942; + coeff_dzz_p1_m1_0 = t1896+t1904+t1942; + coeff_dzz_p2_m1_0 = t1919+t1901+t1942; + coeff_dzz_m1_0_0 = t1893+t1908+t1943; + coeff_dzz_0_0_0 = t1892+t1904+t1948; + coeff_dzz_p1_0_0 = t1950+t1955; + coeff_dzz_p2_0_0 = t1940+t1941; + coeff_dzz_m1_p1_0 = t1904+t1895+t1943; + coeff_dzz_0_p1_0 = t1947+t1950; + coeff_dzz_p1_p1_0 = t1938+t1940; + coeff_dzz_p2_p1_0 = t1953+t1954; + coeff_dzz_m1_p2_0 = t1899+t1919+t1943; + coeff_dzz_0_p2_0 = t1937+t1940; + coeff_dzz_p1_p2_0 = t1951+t1954; + coeff_dzz_p2_p2_0 = t1892+RATIONAL(-7.0,80.0)+t1956; + coeff_dzz_m1_m1_p1 = t1902+t1928; + coeff_dzz_0_m1_p1 = t1902+t1932; + coeff_dzz_p1_m1_p1 = t1902+t1924; + coeff_dzz_p2_m1_p1 = t1902+t1923; + coeff_dzz_m1_0_p1 = t1902+t1931; + coeff_dzz_0_0_p1 = t1902+t1935; + coeff_dzz_p1_0_p1 = t1902+t1933; + coeff_dzz_p2_0_p1 = t1902+t1930; + coeff_dzz_m1_p1_p1 = t1902+t1934; + coeff_dzz_0_p1_p1 = t1902+t1929; + coeff_dzz_p1_p1_p1 = t1902+t1921; + coeff_dzz_p2_p1_p1 = t1902+t1925; + coeff_dzz_m1_p2_p1 = t1902+t1927; + coeff_dzz_0_p2_p1 = t1902+t1920; + coeff_dzz_p1_p2_p1 = t1902+t1926; + coeff_dzz_p2_p2_p1 = t1902+t1922; + coeff_dzz_m1_m1_p2 = t1917+t1901+t1936; + coeff_dzz_0_m1_p2 = t1907+t1896+t1936; + coeff_dzz_p1_m1_p2 = t1897+t1916+t1936; + coeff_dzz_p2_m1_p2 = t1898+t1936; + coeff_dzz_m1_0_p2 = t1903+t1907+t1953; + coeff_dzz_0_0_p2 = t1938+t1952; + coeff_dzz_p1_0_p2 = t1895+t1944; + coeff_dzz_p2_0_p2 = t1895+t1915+t1945; + coeff_dzz_m1_p1_p2 = t1941+t1952; + coeff_dzz_0_p1_p2 = t1903+t1955; + coeff_dzz_p1_p1_p2 = t1893+t1915+t1944; + coeff_dzz_p2_p1_p2 = t1893+t1918+t1945; + coeff_dzz_m1_p2_p2 = t1903+t1946; + coeff_dzz_0_p2_p2 = t1915+t1903+t1949; + coeff_dzz_p1_p2_p2 = t1918+t1894+t1944; + coeff_dzz_p2_p2_p2 = t1903+t1914+t1939; |