diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0')
6 files changed, 874 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..51bec8d --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c @@ -0,0 +1,239 @@ +fp t411; +fp t412; +fp t464; +fp t414; +fp t404; +fp t356; +fp t413; +fp t360; +fp t433; +fp t442; +fp t463; +fp t407; +fp t365; +fp t364; +fp t436; +fp t462; +fp t435; +fp t461; +fp t438; +fp t389; +fp t398; +fp t406; +fp t362; +fp t460; +fp t401; +fp t368; +fp t395; +fp t397; +fp t459; +fp t439; +fp t359; +fp t458; +fp t370; +fp t291; +fp t396; +fp t440; +fp t457; +fp t405; +fp t322; +fp t409; +fp t456; +fp t402; +fp t400; +fp t327; +fp t455; +fp t383; +fp t454; +fp t311; +fp t453; +fp t390; +fp t307; +fp t452; +fp t451; +fp t384; +fp t408; +fp t318; +fp t450; +fp t399; +fp t296; +fp t449; +fp t448; +fp t387; +fp t403; +fp t330; +fp t447; +fp t416; +fp t354; +fp t415; +fp t355; +fp t314; +fp t446; +fp t369; +fp t358; +fp t445; +fp t376; +fp t379; +fp t444; +fp t443; +fp t441; +fp t437; +fp t434; +fp t410; +fp t366; +fp t394; +fp t372; +fp t432; +fp t371; +fp t431; +fp t391; +fp t393; +fp t430; +fp t429; +fp t428; +fp t306; +fp t427; +fp t426; +fp t316; +fp t425; +fp t424; +fp t423; +fp t422; +fp t421; +fp t392; +fp t386; +fp t385; +fp t382; + t411 = x*x; + t412 = y*y; + t464 = t411+t412; + t414 = t411*x; + t404 = RATIONAL(-1.0,60.0); + t356 = t404*t414; + t413 = y*t412; + t360 = t404*t413; + t433 = y*t356+x*t360; + t442 = x*y; + t463 = RATIONAL(11.0,150.0)*t442+t433; + t407 = RATIONAL(1.0,60.0); + t365 = t407*t414; + t364 = t407*t413; + t436 = x*t364+y*t365; + t462 = RATIONAL(-11.0,150.0)*t442+t436; + t435 = t411*t412; + t461 = RATIONAL(-99.0,1225.0)+RATIONAL(-1.0,98.0)*t435; + t438 = y*t411; + t389 = RATIONAL(31.0,210.0); + t398 = RATIONAL(-1.0,140.0); + t406 = RATIONAL(-1.0,30.0); + t362 = t406*t413; + t460 = t389*y+t398*t438+t362; + t401 = RATIONAL(1.0,30.0); + t368 = t401*t413; + t395 = RATIONAL(-31.0,210.0); + t397 = RATIONAL(1.0,140.0); + t459 = t368+t395*y+t397*t438; + t439 = x*t412; + t359 = t401*t414; + t458 = t395*x+t397*t439+t359; + t370 = t406*t414; + t291 = RATIONAL(21.0,200.0)*t442; + t396 = RATIONAL(1.0,120.0); + t440 = x*t413; + t457 = t291+t396*t440+y*t370; + t405 = RATIONAL(-1.0,70.0); + t322 = t405*t438; + t409 = RATIONAL(1.0,84.0); + t456 = t322+t409*y+t364; + t402 = RATIONAL(-1.0,84.0); + t400 = RATIONAL(1.0,70.0); + t327 = t400*t439; + t455 = t327+t356+t402*x; + t383 = RATIONAL(-11.0,105.0); + t454 = t368+t322+t383*y; + t311 = t405*t439; + t453 = t409*x+t311+t365; + t390 = RATIONAL(11.0,105.0); + t307 = t400*t438; + t452 = t307+t390*y+t362; + t451 = t307+t402*y+t360; + t384 = RATIONAL(31.0,420.0); + t408 = RATIONAL(-1.0,35.0); + t318 = t408*t439; + t450 = t384*x+t356+t318; + t399 = RATIONAL(-1.0,120.0); + t296 = RATIONAL(-21.0,200.0)*t442; + t449 = y*t359+t399*t440+t296; + t448 = t389*x+t370+t398*t439; + t387 = RATIONAL(-31.0,420.0); + t403 = RATIONAL(1.0,35.0); + t330 = t403*t439; + t447 = t330+t387*x+t365; + t416 = t412*t412; + t354 = t396*t416; + t415 = t411*t411; + t355 = t396*t415; + t314 = RATIONAL(1.0,49.0)*t435; + t446 = RATIONAL(51.0,1225.0)+t314+t354+t355+t464*RATIONAL(-289.0,5880.0); + t369 = t406*t416; + t358 = t406*t415; + t445 = RATIONAL(1.0,196.0)*t435+t369+RATIONAL(-24.0,1225.0)+t358+t464* +RATIONAL(181.0,1470.0); + t376 = RATIONAL(71.0,5880.0); + t379 = RATIONAL(113.0,735.0); + t444 = t376*t412+t379*t411+t354+t358+t461; + t443 = t369+t376*t411+t379*t412+t355+t461; + t441 = t414*y; + t437 = RATIONAL(-1.0,49.0)*t435+RATIONAL(96.0,1225.0); + t434 = RATIONAL(1.0,98.0)*t435+RATIONAL(246.0,1225.0); + t410 = RATIONAL(1.0,20.0); + t366 = t410*t416; + t394 = RATIONAL(-41.0,196.0); + t372 = RATIONAL(191.0,5880.0); + t432 = t366+t372*t411+t355+t394*t412+t437; + t371 = t410*t415; + t431 = t371+t394*t411+t372*t412+t354+t437; + t391 = RATIONAL(83.0,735.0); + t393 = RATIONAL(-53.0,196.0); + t430 = t369+t391*t412+t371+t393*t411+t434; + t429 = t366+t393*t412+t391*t411+t358+t434; + t428 = RATIONAL(-37.0,300.0)*t442+t436+t445; + t306 = t403*t438; + t427 = t306+t387*y+t364+t446; + t426 = RATIONAL(37.0,300.0)*t442+t433+t445; + t316 = t408*t438; + t425 = t384*y+t360+t316+t446; + t424 = t390*x+t370+t327+t444; + t423 = t383*x+t311+t359+t444; + t422 = x*t368+t296+t399*t441+t443; + t421 = t396*t441+t291+x*t362+t443; + t392 = RATIONAL(-17.0,105.0); + t386 = RATIONAL(17.0,105.0); + t385 = RATIONAL(17.0,420.0); + t382 = RATIONAL(-17.0,420.0); + coeffs_I->coeff_m2_m2 = t425+t450+t462; + coeffs_I->coeff_m1_m2 = t423+t451+t457; + coeffs_I->coeff_0_m2 = t306+t382*y+t360+t431; + coeffs_I->coeff_p1_m2 = t424+t449+t451; + coeffs_I->coeff_p2_m2 = t425+t447+t463; + coeffs_I->coeff_m2_m1 = t421+t454+t455; + coeffs_I->coeff_m1_m1 = t426+t458+t459; + coeffs_I->coeff_0_m1 = t392*y+t368+t307+t430; + coeffs_I->coeff_p1_m1 = t428+t448+t459; + coeffs_I->coeff_p2_m1 = t422+t453+t454; + coeffs_I->coeff_m2_0 = t356+t330+t382*x+t432; + coeffs_I->coeff_m1_0 = t359+t392*x+t327+t429; + coeffs_I->coeff_0_0 = RATIONAL(541.0,1225.0)+t366+t314+t371+t464*RATIONAL +(-57.0,196.0); + coeffs_I->coeff_p1_0 = t311+t370+t386*x+t429; + coeffs_I->coeff_p2_0 = t318+t365+t385*x+t432; + coeffs_I->coeff_m2_p1 = t422+t452+t455; + coeffs_I->coeff_m1_p1 = t428+t458+t460; + coeffs_I->coeff_0_p1 = t386*y+t362+t322+t430; + coeffs_I->coeff_p1_p1 = t426+t448+t460; + coeffs_I->coeff_p2_p1 = t421+t452+t453; + coeffs_I->coeff_m2_p2 = t427+t450+t463; + coeffs_I->coeff_m1_p2 = t423+t449+t456; + coeffs_I->coeff_0_p2 = t316+t364+t385*y+t431; + coeffs_I->coeff_p1_p2 = t424+t456+t457; + coeffs_I->coeff_p2_p2 = t427+t447+t462; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..295e9c5 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c @@ -0,0 +1,155 @@ +fp t532; +fp t503; +fp t533; +fp t534; +fp t551; +fp t569; +fp t511; +fp t550; +fp t568; +fp t557; +fp t567; +fp t522; +fp t535; +fp t507; +fp t566; +fp t481; +fp t565; +fp t526; +fp t512; +fp t564; +fp t505; +fp t491; +fp t563; +fp t506; +fp t490; +fp t562; +fp t556; +fp t561; +fp t524; +fp t501; +fp t560; +fp t559; +fp t558; +fp t508; +fp t555; +fp t554; +fp t482; +fp t525; +fp t509; +fp t553; +fp t521; +fp t504; +fp t552; +fp t475; +fp t549; +fp t467; +fp t548; +fp t514; +fp t547; +fp t546; +fp t545; +fp t544; +fp t486; +fp t543; +fp t542; +fp t541; +fp t540; +fp t539; +fp t538; +fp t483; +fp t479; +fp t473; +fp t469; +fp t466; +fp t465; + t532 = x*x; + t503 = RATIONAL(1.0,20.0)*t532; + t533 = y*y; + t534 = y*t533; + t551 = y*t503+RATIONAL(1.0,60.0)*t534; + t569 = t551+RATIONAL(-37.0,300.0)*y; + t511 = RATIONAL(-1.0,20.0)*t532; + t550 = RATIONAL(-1.0,60.0)*t534+y*t511; + t568 = t550+RATIONAL(37.0,300.0)*y; + t557 = x*t533; + t567 = RATIONAL(1.0,98.0)*t557+RATIONAL(181.0,735.0)*x; + t522 = RATIONAL(1.0,30.0); + t535 = x*t532; + t507 = t522*t535; + t566 = RATIONAL(71.0,2940.0)*x+t507; + t481 = RATIONAL(2.0,49.0)*t557; + t565 = RATIONAL(-289.0,2940.0)*x+t481+t507; + t526 = RATIONAL(-1.0,35.0); + t512 = t526*t533; + t564 = t511+t512+RATIONAL(31.0,420.0); + t505 = RATIONAL(-1.0,10.0)*t532; + t491 = RATIONAL(21.0,200.0)*y; + t563 = y*t505+t491+RATIONAL(1.0,120.0)*t534; + t506 = RATIONAL(1.0,10.0)*t532; + t490 = RATIONAL(-21.0,200.0)*y; + t562 = t490+RATIONAL(-1.0,120.0)*t534+y*t506; + t556 = y*t532; + t561 = RATIONAL(1.0,40.0)*t556+t491+RATIONAL(-1.0,30.0)*t534; + t524 = RATIONAL(1.0,35.0); + t501 = t524*t533; + t560 = t503+RATIONAL(-31.0,420.0)+t501; + t559 = RATIONAL(-1.0,40.0)*t556+t490+t522*t534; + t558 = x*y; + t508 = RATIONAL(-2.0,15.0)*t535; + t555 = t506+t508; + t554 = t505+t508; + t482 = RATIONAL(-1.0,49.0)*t557; + t525 = RATIONAL(-1.0,70.0); + t509 = t525*t533; + t553 = t482+t509; + t521 = RATIONAL(1.0,70.0); + t504 = t521*t533; + t552 = t504+t482; + t475 = RATIONAL(-2.0,49.0)*t557; + t549 = t475+RATIONAL(191.0,2940.0)*x+t507; + t467 = RATIONAL(1.0,49.0)*t557; + t548 = t508+RATIONAL(166.0,735.0)*x+t467; + t514 = RATIONAL(1.0,5.0)*t535; + t547 = t475+t514+RATIONAL(-41.0,98.0)*x; + t546 = t467+RATIONAL(-53.0,98.0)*x+t514; + t545 = t511+RATIONAL(-1.0,84.0)+t552+t566; + t544 = t503+RATIONAL(1.0,84.0)+t553+t566; + t486 = RATIONAL(226.0,735.0)*x; + t543 = t486+RATIONAL(-11.0,105.0)+t553+t555; + t542 = RATIONAL(1.0,140.0)*t533+RATIONAL(-31.0,210.0)+t555+t567; + t541 = t486+RATIONAL(11.0,105.0)+t552+t554; + t540 = RATIONAL(-1.0,140.0)*t533+RATIONAL(31.0,210.0)+t554+t567; + t539 = RATIONAL(-11.0,150.0)*y+t551+t565; + t538 = RATIONAL(11.0,150.0)*y+t550+t565; + t483 = t524*t558; + t479 = RATIONAL(-2.0,35.0)*t558; + t473 = RATIONAL(2.0,35.0)*t558; + t469 = t526*t558; + t466 = t521*t558; + t465 = t525*t558; + coeffs_dx->coeff_m2_m2 = t479+t539+t564; + coeffs_dx->coeff_m1_m2 = t483+t543+t563; + coeffs_dx->coeff_0_m2 = t473+t547; + coeffs_dx->coeff_p1_m2 = t483+t541+t562; + coeffs_dx->coeff_p2_m2 = t479+t538+t560; + coeffs_dx->coeff_m2_m1 = t469+t545+t561; + coeffs_dx->coeff_m1_m1 = t466+t542+t568; + coeffs_dx->coeff_0_m1 = t483+t546; + coeffs_dx->coeff_p1_m1 = t466+t540+t569; + coeffs_dx->coeff_p2_m1 = t469+t544+t559; + coeffs_dx->coeff_m2_0 = t501+t511+RATIONAL(-17.0,420.0)+t549; + coeffs_dx->coeff_m1_0 = t504+RATIONAL(-17.0,105.0)+t506+t548; + coeffs_dx->coeff_0_0 = t514+RATIONAL(-57.0,98.0)*x+t481; + coeffs_dx->coeff_p1_0 = t509+RATIONAL(17.0,105.0)+t505+t548; + coeffs_dx->coeff_p2_0 = RATIONAL(17.0,420.0)+t512+t503+t549; + coeffs_dx->coeff_m2_p1 = t483+t545+t559; + coeffs_dx->coeff_m1_p1 = t465+t542+t569; + coeffs_dx->coeff_0_p1 = t469+t546; + coeffs_dx->coeff_p1_p1 = t465+t540+t568; + coeffs_dx->coeff_p2_p1 = t483+t544+t561; + coeffs_dx->coeff_m2_p2 = t473+t538+t564; + coeffs_dx->coeff_m1_p2 = t469+t543+t562; + coeffs_dx->coeff_0_p2 = t479+t547; + coeffs_dx->coeff_p1_p2 = t469+t541+t563; + coeffs_dx->coeff_p2_p2 = t473+t539+t560; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..c1b3b65 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,107 @@ +fp t704; +fp t708; +fp t685; +fp t709; +fp t688; +fp t726; +fp t683; +fp t695; +fp t725; +fp t724; +fp t723; +fp t722; +fp t679; +fp t692; +fp t694; +fp t687; +fp t681; +fp t721; +fp t693; +fp t720; +fp t686; +fp t719; +fp t718; +fp t682; +fp t717; +fp t678; +fp t716; +fp t684; +fp t715; +fp t714; +fp t713; +fp t712; +fp t689; +fp t711; +fp t710; +fp t677; +fp t676; +fp t675; +fp t674; +fp t673; +fp t672; + t704 = RATIONAL(1.0,10.0); + t708 = x*x; + t685 = t704*t708; + t709 = y*y; + t688 = RATIONAL(2.0,49.0)*t709; + t726 = t685+t688+RATIONAL(-289.0,2940.0); + t683 = RATIONAL(-1.0,49.0)*t709; + t695 = RATIONAL(-2.0,5.0)*t708; + t725 = t683+t695+RATIONAL(226.0,735.0); + t724 = RATIONAL(1.0,98.0)*t709+t695+RATIONAL(181.0,735.0); + t723 = t683+t685+RATIONAL(71.0,2940.0); + t722 = x*y; + t679 = t704*x; + t692 = RATIONAL(-1.0,5.0)*x; + t694 = RATIONAL(1.0,5.0)*x; + t687 = RATIONAL(-1.0,10.0)*x; + t681 = RATIONAL(-2.0,49.0)*t709; + t721 = t681+t685+RATIONAL(191.0,2940.0); + t693 = RATIONAL(3.0,5.0)*t708; + t720 = t681+t693+RATIONAL(-41.0,98.0); + t686 = RATIONAL(1.0,49.0)*t709; + t719 = t686+t695+RATIONAL(166.0,735.0); + t718 = t686+t693+RATIONAL(-53.0,98.0); + t682 = RATIONAL(2.0,35.0)*y; + t717 = t682+t726; + t678 = RATIONAL(-1.0,35.0)*y; + t716 = t678+t725; + t684 = RATIONAL(1.0,35.0)*y; + t715 = t684+t725; + t714 = t679+t723; + t713 = t687+t723; + t712 = RATIONAL(-1.0,70.0)*y+t724; + t689 = RATIONAL(-2.0,35.0)*y; + t711 = t689+t726; + t710 = RATIONAL(1.0,70.0)*y+t724; + t677 = y*t692; + t676 = y*t694; + t675 = RATIONAL(1.0,20.0)*t722; + t674 = y*t687; + t673 = RATIONAL(-1.0,20.0)*t722; + t672 = y*t679; + coeffs_dxx->coeff_m2_m2 = t672+t687+t711; + coeffs_dxx->coeff_m1_m2 = t677+t694+t715; + coeffs_dxx->coeff_0_m2 = t682+t720; + coeffs_dxx->coeff_p1_m2 = t676+t692+t715; + coeffs_dxx->coeff_p2_m2 = t674+t679+t711; + coeffs_dxx->coeff_m2_m1 = t675+t678+t713; + coeffs_dxx->coeff_m1_m1 = t694+t674+t710; + coeffs_dxx->coeff_0_m1 = t684+t718; + coeffs_dxx->coeff_p1_m1 = t692+t672+t710; + coeffs_dxx->coeff_p2_m1 = t678+t673+t714; + coeffs_dxx->coeff_m2_0 = t687+t721; + coeffs_dxx->coeff_m1_0 = t694+t719; + coeffs_dxx->coeff_0_0 = t693+RATIONAL(-57.0,98.0)+t688; + coeffs_dxx->coeff_p1_0 = t692+t719; + coeffs_dxx->coeff_p2_0 = t679+t721; + coeffs_dxx->coeff_m2_p1 = t684+t673+t713; + coeffs_dxx->coeff_m1_p1 = t672+t694+t712; + coeffs_dxx->coeff_0_p1 = t678+t718; + coeffs_dxx->coeff_p1_p1 = t692+t674+t712; + coeffs_dxx->coeff_p2_p1 = t684+t675+t714; + coeffs_dxx->coeff_m2_p2 = t687+t674+t717; + coeffs_dxx->coeff_m1_p2 = t694+t676+t716; + coeffs_dxx->coeff_0_p2 = t689+t720; + coeffs_dxx->coeff_p1_p2 = t677+t692+t716; + coeffs_dxx->coeff_p2_p2 = t679+t672+t717; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c new file mode 100644 index 0000000..dd9a747 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c @@ -0,0 +1,119 @@ +fp t774; +fp t775; +fp t793; +fp t763; +fp t747; +fp t765; +fp t732; +fp t792; +fp t740; +fp t753; +fp t791; +fp t767; +fp t771; +fp t790; +fp t789; +fp t788; +fp t787; +fp t786; +fp t731; +fp t785; +fp t784; +fp t728; +fp t783; +fp t730; +fp t782; +fp t768; +fp t773; +fp t781; +fp t780; +fp t779; +fp t778; +fp t777; +fp t776; +fp t770; +fp t766; +fp t764; +fp t762; +fp t755; +fp t754; +fp t751; +fp t746; +fp t745; +fp t742; +fp t736; +fp t733; +fp t729; +fp t727; + t774 = x*x; + t775 = y*y; + t793 = t774+t775; + t763 = RATIONAL(1.0,35.0); + t747 = t763*x; + t765 = RATIONAL(-1.0,35.0); + t732 = t765*y; + t792 = t747+t732; + t740 = t765*x; + t753 = t763*y; + t791 = t740+t753; + t767 = RATIONAL(1.0,10.0); + t771 = RATIONAL(-1.0,40.0); + t790 = t767*t775+t771*t774; + t789 = t767*t774+t771*t775; + t788 = x*y; + t787 = t793*RATIONAL(-1.0,20.0); + t786 = t793*RATIONAL(1.0,20.0); + t731 = RATIONAL(-2.0,49.0)*t788; + t785 = t731+RATIONAL(21.0,200.0); + t784 = t731+RATIONAL(-21.0,200.0); + t728 = RATIONAL(1.0,49.0)*t788; + t783 = t728+RATIONAL(-37.0,300.0)+t786; + t730 = RATIONAL(4.0,49.0)*t788; + t782 = t730+RATIONAL(-11.0,150.0)+t786; + t768 = RATIONAL(-1.0,10.0); + t773 = RATIONAL(1.0,40.0); + t781 = t773*t775+t768*t774+t785; + t780 = t768*t775+t773*t774+t785; + t779 = t740+t732+t784; + t778 = t730+RATIONAL(11.0,150.0)+t787; + t777 = t747+t753+t784; + t776 = t728+RATIONAL(37.0,300.0)+t787; + t770 = RATIONAL(-2.0,35.0); + t766 = RATIONAL(2.0,35.0); + t764 = RATIONAL(-1.0,70.0); + t762 = RATIONAL(1.0,70.0); + t755 = t762*y; + t754 = t770*y; + t751 = t770*x; + t746 = t766*y; + t745 = t766*x; + t742 = t764*y; + t736 = t762*x; + t733 = t764*x; + t729 = RATIONAL(-4.0,49.0)*t788; + t727 = RATIONAL(2.0,49.0)*t788; + coeffs_dxy->coeff_m2_m2 = t754+t751+t782; + coeffs_dxy->coeff_m1_m2 = t781+t792; + coeffs_dxy->coeff_0_m2 = t729+t745; + coeffs_dxy->coeff_p1_m2 = t777+t789; + coeffs_dxy->coeff_p2_m2 = t746+t751+t778; + coeffs_dxy->coeff_m2_m1 = t780+t791; + coeffs_dxy->coeff_m1_m1 = t755+t736+t776; + coeffs_dxy->coeff_0_m1 = t727+t747; + coeffs_dxy->coeff_p1_m1 = t742+t736+t783; + coeffs_dxy->coeff_p2_m1 = t779+t790; + coeffs_dxy->coeff_m2_0 = t746+t729; + coeffs_dxy->coeff_m1_0 = t727+t753; + coeffs_dxy->coeff_0_0 = t730; + coeffs_dxy->coeff_p1_0 = t732+t727; + coeffs_dxy->coeff_p2_0 = t754+t729; + coeffs_dxy->coeff_m2_p1 = t777+t790; + coeffs_dxy->coeff_m1_p1 = t733+t755+t783; + coeffs_dxy->coeff_0_p1 = t727+t740; + coeffs_dxy->coeff_p1_p1 = t733+t742+t776; + coeffs_dxy->coeff_p2_p1 = t780+t792; + coeffs_dxy->coeff_m2_p2 = t745+t754+t778; + coeffs_dxy->coeff_m1_p2 = t779+t789; + coeffs_dxy->coeff_0_p2 = t729+t751; + coeffs_dxy->coeff_p1_p2 = t781+t791; + coeffs_dxy->coeff_p2_p2 = t745+t746+t782; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c new file mode 100644 index 0000000..a1b62b5 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c @@ -0,0 +1,147 @@ +fp t638; +fp t617; +fp t637; +fp t640; +fp t655; +fp t671; +fp t613; +fp t656; +fp t670; +fp t658; +fp t586; +fp t627; +fp t639; +fp t614; +fp t669; +fp t587; +fp t608; +fp t668; +fp t667; +fp t666; +fp t596; +fp t657; +fp t665; +fp t631; +fp t618; +fp t664; +fp t629; +fp t611; +fp t663; +fp t597; +fp t662; +fp t607; +fp t630; +fp t606; +fp t661; +fp t626; +fp t612; +fp t605; +fp t660; +fp t659; +fp t572; +fp t654; +fp t581; +fp t653; +fp t619; +fp t652; +fp t651; +fp t650; +fp t649; +fp t648; +fp t647; +fp t646; +fp t645; +fp t644; +fp t643; +fp t584; +fp t582; +fp t577; +fp t576; +fp t573; +fp t571; + t638 = y*y; + t617 = RATIONAL(-1.0,20.0)*t638; + t637 = x*x; + t640 = x*t637; + t655 = x*t617+RATIONAL(-1.0,60.0)*t640; + t671 = t655+RATIONAL(37.0,300.0)*x; + t613 = RATIONAL(1.0,20.0)*t638; + t656 = x*t613+RATIONAL(1.0,60.0)*t640; + t670 = RATIONAL(-37.0,300.0)*x+t656; + t658 = y*t637; + t586 = RATIONAL(2.0,49.0)*t658; + t627 = RATIONAL(1.0,30.0); + t639 = y*t638; + t614 = t627*t639; + t669 = RATIONAL(-289.0,2940.0)*y+t586+t614; + t587 = RATIONAL(-1.0,49.0)*t658; + t608 = RATIONAL(-2.0,15.0)*t639; + t668 = t587+RATIONAL(226.0,735.0)*y+t608; + t667 = t587+t614+RATIONAL(71.0,2940.0)*y; + t666 = RATIONAL(1.0,98.0)*t658+RATIONAL(181.0,735.0)*y+t608; + t596 = RATIONAL(-21.0,200.0)*x; + t657 = x*t638; + t665 = t627*t640+RATIONAL(-1.0,40.0)*t657+t596; + t631 = RATIONAL(-1.0,35.0); + t618 = t631*t637; + t664 = t617+t618+RATIONAL(31.0,420.0); + t629 = RATIONAL(1.0,35.0); + t611 = t629*t637; + t663 = t611+t613+RATIONAL(-31.0,420.0); + t597 = RATIONAL(21.0,200.0)*x; + t662 = RATIONAL(1.0,40.0)*t657+RATIONAL(-1.0,30.0)*t640+t597; + t607 = RATIONAL(1.0,10.0)*t638; + t630 = RATIONAL(-1.0,70.0); + t606 = t630*t637; + t661 = RATIONAL(-11.0,105.0)+t607+t606; + t626 = RATIONAL(1.0,70.0); + t612 = t626*t637; + t605 = RATIONAL(-1.0,10.0)*t638; + t660 = t612+t605+RATIONAL(11.0,105.0); + t659 = x*y; + t572 = RATIONAL(1.0,49.0)*t658; + t654 = t572+RATIONAL(166.0,735.0)*y+t608; + t581 = RATIONAL(-2.0,49.0)*t658; + t653 = RATIONAL(191.0,2940.0)*y+t581+t614; + t619 = RATIONAL(1.0,5.0)*t639; + t652 = t581+t619+RATIONAL(-41.0,98.0)*y; + t651 = t572+t619+RATIONAL(-53.0,98.0)*y; + t650 = RATIONAL(-11.0,150.0)*x+t656+t669; + t649 = RATIONAL(1.0,140.0)*t637+t607+RATIONAL(-31.0,210.0)+t666; + t648 = x*t607+t596+RATIONAL(-1.0,120.0)*t640+t668; + t647 = x*t605+t597+RATIONAL(1.0,120.0)*t640+t668; + t646 = RATIONAL(-1.0,140.0)*t637+t605+RATIONAL(31.0,210.0)+t666; + t645 = t612+t617+RATIONAL(-1.0,84.0)+t667; + t644 = t613+t606+RATIONAL(1.0,84.0)+t667; + t643 = RATIONAL(11.0,150.0)*x+t655+t669; + t584 = RATIONAL(-2.0,35.0)*t659; + t582 = t631*t659; + t577 = t629*t659; + t576 = RATIONAL(2.0,35.0)*t659; + t573 = t626*t659; + t571 = t630*t659; + coeffs_dy->coeff_m2_m2 = t584+t650+t664; + coeffs_dy->coeff_m1_m2 = t582+t645+t662; + coeffs_dy->coeff_0_m2 = t611+t617+RATIONAL(-17.0,420.0)+t653; + coeffs_dy->coeff_p1_m2 = t577+t645+t665; + coeffs_dy->coeff_p2_m2 = t576+t643+t664; + coeffs_dy->coeff_m2_m1 = t577+t647+t661; + coeffs_dy->coeff_m1_m1 = t573+t649+t671; + coeffs_dy->coeff_0_m1 = RATIONAL(-17.0,105.0)+t607+t612+t654; + coeffs_dy->coeff_p1_m1 = t571+t649+t670; + coeffs_dy->coeff_p2_m1 = t582+t648+t661; + coeffs_dy->coeff_m2_0 = t576+t652; + coeffs_dy->coeff_m1_0 = t577+t651; + coeffs_dy->coeff_0_0 = RATIONAL(-57.0,98.0)*y+t619+t586; + coeffs_dy->coeff_p1_0 = t582+t651; + coeffs_dy->coeff_p2_0 = t584+t652; + coeffs_dy->coeff_m2_p1 = t577+t648+t660; + coeffs_dy->coeff_m1_p1 = t573+t646+t670; + coeffs_dy->coeff_0_p1 = RATIONAL(17.0,105.0)+t606+t605+t654; + coeffs_dy->coeff_p1_p1 = t571+t646+t671; + coeffs_dy->coeff_p2_p1 = t582+t647+t660; + coeffs_dy->coeff_m2_p2 = t584+t643+t663; + coeffs_dy->coeff_m1_p2 = t582+t644+t665; + coeffs_dy->coeff_0_p2 = t618+RATIONAL(17.0,420.0)+t613+t653; + coeffs_dy->coeff_p1_p2 = t577+t644+t662; + coeffs_dy->coeff_p2_p2 = t576+t650+t663; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c new file mode 100644 index 0000000..29eef2a --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-maximum-degree/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c @@ -0,0 +1,107 @@ +fp t831; +fp t817; +fp t830; +fp t848; +fp t808; +fp t826; +fp t811; +fp t847; +fp t846; +fp t813; +fp t845; +fp t844; +fp t814; +fp t816; +fp t800; +fp t812; +fp t807; +fp t843; +fp t815; +fp t842; +fp t801; +fp t841; +fp t840; +fp t839; +fp t838; +fp t806; +fp t837; +fp t809; +fp t836; +fp t835; +fp t834; +fp t833; +fp t832; +fp t810; +fp t805; +fp t799; +fp t798; +fp t797; +fp t796; +fp t795; +fp t794; + t831 = y*y; + t817 = RATIONAL(-2.0,5.0)*t831; + t830 = x*x; + t848 = RATIONAL(1.0,98.0)*t830+t817+RATIONAL(181.0,735.0); + t808 = RATIONAL(-1.0,49.0)*t830; + t826 = RATIONAL(1.0,10.0); + t811 = t826*t831; + t847 = t808+t811+RATIONAL(71.0,2940.0); + t846 = t808+t817+RATIONAL(226.0,735.0); + t813 = RATIONAL(2.0,49.0)*t830; + t845 = t811+t813+RATIONAL(-289.0,2940.0); + t844 = x*y; + t814 = RATIONAL(-1.0,5.0)*y; + t816 = RATIONAL(1.0,5.0)*y; + t800 = t826*y; + t812 = RATIONAL(-1.0,10.0)*y; + t807 = RATIONAL(-2.0,49.0)*t830; + t843 = t807+t811+RATIONAL(191.0,2940.0); + t815 = RATIONAL(3.0,5.0)*t831; + t842 = t807+t815+RATIONAL(-41.0,98.0); + t801 = RATIONAL(1.0,49.0)*t830; + t841 = t801+t815+RATIONAL(-53.0,98.0); + t840 = t801+RATIONAL(166.0,735.0)+t817; + t839 = t800+t845; + t838 = t812+t845; + t806 = RATIONAL(1.0,35.0)*x; + t837 = t806+t847; + t809 = RATIONAL(-1.0,35.0)*x; + t836 = t809+t847; + t835 = t814+t846; + t834 = t816+t846; + t833 = RATIONAL(1.0,70.0)*x+t848; + t832 = RATIONAL(-1.0,70.0)*x+t848; + t810 = RATIONAL(-2.0,35.0)*x; + t805 = RATIONAL(2.0,35.0)*x; + t799 = x*t814; + t798 = x*t816; + t797 = RATIONAL(1.0,20.0)*t844; + t796 = x*t812; + t795 = RATIONAL(-1.0,20.0)*t844; + t794 = x*t800; + coeffs_dyy->coeff_m2_m2 = t810+t794+t838; + coeffs_dyy->coeff_m1_m2 = t797+t812+t836; + coeffs_dyy->coeff_0_m2 = t812+t843; + coeffs_dyy->coeff_p1_m2 = t812+t795+t837; + coeffs_dyy->coeff_p2_m2 = t796+t805+t838; + coeffs_dyy->coeff_m2_m1 = t799+t806+t834; + coeffs_dyy->coeff_m1_m1 = t796+t816+t833; + coeffs_dyy->coeff_0_m1 = t816+t840; + coeffs_dyy->coeff_p1_m1 = t794+t816+t832; + coeffs_dyy->coeff_p2_m1 = t809+t798+t834; + coeffs_dyy->coeff_m2_0 = t805+t842; + coeffs_dyy->coeff_m1_0 = t806+t841; + coeffs_dyy->coeff_0_0 = t813+t815+RATIONAL(-57.0,98.0); + coeffs_dyy->coeff_p1_0 = t809+t841; + coeffs_dyy->coeff_p2_0 = t810+t842; + coeffs_dyy->coeff_m2_p1 = t806+t798+t835; + coeffs_dyy->coeff_m1_p1 = t794+t814+t833; + coeffs_dyy->coeff_0_p1 = t814+t840; + coeffs_dyy->coeff_p1_p1 = t814+t796+t832; + coeffs_dyy->coeff_p2_p1 = t799+t809+t835; + coeffs_dyy->coeff_m2_p2 = t796+t810+t839; + coeffs_dyy->coeff_m1_p2 = t795+t800+t836; + coeffs_dyy->coeff_0_p2 = t800+t843; + coeffs_dyy->coeff_p1_p2 = t800+t797+t837; + coeffs_dyy->coeff_p2_p2 = t794+t805+t839; |