diff options
author | jthorn <jthorn@0f49ee68-0e4f-0410-9b9c-b2c123ded7ef> | 2003-07-06 11:16:19 +0000 |
---|---|---|
committer | jthorn <jthorn@0f49ee68-0e4f-0410-9b9c-b2c123ded7ef> | 2003-07-06 11:16:19 +0000 |
commit | a07489dec7a4e1153624e158a2c5f2837b9247de (patch) | |
tree | 83be503e3cdaf39c578202c0fcdcf71337845e42 /src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0 | |
parent | 892b8a2d121db4c1e436177cb19baa06eb8d0e4a (diff) |
This commit was generated by cvs2svn to compensate for changes in r2,
which included commits to RCS files with non-trunk default branches.
git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/AEILocalInterp/trunk@3 0f49ee68-0e4f-0410-9b9c-b2c123ded7ef
Diffstat (limited to 'src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0')
10 files changed, 3648 insertions, 0 deletions
diff --git a/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-I.compute.c b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..cb38b59 --- /dev/null +++ b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-I.compute.c @@ -0,0 +1,938 @@ +fp t652; +fp t401; +fp t470; +fp t748; +fp t509; +fp t489; +fp t341; +fp t260; +fp t510; +fp t506; +fp t498; +fp t378; +fp t490; +fp t511; +fp t346; +fp t233; +fp t747; +fp t642; +fp t746; +fp t448; +fp t451; +fp t638; +fp t745; +fp t493; +fp t607; +fp t274; +fp t508; +fp t497; +fp t365; +fp t744; +fp t612; +fp t265; +fp t375; +fp t743; +fp t350; +fp t237; +fp t502; +fp t364; +fp t742; +fp t407; +fp t482; +fp t339; +fp t483; +fp t631; +fp t741; +fp t349; +fp t740; +fp t433; +fp t495; +fp t466; +fp t649; +fp t739; +fp t738; +fp t501; +fp t380; +fp t654; +fp t737; +fp t384; +fp t226; +fp t601; +fp t736; +fp t611; +fp t256; +fp t735; +fp t481; +fp t352; +fp t276; +fp t734; +fp t353; +fp t232; +fp t733; +fp t253; +fp t351; +fp t230; +fp t507; +fp t363; +fp t732; +fp t347; +fp t249; +fp t586; +fp t731; +fp t467; +fp t206; +fp t160; +fp t730; +fp t651; +fp t616; +fp t164; +fp t729; +fp t157; +fp t728; +fp t505; +fp t374; +fp t727; +fp t726; +fp t392; +fp t725; +fp t245; +fp t724; +fp t438; +fp t723; +fp t418; +fp t462; +fp t722; +fp t376; +fp t721; +fp t261; +fp t610; +fp t720; +fp t719; +fp t718; +fp t336; +fp t251; +fp t717; +fp t453; +fp t476; +fp t639; +fp t214; +fp t163; +fp t716; +fp t429; +fp t405; +fp t715; +fp t258; +fp t590; +fp t714; +fp t367; +fp t345; +fp t243; +fp t713; +fp t712; +fp t469; +fp t473; +fp t711; +fp t248; +fp t278; +fp t604; +fp t389; +fp t710; +fp t262; +fp t263; +fp t709; +fp t707; +fp t706; +fp t270; +fp t705; +fp t161; +fp t704; +fp t463; +fp t194; +fp t158; +fp t703; +fp t640; +fp t702; +fp t641; +fp t216; +fp t162; +fp t701; +fp t630; +fp t700; +fp t159; +fp t699; +fp t577; +fp t698; +fp t697; +fp t488; +fp t628; +fp t396; +fp t696; +fp t695; +fp t694; +fp t693; +fp t264; +fp t582; +fp t354; +fp t692; +fp t472; +fp t643; +fp t397; +fp t691; +fp t447; +fp t168; +fp t480; +fp t494; +fp t690; +fp t404; +fp t689; +fp t503; +fp t382; +fp t443; +fp t459; +fp t619; +fp t688; +fp t632; +fp t195; +fp t687; +fp t446; +fp t622; +fp t182; +fp t458; +fp t499; +fp t379; +fp t686; +fp t492; +fp t393; +fp t685; +fp t636; +fp t684; +fp t383; +fp t454; +fp t436; +fp t634; +fp t683; +fp t464; +fp t330; +fp t478; +fp t621; +fp t402; +fp t682; +fp t414; +fp t372; +fp t487; +fp t681; +fp t440; +fp t449; +fp t626; +fp t496; +fp t680; +fp t427; +fp t633; +fp t403; +fp t679; +fp t678; +fp t391; +fp t624; +fp t677; +fp t676; +fp t355; +fp t675; +fp t387; +fp t185; +fp t674; +fp t183; +fp t362; +fp t673; +fp t672; +fp t366; +fp t671; +fp t670; +fp t629; +fp t669; +fp t668; +fp t667; +fp t500; +fp t360; +fp t644; +fp t666; +fp t333; +fp t665; +fp t664; +fp t356; +fp t663; +fp t662; +fp t661; +fp t332; +fp t660; +fp t659; +fp t450; +fp t658; +fp t439; +fp t484; +fp t348; +fp t657; +fp t406; +fp t340; +fp t656; +fp t409; +fp t655; +fp t653; +fp t431; +fp t650; +fp t435; +fp t648; +fp t434; +fp t647; +fp t491; +fp t646; +fp t475; +fp t645; +fp t417; +fp t637; +fp t428; +fp t635; +fp t627; +fp t442; +fp t625; +fp t411; +fp t623; +fp t441; +fp t620; +fp t479; +fp t618; +fp t452; +fp t617; +fp t239; +fp t614; +fp t342; +fp t343; +fp t609; +fp t605; +fp t603; +fp t602; +fp t600; +fp t369; +fp t599; +fp t598; +fp t597; +fp t596; +fp t595; +fp t594; +fp t593; +fp t592; +fp t588; +fp t587; +fp t585; +fp t584; +fp t381; +fp t583; +fp t581; +fp t579; +fp t578; +fp t575; +fp t574; +fp t573; +fp t572; +fp t571; +fp t570; +fp t569; +fp t568; +fp t567; +fp t566; +fp t565; +fp t564; +fp t563; +fp t562; +fp t561; +fp t560; +fp t559; +fp t558; +fp t486; +fp t557; +fp t556; +fp t555; +fp t554; +fp t553; +fp t552; +fp t551; +fp t550; +fp t549; +fp t548; +fp t547; +fp t546; +fp t545; +fp t544; +fp t432; +fp t197; +fp t543; +fp t542; +fp t541; +fp t540; +fp t539; +fp t538; +fp t537; +fp t536; +fp t474; +fp t535; +fp t534; +fp t533; +fp t532; +fp t457; +fp t531; +fp t530; +fp t529; +fp t528; +fp t421; +fp t527; +fp t526; +fp t525; +fp t394; +fp t524; +fp t523; +fp t522; +fp t521; +fp t520; +fp t444; +fp t519; +fp t518; +fp t517; +fp t516; +fp t515; +fp t504; +fp t485; +fp t477; +fp t471; +fp t468; +fp t465; +fp t456; +fp t455; +fp t430; +fp t424; +fp t422; +fp t420; +fp t419; +fp t416; +fp t413; +fp t412; +fp t408; +fp t400; +fp t399; +fp t398; +fp t395; +fp t390; +fp t388; +fp t386; +fp t199; +fp t181; +fp t179; + t652 = y*z; + t401 = RATIONAL(-607.0,12000.0); + t470 = RATIONAL(-1.0,2000.0); + t748 = t401*x+t470*t652; + t509 = z*z; + t489 = RATIONAL(1.0,160.0); + t341 = t489*t509; + t260 = x*t341; + t510 = x*x; + t506 = t510*x; + t498 = RATIONAL(-1.0,32.0); + t378 = t498*t506; + t490 = RATIONAL(-1.0,160.0); + t511 = y*y; + t346 = t490*t511; + t233 = x*t346; + t747 = t260+t378+t233; + t642 = x*t470; + t746 = t401*y+z*t642; + t448 = RATIONAL(-97.0,3000.0); + t451 = RATIONAL(-37.0,1000.0); + t638 = y*t451; + t745 = t448*x+z*t638; + t493 = RATIONAL(-3.0,160.0); + t607 = t493*t510; + t274 = y*t607; + t508 = t511*y; + t497 = RATIONAL(1.0,96.0); + t365 = t497*t508; + t744 = t274+t365; + t612 = t511*t493; + t265 = x*t612; + t375 = t497*t506; + t743 = t265+t375; + t350 = t490*t509; + t237 = y*t350; + t502 = RATIONAL(1.0,32.0); + t364 = t502*t508; + t742 = t237+t364; + t407 = RATIONAL(-137.0,4000.0); + t482 = RATIONAL(-9.0,160.0); + t339 = t482*t511; + t483 = RATIONAL(2.0,125.0); + t631 = z*t483; + t741 = x*t631+t339+t407*y; + t349 = t482*t510; + t740 = t349+y*t631+t407*x; + t433 = RATIONAL(-49.0,2000.0); + t495 = RATIONAL(-3.0,40.0); + t466 = RATIONAL(-7.0,2000.0); + t649 = y*t466; + t739 = t495*t509+t433*z+x*t649; + t738 = t433*x+z*t649+t495*t510; + t501 = RATIONAL(3.0,80.0); + t380 = t501*t511; + t654 = x*z; + t737 = t451*t654+t448*y+t380; + t384 = t502*t506; + t226 = x*t350; + t601 = t233+t226; + t736 = t384+t601; + t611 = t493*t509; + t256 = x*t611; + t735 = t256+t743; + t481 = RATIONAL(3.0,160.0); + t352 = t481*t509; + t276 = y*t352; + t734 = t276+t744; + t353 = t489*t510; + t232 = y*t353; + t733 = t232+t742; + t253 = z*t353; + t351 = t489*t511; + t230 = z*t351; + t507 = t509*z; + t363 = t498*t507; + t732 = t253+t230+t363; + t347 = t481*t510; + t249 = y*t347; + t586 = t249+t276; + t731 = t365+t586; + t467 = RATIONAL(-3.0,1000.0); + t206 = t467*t654; + t160 = y*t206; + t730 = t160+RATIONAL(49.0,4000.0); + t651 = x*y; + t616 = z*t651; + t164 = RATIONAL(1.0,1000.0)*t616; + t729 = RATIONAL(397.0,4000.0)+t164; + t157 = RATIONAL(27.0,1000.0)*t616; + t728 = t157+RATIONAL(-59.0,1000.0); + t505 = RATIONAL(-1.0,96.0); + t374 = t505*t508; + t727 = t274+t374; + t726 = t260+t384; + t392 = t505*t506; + t725 = t256+t392; + t245 = y*t341; + t724 = t245+t747; + t438 = RATIONAL(-87.0,2000.0); + t723 = t438*y+t206; + t418 = RATIONAL(183.0,4000.0); + t462 = RATIONAL(37.0,2000.0); + t722 = t462*t654+t418*y; + t376 = t502*t507; + t721 = t253+t376; + t261 = x*t351; + t610 = t260+t261; + t720 = t610+t378; + t719 = t265+t392; + t718 = t384+t261+t226; + t336 = t481*t511; + t251 = x*t336; + t717 = t251+t375; + t453 = RATIONAL(49.0,2000.0); + t476 = RATIONAL(-1.0,1000.0); + t639 = x*t476; + t214 = z*t639; + t163 = y*t214; + t716 = t453+t163; + t429 = RATIONAL(-31.0,1000.0); + t405 = RATIONAL(141.0,2000.0); + t715 = t429*t651+t405*z; + t258 = x*t352; + t590 = t251+t258; + t714 = t392+t590; + t367 = t498*t508; + t345 = t490*t510; + t243 = y*t345; + t713 = t367+t243; + t712 = t365+t249; + t469 = RATIONAL(-9.0,1000.0); + t473 = RATIONAL(-41.0,750.0); + t711 = t473*x+t469*t652; + t248 = z*t347; + t278 = z*t612; + t604 = t248+t278; + t389 = t505*t507; + t710 = t604+t389; + t262 = z*t345; + t263 = z*t346; + t709 = t262+t263+t376; + t707 = t245+t364; + t706 = t376+t230; + t270 = y*t611; + t705 = t249+t270+t374; + t161 = RATIONAL(9.0,1000.0)*t616; + t704 = RATIONAL(73.0,4000.0)+t161; + t463 = RATIONAL(-27.0,1000.0); + t194 = t463*t654; + t158 = y*t194; + t703 = t158+RATIONAL(-59.0,4000.0); + t640 = z*t429; + t702 = t405*x+y*t640; + t641 = x*t469; + t216 = z*t641; + t162 = y*t216; + t701 = RATIONAL(-59.0,2000.0)+t162; + t630 = y*t462; + t700 = t418*x+z*t630; + t159 = RATIONAL(3.0,1000.0)*t616; + t699 = t159+RATIONAL(49.0,1000.0); + t577 = t274+t270; + t698 = t365+t577; + t697 = y*t641+t473*z; + t488 = RATIONAL(-9.0,250.0); + t628 = z*t488; + t396 = RATIONAL(281.0,12000.0); + t696 = t396*y+x*t628; + t695 = t375+t258; + t694 = t263+t363; + t693 = t262+t226; + t264 = z*t607; + t582 = t278+t264; + t354 = t497*t507; + t692 = t582+t354; + t472 = RATIONAL(11.0,500.0); + t643 = x*t472; + t397 = RATIONAL(-167.0,12000.0); + t691 = t397*y+z*t643+t346; + t447 = RATIONAL(31.0,2000.0); + t168 = t447*t654; + t480 = RATIONAL(-4.0,375.0); + t494 = RATIONAL(-1.0,80.0); + t690 = t168+t480*y+t494*t511; + t404 = RATIONAL(53.0,12000.0); + t689 = t453*t652+t404*x+t347; + t503 = RATIONAL(1.0,40.0); + t382 = t503*t509; + t443 = RATIONAL(-29.0,6000.0); + t459 = RATIONAL(13.0,1000.0); + t619 = y*t459; + t688 = t382+x*t619+t443*z; + t632 = y*t447; + t195 = z*t632; + t687 = t480*x+t195+t494*t510; + t446 = RATIONAL(-19.0,2000.0); + t622 = y*t446; + t182 = z*t622; + t458 = RATIONAL(-17.0,3000.0); + t499 = RATIONAL(1.0,80.0); + t379 = t499*t510; + t686 = t182+t458*x+t379; + t492 = RATIONAL(-7.0,250.0); + t393 = RATIONAL(-137.0,12000.0); + t685 = t341+t393*z+t492*t651; + t636 = x*t453; + t684 = t336+z*t636+t404*y; + t383 = t501*t509; + t454 = RATIONAL(53.0,1000.0); + t436 = RATIONAL(19.0,1000.0); + t634 = y*t436; + t683 = t383+x*t634+t454*z; + t464 = RATIONAL(-11.0,160.0); + t330 = t464*t511; + t478 = RATIONAL(-17.0,500.0); + t621 = z*t478; + t402 = RATIONAL(-107.0,4000.0); + t682 = t402*y+t330+x*t621; + t414 = RATIONAL(299.0,4000.0); + t372 = t502*t509; + t487 = RATIONAL(-1.0,500.0); + t681 = t372+t414*z+t487*t651; + t440 = RATIONAL(-79.0,2000.0); + t449 = RATIONAL(43.0,2000.0); + t626 = y*t449; + t496 = RATIONAL(-1.0,20.0); + t680 = z*t626+t496*t510+t440*x; + t427 = RATIONAL(-13.0,2000.0); + t633 = x*t427; + t403 = RATIONAL(323.0,4000.0); + t679 = t352+t403*z+y*t633; + t678 = t454*y+t436*t654+t380; + t391 = t502*t510; + t624 = z*t487; + t677 = t414*x+t391+y*t624; + t676 = t352+y*t636+t404*z; + t355 = t503*t510; + t675 = t355+t443*x+z*t619; + t387 = t499*t509; + t185 = x*t622; + t674 = t387+t458*z+t185; + t183 = t446*t654; + t362 = t499*t511; + t673 = t458*y+t183+t362; + t672 = t496*t511+t440*y+t449*t654; + t366 = t501*t510; + t671 = t366+t454*x+z*t634; + t670 = t427*t652+t403*x+t347; + t629 = z*t492; + t669 = t393*x+y*t629+t353; + t668 = t336+t403*y+z*t633; + t667 = x*t626+t496*t509+t440*z; + t500 = RATIONAL(-1.0,16.0); + t360 = t500*t509; + t644 = y*t467; + t666 = t360+t438*z+x*t644; + t333 = t464*t509; + t665 = t402*z+t478*t651+t333; + t664 = t351+x*t629+t393*y; + t356 = t503*t511; + t663 = t443*y+t356+t459*t654; + t662 = t396*x+t391+y*t628; + t661 = t472*t652+t345+t397*x; + t332 = t464*t510; + t660 = t402*x+y*t621+t332; + t659 = t382+t355+t356+RATIONAL(-61.0,4000.0); + t450 = RATIONAL(-97.0,4000.0); + t658 = t448*z+t450+t161+x*t638; + t439 = RATIONAL(71.0,2000.0); + t484 = RATIONAL(7.0,160.0); + t348 = t484*t509; + t657 = y*t642+t401*z+t439+t348+t162; + t406 = RATIONAL(143.0,2000.0); + t340 = t482*t509; + t656 = t406+t159+t340+t407*z+t483*t651; + t409 = RATIONAL(469.0,4000.0); + t655 = t466*t654+t409+t160+t495*t511+t433*y; + t653 = y+z; + t431 = RATIONAL(-17.0,2000.0); + t650 = y*t431; + t435 = RATIONAL(-23.0,2000.0); + t648 = x*t435; + t434 = RATIONAL(-39.0,2000.0); + t647 = y*t434; + t491 = RATIONAL(-9.0,500.0); + t646 = y*t491; + t475 = RATIONAL(-3.0,2000.0); + t645 = t475*z; + t417 = RATIONAL(111.0,2000.0); + t637 = y*t417; + t428 = RATIONAL(33.0,2000.0); + t635 = x*t428; + t627 = y*t463; + t442 = RATIONAL(93.0,2000.0); + t625 = x*t442; + t411 = RATIONAL(129.0,2000.0); + t623 = x*t411; + t441 = RATIONAL(-21.0,2000.0); + t620 = t441*y; + t479 = RATIONAL(-1.0,250.0); + t618 = y*t479; + t452 = RATIONAL(27.0,2000.0); + t617 = x*t452; + t239 = z*t336; + t614 = t239+t264; + t342 = t484*t510; + t343 = t484*t511; + t609 = t342+t343; + t605 = t232+t245; + t603 = t237+t243; + t602 = t239+t265; + t600 = t256+t276; + t369 = t502*t511; + t599 = t369+t367; + t598 = t261+t263; + t597 = t369+t374; + t596 = t262+t232; + t595 = t372+t389; + t594 = t270+t258; + t593 = t253+t243; + t592 = t230+t233; + t588 = t256+t270; + t587 = t248+t239; + t585 = t384+t500*t510; + t584 = t253+t232; + t381 = t500*t511; + t583 = t364+t381; + t581 = t263+t233; + t579 = t248+t243; + t578 = t251+t278; + t575 = t383+t366+t367; + t574 = t160+t374+t586; + t573 = t161+t732; + t572 = t159+t692; + t571 = t712+t719; + t570 = t367+t232+t743; + t569 = t262+t363+t726; + t568 = t248+t354+t725; + t567 = t706+t707; + t566 = t717+t727; + t565 = t162+t709; + t564 = t378+t226+t721; + t563 = t264+t389+t695; + t562 = t719+t727; + t561 = t694+t707; + t560 = t237+t367+t706; + t559 = t261+t378+t744; + t558 = t378+t363+t693; + t486 = RATIONAL(7.0,250.0); + t557 = t350+t397*z+t486+y*t643+t354+t162; + t556 = t571+t594; + t555 = t569+t592; + t554 = t354+t614+t705; + t553 = t710+t734; + t552 = t568+t602; + t551 = t251+t389+t614+t725; + t550 = t560+t596; + t549 = t718+t733; + t548 = t564+t598; + t547 = t563+t602; + t546 = t354+t375+t587+t590; + t545 = t262+t243+t567; + t544 = t564+t592; + t432 = RATIONAL(-67.0,4000.0); + t197 = x*t632; + t543 = t354+t494*t509+t480*z+t161+t432+t197; + t542 = t245+t713+t718; + t541 = t588+t712+t717; + t540 = t561+t593; + t539 = t733+t747; + t538 = t705+t710; + t537 = t558+t581; + t536 = t276+t258+t562; + t474 = RATIONAL(-23.0,500.0); + t535 = t474+t159+x*t630+t418*z+t363+t348; + t534 = t581+t721+t726; + t533 = t566+t600; + t532 = t568+t578; + t457 = RATIONAL(79.0,3000.0); + t531 = t457*z+y*t639+t692+t730; + t530 = t457*y+t214+t698+t730; + t529 = t473*y+t216+t703+t731; + t528 = t457*x+t476*t652+t730+t735; + t421 = RATIONAL(191.0,4000.0); + t527 = t372+t421*z+x*t646+t701+t732; + t526 = t438*x+z*t644+t585+t610+t729; + t525 = x*t640+t160+t367+t405*y+t659; + t394 = RATIONAL(-691.0,12000.0); + t524 = t479*t654+t394*y+t586+t597+t699; + t523 = t491*t654+t421*y+t599+t605+t701; + t522 = t394*x+t391+z*t618+t699+t714; + t521 = x*t618+t394*z+t587+t595+t699; + t520 = t396*z+t488*t651+t582+t595+t728; + t444 = RATIONAL(-33.0,2000.0); + t519 = t360+t444*z+x*t627+t704+t709; + t518 = t444*x+z*t627+t585+t601+t704; + t517 = z*t646+t391+t421*x+t701+t720; + t516 = t194+t444*y+t583+t603+t704; + t515 = t414*y+x*t624+t599+t603+t716; + t504 = RATIONAL(1.0,20.0); + t485 = RATIONAL(-7.0,160.0); + t477 = RATIONAL(-13.0,160.0); + t471 = RATIONAL(21.0,250.0); + t468 = RATIONAL(13.0,500.0); + t465 = RATIONAL(-17.0,375.0); + t456 = RATIONAL(-39.0,4000.0); + t455 = RATIONAL(-47.0,1500.0); + t430 = RATIONAL(29.0,2000.0); + t424 = RATIONAL(-57.0,2000.0); + t422 = RATIONAL(133.0,2000.0); + t420 = RATIONAL(147.0,2000.0); + t419 = RATIONAL(307.0,6000.0); + t416 = RATIONAL(157.0,2000.0); + t413 = RATIONAL(-129.0,4000.0); + t412 = RATIONAL(-191.0,4000.0); + t408 = RATIONAL(-437.0,6000.0); + t400 = RATIONAL(319.0,12000.0); + t399 = RATIONAL(-539.0,12000.0); + t398 = RATIONAL(-449.0,12000.0); + t395 = RATIONAL(349.0,12000.0); + t390 = t504*t510; + t388 = t504*t511; + t386 = t504*t509; + t199 = t441*t654; + t181 = x*t620; + t179 = z*t620; + coeffs_I->coeff_m1_m1_m1 = RATIONAL(-269.0,4000.0)+t390+t158+t386+t389+ +t455*z+t388+(t420*z+t455)*y+(t455+t653*t420)*x+t562+t582+t588; + coeffs_I->coeff_0_m1_m1 = t477*t510+t456*x+t411*t652+t538+t343+t657+t736+ +t746; + coeffs_I->coeff_p1_m1_m1 = t383+t416*x+z*t637+t538+t379+t658+t720+t737; + coeffs_I->coeff_p2_m1_m1 = t442*t652+t369+t398*x+t520+t347+t566+t594+t696 +; + coeffs_I->coeff_m1_0_m1 = t477*t511+t456*y+z*t623+t342+t551+t364+t603+ +t657+t748; + coeffs_I->coeff_0_0_m1 = t383+t408*z+t549+t389+x*t650+t587+t655+t738; + coeffs_I->coeff_p1_0_m1 = t521+t539+t670+t682; + coeffs_I->coeff_p2_0_m1 = t516+t547+t686+t688; + coeffs_I->coeff_m1_p1_m1 = t416*y+t417*t654+t362+t551+t575+t605+t658+t745 +; + coeffs_I->coeff_0_p1_m1 = t521+t542+t660+t668; + coeffs_I->coeff_p1_p1_m1 = t239+t525+y*t635+t389+t465*z+t579+t702+t724; + coeffs_I->coeff_p2_p1_m1 = t523+t547+t669+t676; + coeffs_I->coeff_m1_p2_m1 = z*t625+t398*y+t336+t520+t571+t600+t662; + coeffs_I->coeff_0_p2_m1 = t518+t553+t673+t688; + coeffs_I->coeff_p1_p2_m1 = t517+t553+t664+t676; + coeffs_I->coeff_p2_p2_m1 = t387+t475*t651+t419*z+t529+t563+t578+t711; + coeffs_I->coeff_m1_m1_0 = t477*t509+t536+t456*z+y*t623+t439+t565+t609+ +t746+t748; + coeffs_I->coeff_0_m1_0 = t534+t408*y+t409+t380+t431*t654+t574+t738+t739; + coeffs_I->coeff_p1_m1_0 = t524+t548+t665+t670; + coeffs_I->coeff_p2_m1_0 = t519+t533+t663+t686; + coeffs_I->coeff_m1_0_0 = t545+t366+t408*x+z*t650+t655+t714+t739; + coeffs_I->coeff_0_0_0 = RATIONAL(177.0,1000.0)+t163+t384+t183+t333+t332+ +t185+t182+t330+(x+t653)*RATIONAL(-161.0,4000.0)+t567+t584+t610; + coeffs_I->coeff_p1_0_0 = t179+t355+t544+t471*x+t583+t605+t666+t723+t729; + coeffs_I->coeff_p2_0_0 = t545+t353+t435*t652+t395*x+t656+t735+t741; + coeffs_I->coeff_m1_p1_0 = t522+t550+t665+t668; + coeffs_I->coeff_0_p1_0 = t356+t526+t471*y+t199+t560+t593+t666; + coeffs_I->coeff_p1_p1_0 = t412*z+t340+t515+t197+t544+t677; + coeffs_I->coeff_p2_p1_0 = t528+t550+t667+t678; + coeffs_I->coeff_m1_p2_0 = t519+t556+t673+t675; + coeffs_I->coeff_0_p2_0 = t534+t351+z*t648+t395*y+t656+t698+t740; + coeffs_I->coeff_p1_p2_0 = t530+t548+t667+t671; + coeffs_I->coeff_p2_p2_0 = t413*z+t485*t509+t181+t486+t541+t565+t661+t691; + coeffs_I->coeff_m1_m1_p1 = t366+t450+t536+t387+t416*z+x*t637+t573+t737+ +t745; + coeffs_I->coeff_0_m1_p1 = t524+t555+t660+t679; + coeffs_I->coeff_p1_m1_p1 = t261+t558+z*t635+t230+t465*y+t574+t659+t702+ +t715; + coeffs_I->coeff_p2_m1_p1 = t527+t533+t669+t684; + coeffs_I->coeff_m1_0_p1 = t522+t540+t679+t682; + coeffs_I->coeff_0_0_p1 = t382+t181+t526+t381+t471*z+t561+t596+t723; + coeffs_I->coeff_p1_0_p1 = t168+t412*y+t537+t339+t364+t605+t677+t681+t716; + coeffs_I->coeff_p2_0_p1 = t528+t540+t672+t683; + coeffs_I->coeff_m1_p1_p1 = t237+t525+t428*t652+t465*x+t584+t694+t714+t715 +; + coeffs_I->coeff_0_p1_p1 = t349+t412*x+t515+t195+t569+t598+t681; + coeffs_I->coeff_p1_p1_p1 = RATIONAL(-193.0,4000.0)+t422*z+t164+t380+(t430 +*z+t422)*y+(t422+t653*t430)*x+t537+t575+t603; + coeffs_I->coeff_p2_p1_p1 = t237+t263+t345+t253+t535+t256+t452*t652+t343+ +t400*x+t570+t722; + coeffs_I->coeff_m1_p2_p1 = t527+t556+t664+t689; + coeffs_I->coeff_0_p2_p1 = t530+t555+t680+t683; + coeffs_I->coeff_p1_p2_p1 = t535+t270+z*t617+t342+t400*y+t346+t230+t559+ +t693+t700; + coeffs_I->coeff_p2_p2_p1 = t386+t468*z+x*t647+t432+t541+t573+t687+t690; + coeffs_I->coeff_m1_m1_p2 = t352+y*t625+t398*z+t552+t577+t597+t662+t696+ +t728; + coeffs_I->coeff_0_m1_p2 = t518+t554+t663+t674; + coeffs_I->coeff_p1_m1_p2 = t517+t554+t684+t685; + coeffs_I->coeff_p2_m1_p2 = t419*y+x*t645+t546+t374+t362+t577+t697+t703+ +t711; + coeffs_I->coeff_m1_0_p2 = t516+t532+t674+t675; + coeffs_I->coeff_0_0_p2 = t406+t341+t549+y*t648+t395*z+t572+t740+t741; + coeffs_I->coeff_p1_0_p2 = t531+t539+t671+t672; + coeffs_I->coeff_p2_0_p2 = t557+t413*y+t485*t511+t265+t278+t199+t579+t661+ +t695+t742; + coeffs_I->coeff_m1_p1_p2 = t523+t532+t685+t689; + coeffs_I->coeff_0_p1_p2 = t531+t542+t678+t680; + coeffs_I->coeff_p1_p1_p2 = t474+t350+y*t617+t400*z+t572+t609+t700+t713+ +t722+t724; + coeffs_I->coeff_p2_p1_p2 = t543+t468*y+t388+t434*t654+t258+t245+t570+t604 ++t687; + coeffs_I->coeff_m1_p2_p2 = t419*x+y*t645+t379+t529+t552+t697; + coeffs_I->coeff_0_p2_p2 = t179+t413*x+t485*t510+t557+t614+t691+t734+t736; + coeffs_I->coeff_p1_p2_p2 = t260+t543+t390+t468*x+z*t647+t276+t559+t614+ +t690; + coeffs_I->coeff_p2_p2_p2 = RATIONAL(77.0,2000.0)+t399*z+t157+(t424*z+t399 +)*y+(t509+t510+t511)*t493+(t399+t653*t424)*x+t546+t731; diff --git a/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dx.compute.c b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..1753d9b --- /dev/null +++ b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dx.compute.c @@ -0,0 +1,518 @@ +fp t911; +fp t898; +fp t808; +fp t913; +fp t809; +fp t1006; +fp t912; +fp t826; +fp t906; +fp t804; +fp t962; +fp t1005; +fp t840; +fp t1004; +fp t888; +fp t797; +fp t755; +fp t1003; +fp t829; +fp t1002; +fp t903; +fp t812; +fp t803; +fp t837; +fp t1001; +fp t834; +fp t805; +fp t963; +fp t1000; +fp t999; +fp t998; +fp t997; +fp t833; +fp t996; +fp t835; +fp t995; +fp t832; +fp t759; +fp t994; +fp t993; +fp t902; +fp t816; +fp t958; +fp t813; +fp t992; +fp t838; +fp t961; +fp t991; +fp t827; +fp t990; +fp t989; +fp t822; +fp t988; +fp t987; +fp t986; +fp t985; +fp t984; +fp t983; +fp t982; +fp t886; +fp t783; +fp t889; +fp t799; +fp t894; +fp t793; +fp t758; +fp t760; +fp t981; +fp t757; +fp t762; +fp t980; +fp t979; +fp t978; +fp t761; +fp t764; +fp t976; +fp t763; +fp t975; +fp t974; +fp t973; +fp t971; +fp t970; +fp t969; +fp t967; +fp t966; +fp t965; +fp t964; +fp t960; +fp t957; +fp t955; +fp t954; +fp t751; +fp t953; +fp t756; +fp t952; +fp t951; +fp t950; +fp t830; +fp t949; +fp t752; +fp t948; +fp t947; +fp t946; +fp t836; +fp t945; +fp t944; +fp t943; +fp t753; +fp t942; +fp t941; +fp t940; +fp t939; +fp t938; +fp t749; +fp t937; +fp t754; +fp t936; +fp t935; +fp t934; +fp t933; +fp t932; +fp t750; +fp t931; +fp t930; +fp t929; +fp t928; +fp t927; +fp t926; +fp t925; +fp t924; +fp t923; +fp t922; +fp t921; +fp t920; +fp t841; +fp t919; +fp t918; +fp t824; +fp t917; +fp t916; +fp t915; +fp t914; +fp t905; +fp t904; +fp t901; +fp t900; +fp t899; +fp t896; +fp t895; +fp t893; +fp t891; +fp t890; +fp t887; +fp t885; +fp t884; +fp t880; +fp t879; +fp t878; +fp t877; +fp t875; +fp t874; +fp t871; +fp t870; +fp t867; +fp t866; +fp t865; +fp t863; +fp t861; +fp t860; +fp t859; +fp t855; +fp t854; +fp t852; +fp t849; +fp t848; +fp t846; +fp t823; +fp t821; +fp t820; +fp t819; +fp t818; +fp t817; +fp t815; +fp t814; +fp t811; +fp t810; +fp t807; +fp t806; +fp t802; +fp t801; +fp t800; +fp t798; +fp t796; +fp t795; +fp t794; +fp t792; +fp t791; +fp t790; +fp t789; +fp t788; +fp t787; +fp t786; +fp t785; +fp t784; +fp t782; +fp t781; +fp t780; +fp t779; +fp t778; +fp t777; +fp t776; +fp t775; +fp t774; +fp t773; +fp t772; +fp t771; +fp t770; +fp t769; +fp t768; +fp t767; +fp t766; +fp t765; + t911 = z*z; + t898 = RATIONAL(3.0,160.0); + t808 = t898*t911; + t913 = y*y; + t809 = t898*t913; + t1006 = t808+t809; + t912 = x*x; + t826 = RATIONAL(-1.0,32.0)*t912; + t906 = RATIONAL(-3.0,160.0); + t804 = t906*t913; + t962 = t808+t804; + t1005 = t826+t962; + t840 = RATIONAL(3.0,40.0)*x; + t1004 = t840+RATIONAL(-97.0,3000.0); + t888 = RATIONAL(-3.0,1000.0); + t797 = t888*y; + t755 = z*t797; + t1003 = RATIONAL(53.0,1000.0)+t840+t755; + t829 = RATIONAL(-1.0,80.0)*x; + t1002 = RATIONAL(-167.0,12000.0)+t829; + t903 = RATIONAL(-1.0,160.0); + t812 = t903*t913; + t803 = t903*t911; + t837 = RATIONAL(-3.0,32.0)*t912; + t1001 = t812+t837+t803; + t834 = RATIONAL(1.0,32.0)*t912; + t805 = t906*t911; + t963 = t804+t805; + t1000 = t834+t963; + t999 = t834+t1006; + t998 = RATIONAL(-79.0,2000.0)+RATIONAL(-1.0,10.0)*x; + t997 = t834+t962; + t833 = RATIONAL(-9.0,80.0)*x; + t996 = RATIONAL(-137.0,4000.0)+t833; + t835 = RATIONAL(1.0,40.0)*x; + t995 = RATIONAL(-17.0,3000.0)+t835; + t832 = RATIONAL(3.0,80.0)*x; + t759 = y*t832; + t994 = t759+t832; + t993 = RATIONAL(-49.0,2000.0)+RATIONAL(-3.0,20.0)*x; + t902 = RATIONAL(1.0,160.0); + t816 = t902*t913; + t958 = t816+t837; + t813 = t902*t911; + t992 = t958+t813; + t838 = RATIONAL(3.0,32.0)*t912; + t961 = t803+t838; + t991 = t812+t961; + t827 = RATIONAL(1.0,80.0)*x; + t990 = RATIONAL(-137.0,12000.0)+t827; + t989 = t826+t1006; + t822 = RATIONAL(-11.0,80.0)*x; + t988 = RATIONAL(-107.0,4000.0)+t822; + t987 = RATIONAL(-4.0,375.0)+RATIONAL(-1.0,40.0)*x; + t986 = t816+t813+t838; + t985 = t826+t963; + t984 = y+z; + t983 = y*z; + t982 = RATIONAL(-3.0,80.0)*x; + t886 = RATIONAL(-27.0,1000.0); + t783 = t886*y; + t889 = RATIONAL(-9.0,1000.0); + t799 = t889*y; + t894 = RATIONAL(-1.0,1000.0); + t793 = t894*z; + t758 = z*t832; + t760 = y*t829; + t981 = t758+t760; + t757 = y*t982; + t762 = z*t827; + t980 = t757+t762; + t979 = t816+t803; + t978 = t809+t805; + t761 = y*t827; + t764 = z*t982; + t976 = t761+t764; + t763 = z*t829; + t975 = t760+t763; + t974 = t758+t761; + t973 = t757+t764; + t971 = t757+t763; + t970 = t759+t763; + t969 = t759+t762; + t967 = t760+t764; + t966 = t761+t763; + t965 = t758+t759; + t964 = t760+t762; + t960 = t759+t764; + t957 = t761+t762; + t955 = t757+t758; + t954 = t812+t813; + t751 = z*t799; + t953 = t751+t978; + t756 = RATIONAL(3.0,1000.0)*t983; + t952 = t756+t954; + t951 = RATIONAL(1.0,10.0)*x+t973; + t950 = t832+t974; + t830 = RATIONAL(7.0,80.0)*x; + t949 = t830+t971; + t752 = RATIONAL(9.0,1000.0)*t983; + t948 = t752+t978; + t947 = t840+t975; + t946 = t756+t979; + t836 = RATIONAL(1.0,20.0)*x; + t945 = t836+t957; + t944 = t835+t965; + t943 = t751+t991; + t753 = y*t793; + t942 = t753+t986; + t941 = t838+t952; + t940 = t755+t989; + t939 = t762+t803+t958; + t938 = t834+t953; + t749 = z*t783; + t937 = t749+t985; + t754 = RATIONAL(1.0,1000.0)*t983; + t936 = t754+t1001; + t935 = t764+t997; + t934 = t762+t1005; + t933 = t758+t997; + t932 = t755+t838+t954; + t750 = RATIONAL(27.0,1000.0)*t983; + t931 = t750+t999; + t930 = t755+t816+t961; + t929 = t756+t1000; + t928 = t752+t826+t836+RATIONAL(-29.0,6000.0); + t927 = t755+RATIONAL(141.0,2000.0)+t836+t837; + t926 = t764+t837+t954; + t925 = t826+t953; + t924 = t834+t948; + t923 = t752+t992; + t922 = t838+t946; + t921 = t755+RATIONAL(79.0,3000.0)+t1000; + t920 = RATIONAL(-41.0,750.0)+t749+t999; + t841 = RATIONAL(-1.0,8.0)*x; + t919 = t752+RATIONAL(-33.0,2000.0)+t841+t991; + t918 = t841+RATIONAL(-87.0,2000.0)+t754+t986; + t824 = RATIONAL(1.0,16.0)*x; + t917 = t756+RATIONAL(-691.0,12000.0)+t824+t989; + t916 = t753+RATIONAL(299.0,4000.0)+t824+t1001; + t915 = RATIONAL(281.0,12000.0)+t750+t824+t985; + t914 = t751+t824+RATIONAL(191.0,4000.0)+t992; + t905 = RATIONAL(-7.0,250.0); + t904 = RATIONAL(-9.0,500.0); + t901 = RATIONAL(-9.0,250.0); + t900 = RATIONAL(-1.0,500.0); + t899 = RATIONAL(2.0,125.0); + t896 = RATIONAL(-1.0,250.0); + t895 = RATIONAL(-17.0,500.0); + t893 = RATIONAL(-3.0,2000.0); + t891 = RATIONAL(11.0,500.0); + t890 = RATIONAL(-1.0,2000.0); + t887 = RATIONAL(-7.0,2000.0); + t885 = RATIONAL(37.0,2000.0); + t884 = RATIONAL(13.0,1000.0); + t880 = RATIONAL(49.0,2000.0); + t879 = RATIONAL(27.0,2000.0); + t878 = RATIONAL(-37.0,1000.0); + t877 = RATIONAL(43.0,2000.0); + t875 = RATIONAL(31.0,2000.0); + t874 = RATIONAL(-19.0,2000.0); + t871 = RATIONAL(93.0,2000.0); + t870 = RATIONAL(-21.0,2000.0); + t867 = RATIONAL(19.0,1000.0); + t866 = RATIONAL(-23.0,2000.0); + t865 = RATIONAL(-39.0,2000.0); + t863 = RATIONAL(-17.0,2000.0); + t861 = RATIONAL(-31.0,1000.0); + t860 = RATIONAL(33.0,2000.0); + t859 = RATIONAL(-13.0,2000.0); + t855 = RATIONAL(183.0,4000.0); + t854 = RATIONAL(111.0,2000.0); + t852 = RATIONAL(129.0,2000.0); + t849 = RATIONAL(53.0,12000.0); + t848 = RATIONAL(323.0,4000.0); + t846 = RATIONAL(-607.0,12000.0); + t823 = t900*y; + t821 = t896*z; + t820 = t905*y; + t819 = t904*y; + t818 = t899*z; + t817 = t899*y; + t815 = t905*z; + t814 = t900*z; + t811 = t901*y; + t810 = t901*z; + t807 = t904*z; + t806 = t896*y; + t802 = t887*y; + t801 = t895*y; + t800 = t890*z; + t798 = t888*z; + t796 = t894*y; + t795 = t891*y; + t794 = t889*z; + t792 = t891*z; + t791 = t887*z; + t790 = t895*z; + t789 = t890*y; + t788 = t875*z; + t787 = t885*z; + t786 = t859*z; + t785 = t885*y; + t784 = t877*z; + t782 = t884*z; + t781 = t861*y; + t780 = t874*y; + t779 = t874*z; + t778 = t867*y; + t777 = t861*z; + t776 = t880*z; + t775 = t870*z; + t774 = t878*y; + t773 = t875*y; + t772 = t877*y; + t771 = t878*z; + t770 = t867*z; + t769 = t870*y; + t768 = t880*y; + t767 = t884*y; + t766 = t886*z; + t765 = t859*y; + coeffs_dx->coeff_m1_m1_m1 = RATIONAL(-47.0,1500.0)+t984*RATIONAL(147.0, +2000.0)+t937+t951; + coeffs_dx->coeff_0_m1_m1 = RATIONAL(-39.0,4000.0)+t789+t800+RATIONAL( +-13.0,80.0)*x+t943+t965; + coeffs_dx->coeff_p1_m1_m1 = RATIONAL(157.0,2000.0)+t774+t771+t923+t944; + coeffs_dx->coeff_p2_m1_m1 = t811+t832+t810+RATIONAL(-449.0,12000.0)+t931+ +t973; + coeffs_dx->coeff_m1_0_m1 = t830+t846+t852*z+t789+t925+t967; + coeffs_dx->coeff_0_0_m1 = t863*y+t791+t930+t974+t993; + coeffs_dx->coeff_p1_0_m1 = t837+t790+t806+t848+t950+t952; + coeffs_dx->coeff_p2_0_m1 = t766+t760+t767+t752+t935+t995; + coeffs_dx->coeff_m1_p1_m1 = t854*z+t774+t826+t948+t976+t1004; + coeffs_dx->coeff_0_p1_m1 = t786+t806+t922+t981+t988; + coeffs_dx->coeff_p1_p1_m1 = t777+t860*y+t927+t954+t981; + coeffs_dx->coeff_p2_p1_m1 = t761+t807+t751+t768+t935+t990; + coeffs_dx->coeff_m1_p2_m1 = t811+t871*z+t915+t960; + coeffs_dx->coeff_0_p2_m1 = t767+t779+t919+t955; + coeffs_dx->coeff_p1_p2_m1 = t815+t768+t914+t955; + coeffs_dx->coeff_p2_p2_m1 = t794+t893*y+t920+t960; + coeffs_dx->coeff_m1_m1_0 = t751+t852*y+t846+t800+t949+t1005; + coeffs_dx->coeff_0_m1_0 = t802+t863*z+t932+t969+t993; + coeffs_dx->coeff_p1_m1_0 = t821+t756+t801+t848+t939+t994; + coeffs_dx->coeff_p2_m1_0 = t783+t782+t924+t971+t995; + coeffs_dx->coeff_m1_0_0 = t791+RATIONAL(-437.0,6000.0)+t802+t940+t947; + coeffs_dx->coeff_0_0_0 = RATIONAL(-161.0,4000.0)+t822+t779+t780+t942+t957 +; + coeffs_dx->coeff_p1_0_0 = RATIONAL(21.0,250.0)+t797+t798+t936+t945; + coeffs_dx->coeff_p2_0_0 = t827+t817+RATIONAL(349.0,12000.0)+t818+t929+ +t975; + coeffs_dx->coeff_m1_p1_0 = t786+t801+t917+t966; + coeffs_dx->coeff_0_p1_0 = t775+t797+t918+t964; + coeffs_dx->coeff_p1_p1_0 = t773+t814+t916+t964; + coeffs_dx->coeff_p2_p1_0 = t772+t770+t921+t966; + coeffs_dx->coeff_m1_p2_0 = t783+t779+t928+t962+t970; + coeffs_dx->coeff_0_p2_0 = t817+t866*z+t941+t980+t996; + coeffs_dx->coeff_p1_p2_0 = t772+t757+t793+t939+t1003; + coeffs_dx->coeff_p2_p2_0 = t792+t769+t938+t970+t1002; + coeffs_dx->coeff_m1_m1_p1 = t752+t757+t854*y+t771+t934+t1004; + coeffs_dx->coeff_0_m1_p1 = t765+t821+t941+t970+t988; + coeffs_dx->coeff_p1_m1_p1 = t860*z+t781+t927+t970+t979; + coeffs_dx->coeff_p2_m1_p1 = t819+t776+t938+t980+t990; + coeffs_dx->coeff_m1_0_p1 = t790+t765+t917+t964; + coeffs_dx->coeff_0_0_p1 = t769+t798+t918+t966; + coeffs_dx->coeff_p1_0_p1 = t823+t788+t916+t966; + coeffs_dx->coeff_p2_0_p1 = t778+t784+t921+t964; + coeffs_dx->coeff_m1_p1_p1 = RATIONAL(-17.0,375.0)+t777+t781+t940+t945; + coeffs_dx->coeff_0_p1_p1 = RATIONAL(-191.0,4000.0)+t833+t814+t823+t942+ +t975; + coeffs_dx->coeff_p1_p1_p1 = RATIONAL(133.0,2000.0)+t984*RATIONAL(29.0, +2000.0)+t936+t947; + coeffs_dx->coeff_p2_p1_p1 = RATIONAL(319.0,12000.0)+t829+t787+t785+t929+ +t957; + coeffs_dx->coeff_m1_p2_p1 = t849+t819+t751+t815+t934+t994; + coeffs_dx->coeff_0_p2_p1 = t793+t778+t932+t971+t998; + coeffs_dx->coeff_p1_p2_p1 = t855+t837+t879*z+t785+t946+t949; + coeffs_dx->coeff_p2_p2_p1 = t865*y+t788+t924+t969+t987; + coeffs_dx->coeff_m1_m1_p2 = t810+t871*y+t915+t955; + coeffs_dx->coeff_0_m1_p2 = t780+t782+t919+t960; + coeffs_dx->coeff_p1_m1_p2 = t776+t820+t914+t960; + coeffs_dx->coeff_p2_m1_p2 = t799+t893*z+t920+t955; + coeffs_dx->coeff_m1_0_p2 = t780+t766+t928+t978+t981; + coeffs_dx->coeff_0_0_p2 = t866*y+t818+t922+t976+t996; + coeffs_dx->coeff_p1_0_p2 = t796+t761+t784+t926+t1003; + coeffs_dx->coeff_p2_0_p2 = t775+t795+t760+t751+t933+t1002; + coeffs_dx->coeff_m1_p1_p2 = t849+t820+t807+t925+t950; + coeffs_dx->coeff_0_p1_p2 = t770+t796+t930+t967+t998; + coeffs_dx->coeff_p1_p1_p2 = t879*y+t760+t855+t830+t756+t787+t926; + coeffs_dx->coeff_p2_p1_p2 = t761+t752+t865*z+t773+t933+t987; + coeffs_dx->coeff_m1_p2_p2 = t794+RATIONAL(307.0,6000.0)+t799+t937+t944; + coeffs_dx->coeff_0_p2_p2 = t792+t795+RATIONAL(-7.0,80.0)*x+RATIONAL( +-129.0,4000.0)+t943+t973; + coeffs_dx->coeff_p1_p2_p2 = RATIONAL(13.0,500.0)+t788+t773+t923+t951; + coeffs_dx->coeff_p2_p2_p2 = RATIONAL(-539.0,12000.0)+t982+t984*RATIONAL( +-57.0,2000.0)+t931+t965; diff --git a/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxx.compute.c b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..cd2e06d --- /dev/null +++ b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,182 @@ +fp t1543; +fp t1525; +fp t1539; +fp t1578; +fp t1531; +fp t1519; +fp t1529; +fp t1577; +fp t1518; +fp t1538; +fp t1527; +fp t1576; +fp t1524; +fp t1528; +fp t1575; +fp t1574; +fp t1573; +fp t1521; +fp t1572; +fp t1571; +fp t1540; +fp t1520; +fp t1570; +fp t1569; +fp t1568; +fp t1567; +fp t1526; +fp t1566; +fp t1565; +fp t1564; +fp t1563; +fp t1562; +fp t1533; +fp t1561; +fp t1544; +fp t1560; +fp t1541; +fp t1559; +fp t1558; +fp t1557; +fp t1556; +fp t1555; +fp t1554; +fp t1553; +fp t1552; +fp t1551; +fp t1550; +fp t1535; +fp t1549; +fp t1548; +fp t1547; +fp t1546; +fp t1545; +fp t1542; +fp t1537; +fp t1534; +fp t1530; +fp t1523; +fp t1522; + t1543 = RATIONAL(1.0,16.0); + t1525 = t1543*x; + t1539 = RATIONAL(1.0,80.0); + t1578 = t1525+t1539; + t1531 = RATIONAL(-3.0,80.0); + t1519 = t1531*y; + t1529 = t1531*z; + t1577 = t1519+t1529; + t1518 = t1539*y; + t1538 = RATIONAL(-1.0,80.0); + t1527 = t1538*z; + t1576 = t1518+t1527; + t1524 = t1538*y; + t1528 = t1539*z; + t1575 = t1524+t1528; + t1574 = t1518+t1528; + t1573 = t1518+t1529; + t1521 = RATIONAL(3.0,16.0)*x; + t1572 = t1518+t1521; + t1571 = t1521+t1528; + t1540 = RATIONAL(3.0,80.0); + t1520 = t1540*y; + t1570 = t1520+t1525; + t1569 = t1524+t1529; + t1568 = t1521+RATIONAL(-1.0,10.0); + t1567 = t1519+t1527; + t1526 = t1540*z; + t1566 = t1525+t1526; + t1565 = t1520+t1521; + t1564 = t1524+t1527; + t1563 = t1521+t1526; + t1562 = t1519+t1525; + t1533 = RATIONAL(3.0,40.0); + t1561 = t1519+t1528+t1533; + t1544 = RATIONAL(1.0,20.0); + t1560 = t1544+t1574; + t1541 = RATIONAL(1.0,40.0); + t1559 = t1520+t1526+t1541; + t1558 = t1520+t1527+t1544; + t1557 = RATIONAL(1.0,10.0)+t1577; + t1556 = t1524+t1526+t1544; + t1555 = t1519+t1526+t1543; + t1554 = t1533+t1573; + t1553 = t1518+t1526+t1540; + t1552 = t1520+t1528+t1540; + t1551 = t1520+t1529+t1543; + t1550 = t1543+t1575; + t1535 = RATIONAL(7.0,80.0); + t1549 = t1535+t1569; + t1548 = t1533+t1564; + t1547 = t1543+t1576; + t1546 = t1535+t1567; + t1545 = RATIONAL(-1.0,8.0); + t1542 = RATIONAL(-1.0,40.0); + t1537 = RATIONAL(-9.0,80.0); + t1534 = RATIONAL(-3.0,20.0); + t1530 = RATIONAL(-11.0,80.0); + t1523 = RATIONAL(-3.0,16.0)*x; + t1522 = RATIONAL(-1.0,16.0)*x; + coeffs_dxx->coeff_m1_m1_m1 = t1522+t1557; + coeffs_dxx->coeff_0_m1_m1 = RATIONAL(-13.0,80.0)+t1520+t1563; + coeffs_dxx->coeff_p1_m1_m1 = t1523+t1559; + coeffs_dxx->coeff_p2_m1_m1 = t1540+t1529+t1562; + coeffs_dxx->coeff_m1_0_m1 = t1522+t1549; + coeffs_dxx->coeff_0_0_m1 = t1534+t1518+t1563; + coeffs_dxx->coeff_p1_0_m1 = t1523+t1553; + coeffs_dxx->coeff_p2_0_m1 = t1525+t1541+t1569; + coeffs_dxx->coeff_m1_p1_m1 = t1522+t1554; + coeffs_dxx->coeff_0_p1_m1 = t1524+t1530+t1563; + coeffs_dxx->coeff_p1_p1_m1 = t1523+t1556; + coeffs_dxx->coeff_p2_p1_m1 = t1573+t1578; + coeffs_dxx->coeff_m1_p2_m1 = t1522+t1551; + coeffs_dxx->coeff_0_p2_m1 = t1545+t1519+t1563; + coeffs_dxx->coeff_p1_p2_m1 = t1523+t1555; + coeffs_dxx->coeff_p2_p2_m1 = t1529+t1570; + coeffs_dxx->coeff_m1_m1_0 = t1522+t1546; + coeffs_dxx->coeff_0_m1_0 = t1528+t1534+t1565; + coeffs_dxx->coeff_p1_m1_0 = t1523+t1552; + coeffs_dxx->coeff_p2_m1_0 = t1541+t1527+t1562; + coeffs_dxx->coeff_m1_0_0 = t1522+t1548; + coeffs_dxx->coeff_0_0_0 = t1530+t1518+t1571; + coeffs_dxx->coeff_p1_0_0 = t1523+t1560; + coeffs_dxx->coeff_p2_0_0 = t1564+t1578; + coeffs_dxx->coeff_m1_p1_0 = t1522+t1547; + coeffs_dxx->coeff_0_p1_0 = t1524+t1545+t1571; + coeffs_dxx->coeff_p1_p1_0 = t1523+t1550; + coeffs_dxx->coeff_p2_p1_0 = t1525+t1576; + coeffs_dxx->coeff_m1_p2_0 = t1522+t1558; + coeffs_dxx->coeff_0_p2_0 = t1537+t1519+t1571; + coeffs_dxx->coeff_p1_p2_0 = t1523+t1561; + coeffs_dxx->coeff_p2_p2_0 = t1527+t1538+t1570; + coeffs_dxx->coeff_m1_m1_p1 = t1522+t1561; + coeffs_dxx->coeff_0_m1_p1 = t1527+t1530+t1565; + coeffs_dxx->coeff_p1_m1_p1 = t1523+t1558; + coeffs_dxx->coeff_p2_m1_p1 = t1528+t1539+t1562; + coeffs_dxx->coeff_m1_0_p1 = t1522+t1550; + coeffs_dxx->coeff_0_0_p1 = t1545+t1527+t1572; + coeffs_dxx->coeff_p1_0_p1 = t1523+t1547; + coeffs_dxx->coeff_p2_0_p1 = t1525+t1575; + coeffs_dxx->coeff_m1_p1_p1 = t1522+t1560; + coeffs_dxx->coeff_0_p1_p1 = t1537+t1521+t1564; + coeffs_dxx->coeff_p1_p1_p1 = t1523+t1548; + coeffs_dxx->coeff_p2_p1_p1 = t1538+t1525+t1574; + coeffs_dxx->coeff_m1_p2_p1 = t1522+t1552; + coeffs_dxx->coeff_0_p2_p1 = t1567+t1568; + coeffs_dxx->coeff_p1_p2_p1 = t1523+t1546; + coeffs_dxx->coeff_p2_p2_p1 = t1528+t1542+t1570; + coeffs_dxx->coeff_m1_m1_p2 = t1522+t1555; + coeffs_dxx->coeff_0_m1_p2 = t1545+t1529+t1565; + coeffs_dxx->coeff_p1_m1_p2 = t1523+t1551; + coeffs_dxx->coeff_p2_m1_p2 = t1526+t1562; + coeffs_dxx->coeff_m1_0_p2 = t1522+t1556; + coeffs_dxx->coeff_0_0_p2 = t1529+t1537+t1572; + coeffs_dxx->coeff_p1_0_p2 = t1523+t1554; + coeffs_dxx->coeff_p2_0_p2 = t1538+t1524+t1566; + coeffs_dxx->coeff_m1_p1_p2 = t1522+t1553; + coeffs_dxx->coeff_0_p1_p2 = t1568+t1569; + coeffs_dxx->coeff_p1_p1_p2 = t1523+t1549; + coeffs_dxx->coeff_p2_p1_p2 = t1542+t1518+t1566; + coeffs_dxx->coeff_m1_p2_p2 = t1522+t1559; + coeffs_dxx->coeff_0_p2_p2 = RATIONAL(-7.0,80.0)+t1521+t1577; + coeffs_dxx->coeff_p1_p2_p2 = t1523+t1557; + coeffs_dxx->coeff_p2_p2_p2 = t1531+t1520+t1566; diff --git a/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxy.compute.c b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxy.compute.c new file mode 100644 index 0000000..e7d150e --- /dev/null +++ b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxy.compute.c @@ -0,0 +1,208 @@ +fp t1622; +fp t1588; +fp t1590; +fp t1660; +fp t1611; +fp t1584; +fp t1659; +fp t1621; +fp t1587; +fp t1658; +fp t1583; +fp t1657; +fp t1613; +fp t1585; +fp t1656; +fp t1586; +fp t1599; +fp t1655; +fp t1581; +fp t1593; +fp t1654; +fp t1653; +fp t1623; +fp t1594; +fp t1652; +fp t1651; +fp t1650; +fp t1649; +fp t1648; +fp t1647; +fp t1603; +fp t1580; +fp t1646; +fp t1645; +fp t1589; +fp t1644; +fp t1643; +fp t1642; +fp t1620; +fp t1591; +fp t1592; +fp t1641; +fp t1640; +fp t1639; +fp t1638; +fp t1637; +fp t1636; +fp t1635; +fp t1634; +fp t1633; +fp t1632; +fp t1579; +fp t1631; +fp t1630; +fp t1629; +fp t1628; +fp t1627; +fp t1607; +fp t1582; +fp t1626; +fp t1625; +fp t1624; +fp t1618; +fp t1616; +fp t1612; +fp t1609; +fp t1604; +fp t1602; +fp t1601; +fp t1598; +fp t1597; +fp t1595; + t1622 = RATIONAL(1.0,80.0); + t1588 = t1622*y; + t1590 = t1622*x; + t1660 = t1588+t1590; + t1611 = RATIONAL(-3.0,1000.0); + t1584 = t1611*z; + t1659 = t1584+RATIONAL(-31.0,1000.0); + t1621 = RATIONAL(-1.0,80.0); + t1587 = t1621*x; + t1658 = t1587+t1588; + t1583 = RATIONAL(3.0,1000.0)*z; + t1657 = t1583+RATIONAL(37.0,2000.0); + t1613 = RATIONAL(-9.0,1000.0); + t1585 = t1613*z; + t1656 = t1585+RATIONAL(49.0,2000.0); + t1586 = RATIONAL(9.0,1000.0)*z; + t1599 = RATIONAL(-19.0,2000.0); + t1655 = t1586+t1599; + t1581 = RATIONAL(1.0,1000.0)*z; + t1593 = t1621*y; + t1654 = t1581+t1593; + t1653 = t1583+RATIONAL(-17.0,500.0); + t1623 = RATIONAL(3.0,80.0); + t1594 = t1623*x; + t1652 = t1594+t1593; + t1651 = t1583+RATIONAL(-1.0,250.0); + t1650 = t1584+t1593; + t1649 = t1584+t1587; + t1648 = t1583+RATIONAL(2.0,125.0); + t1647 = t1581+t1588; + t1603 = RATIONAL(-27.0,1000.0); + t1580 = t1603*z; + t1646 = t1580+t1613; + t1645 = t1585+RATIONAL(-9.0,500.0); + t1589 = t1623*y; + t1644 = t1589+t1594; + t1643 = t1585+RATIONAL(-7.0,250.0); + t1642 = t1585+t1594; + t1620 = RATIONAL(-3.0,80.0); + t1591 = t1620*x; + t1592 = t1620*y; + t1641 = t1591+t1592; + t1640 = t1589+t1590; + t1639 = t1586+t1589; + t1638 = t1588+t1594; + t1637 = t1592+t1594; + t1636 = t1590+t1592; + t1635 = t1588+t1591; + t1634 = t1591+t1593; + t1633 = t1584+t1591; + t1632 = t1589+t1591; + t1579 = RATIONAL(27.0,1000.0)*z; + t1631 = t1592+t1579; + t1630 = t1587+t1592; + t1629 = t1586+t1592; + t1628 = t1585+RATIONAL(11.0,500.0); + t1627 = t1586+t1591; + t1607 = RATIONAL(-1.0,1000.0); + t1582 = t1607*z; + t1626 = t1582+t1593; + t1625 = t1583+RATIONAL(-13.0,2000.0); + t1624 = t1587+t1589; + t1618 = RATIONAL(-1.0,500.0); + t1616 = RATIONAL(-9.0,250.0); + t1612 = RATIONAL(-1.0,2000.0); + t1609 = RATIONAL(-7.0,2000.0); + t1604 = RATIONAL(19.0,1000.0); + t1602 = RATIONAL(-37.0,1000.0); + t1601 = RATIONAL(13.0,1000.0); + t1598 = RATIONAL(43.0,2000.0); + t1597 = RATIONAL(31.0,2000.0); + t1595 = RATIONAL(-21.0,2000.0); + coeffs_dxy->coeff_m1_m1_m1 = t1580+RATIONAL(147.0,2000.0)+t1641; + coeffs_dxy->coeff_0_m1_m1 = t1593+t1612+t1642; + coeffs_dxy->coeff_p1_m1_m1 = t1586+t1602+t1638; + coeffs_dxy->coeff_p2_m1_m1 = t1579+t1616+t1632; + coeffs_dxy->coeff_m1_0_m1 = t1585+t1612+t1624; + coeffs_dxy->coeff_0_0_m1 = t1584+RATIONAL(-17.0,2000.0)+t1660; + coeffs_dxy->coeff_p1_0_m1 = t1593+t1590+t1651; + coeffs_dxy->coeff_p2_0_m1 = t1587+t1601+t1629; + coeffs_dxy->coeff_m1_p1_m1 = t1602+t1590+t1639; + coeffs_dxy->coeff_0_p1_m1 = t1651+t1658; + coeffs_dxy->coeff_p1_p1_m1 = RATIONAL(33.0,2000.0)+t1593+t1649; + coeffs_dxy->coeff_p2_p1_m1 = t1636+t1656; + coeffs_dxy->coeff_m1_p2_m1 = t1594+t1616+t1631; + coeffs_dxy->coeff_0_p2_m1 = t1601+t1593+t1627; + coeffs_dxy->coeff_p1_p2_m1 = t1635+t1656; + coeffs_dxy->coeff_p2_p2_m1 = RATIONAL(-3.0,2000.0)+t1580+t1644; + coeffs_dxy->coeff_m1_m1_0 = RATIONAL(129.0,2000.0)+t1585+t1641; + coeffs_dxy->coeff_0_m1_0 = t1609+t1594+t1650; + coeffs_dxy->coeff_p1_m1_0 = t1638+t1653; + coeffs_dxy->coeff_p2_m1_0 = t1589+t1603+t1627; + coeffs_dxy->coeff_m1_0_0 = t1609+t1584+t1624; + coeffs_dxy->coeff_0_0_0 = t1599+t1582+t1660; + coeffs_dxy->coeff_p1_0_0 = t1611+t1590+t1654; + coeffs_dxy->coeff_p2_0_0 = t1630+t1648; + coeffs_dxy->coeff_m1_p1_0 = t1640+t1653; + coeffs_dxy->coeff_0_p1_0 = t1587+t1611+t1647; + coeffs_dxy->coeff_p1_p1_0 = t1587+t1597+t1626; + coeffs_dxy->coeff_p2_p1_0 = t1598+t1584+t1636; + coeffs_dxy->coeff_m1_p2_0 = t1603+t1594+t1629; + coeffs_dxy->coeff_0_p2_0 = t1634+t1648; + coeffs_dxy->coeff_p1_p2_0 = t1588+t1598+t1633; + coeffs_dxy->coeff_p2_p2_0 = t1595+t1589+t1642; + coeffs_dxy->coeff_m1_m1_p1 = t1592+RATIONAL(111.0,2000.0)+t1627; + coeffs_dxy->coeff_0_m1_p1 = t1625+t1652; + coeffs_dxy->coeff_p1_m1_p1 = t1638+t1659; + coeffs_dxy->coeff_p2_m1_p1 = t1632+t1645; + coeffs_dxy->coeff_m1_0_p1 = t1624+t1625; + coeffs_dxy->coeff_0_0_p1 = t1590+t1595+t1647; + coeffs_dxy->coeff_p1_0_p1 = t1618+t1590+t1626; + coeffs_dxy->coeff_p2_0_p1 = t1584+t1604+t1630; + coeffs_dxy->coeff_m1_p1_p1 = t1640+t1659; + coeffs_dxy->coeff_0_p1_p1 = t1582+t1618+t1658; + coeffs_dxy->coeff_p1_p1_p1 = t1587+RATIONAL(29.0,2000.0)+t1654; + coeffs_dxy->coeff_p2_p1_p1 = t1636+t1657; + coeffs_dxy->coeff_m1_p2_p1 = t1637+t1645; + coeffs_dxy->coeff_0_p2_p1 = t1604+t1593+t1633; + coeffs_dxy->coeff_p1_p2_p1 = t1635+t1657; + coeffs_dxy->coeff_p2_p2_p1 = t1594+RATIONAL(-39.0,2000.0)+t1639; + coeffs_dxy->coeff_m1_m1_p2 = t1591+RATIONAL(93.0,2000.0)+t1631; + coeffs_dxy->coeff_0_m1_p2 = t1652+t1655; + coeffs_dxy->coeff_p1_m1_p2 = t1638+t1643; + coeffs_dxy->coeff_p2_m1_p2 = t1632+t1646; + coeffs_dxy->coeff_m1_0_p2 = t1624+t1655; + coeffs_dxy->coeff_0_0_p2 = RATIONAL(-23.0,2000.0)+t1583+t1660; + coeffs_dxy->coeff_p1_0_p2 = t1590+t1607+t1650; + coeffs_dxy->coeff_p2_0_p2 = t1628+t1630; + coeffs_dxy->coeff_m1_p1_p2 = t1640+t1643; + coeffs_dxy->coeff_0_p1_p2 = t1607+t1588+t1649; + coeffs_dxy->coeff_p1_p1_p2 = t1593+RATIONAL(27.0,2000.0)+t1583+t1587; + coeffs_dxy->coeff_p2_p1_p2 = t1590+t1597+t1629; + coeffs_dxy->coeff_m1_p2_p2 = t1637+t1646; + coeffs_dxy->coeff_0_p2_p2 = t1628+t1634; + coeffs_dxy->coeff_p1_p2_p2 = t1597+t1588+t1627; + coeffs_dxy->coeff_p2_p2_p2 = t1579+RATIONAL(-57.0,2000.0)+t1644; diff --git a/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxz.compute.c b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxz.compute.c new file mode 100644 index 0000000..1f44ca7 --- /dev/null +++ b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxz.compute.c @@ -0,0 +1,202 @@ +fp t1695; +fp t1666; +fp t1741; +fp t1685; +fp t1662; +fp t1740; +fp t1739; +fp t1661; +fp t1738; +fp t1737; +fp t1689; +fp t1664; +fp t1703; +fp t1674; +fp t1736; +fp t1705; +fp t1670; +fp t1671; +fp t1735; +fp t1734; +fp t1702; +fp t1676; +fp t1673; +fp t1733; +fp t1667; +fp t1732; +fp t1731; +fp t1730; +fp t1665; +fp t1729; +fp t1675; +fp t1728; +fp t1663; +fp t1704; +fp t1672; +fp t1727; +fp t1693; +fp t1668; +fp t1726; +fp t1725; +fp t1724; +fp t1723; +fp t1722; +fp t1679; +fp t1721; +fp t1669; +fp t1720; +fp t1719; +fp t1718; +fp t1717; +fp t1716; +fp t1715; +fp t1714; +fp t1713; +fp t1712; +fp t1711; +fp t1710; +fp t1709; +fp t1708; +fp t1707; +fp t1706; +fp t1700; +fp t1697; +fp t1691; +fp t1688; +fp t1686; +fp t1681; +fp t1678; +fp t1677; + t1695 = RATIONAL(-9.0,1000.0); + t1666 = t1695*y; + t1741 = t1666+RATIONAL(11.0,500.0); + t1685 = RATIONAL(-27.0,1000.0); + t1662 = t1685*y; + t1740 = t1695+t1662; + t1739 = t1666+RATIONAL(-7.0,250.0); + t1661 = RATIONAL(27.0,1000.0)*y; + t1738 = t1661+RATIONAL(-9.0,250.0); + t1737 = t1666+RATIONAL(-9.0,500.0); + t1689 = RATIONAL(-1.0,1000.0); + t1664 = t1689*y; + t1703 = RATIONAL(-1.0,80.0); + t1674 = t1703*z; + t1736 = t1664+t1674; + t1705 = RATIONAL(3.0,80.0); + t1670 = t1705*z; + t1671 = t1703*x; + t1735 = t1670+t1671; + t1734 = t1666+RATIONAL(-1.0,2000.0); + t1702 = RATIONAL(-3.0,80.0); + t1676 = t1702*z; + t1673 = t1705*x; + t1733 = t1676+t1673; + t1667 = RATIONAL(3.0,1000.0)*y; + t1732 = t1667+t1674; + t1731 = t1666+RATIONAL(49.0,2000.0); + t1730 = t1670+t1673; + t1665 = RATIONAL(9.0,1000.0)*y; + t1729 = t1665+t1673; + t1675 = t1702*x; + t1728 = t1670+t1675; + t1663 = RATIONAL(1.0,1000.0)*y; + t1704 = RATIONAL(1.0,80.0); + t1672 = t1704*z; + t1727 = t1663+t1672; + t1693 = RATIONAL(-3.0,1000.0); + t1668 = t1693*y; + t1726 = RATIONAL(43.0,2000.0)+t1668; + t1725 = t1664+t1672; + t1724 = RATIONAL(-37.0,1000.0)+t1665; + t1723 = t1672+t1675; + t1722 = t1663+t1674; + t1679 = RATIONAL(31.0,2000.0); + t1721 = t1679+t1665; + t1669 = t1704*x; + t1720 = t1669+t1670; + t1719 = t1667+t1672; + t1718 = t1672+t1673; + t1717 = t1667+t1669; + t1716 = RATIONAL(-31.0,1000.0)+t1668; + t1715 = t1676+t1675; + t1714 = t1667+RATIONAL(2.0,125.0); + t1713 = t1674+t1675; + t1712 = RATIONAL(13.0,1000.0)+t1665; + t1711 = t1676+t1671; + t1710 = t1667+RATIONAL(-17.0,500.0); + t1709 = t1668+t1674; + t1708 = t1668+t1671; + t1707 = t1676+t1669; + t1706 = t1685+t1665; + t1700 = RATIONAL(-1.0,500.0); + t1697 = RATIONAL(-1.0,250.0); + t1691 = RATIONAL(-7.0,2000.0); + t1688 = RATIONAL(37.0,2000.0); + t1686 = RATIONAL(19.0,1000.0); + t1681 = RATIONAL(-19.0,2000.0); + t1678 = RATIONAL(-13.0,2000.0); + t1677 = RATIONAL(-21.0,2000.0); + coeffs_dxz->coeff_m1_m1_m1 = RATIONAL(147.0,2000.0)+t1662+t1715; + coeffs_dxz->coeff_0_m1_m1 = t1673+t1674+t1734; + coeffs_dxz->coeff_p1_m1_m1 = t1718+t1724; + coeffs_dxz->coeff_p2_m1_m1 = t1728+t1738; + coeffs_dxz->coeff_m1_0_m1 = RATIONAL(129.0,2000.0)+t1666+t1715; + coeffs_dxz->coeff_0_0_m1 = t1691+t1673+t1709; + coeffs_dxz->coeff_p1_0_m1 = t1710+t1718; + coeffs_dxz->coeff_p2_0_m1 = t1706+t1728; + coeffs_dxz->coeff_m1_p1_m1 = RATIONAL(111.0,2000.0)+t1665+t1715; + coeffs_dxz->coeff_0_p1_m1 = t1673+t1678+t1732; + coeffs_dxz->coeff_p1_p1_m1 = t1716+t1718; + coeffs_dxz->coeff_p2_p1_m1 = t1728+t1737; + coeffs_dxz->coeff_m1_p2_m1 = RATIONAL(93.0,2000.0)+t1661+t1715; + coeffs_dxz->coeff_0_p2_m1 = t1674+t1681+t1729; + coeffs_dxz->coeff_p1_p2_m1 = t1718+t1739; + coeffs_dxz->coeff_p2_p2_m1 = t1728+t1740; + coeffs_dxz->coeff_m1_m1_0 = t1734+t1735; + coeffs_dxz->coeff_0_m1_0 = t1669+t1668+t1672+RATIONAL(-17.0,2000.0); + coeffs_dxz->coeff_p1_m1_0 = t1674+t1697+t1717; + coeffs_dxz->coeff_p2_m1_0 = t1711+t1712; + coeffs_dxz->coeff_m1_0_0 = t1691+t1670+t1708; + coeffs_dxz->coeff_0_0_0 = t1669+t1681+t1725; + coeffs_dxz->coeff_p1_0_0 = t1693+t1669+t1722; + coeffs_dxz->coeff_p2_0_0 = t1711+t1714; + coeffs_dxz->coeff_m1_p1_0 = t1678+t1667+t1735; + coeffs_dxz->coeff_0_p1_0 = t1669+t1677+t1727; + coeffs_dxz->coeff_p1_p1_0 = t1700+t1669+t1736; + coeffs_dxz->coeff_p2_p1_0 = t1676+t1686+t1708; + coeffs_dxz->coeff_m1_p2_0 = t1665+t1681+t1735; + coeffs_dxz->coeff_0_p2_0 = t1672+RATIONAL(-23.0,2000.0)+t1717; + coeffs_dxz->coeff_p1_p2_0 = t1689+t1669+t1709; + coeffs_dxz->coeff_p2_p2_0 = t1711+t1741; + coeffs_dxz->coeff_m1_m1_p1 = t1720+t1724; + coeffs_dxz->coeff_0_m1_p1 = t1671+t1697+t1719; + coeffs_dxz->coeff_p1_m1_p1 = RATIONAL(33.0,2000.0)+t1674+t1708; + coeffs_dxz->coeff_p2_m1_p1 = t1707+t1731; + coeffs_dxz->coeff_m1_0_p1 = t1710+t1720; + coeffs_dxz->coeff_0_0_p1 = t1671+t1693+t1727; + coeffs_dxz->coeff_p1_0_p1 = t1671+t1679+t1736; + coeffs_dxz->coeff_p2_0_p1 = t1707+t1726; + coeffs_dxz->coeff_m1_p1_p1 = t1716+t1720; + coeffs_dxz->coeff_0_p1_p1 = t1700+t1671+t1725; + coeffs_dxz->coeff_p1_p1_p1 = RATIONAL(29.0,2000.0)+t1671+t1722; + coeffs_dxz->coeff_p2_p1_p1 = t1667+t1688+t1707; + coeffs_dxz->coeff_m1_p2_p1 = t1720+t1739; + coeffs_dxz->coeff_0_p2_p1 = t1689+t1672+t1708; + coeffs_dxz->coeff_p1_p2_p1 = RATIONAL(27.0,2000.0)+t1671+t1732; + coeffs_dxz->coeff_p2_p2_p1 = t1707+t1721; + coeffs_dxz->coeff_m1_m1_p2 = t1733+t1738; + coeffs_dxz->coeff_0_m1_p2 = t1712+t1713; + coeffs_dxz->coeff_p1_m1_p2 = t1723+t1731; + coeffs_dxz->coeff_p2_m1_p2 = t1662+RATIONAL(-3.0,2000.0)+t1730; + coeffs_dxz->coeff_m1_0_p2 = t1706+t1733; + coeffs_dxz->coeff_0_0_p2 = t1713+t1714; + coeffs_dxz->coeff_p1_0_p2 = t1723+t1726; + coeffs_dxz->coeff_p2_0_p2 = t1666+t1677+t1730; + coeffs_dxz->coeff_m1_p1_p2 = t1733+t1737; + coeffs_dxz->coeff_0_p1_p2 = t1686+t1675+t1709; + coeffs_dxz->coeff_p1_p1_p2 = t1688+t1675+t1719; + coeffs_dxz->coeff_p2_p1_p2 = RATIONAL(-39.0,2000.0)+t1670+t1729; + coeffs_dxz->coeff_m1_p2_p2 = t1733+t1740; + coeffs_dxz->coeff_0_p2_p2 = t1713+t1741; + coeffs_dxz->coeff_p1_p2_p2 = t1721+t1723; + coeffs_dxz->coeff_p2_p2_p2 = RATIONAL(-57.0,2000.0)+t1661+t1730; diff --git a/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dy.compute.c b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dy.compute.c new file mode 100644 index 0000000..c483dac --- /dev/null +++ b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dy.compute.c @@ -0,0 +1,510 @@ +fp t1156; +fp t1170; +fp t1069; +fp t1164; +fp t1169; +fp t1076; +fp t1237; +fp t1078; +fp t1073; +fp t1231; +fp t1171; +fp t1086; +fp t1160; +fp t1065; +fp t1265; +fp t1063; +fp t1264; +fp t1085; +fp t1263; +fp t1094; +fp t1262; +fp t1091; +fp t1261; +fp t1216; +fp t1260; +fp t1259; +fp t1080; +fp t1258; +fp t1161; +fp t1068; +fp t1070; +fp t1230; +fp t1257; +fp t1154; +fp t1017; +fp t1256; +fp t1255; +fp t1254; +fp t1082; +fp t1253; +fp t1162; +fp t1147; +fp t1051; +fp t1014; +fp t1252; +fp t1098; +fp t1251; +fp t1226; +fp t1250; +fp t1093; +fp t1249; +fp t1248; +fp t1096; +fp t1247; +fp t1246; +fp t1088; +fp t1245; +fp t1083; +fp t1244; +fp t1240; +fp t1013; +fp t1243; +fp t1242; +fp t1241; +fp t1144; +fp t1027; +fp t1146; +fp t1050; +fp t1239; +fp t1152; +fp t1057; +fp t1021; +fp t1016; +fp t1238; +fp t1022; +fp t1019; +fp t1236; +fp t1020; +fp t1235; +fp t1015; +fp t1234; +fp t1018; +fp t1233; +fp t1232; +fp t1229; +fp t1228; +fp t1227; +fp t1225; +fp t1224; +fp t1223; +fp t1217; +fp t1090; +fp t1215; +fp t1214; +fp t1089; +fp t1213; +fp t1212; +fp t1211; +fp t1210; +fp t1209; +fp t1208; +fp t1012; +fp t1207; +fp t1097; +fp t1206; +fp t1205; +fp t1011; +fp t1204; +fp t1203; +fp t1202; +fp t1201; +fp t1008; +fp t1200; +fp t1199; +fp t1009; +fp t1198; +fp t1197; +fp t1196; +fp t1007; +fp t1195; +fp t1194; +fp t1010; +fp t1193; +fp t1192; +fp t1191; +fp t1190; +fp t1189; +fp t1188; +fp t1187; +fp t1186; +fp t1185; +fp t1184; +fp t1183; +fp t1182; +fp t1181; +fp t1180; +fp t1179; +fp t1178; +fp t1177; +fp t1099; +fp t1176; +fp t1175; +fp t1174; +fp t1173; +fp t1172; +fp t1163; +fp t1159; +fp t1158; +fp t1157; +fp t1153; +fp t1151; +fp t1149; +fp t1148; +fp t1145; +fp t1143; +fp t1142; +fp t1138; +fp t1137; +fp t1136; +fp t1135; +fp t1133; +fp t1132; +fp t1129; +fp t1128; +fp t1125; +fp t1124; +fp t1123; +fp t1121; +fp t1119; +fp t1118; +fp t1117; +fp t1112; +fp t1110; +fp t1081; +fp t1079; +fp t1077; +fp t1075; +fp t1074; +fp t1072; +fp t1067; +fp t1064; +fp t1062; +fp t1060; +fp t1059; +fp t1058; +fp t1056; +fp t1055; +fp t1054; +fp t1053; +fp t1052; +fp t1049; +fp t1048; +fp t1047; +fp t1046; +fp t1045; +fp t1044; +fp t1043; +fp t1042; +fp t1041; +fp t1040; +fp t1039; +fp t1038; +fp t1037; +fp t1036; +fp t1035; +fp t1034; +fp t1033; +fp t1032; +fp t1031; +fp t1030; +fp t1029; +fp t1028; +fp t1026; +fp t1025; +fp t1024; +fp t1023; + t1156 = RATIONAL(3.0,160.0); + t1170 = x*x; + t1069 = t1156*t1170; + t1164 = RATIONAL(-3.0,160.0); + t1169 = z*z; + t1076 = t1164*t1169; + t1237 = t1069+t1076; + t1078 = t1164*t1170; + t1073 = t1156*t1169; + t1231 = t1078+t1073; + t1171 = y*y; + t1086 = RATIONAL(3.0,32.0)*t1171; + t1160 = RATIONAL(1.0,160.0); + t1065 = t1160*t1170; + t1265 = t1086+t1065; + t1063 = RATIONAL(-11.0,80.0)*y; + t1264 = t1063+RATIONAL(-107.0,4000.0); + t1085 = RATIONAL(7.0,80.0)*y; + t1263 = RATIONAL(183.0,4000.0)+t1085; + t1094 = RATIONAL(3.0,80.0)*y; + t1262 = RATIONAL(53.0,12000.0)+t1094; + t1091 = RATIONAL(1.0,32.0)*t1171; + t1261 = t1091+t1237; + t1216 = t1078+t1076; + t1260 = t1216+t1091; + t1259 = RATIONAL(-79.0,2000.0)+RATIONAL(-1.0,10.0)*y; + t1080 = t1160*t1169; + t1258 = t1080+t1265; + t1161 = RATIONAL(-1.0,160.0); + t1068 = t1161*t1170; + t1070 = t1161*t1169; + t1230 = t1068+t1070; + t1257 = t1086+t1230; + t1154 = RATIONAL(-1.0,250.0); + t1017 = z*t1094; + t1256 = t1154*x+t1017; + t1255 = t1070+t1265; + t1254 = RATIONAL(323.0,4000.0)+t1094; + t1082 = RATIONAL(-9.0,80.0)*y; + t1253 = RATIONAL(-137.0,4000.0)+t1082; + t1162 = RATIONAL(-9.0,500.0); + t1147 = RATIONAL(-9.0,1000.0); + t1051 = t1147*x; + t1014 = z*t1051; + t1252 = t1162*x+t1014; + t1098 = RATIONAL(3.0,40.0)*y; + t1251 = t1098+RATIONAL(53.0,1000.0); + t1226 = t1069+t1073; + t1250 = t1226+t1091; + t1093 = RATIONAL(1.0,20.0)*y; + t1249 = t1093+RATIONAL(141.0,2000.0); + t1248 = RATIONAL(-3.0,20.0)*y+RATIONAL(-49.0,2000.0); + t1096 = RATIONAL(-3.0,32.0)*t1171; + t1247 = t1096+t1065+t1070; + t1246 = t1091+t1231; + t1088 = RATIONAL(1.0,80.0)*y; + t1245 = t1088+RATIONAL(-137.0,12000.0); + t1083 = RATIONAL(-1.0,80.0)*y; + t1244 = RATIONAL(-167.0,12000.0)+t1014+t1083; + t1240 = x*z; + t1013 = RATIONAL(9.0,1000.0)*t1240; + t1243 = t1093+t1013+RATIONAL(-29.0,6000.0); + t1242 = RATIONAL(-1.0,40.0)*y+t1013+RATIONAL(-4.0,375.0); + t1241 = z+x; + t1144 = RATIONAL(-27.0,1000.0); + t1027 = t1144*x; + t1146 = RATIONAL(-3.0,1000.0); + t1050 = t1146*x; + t1239 = RATIONAL(-3.0,80.0)*y; + t1152 = RATIONAL(-1.0,1000.0); + t1057 = t1152*z; + t1021 = x*t1094; + t1016 = z*t1088; + t1238 = t1021+t1016; + t1022 = x*t1083; + t1019 = z*t1239; + t1236 = t1022+t1019; + t1020 = x*t1239; + t1235 = t1019+t1020; + t1015 = x*t1088; + t1234 = t1015+t1016; + t1018 = z*t1083; + t1233 = t1018+t1020; + t1232 = t1015+t1019; + t1229 = t1017+t1020; + t1228 = t1021+t1018; + t1227 = t1065+t1080; + t1225 = t1068+t1080; + t1224 = t1021+t1019; + t1223 = t1021+t1017; + t1217 = t1016+t1022; + t1090 = RATIONAL(-1.0,32.0)*t1171; + t1215 = t1069+t1090; + t1214 = t1016+t1020; + t1089 = RATIONAL(1.0,16.0)*y; + t1213 = t1096+t1089; + t1212 = t1015+t1018; + t1211 = t1022+t1018; + t1210 = t1089+t1090; + t1209 = t1014+t1015+t1162*z; + t1208 = RATIONAL(1.0,10.0)*y+t1235; + t1012 = RATIONAL(3.0,1000.0)*t1240; + t1207 = t1012+t1225; + t1097 = RATIONAL(1.0,40.0)*y; + t1206 = t1097+t1223; + t1205 = t1093+t1234; + t1011 = z*t1050; + t1204 = t1011+t1225; + t1203 = t1098+t1211; + t1202 = t1015+t1255; + t1201 = t1098+RATIONAL(-97.0,3000.0)+t1013+t1090; + t1008 = RATIONAL(27.0,1000.0)*t1240; + t1200 = t1008+t1250; + t1199 = t1017+t1246; + t1009 = RATIONAL(1.0,1000.0)*t1240; + t1198 = t1009+t1096+t1230; + t1197 = t1096+t1013+t1227; + t1196 = t1096+t1207; + t1007 = z*t1027; + t1195 = t1090+t1007+t1216; + t1194 = t1086+t1207; + t1010 = x*t1057; + t1193 = t1010+t1258; + t1192 = t1011+t1073+t1215; + t1191 = t1022+t1255; + t1190 = t1016+t1261; + t1189 = t1096+t1204; + t1188 = t1017+t1076+t1215; + t1187 = t1021+t1090+t1231; + t1186 = t1011+t1247; + t1185 = t1012+t1247; + t1184 = t1012+t1260; + t1183 = RATIONAL(-17.0,3000.0)+t1097+t1013+t1091; + t1182 = t1086+t1204; + t1181 = t1014+t1257; + t1180 = RATIONAL(-607.0,12000.0)+t1014+t1090+t1085; + t1179 = RATIONAL(-41.0,750.0)+t1007+t1250; + t1178 = RATIONAL(79.0,3000.0)+t1011+t1260; + t1177 = t1014+RATIONAL(191.0,4000.0)+t1213+t1227; + t1099 = RATIONAL(-1.0,8.0)*y; + t1176 = t1009+t1099+RATIONAL(-87.0,2000.0)+t1258; + t1175 = t1008+RATIONAL(281.0,12000.0)+t1210+t1216; + t1174 = t1099+t1013+RATIONAL(-33.0,2000.0)+t1257; + t1173 = RATIONAL(-691.0,12000.0)+t1012+t1210+t1226; + t1172 = t1010+RATIONAL(299.0,4000.0)+t1213+t1230; + t1163 = RATIONAL(-7.0,250.0); + t1159 = RATIONAL(-9.0,250.0); + t1158 = RATIONAL(-1.0,500.0); + t1157 = RATIONAL(2.0,125.0); + t1153 = RATIONAL(-17.0,500.0); + t1151 = RATIONAL(-3.0,2000.0); + t1149 = RATIONAL(11.0,500.0); + t1148 = RATIONAL(-1.0,2000.0); + t1145 = RATIONAL(-7.0,2000.0); + t1143 = RATIONAL(37.0,2000.0); + t1142 = RATIONAL(13.0,1000.0); + t1138 = RATIONAL(49.0,2000.0); + t1137 = RATIONAL(27.0,2000.0); + t1136 = RATIONAL(-37.0,1000.0); + t1135 = RATIONAL(43.0,2000.0); + t1133 = RATIONAL(31.0,2000.0); + t1132 = RATIONAL(-19.0,2000.0); + t1129 = RATIONAL(93.0,2000.0); + t1128 = RATIONAL(-21.0,2000.0); + t1125 = RATIONAL(19.0,1000.0); + t1124 = RATIONAL(-23.0,2000.0); + t1123 = RATIONAL(-39.0,2000.0); + t1121 = RATIONAL(-17.0,2000.0); + t1119 = RATIONAL(-31.0,1000.0); + t1118 = RATIONAL(33.0,2000.0); + t1117 = RATIONAL(-13.0,2000.0); + t1112 = RATIONAL(111.0,2000.0); + t1110 = RATIONAL(129.0,2000.0); + t1081 = t1158*x; + t1079 = t1157*z; + t1077 = t1158*z; + t1075 = t1163*z; + t1074 = t1159*z; + t1072 = t1159*x; + t1067 = t1163*x; + t1064 = t1157*x; + t1062 = t1154*z; + t1060 = t1153*z; + t1059 = t1147*z; + t1058 = t1148*x; + t1056 = t1145*x; + t1055 = t1153*x; + t1054 = t1148*z; + t1053 = t1149*z; + t1052 = t1145*z; + t1049 = t1152*x; + t1048 = t1149*x; + t1047 = t1146*z; + t1046 = t1117*z; + t1045 = t1132*z; + t1044 = t1142*x; + t1043 = t1136*z; + t1042 = t1125*x; + t1041 = t1119*z; + t1040 = t1133*z; + t1039 = t1138*x; + t1038 = t1117*x; + t1037 = t1133*x; + t1036 = t1125*z; + t1035 = t1143*x; + t1034 = t1144*z; + t1033 = t1138*z; + t1032 = t1135*z; + t1031 = t1143*z; + t1030 = t1119*x; + t1029 = t1135*x; + t1028 = t1142*z; + t1026 = t1128*x; + t1025 = t1132*x; + t1024 = t1136*x; + t1023 = t1128*z; + coeffs_dy->coeff_m1_m1_m1 = RATIONAL(-47.0,1500.0)+t1241*RATIONAL(147.0, +2000.0)+t1195+t1208; + coeffs_dy->coeff_0_m1_m1 = t1058+t1110*z+t1180+t1236+t1237; + coeffs_dy->coeff_p1_m1_m1 = t1024+t1112*z+t1201+t1232+t1237; + coeffs_dy->coeff_p2_m1_m1 = t1129*z+t1072+t1175+t1224; + coeffs_dy->coeff_m1_0_m1 = RATIONAL(-13.0,80.0)*y+t1058+t1054+RATIONAL( +-39.0,4000.0)+t1181+t1223; + coeffs_dy->coeff_0_0_m1 = t1017+t1121*x+t1052+t1011+t1202+t1248; + coeffs_dy->coeff_p1_0_m1 = t1046+t1012+t1191+t1256+t1264; + coeffs_dy->coeff_p2_0_m1 = t1045+t1044+t1174+t1229; + coeffs_dy->coeff_m1_p1_m1 = t1043+t1024+RATIONAL(157.0,2000.0)+t1197+ +t1206; + coeffs_dy->coeff_0_p1_m1 = t1060+t1015+t1196+t1254+t1256; + coeffs_dy->coeff_p1_p1_m1 = t1022+t1017+t1041+t1118*x+t1189+t1249; + coeffs_dy->coeff_p2_p1_m1 = t1039+t1075+t1177+t1229; + coeffs_dy->coeff_m1_p2_m1 = t1074+t1072+RATIONAL(-449.0,12000.0)+t1094+ +t1200+t1235; + coeffs_dy->coeff_0_p2_m1 = t1044+t1034+t1183+t1231+t1236; + coeffs_dy->coeff_p1_p2_m1 = t1019+t1039+t1209+t1245+t1246; + coeffs_dy->coeff_p2_p2_m1 = t1059+t1151*x+t1179+t1224; + coeffs_dy->coeff_m1_m1_0 = t1054+t1110*x+t1180+t1231+t1233; + coeffs_dy->coeff_0_m1_0 = t1052+RATIONAL(-437.0,6000.0)+t1056+t1192+t1203 +; + coeffs_dy->coeff_p1_m1_0 = t1055+t1046+t1173+t1212; + coeffs_dy->coeff_p2_m1_0 = t1027+t1045+t1018+t1187+t1243; + coeffs_dy->coeff_m1_0_0 = t1121*z+t1056+t1182+t1238+t1248; + coeffs_dy->coeff_0_0_0 = t1025+t1063+t1045+RATIONAL(-161.0,4000.0)+t1193+ +t1234; + coeffs_dy->coeff_p1_0_0 = t1050+t1023+t1176+t1217; + coeffs_dy->coeff_p2_0_0 = t1124*z+t1064+t1194+t1214+t1253; + coeffs_dy->coeff_m1_p1_0 = t1055+t1062+t1185+t1238+t1254; + coeffs_dy->coeff_0_p1_0 = t1050+t1047+RATIONAL(21.0,250.0)+t1198+t1205; + coeffs_dy->coeff_p1_p1_0 = t1077+t1037+t1172+t1217; + coeffs_dy->coeff_p2_p1_0 = t1057+t1029+t1186+t1214+t1251; + coeffs_dy->coeff_m1_p2_0 = t1028+t1027+t1183+t1233+t1237; + coeffs_dy->coeff_0_p2_0 = t1064+t1088+t1079+RATIONAL(349.0,12000.0)+t1184 ++t1211; + coeffs_dy->coeff_p1_p2_0 = t1036+t1029+t1178+t1212; + coeffs_dy->coeff_p2_p2_0 = t1026+t1053+t1228+t1244+t1261; + coeffs_dy->coeff_m1_m1_p1 = t1043+t1112*x+t1201+t1214+t1231; + coeffs_dy->coeff_0_m1_p1 = t1060+t1038+t1173+t1217; + coeffs_dy->coeff_p1_m1_p1 = RATIONAL(-17.0,375.0)+t1041+t1030+t1192+t1205 +; + coeffs_dy->coeff_p2_m1_p1 = t1075+t1016+t1187+t1252+t1262; + coeffs_dy->coeff_m1_0_p1 = t1062+t1038+t1194+t1228+t1264; + coeffs_dy->coeff_0_0_p1 = t1026+t1047+t1176+t1212; + coeffs_dy->coeff_p1_0_p1 = t1077+t1082+RATIONAL(-191.0,4000.0)+t1081+ +t1193+t1211; + coeffs_dy->coeff_p2_0_p1 = t1057+t1042+t1182+t1233+t1259; + coeffs_dy->coeff_m1_p1_p1 = t1030+t1118*z+t1186+t1228+t1249; + coeffs_dy->coeff_0_p1_p1 = t1081+t1040+t1172+t1212; + coeffs_dy->coeff_p1_p1_p1 = RATIONAL(133.0,2000.0)+t1241*RATIONAL(29.0, +2000.0)+t1198+t1203; + coeffs_dy->coeff_p2_p1_p1 = t1137*z+t1035+t1185+t1233+t1263; + coeffs_dy->coeff_m1_p2_p1 = t1020+t1033+t1190+t1245+t1252; + coeffs_dy->coeff_0_p2_p1 = t1042+t1032+t1178+t1217; + coeffs_dy->coeff_p1_p2_p1 = t1035+RATIONAL(319.0,12000.0)+t1031+t1083+ +t1184+t1234; + coeffs_dy->coeff_p2_p2_p1 = t1021+t1040+t1123*x+t1190+t1242; + coeffs_dy->coeff_m1_m1_p2 = t1074+t1129*x+t1175+t1229; + coeffs_dy->coeff_0_m1_p2 = t1025+t1034+t1022+t1188+t1243; + coeffs_dy->coeff_p1_m1_p2 = t1067+t1188+t1209+t1262; + coeffs_dy->coeff_p2_m1_p2 = t1059+t1051+RATIONAL(307.0,6000.0)+t1195+ +t1206; + coeffs_dy->coeff_m1_0_p2 = t1025+t1028+t1174+t1224; + coeffs_dy->coeff_0_0_p2 = t1079+t1124*x+t1019+t1012+t1202+t1253; + coeffs_dy->coeff_p1_0_p2 = t1011+t1036+t1049+t1019+t1191+t1259; + coeffs_dy->coeff_p2_0_p2 = t1048+t1053+RATIONAL(-129.0,4000.0)+RATIONAL( +-7.0,80.0)*y+t1181+t1235; + coeffs_dy->coeff_m1_p1_p2 = t1033+t1067+t1177+t1224; + coeffs_dy->coeff_0_p1_p2 = t1049+t1032+t1189+t1232+t1251; + coeffs_dy->coeff_p1_p1_p2 = t1031+t1137*x+t1196+t1236+t1263; + coeffs_dy->coeff_p2_p1_p2 = t1040+t1037+RATIONAL(13.0,500.0)+t1197+t1208; + coeffs_dy->coeff_m1_p2_p2 = t1151*z+t1051+t1179+t1229; + coeffs_dy->coeff_0_p2_p2 = t1023+t1022+t1048+t1199+t1244; + coeffs_dy->coeff_p1_p2_p2 = t1037+t1123*z+t1015+t1199+t1242; + coeffs_dy->coeff_p2_p2_p2 = t1239+RATIONAL(-539.0,12000.0)+t1241*RATIONAL +(-57.0,2000.0)+t1200+t1223; diff --git a/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dyy.compute.c b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dyy.compute.c new file mode 100644 index 0000000..52494e5 --- /dev/null +++ b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dyy.compute.c @@ -0,0 +1,186 @@ +fp t1761; +fp t1745; +fp t1755; +fp t1751; +fp t1804; +fp t1756; +fp t1744; +fp t1752; +fp t1803; +fp t1753; +fp t1802; +fp t1801; +fp t1760; +fp t1749; +fp t1800; +fp t1769; +fp t1746; +fp t1799; +fp t1798; +fp t1747; +fp t1797; +fp t1796; +fp t1795; +fp t1794; +fp t1743; +fp t1793; +fp t1792; +fp t1791; +fp t1762; +fp t1748; +fp t1790; +fp t1789; +fp t1788; +fp t1787; +fp t1786; +fp t1785; +fp t1784; +fp t1783; +fp t1768; +fp t1782; +fp t1763; +fp t1781; +fp t1780; +fp t1779; +fp t1778; +fp t1777; +fp t1776; +fp t1775; +fp t1765; +fp t1774; +fp t1773; +fp t1772; +fp t1766; +fp t1771; +fp t1770; +fp t1764; +fp t1758; +fp t1757; +fp t1754; +fp t1750; +fp t1742; + t1761 = RATIONAL(3.0,80.0); + t1745 = t1761*z; + t1755 = RATIONAL(-3.0,80.0); + t1751 = t1755*x; + t1804 = t1745+t1751; + t1756 = RATIONAL(-1.0,80.0); + t1744 = t1756*z; + t1752 = t1756*x; + t1803 = t1744+t1752; + t1753 = t1755*z; + t1802 = t1752+t1753; + t1801 = t1751+t1753; + t1760 = RATIONAL(1.0,80.0); + t1749 = t1760*x; + t1800 = t1745+t1749; + t1769 = RATIONAL(-1.0,8.0); + t1746 = RATIONAL(3.0,16.0)*y; + t1799 = t1769+t1746; + t1798 = t1744+t1751; + t1747 = t1761*x; + t1797 = t1746+t1747; + t1796 = t1744+t1749; + t1795 = t1745+t1752; + t1794 = t1746+t1752; + t1743 = t1760*z; + t1793 = t1743+t1746; + t1792 = t1743+t1752; + t1791 = t1749+t1753; + t1762 = RATIONAL(1.0,16.0); + t1748 = t1762*y; + t1790 = t1748+t1749; + t1789 = t1747+t1748; + t1788 = t1746+RATIONAL(-1.0,10.0); + t1787 = t1748+t1751; + t1786 = t1743+t1748; + t1785 = t1743+t1747+t1761; + t1784 = RATIONAL(1.0,10.0)+t1801; + t1783 = t1761+t1800; + t1768 = RATIONAL(1.0,40.0); + t1782 = t1768+t1745+t1747; + t1763 = RATIONAL(3.0,40.0); + t1781 = t1763+t1791; + t1780 = t1747+t1753+t1762; + t1779 = t1762+t1804; + t1778 = t1762+t1792; + t1777 = t1762+t1796; + t1776 = t1763+t1803; + t1775 = t1743+t1751+t1763; + t1765 = RATIONAL(1.0,20.0); + t1774 = t1765+t1795; + t1773 = t1743+t1749+t1765; + t1772 = t1744+t1747+t1765; + t1766 = RATIONAL(7.0,80.0); + t1771 = t1766+t1802; + t1770 = t1766+t1798; + t1764 = RATIONAL(-9.0,80.0); + t1758 = RATIONAL(-3.0,20.0); + t1757 = RATIONAL(-1.0,40.0); + t1754 = RATIONAL(-11.0,80.0); + t1750 = RATIONAL(-1.0,16.0)*y; + t1742 = RATIONAL(-3.0,16.0)*y; + coeffs_dyy->coeff_m1_m1_m1 = t1750+t1784; + coeffs_dyy->coeff_0_m1_m1 = t1750+t1771; + coeffs_dyy->coeff_p1_m1_m1 = t1750+t1781; + coeffs_dyy->coeff_p2_m1_m1 = t1750+t1780; + coeffs_dyy->coeff_m1_0_m1 = RATIONAL(-13.0,80.0)+t1745+t1797; + coeffs_dyy->coeff_0_0_m1 = t1746+t1758+t1800; + coeffs_dyy->coeff_p1_0_m1 = t1754+t1745+t1794; + coeffs_dyy->coeff_p2_0_m1 = t1799+t1804; + coeffs_dyy->coeff_m1_p1_m1 = t1742+t1782; + coeffs_dyy->coeff_0_p1_m1 = t1742+t1783; + coeffs_dyy->coeff_p1_p1_m1 = t1742+t1774; + coeffs_dyy->coeff_p2_p1_m1 = t1742+t1779; + coeffs_dyy->coeff_m1_p2_m1 = t1753+t1761+t1787; + coeffs_dyy->coeff_0_p2_m1 = t1748+t1768+t1802; + coeffs_dyy->coeff_p1_p2_m1 = t1753+t1760+t1790; + coeffs_dyy->coeff_p2_p2_m1 = t1753+t1789; + coeffs_dyy->coeff_m1_m1_0 = t1750+t1770; + coeffs_dyy->coeff_0_m1_0 = t1750+t1776; + coeffs_dyy->coeff_p1_m1_0 = t1750+t1777; + coeffs_dyy->coeff_p2_m1_0 = t1750+t1772; + coeffs_dyy->coeff_m1_0_0 = t1747+t1758+t1793; + coeffs_dyy->coeff_0_0_0 = t1749+t1754+t1793; + coeffs_dyy->coeff_p1_0_0 = t1792+t1799; + coeffs_dyy->coeff_p2_0_0 = t1751+t1764+t1793; + coeffs_dyy->coeff_m1_p1_0 = t1742+t1785; + coeffs_dyy->coeff_0_p1_0 = t1742+t1773; + coeffs_dyy->coeff_p1_p1_0 = t1742+t1778; + coeffs_dyy->coeff_p2_p1_0 = t1742+t1775; + coeffs_dyy->coeff_m1_p2_0 = t1744+t1768+t1787; + coeffs_dyy->coeff_0_p2_0 = t1748+t1760+t1803; + coeffs_dyy->coeff_p1_p2_0 = t1744+t1790; + coeffs_dyy->coeff_p2_p2_0 = t1744+t1756+t1789; + coeffs_dyy->coeff_m1_m1_p1 = t1750+t1775; + coeffs_dyy->coeff_0_m1_p1 = t1750+t1778; + coeffs_dyy->coeff_p1_m1_p1 = t1750+t1773; + coeffs_dyy->coeff_p2_m1_p1 = t1750+t1785; + coeffs_dyy->coeff_m1_0_p1 = t1744+t1754+t1797; + coeffs_dyy->coeff_0_0_p1 = t1796+t1799; + coeffs_dyy->coeff_p1_0_p1 = t1744+t1764+t1794; + coeffs_dyy->coeff_p2_0_p1 = t1788+t1798; + coeffs_dyy->coeff_m1_p1_p1 = t1742+t1772; + coeffs_dyy->coeff_0_p1_p1 = t1742+t1777; + coeffs_dyy->coeff_p1_p1_p1 = t1742+t1776; + coeffs_dyy->coeff_p2_p1_p1 = t1742+t1770; + coeffs_dyy->coeff_m1_p2_p1 = t1760+t1751+t1786; + coeffs_dyy->coeff_0_p2_p1 = t1752+t1786; + coeffs_dyy->coeff_p1_p2_p1 = t1749+t1756+t1786; + coeffs_dyy->coeff_p2_p2_p1 = t1757+t1747+t1786; + coeffs_dyy->coeff_m1_m1_p2 = t1750+t1779; + coeffs_dyy->coeff_0_m1_p2 = t1750+t1774; + coeffs_dyy->coeff_p1_m1_p2 = t1750+t1783; + coeffs_dyy->coeff_p2_m1_p2 = t1750+t1782; + coeffs_dyy->coeff_m1_0_p2 = t1769+t1753+t1797; + coeffs_dyy->coeff_0_0_p2 = t1746+t1764+t1791; + coeffs_dyy->coeff_p1_0_p2 = t1788+t1802; + coeffs_dyy->coeff_p2_0_p2 = RATIONAL(-7.0,80.0)+t1746+t1801; + coeffs_dyy->coeff_m1_p1_p2 = t1742+t1780; + coeffs_dyy->coeff_0_p1_p2 = t1742+t1781; + coeffs_dyy->coeff_p1_p1_p2 = t1742+t1771; + coeffs_dyy->coeff_p2_p1_p2 = t1742+t1784; + coeffs_dyy->coeff_m1_p2_p2 = t1745+t1787; + coeffs_dyy->coeff_0_p2_p2 = t1756+t1748+t1795; + coeffs_dyy->coeff_p1_p2_p2 = t1745+t1757+t1790; + coeffs_dyy->coeff_p2_p2_p2 = t1755+t1745+t1789; diff --git a/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dyz.compute.c b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dyz.compute.c new file mode 100644 index 0000000..8af310e --- /dev/null +++ b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dyz.compute.c @@ -0,0 +1,210 @@ +fp t1848; +fp t1813; +fp t1847; +fp t1815; +fp t1883; +fp t1809; +fp t1882; +fp t1846; +fp t1819; +fp t1820; +fp t1881; +fp t1849; +fp t1816; +fp t1880; +fp t1879; +fp t1807; +fp t1834; +fp t1878; +fp t1814; +fp t1877; +fp t1822; +fp t1811; +fp t1876; +fp t1838; +fp t1824; +fp t1806; +fp t1875; +fp t1817; +fp t1874; +fp t1818; +fp t1873; +fp t1833; +fp t1808; +fp t1872; +fp t1871; +fp t1870; +fp t1869; +fp t1868; +fp t1812; +fp t1867; +fp t1866; +fp t1865; +fp t1864; +fp t1863; +fp t1862; +fp t1861; +fp t1860; +fp t1810; +fp t1859; +fp t1858; +fp t1857; +fp t1856; +fp t1855; +fp t1854; +fp t1805; +fp t1853; +fp t1852; +fp t1851; +fp t1850; +fp t1843; +fp t1841; +fp t1840; +fp t1839; +fp t1837; +fp t1836; +fp t1835; +fp t1832; +fp t1831; +fp t1829; +fp t1828; +fp t1827; +fp t1823; +fp t1821; + t1848 = RATIONAL(1.0,80.0); + t1813 = t1848*y; + t1847 = RATIONAL(-1.0,80.0); + t1815 = t1847*z; + t1883 = t1813+t1815; + t1809 = RATIONAL(3.0,1000.0)*x; + t1882 = RATIONAL(-1.0,250.0)+t1809; + t1846 = RATIONAL(-3.0,80.0); + t1819 = t1846*z; + t1820 = t1846*y; + t1881 = t1819+t1820; + t1849 = RATIONAL(3.0,80.0); + t1816 = t1849*y; + t1880 = t1819+t1816; + t1879 = RATIONAL(2.0,125.0)+t1809; + t1807 = RATIONAL(1.0,1000.0)*x; + t1834 = RATIONAL(-3.0,1000.0); + t1878 = t1807+t1834; + t1814 = t1847*y; + t1877 = t1819+t1814; + t1822 = RATIONAL(-19.0,2000.0); + t1811 = RATIONAL(9.0,1000.0)*x; + t1876 = t1822+t1811; + t1838 = RATIONAL(-9.0,1000.0); + t1824 = RATIONAL(-27.0,1000.0); + t1806 = t1824*x; + t1875 = t1838+t1806; + t1817 = t1849*z; + t1874 = t1817+t1816; + t1818 = t1848*z; + t1873 = t1818+t1814; + t1833 = RATIONAL(-1.0,1000.0); + t1808 = t1833*x; + t1872 = RATIONAL(-1.0,500.0)+t1808; + t1871 = t1817+t1809; + t1870 = t1819+t1813; + t1869 = RATIONAL(37.0,2000.0)+t1809; + t1868 = t1809+t1816; + t1812 = t1838*x; + t1867 = t1812+t1816; + t1866 = t1818+t1820; + t1865 = t1814+t1815; + t1864 = t1811+t1813; + t1863 = t1824+t1811; + t1862 = t1817+t1820; + t1861 = RATIONAL(13.0,1000.0)+t1811; + t1860 = t1820+t1815; + t1810 = t1834*x; + t1859 = t1810+t1815; + t1858 = t1820+t1812; + t1857 = t1817+t1812; + t1856 = RATIONAL(49.0,2000.0)+t1812; + t1855 = t1810+t1814; + t1854 = t1818+t1813; + t1805 = RATIONAL(27.0,1000.0)*x; + t1853 = t1819+t1805; + t1852 = t1810+t1813; + t1851 = t1818+t1810; + t1850 = t1817+t1805; + t1843 = RATIONAL(-7.0,250.0); + t1841 = RATIONAL(-9.0,250.0); + t1840 = RATIONAL(-9.0,500.0); + t1839 = RATIONAL(-17.0,500.0); + t1837 = RATIONAL(11.0,500.0); + t1836 = RATIONAL(-1.0,2000.0); + t1835 = RATIONAL(-7.0,2000.0); + t1832 = RATIONAL(-37.0,1000.0); + t1831 = RATIONAL(-31.0,1000.0); + t1829 = RATIONAL(19.0,1000.0); + t1828 = RATIONAL(31.0,2000.0); + t1827 = RATIONAL(43.0,2000.0); + t1823 = RATIONAL(-13.0,2000.0); + t1821 = RATIONAL(-21.0,2000.0); + coeffs_dyz->coeff_m1_m1_m1 = t1806+RATIONAL(147.0,2000.0)+t1881; + coeffs_dyz->coeff_0_m1_m1 = RATIONAL(129.0,2000.0)+t1819+t1858; + coeffs_dyz->coeff_p1_m1_m1 = t1811+RATIONAL(111.0,2000.0)+t1881; + coeffs_dyz->coeff_p2_m1_m1 = t1820+RATIONAL(93.0,2000.0)+t1853; + coeffs_dyz->coeff_m1_0_m1 = t1836+t1815+t1867; + coeffs_dyz->coeff_0_0_m1 = t1835+t1816+t1859; + coeffs_dyz->coeff_p1_0_m1 = t1815+t1823+t1868; + coeffs_dyz->coeff_p2_0_m1 = t1815+t1816+t1876; + coeffs_dyz->coeff_m1_p1_m1 = t1811+t1832+t1818+t1816; + coeffs_dyz->coeff_0_p1_m1 = t1818+t1839+t1868; + coeffs_dyz->coeff_p1_p1_m1 = t1816+t1831+t1851; + coeffs_dyz->coeff_p2_p1_m1 = t1843+t1818+t1867; + coeffs_dyz->coeff_m1_p2_m1 = t1820+t1841+t1850; + coeffs_dyz->coeff_0_p2_m1 = t1862+t1863; + coeffs_dyz->coeff_p1_p2_m1 = t1840+t1820+t1857; + coeffs_dyz->coeff_p2_p2_m1 = t1862+t1875; + coeffs_dyz->coeff_m1_m1_0 = t1814+t1836+t1857; + coeffs_dyz->coeff_0_m1_0 = t1817+t1835+t1855; + coeffs_dyz->coeff_p1_m1_0 = t1823+t1814+t1871; + coeffs_dyz->coeff_p2_m1_0 = t1814+t1817+t1876; + coeffs_dyz->coeff_m1_0_0 = t1813+RATIONAL(-17.0,2000.0)+t1851; + coeffs_dyz->coeff_0_0_0 = t1822+t1808+t1854; + coeffs_dyz->coeff_p1_0_0 = t1821+t1807+t1854; + coeffs_dyz->coeff_p2_0_0 = t1809+RATIONAL(-23.0,2000.0)+t1854; + coeffs_dyz->coeff_m1_p1_0 = t1882+t1883; + coeffs_dyz->coeff_0_p1_0 = t1878+t1883; + coeffs_dyz->coeff_p1_p1_0 = t1872+t1883; + coeffs_dyz->coeff_p2_p1_0 = t1815+t1833+t1852; + coeffs_dyz->coeff_m1_p2_0 = t1861+t1877; + coeffs_dyz->coeff_0_p2_0 = t1877+t1879; + coeffs_dyz->coeff_p1_p2_0 = t1829+t1819+t1855; + coeffs_dyz->coeff_p2_p2_0 = t1812+t1837+t1877; + coeffs_dyz->coeff_m1_m1_p1 = t1817+t1832+t1864; + coeffs_dyz->coeff_0_m1_p1 = t1813+t1839+t1871; + coeffs_dyz->coeff_p1_m1_p1 = t1817+t1831+t1852; + coeffs_dyz->coeff_p2_m1_p1 = t1843+t1813+t1857; + coeffs_dyz->coeff_m1_0_p1 = t1873+t1882; + coeffs_dyz->coeff_0_0_p1 = t1873+t1878; + coeffs_dyz->coeff_p1_0_p1 = t1872+t1873; + coeffs_dyz->coeff_p2_0_p1 = t1833+t1814+t1851; + coeffs_dyz->coeff_m1_p1_p1 = t1815+RATIONAL(33.0,2000.0)+t1855; + coeffs_dyz->coeff_0_p1_p1 = t1828+t1808+t1865; + coeffs_dyz->coeff_p1_p1_p1 = RATIONAL(29.0,2000.0)+t1807+t1865; + coeffs_dyz->coeff_p2_p1_p1 = RATIONAL(27.0,2000.0)+t1809+t1865; + coeffs_dyz->coeff_m1_p2_p1 = t1856+t1870; + coeffs_dyz->coeff_0_p2_p1 = t1819+t1827+t1852; + coeffs_dyz->coeff_p1_p2_p1 = t1869+t1870; + coeffs_dyz->coeff_p2_p2_p1 = t1819+t1828+t1864; + coeffs_dyz->coeff_m1_m1_p2 = t1841+t1816+t1853; + coeffs_dyz->coeff_0_m1_p2 = t1863+t1880; + coeffs_dyz->coeff_p1_m1_p2 = t1840+t1819+t1867; + coeffs_dyz->coeff_p2_m1_p2 = t1875+t1880; + coeffs_dyz->coeff_m1_0_p2 = t1860+t1861; + coeffs_dyz->coeff_0_0_p2 = t1860+t1879; + coeffs_dyz->coeff_p1_0_p2 = t1829+t1820+t1859; + coeffs_dyz->coeff_p2_0_p2 = t1837+t1815+t1858; + coeffs_dyz->coeff_m1_p1_p2 = t1856+t1866; + coeffs_dyz->coeff_0_p1_p2 = t1827+t1820+t1851; + coeffs_dyz->coeff_p1_p1_p2 = t1866+t1869; + coeffs_dyz->coeff_p2_p1_p2 = t1828+t1811+t1866; + coeffs_dyz->coeff_m1_p2_p2 = RATIONAL(-3.0,2000.0)+t1806+t1874; + coeffs_dyz->coeff_0_p2_p2 = t1821+t1816+t1857; + coeffs_dyz->coeff_p1_p2_p2 = t1811+RATIONAL(-39.0,2000.0)+t1874; + coeffs_dyz->coeff_p2_p2_p2 = t1816+RATIONAL(-57.0,2000.0)+t1850; diff --git a/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dz.compute.c b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dz.compute.c new file mode 100644 index 0000000..dc79363 --- /dev/null +++ b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dz.compute.c @@ -0,0 +1,502 @@ +fp t1423; +fp t1429; +fp t1339; +fp t1430; +fp t1336; +fp t1477; +fp t1352; +fp t1517; +fp t1516; +fp t1337; +fp t1515; +fp t1428; +fp t1342; +fp t1419; +fp t1320; +fp t1328; +fp t1514; +fp t1348; +fp t1513; +fp t1356; +fp t1420; +fp t1321; +fp t1322; +fp t1512; +fp t1344; +fp t1511; +fp t1403; +fp t1500; +fp t1270; +fp t1510; +fp t1341; +fp t1509; +fp t1351; +fp t1508; +fp t1468; +fp t1507; +fp t1467; +fp t1506; +fp t1405; +fp t1306; +fp t1268; +fp t1505; +fp t1415; +fp t1331; +fp t1335; +fp t1469; +fp t1504; +fp t1347; +fp t1406; +fp t1308; +fp t1269; +fp t1503; +fp t1353; +fp t1502; +fp t1271; +fp t1345; +fp t1501; +fp t1499; +fp t1411; +fp t1319; +fp t1498; +fp t1346; +fp t1297; +fp t1280; +fp t1281; +fp t1497; +fp t1277; +fp t1496; +fp t1275; +fp t1276; +fp t1495; +fp t1279; +fp t1494; +fp t1493; +fp t1491; +fp t1489; +fp t1488; +fp t1278; +fp t1487; +fp t1350; +fp t1354; +fp t1486; +fp t1485; +fp t1484; +fp t1483; +fp t1482; +fp t1481; +fp t1274; +fp t1480; +fp t1479; +fp t1478; +fp t1475; +fp t1474; +fp t1472; +fp t1471; +fp t1470; +fp t1466; +fp t1465; +fp t1464; +fp t1463; +fp t1462; +fp t1461; +fp t1266; +fp t1460; +fp t1273; +fp t1459; +fp t1458; +fp t1457; +fp t1456; +fp t1455; +fp t1454; +fp t1453; +fp t1452; +fp t1267; +fp t1451; +fp t1450; +fp t1449; +fp t1448; +fp t1447; +fp t1446; +fp t1445; +fp t1444; +fp t1443; +fp t1442; +fp t1441; +fp t1272; +fp t1440; +fp t1439; +fp t1438; +fp t1437; +fp t1436; +fp t1435; +fp t1358; +fp t1434; +fp t1433; +fp t1432; +fp t1431; +fp t1422; +fp t1421; +fp t1418; +fp t1417; +fp t1416; +fp t1413; +fp t1412; +fp t1410; +fp t1408; +fp t1407; +fp t1404; +fp t1402; +fp t1401; +fp t1397; +fp t1396; +fp t1395; +fp t1394; +fp t1392; +fp t1391; +fp t1388; +fp t1387; +fp t1384; +fp t1383; +fp t1382; +fp t1380; +fp t1378; +fp t1377; +fp t1376; +fp t1371; +fp t1369; +fp t1340; +fp t1338; +fp t1334; +fp t1333; +fp t1332; +fp t1330; +fp t1329; +fp t1327; +fp t1326; +fp t1325; +fp t1324; +fp t1323; +fp t1318; +fp t1317; +fp t1316; +fp t1315; +fp t1314; +fp t1313; +fp t1312; +fp t1311; +fp t1310; +fp t1309; +fp t1307; +fp t1305; +fp t1304; +fp t1303; +fp t1302; +fp t1301; +fp t1300; +fp t1299; +fp t1298; +fp t1296; +fp t1295; +fp t1294; +fp t1293; +fp t1292; +fp t1291; +fp t1290; +fp t1289; +fp t1288; +fp t1287; +fp t1285; +fp t1284; +fp t1283; +fp t1282; + t1423 = RATIONAL(-3.0,160.0); + t1429 = x*x; + t1339 = t1423*t1429; + t1430 = y*y; + t1336 = t1423*t1430; + t1477 = t1339+t1336; + t1352 = RATIONAL(-9.0,80.0)*z; + t1517 = t1352+RATIONAL(-137.0,4000.0); + t1516 = RATIONAL(-4.0,375.0)+RATIONAL(-1.0,40.0)*z; + t1337 = RATIONAL(-11.0,80.0)*z; + t1515 = t1337+RATIONAL(-107.0,4000.0); + t1428 = z*z; + t1342 = RATIONAL(3.0,32.0)*t1428; + t1419 = RATIONAL(1.0,160.0); + t1320 = t1419*t1429; + t1328 = t1419*t1430; + t1514 = t1342+t1328+t1320; + t1348 = RATIONAL(1.0,20.0)*z; + t1513 = RATIONAL(-29.0,6000.0)+t1348; + t1356 = RATIONAL(-3.0,32.0)*t1428; + t1420 = RATIONAL(-1.0,160.0); + t1321 = t1420*t1430; + t1322 = t1420*t1429; + t1512 = t1356+t1321+t1322; + t1344 = RATIONAL(1.0,40.0)*z; + t1511 = RATIONAL(-17.0,3000.0)+t1344; + t1403 = RATIONAL(-27.0,1000.0); + t1500 = x*y; + t1270 = RATIONAL(9.0,1000.0)*t1500; + t1510 = t1403*y+t1270; + t1341 = RATIONAL(1.0,32.0)*t1428; + t1509 = t1341+t1477; + t1351 = RATIONAL(3.0,40.0)*z; + t1508 = RATIONAL(-97.0,3000.0)+t1351; + t1468 = t1342+t1321; + t1507 = t1322+t1468; + t1467 = t1320+t1356; + t1506 = t1467+t1328; + t1405 = RATIONAL(-3.0,1000.0); + t1306 = t1405*y; + t1268 = x*t1306; + t1505 = RATIONAL(53.0,1000.0)+t1351+t1268; + t1415 = RATIONAL(3.0,160.0); + t1331 = t1415*t1430; + t1335 = t1415*t1429; + t1469 = t1341+t1335; + t1504 = t1331+t1469; + t1347 = RATIONAL(1.0,80.0)*z; + t1406 = RATIONAL(-9.0,1000.0); + t1308 = t1406*x; + t1269 = y*t1308; + t1503 = t1347+t1269+RATIONAL(-137.0,12000.0); + t1353 = RATIONAL(7.0,80.0)*z; + t1502 = t1353+t1269+RATIONAL(-607.0,12000.0); + t1271 = RATIONAL(3.0,1000.0)*t1500; + t1345 = RATIONAL(3.0,80.0)*z; + t1501 = t1271+RATIONAL(323.0,4000.0)+t1345; + t1499 = y+x; + t1411 = RATIONAL(-1.0,1000.0); + t1319 = t1411*x; + t1498 = RATIONAL(-3.0,80.0)*z; + t1346 = RATIONAL(-1.0,80.0)*z; + t1297 = t1403*x; + t1280 = x*t1345; + t1281 = y*t1345; + t1497 = t1280+t1281; + t1277 = y*t1346; + t1496 = t1280+t1277; + t1275 = y*t1347; + t1276 = x*t1498; + t1495 = t1275+t1276; + t1279 = x*t1347; + t1494 = t1279+t1275; + t1493 = t1279+t1281; + t1491 = t1320+t1321; + t1489 = t1280+t1275; + t1488 = t1281+t1276; + t1278 = y*t1498; + t1487 = t1279+t1278; + t1350 = RATIONAL(1.0,16.0)*z; + t1354 = RATIONAL(-1.0,32.0)*t1428; + t1486 = t1350+t1354; + t1485 = t1322+t1328; + t1484 = t1280+t1278; + t1483 = t1279+t1277; + t1482 = t1342+t1268; + t1481 = t1276+t1278; + t1274 = x*t1346; + t1480 = t1275+t1274; + t1479 = t1276+t1277; + t1478 = t1335+t1336; + t1475 = t1278+t1274; + t1474 = t1339+t1331; + t1472 = t1331+t1335; + t1471 = t1277+t1274; + t1470 = t1281+t1274; + t1466 = t1344+t1497; + t1465 = t1270+t1478; + t1464 = t1348+t1494; + t1463 = t1270+t1474; + t1462 = t1351+t1471; + t1461 = RATIONAL(1.0,10.0)*z+t1481; + t1266 = RATIONAL(27.0,1000.0)*t1500; + t1460 = t1266+t1504; + t1273 = y*t1319; + t1459 = t1273+t1514; + t1458 = t1275+t1321+t1467; + t1457 = t1342+t1271+t1485; + t1456 = t1271+t1509; + t1455 = t1354+t1269+RATIONAL(53.0,12000.0)+t1345; + t1454 = t1279+t1356+t1485; + t1453 = t1269+t1507; + t1452 = RATIONAL(-3.0,20.0)*z+RATIONAL(-49.0,2000.0)+t1482; + t1267 = y*t1297; + t1451 = t1267+t1354+t1477; + t1450 = t1341+t1465; + t1449 = t1354+t1268+t1472; + t1448 = t1341+t1278+t1474; + t1447 = RATIONAL(-1.0,10.0)*z+RATIONAL(-79.0,2000.0)+t1482; + t1446 = t1270+t1506; + t1445 = t1356+t1268+t1348+RATIONAL(141.0,2000.0); + t1444 = t1354+t1276+t1474; + t1443 = t1346+t1341+t1269+RATIONAL(-167.0,12000.0); + t1442 = t1354+t1274+t1478; + t1441 = t1320+t1271+t1468; + t1272 = RATIONAL(1.0,1000.0)*t1500; + t1440 = t1272+t1512; + t1439 = t1353+t1356+t1271+RATIONAL(183.0,4000.0); + t1438 = RATIONAL(79.0,3000.0)+t1268+t1509; + t1437 = RATIONAL(-41.0,750.0)+t1267+t1504; + t1436 = RATIONAL(-691.0,12000.0)+t1271+t1472+t1486; + t1435 = t1266+RATIONAL(281.0,12000.0)+t1477+t1486; + t1358 = RATIONAL(-1.0,8.0)*z; + t1434 = RATIONAL(-33.0,2000.0)+t1358+t1270+t1507; + t1433 = t1350+RATIONAL(299.0,4000.0)+t1273+t1512; + t1432 = t1358+t1272+RATIONAL(-87.0,2000.0)+t1514; + t1431 = t1350+t1269+RATIONAL(191.0,4000.0)+t1506; + t1422 = RATIONAL(-7.0,250.0); + t1421 = RATIONAL(-9.0,500.0); + t1418 = RATIONAL(-9.0,250.0); + t1417 = RATIONAL(-1.0,500.0); + t1416 = RATIONAL(2.0,125.0); + t1413 = RATIONAL(-1.0,250.0); + t1412 = RATIONAL(-17.0,500.0); + t1410 = RATIONAL(-3.0,2000.0); + t1408 = RATIONAL(11.0,500.0); + t1407 = RATIONAL(-1.0,2000.0); + t1404 = RATIONAL(-7.0,2000.0); + t1402 = RATIONAL(37.0,2000.0); + t1401 = RATIONAL(13.0,1000.0); + t1397 = RATIONAL(49.0,2000.0); + t1396 = RATIONAL(27.0,2000.0); + t1395 = RATIONAL(-37.0,1000.0); + t1394 = RATIONAL(43.0,2000.0); + t1392 = RATIONAL(31.0,2000.0); + t1391 = RATIONAL(-19.0,2000.0); + t1388 = RATIONAL(93.0,2000.0); + t1387 = RATIONAL(-21.0,2000.0); + t1384 = RATIONAL(19.0,1000.0); + t1383 = RATIONAL(-23.0,2000.0); + t1382 = RATIONAL(-39.0,2000.0); + t1380 = RATIONAL(-17.0,2000.0); + t1378 = RATIONAL(-31.0,1000.0); + t1377 = RATIONAL(33.0,2000.0); + t1376 = RATIONAL(-13.0,2000.0); + t1371 = RATIONAL(111.0,2000.0); + t1369 = RATIONAL(129.0,2000.0); + t1340 = t1416*y; + t1338 = t1413*x; + t1334 = t1417*y; + t1333 = t1417*x; + t1332 = t1418*x; + t1330 = t1416*x; + t1329 = t1413*y; + t1327 = t1421*x; + t1326 = t1418*y; + t1325 = t1422*y; + t1324 = t1422*x; + t1323 = t1421*y; + t1318 = t1408*x; + t1317 = t1408*y; + t1316 = t1407*x; + t1315 = t1405*x; + t1314 = t1407*y; + t1313 = t1404*y; + t1312 = t1404*x; + t1311 = t1412*x; + t1310 = t1412*y; + t1309 = t1411*y; + t1307 = t1406*y; + t1305 = t1384*y; + t1304 = t1387*y; + t1303 = t1395*x; + t1302 = t1397*x; + t1301 = t1392*x; + t1300 = t1376*y; + t1299 = t1384*x; + t1298 = t1391*y; + t1296 = t1401*x; + t1295 = t1395*y; + t1294 = t1394*x; + t1293 = t1378*x; + t1292 = t1392*y; + t1291 = t1376*x; + t1290 = t1378*y; + t1289 = t1402*x; + t1288 = t1401*y; + t1287 = t1397*y; + t1285 = t1391*x; + t1284 = t1387*x; + t1283 = t1394*y; + t1282 = t1402*y; + coeffs_dz->coeff_m1_m1_m1 = RATIONAL(-47.0,1500.0)+t1499*RATIONAL(147.0, +2000.0)+t1451+t1461; + coeffs_dz->coeff_0_m1_m1 = t1316+t1278+t1369*y+t1442+t1502; + coeffs_dz->coeff_p1_m1_m1 = t1371*y+t1354+t1303+t1465+t1487+t1508; + coeffs_dz->coeff_p2_m1_m1 = t1388*y+t1332+t1435+t1484; + coeffs_dz->coeff_m1_0_m1 = t1277+t1314+t1369*x+t1444+t1502; + coeffs_dz->coeff_0_0_m1 = t1312+t1313+RATIONAL(-437.0,6000.0)+t1449+t1462 +; + coeffs_dz->coeff_p1_0_m1 = t1311+t1300+t1436+t1483; + coeffs_dz->coeff_p2_0_m1 = t1354+t1297+t1298+t1463+t1496+t1513; + coeffs_dz->coeff_m1_p1_m1 = t1270+t1371*x+t1275+t1295+t1444+t1508; + coeffs_dz->coeff_0_p1_m1 = t1310+t1291+t1436+t1480; + coeffs_dz->coeff_p1_p1_m1 = RATIONAL(-17.0,375.0)+t1293+t1290+t1449+t1464 +; + coeffs_dz->coeff_p2_p1_m1 = t1327+t1325+t1455+t1474+t1489; + coeffs_dz->coeff_m1_p2_m1 = t1326+t1388*x+t1435+t1488; + coeffs_dz->coeff_0_p2_m1 = t1285+t1281+t1442+t1510+t1513; + coeffs_dz->coeff_p1_p2_m1 = t1324+t1323+t1455+t1478+t1493; + coeffs_dz->coeff_p2_p2_m1 = t1308+RATIONAL(307.0,6000.0)+t1307+t1451+ +t1466; + coeffs_dz->coeff_m1_m1_0 = t1316+t1314+RATIONAL(-39.0,4000.0)+RATIONAL( +-13.0,80.0)*z+t1453+t1497; + coeffs_dz->coeff_0_m1_0 = t1380*x+t1313+t1452+t1491+t1493; + coeffs_dz->coeff_p1_m1_0 = t1338+t1300+t1441+t1470+t1515; + coeffs_dz->coeff_p2_m1_0 = t1298+t1296+t1434+t1488; + coeffs_dz->coeff_m1_0_0 = t1380*y+t1312+t1452+t1485+t1489; + coeffs_dz->coeff_0_0_0 = t1337+t1298+RATIONAL(-161.0,4000.0)+t1285+t1459+ +t1494; + coeffs_dz->coeff_p1_0_0 = t1315+t1304+t1432+t1480; + coeffs_dz->coeff_p2_0_0 = t1383*y+t1330+t1457+t1495+t1517; + coeffs_dz->coeff_m1_p1_0 = t1329+t1291+t1457+t1496+t1515; + coeffs_dz->coeff_0_p1_0 = t1284+t1306+t1432+t1483; + coeffs_dz->coeff_p1_p1_0 = t1352+RATIONAL(-191.0,4000.0)+t1334+t1333+ +t1459+t1471; + coeffs_dz->coeff_p2_p1_0 = t1299+t1309+t1447+t1479+t1485; + coeffs_dz->coeff_m1_p2_0 = t1285+t1288+t1434+t1484; + coeffs_dz->coeff_0_p2_0 = t1383*x+t1340+t1441+t1487+t1517; + coeffs_dz->coeff_p1_p2_0 = t1305+t1319+t1447+t1475+t1491; + coeffs_dz->coeff_p2_p2_0 = t1318+RATIONAL(-7.0,80.0)*z+t1317+RATIONAL( +-129.0,4000.0)+t1453+t1481; + coeffs_dz->coeff_m1_m1_p1 = t1303+t1295+RATIONAL(157.0,2000.0)+t1446+ +t1466; + coeffs_dz->coeff_0_m1_p1 = t1281+t1310+t1338+t1454+t1501; + coeffs_dz->coeff_p1_m1_p1 = t1377*x+t1290+t1445+t1470+t1485; + coeffs_dz->coeff_p2_m1_p1 = t1302+t1325+t1431+t1488; + coeffs_dz->coeff_m1_0_p1 = t1280+t1329+t1311+t1458+t1501; + coeffs_dz->coeff_0_0_p1 = RATIONAL(21.0,250.0)+t1315+t1306+t1440+t1464; + coeffs_dz->coeff_p1_0_p1 = t1334+t1301+t1433+t1480; + coeffs_dz->coeff_p2_0_p1 = t1309+t1294+t1276+t1458+t1505; + coeffs_dz->coeff_m1_p1_p1 = t1377*y+t1293+t1445+t1491+t1496; + coeffs_dz->coeff_0_p1_p1 = t1292+t1333+t1433+t1483; + coeffs_dz->coeff_p1_p1_p1 = RATIONAL(133.0,2000.0)+t1499*RATIONAL(29.0, +2000.0)+t1440+t1462; + coeffs_dz->coeff_p2_p1_p1 = t1289+t1396*y+t1439+t1479+t1491; + coeffs_dz->coeff_m1_p2_p1 = t1287+t1324+t1431+t1484; + coeffs_dz->coeff_0_p2_p1 = t1319+t1283+t1278+t1454+t1505; + coeffs_dz->coeff_p1_p2_p1 = t1396*x+t1282+t1439+t1475+t1485; + coeffs_dz->coeff_p2_p2_p1 = t1301+RATIONAL(13.0,500.0)+t1292+t1446+t1461; + coeffs_dz->coeff_m1_m1_p2 = t1345+t1326+t1332+RATIONAL(-449.0,12000.0)+ +t1460+t1481; + coeffs_dz->coeff_0_m1_p2 = t1274+t1296+t1448+t1510+t1511; + coeffs_dz->coeff_p1_m1_p2 = t1302+t1323+t1279+t1448+t1503; + coeffs_dz->coeff_p2_m1_p2 = t1410*x+t1307+t1437+t1484; + coeffs_dz->coeff_m1_0_p2 = t1288+t1297+t1450+t1479+t1511; + coeffs_dz->coeff_0_0_p2 = RATIONAL(349.0,12000.0)+t1347+t1340+t1330+t1456 ++t1471; + coeffs_dz->coeff_p1_0_p2 = t1305+t1294+t1438+t1483; + coeffs_dz->coeff_p2_0_p2 = t1284+t1317+t1443+t1478+t1496; + coeffs_dz->coeff_m1_p1_p2 = t1327+t1336+t1287+t1469+t1495+t1503; + coeffs_dz->coeff_0_p1_p2 = t1283+t1299+t1438+t1480; + coeffs_dz->coeff_p1_p1_p2 = t1282+RATIONAL(319.0,12000.0)+t1346+t1289+ +t1456+t1494; + coeffs_dz->coeff_p2_p1_p2 = t1382*x+t1292+t1450+t1489+t1516; + coeffs_dz->coeff_m1_p2_p2 = t1308+t1410*y+t1437+t1488; + coeffs_dz->coeff_0_p2_p2 = t1304+t1318+t1443+t1470+t1474; + coeffs_dz->coeff_p1_p2_p2 = t1382*y+t1341+t1301+t1463+t1493+t1516; + coeffs_dz->coeff_p2_p2_p2 = t1498+RATIONAL(-539.0,12000.0)+t1499*RATIONAL +(-57.0,2000.0)+t1460+t1497; diff --git a/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dzz.compute.c b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dzz.compute.c new file mode 100644 index 0000000..8de80ca --- /dev/null +++ b/src/Lagrange-maximum-degree/3d.coeffs/3d.cube.order3.smooth0/coeffs-dzz.compute.c @@ -0,0 +1,192 @@ +fp t1886; +fp t1896; +fp t1949; +fp t1897; +fp t1884; +fp t1905; +fp t1892; +fp t1948; +fp t1885; +fp t1900; +fp t1888; +fp t1947; +fp t1906; +fp t1890; +fp t1946; +fp t1895; +fp t1945; +fp t1903; +fp t1889; +fp t1944; +fp t1894; +fp t1943; +fp t1942; +fp t1891; +fp t1941; +fp t1901; +fp t1940; +fp t1939; +fp t1938; +fp t1937; +fp t1936; +fp t1935; +fp t1934; +fp t1933; +fp t1932; +fp t1931; +fp t1930; +fp t1929; +fp t1928; +fp t1907; +fp t1927; +fp t1910; +fp t1926; +fp t1925; +fp t1924; +fp t1923; +fp t1922; +fp t1921; +fp t1920; +fp t1919; +fp t1904; +fp t1918; +fp t1917; +fp t1916; +fp t1915; +fp t1914; +fp t1913; +fp t1912; +fp t1911; +fp t1909; +fp t1899; +fp t1898; +fp t1893; +fp t1887; + t1886 = RATIONAL(3.0,16.0)*z; + t1896 = RATIONAL(-11.0,80.0); + t1949 = t1886+t1896; + t1897 = RATIONAL(-3.0,80.0); + t1884 = t1897*y; + t1905 = RATIONAL(1.0,80.0); + t1892 = t1905*x; + t1948 = t1884+t1892; + t1885 = t1897*x; + t1900 = RATIONAL(-1.0,80.0); + t1888 = t1900*y; + t1947 = t1885+t1888; + t1906 = RATIONAL(3.0,80.0); + t1890 = t1906*x; + t1946 = t1884+t1890; + t1895 = t1905*y; + t1945 = t1890+t1895; + t1903 = RATIONAL(1.0,16.0); + t1889 = t1903*z; + t1944 = t1889+t1900; + t1894 = t1900*x; + t1943 = t1895+t1894; + t1942 = t1888+t1892; + t1891 = t1906*y; + t1941 = t1889+t1891; + t1901 = RATIONAL(1.0,40.0); + t1940 = t1889+t1901; + t1939 = t1885+t1895; + t1938 = t1886+t1885; + t1937 = t1886+RATIONAL(-1.0,10.0); + t1936 = t1888+t1894; + t1935 = t1886+t1894; + t1934 = t1886+t1892; + t1933 = t1884+t1894; + t1932 = t1885+t1884; + t1931 = t1888+t1890; + t1930 = t1889+t1905; + t1929 = t1890+t1891; + t1928 = t1892+t1895; + t1907 = RATIONAL(7.0,80.0); + t1927 = t1907+t1933; + t1910 = RATIONAL(1.0,20.0); + t1926 = t1910+t1891+t1894; + t1925 = t1907+t1947; + t1924 = t1903+t1942; + t1923 = t1901+t1929; + t1922 = RATIONAL(1.0,10.0)+t1932; + t1921 = t1903+t1946; + t1920 = t1885+t1891+t1903; + t1919 = t1903+t1943; + t1904 = RATIONAL(3.0,40.0); + t1918 = t1904+t1939; + t1917 = t1904+t1948; + t1916 = t1904+t1936; + t1915 = t1910+t1931; + t1914 = t1910+t1928; + t1913 = t1891+t1892+t1906; + t1912 = t1906+t1945; + t1911 = RATIONAL(-1.0,8.0); + t1909 = RATIONAL(-9.0,80.0); + t1899 = RATIONAL(-3.0,20.0); + t1898 = RATIONAL(-1.0,40.0); + t1893 = RATIONAL(-3.0,16.0)*z; + t1887 = RATIONAL(-1.0,16.0)*z; + coeffs_dzz->coeff_m1_m1_m1 = t1887+t1922; + coeffs_dzz->coeff_0_m1_m1 = t1887+t1927; + coeffs_dzz->coeff_p1_m1_m1 = t1887+t1917; + coeffs_dzz->coeff_p2_m1_m1 = t1887+t1921; + coeffs_dzz->coeff_m1_0_m1 = t1887+t1925; + coeffs_dzz->coeff_0_0_m1 = t1887+t1916; + coeffs_dzz->coeff_p1_0_m1 = t1887+t1924; + coeffs_dzz->coeff_p2_0_m1 = t1887+t1915; + coeffs_dzz->coeff_m1_p1_m1 = t1887+t1918; + coeffs_dzz->coeff_0_p1_m1 = t1887+t1919; + coeffs_dzz->coeff_p1_p1_m1 = t1887+t1914; + coeffs_dzz->coeff_p2_p1_m1 = t1887+t1912; + coeffs_dzz->coeff_m1_p2_m1 = t1887+t1920; + coeffs_dzz->coeff_0_p2_m1 = t1887+t1926; + coeffs_dzz->coeff_p1_p2_m1 = t1887+t1913; + coeffs_dzz->coeff_p2_p2_m1 = t1887+t1923; + coeffs_dzz->coeff_m1_m1_0 = RATIONAL(-13.0,80.0)+t1886+t1929; + coeffs_dzz->coeff_0_m1_0 = t1899+t1891+t1934; + coeffs_dzz->coeff_p1_m1_0 = t1896+t1891+t1935; + coeffs_dzz->coeff_p2_m1_0 = t1911+t1891+t1938; + coeffs_dzz->coeff_m1_0_0 = t1886+t1899+t1945; + coeffs_dzz->coeff_0_0_0 = t1928+t1949; + coeffs_dzz->coeff_p1_0_0 = t1911+t1895+t1935; + coeffs_dzz->coeff_p2_0_0 = t1909+t1895+t1938; + coeffs_dzz->coeff_m1_p1_0 = t1931+t1949; + coeffs_dzz->coeff_0_p1_0 = t1911+t1888+t1934; + coeffs_dzz->coeff_p1_p1_0 = t1888+t1909+t1935; + coeffs_dzz->coeff_p2_p1_0 = t1937+t1947; + coeffs_dzz->coeff_m1_p2_0 = t1886+t1911+t1946; + coeffs_dzz->coeff_0_p2_0 = t1884+t1909+t1934; + coeffs_dzz->coeff_p1_p2_0 = t1933+t1937; + coeffs_dzz->coeff_p2_p2_0 = t1886+RATIONAL(-7.0,80.0)+t1932; + coeffs_dzz->coeff_m1_m1_p1 = t1893+t1923; + coeffs_dzz->coeff_0_m1_p1 = t1893+t1913; + coeffs_dzz->coeff_p1_m1_p1 = t1893+t1926; + coeffs_dzz->coeff_p2_m1_p1 = t1893+t1920; + coeffs_dzz->coeff_m1_0_p1 = t1893+t1912; + coeffs_dzz->coeff_0_0_p1 = t1893+t1914; + coeffs_dzz->coeff_p1_0_p1 = t1893+t1919; + coeffs_dzz->coeff_p2_0_p1 = t1893+t1918; + coeffs_dzz->coeff_m1_p1_p1 = t1893+t1915; + coeffs_dzz->coeff_0_p1_p1 = t1893+t1924; + coeffs_dzz->coeff_p1_p1_p1 = t1893+t1916; + coeffs_dzz->coeff_p2_p1_p1 = t1893+t1925; + coeffs_dzz->coeff_m1_p2_p1 = t1893+t1921; + coeffs_dzz->coeff_0_p2_p1 = t1893+t1917; + coeffs_dzz->coeff_p1_p2_p1 = t1893+t1927; + coeffs_dzz->coeff_p2_p2_p1 = t1893+t1922; + coeffs_dzz->coeff_m1_m1_p2 = t1889+t1906+t1932; + coeffs_dzz->coeff_0_m1_p2 = t1933+t1940; + coeffs_dzz->coeff_p1_m1_p2 = t1930+t1948; + coeffs_dzz->coeff_p2_m1_p2 = t1889+t1946; + coeffs_dzz->coeff_m1_0_p2 = t1940+t1947; + coeffs_dzz->coeff_0_0_p2 = t1930+t1936; + coeffs_dzz->coeff_p1_0_p2 = t1889+t1942; + coeffs_dzz->coeff_p2_0_p2 = t1931+t1944; + coeffs_dzz->coeff_m1_p1_p2 = t1930+t1939; + coeffs_dzz->coeff_0_p1_p2 = t1889+t1943; + coeffs_dzz->coeff_p1_p1_p2 = t1928+t1944; + coeffs_dzz->coeff_p2_p1_p2 = t1898+t1889+t1945; + coeffs_dzz->coeff_m1_p2_p2 = t1885+t1941; + coeffs_dzz->coeff_0_p2_p2 = t1900+t1894+t1941; + coeffs_dzz->coeff_p1_p2_p2 = t1898+t1892+t1941; + coeffs_dzz->coeff_p2_p2_p2 = t1897+t1889+t1929; |