aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c')
-rw-r--r--src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c581
1 files changed, 294 insertions, 287 deletions
diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c
index 1c237bc..0eefdd1 100644
--- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c
+++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c
@@ -1,287 +1,294 @@
-fp t519;
-fp t520;
-fp t522;
-fp t525;
-fp t616;
-fp t513;
-fp t512;
-fp t615;
-fp t507;
-fp t592;
-fp t521;
-fp t523;
-fp t524;
-fp t492;
-fp t549;
-fp t614;
-fp t490;
-fp t588;
-fp t493;
-fp t570;
-fp t613;
-fp t574;
-fp t488;
-fp t496;
-fp t556;
-fp t526;
-fp t563;
-fp t582;
-fp t508;
-fp t552;
-fp t484;
-fp t612;
-fp t557;
-fp t554;
-fp t611;
-fp t487;
-fp t559;
-fp t590;
-fp t562;
-fp t489;
-fp t579;
-fp t610;
-fp t583;
-fp t565;
-fp t553;
-fp t609;
-fp t581;
-fp t608;
-fp t558;
-fp t517;
-fp t518;
-fp t576;
-fp t580;
-fp t607;
-fp t584;
-fp t606;
-fp t586;
-fp t605;
-fp t561;
-fp t564;
-fp t604;
-fp t578;
-fp t603;
-fp t577;
-fp t602;
-fp t560;
-fp t568;
-fp t551;
-fp t601;
-fp t585;
-fp t550;
-fp t600;
-fp t599;
-fp t589;
-fp t555;
-fp t598;
-fp t597;
-fp t596;
-fp t595;
-fp t594;
-fp t593;
-fp t591;
-fp t587;
-fp t506;
-fp t479;
-fp t477;
-fp t514;
-fp t480;
-fp t483;
-fp t510;
-fp t572;
-fp t571;
-fp t497;
-fp t569;
-fp t567;
-fp t548;
-fp t547;
-fp t495;
-fp t546;
-fp t499;
-fp t502;
-fp t545;
-fp t544;
-fp t543;
-fp t505;
-fp t500;
-fp t503;
-fp t542;
-fp t501;
-fp t504;
-fp t498;
-fp t541;
-fp t540;
-fp t539;
-fp t538;
-fp t537;
-fp t536;
-fp t535;
-fp t534;
-fp t533;
-fp t532;
-fp t531;
-fp t516;
-fp t515;
-fp t511;
-fp t509;
-fp t494;
-fp t491;
- t519 = RATIONAL(-1.0,9.0);
- t520 = RATIONAL(4.0,9.0);
- t522 = y*y;
- t525 = t522*t522;
- t616 = t520*t522+t519*t525;
- t513 = RATIONAL(1.0,6.0);
- t512 = RATIONAL(2.0,3.0);
- t615 = t512+t513*t525;
- t507 = RATIONAL(-1.0,18.0);
- t592 = x*y;
- t521 = x*x;
- t523 = t521*x;
- t524 = y*t522;
- t492 = RATIONAL(-1.0,72.0);
- t549 = t524*t492;
- t614 = t507*t592+t523*t549;
- t490 = RATIONAL(1.0,18.0);
- t588 = x*t490;
- t493 = RATIONAL(1.0,72.0);
- t570 = t523*t493;
- t613 = t524*t570+y*t588;
- t574 = t521*t525;
- t488 = RATIONAL(-1.0,144.0);
- t496 = RATIONAL(-1.0,36.0);
- t556 = t521*t496;
- t526 = t521*t521;
- t563 = t526*t525;
- t582 = t488*t563+t522*t556;
- t508 = RATIONAL(1.0,36.0);
- t552 = t526*t508;
- t484 = RATIONAL(1.0,144.0);
- t612 = t484*t574+t582+t522*t552;
- t557 = t522*t526;
- t554 = t525*t508;
- t611 = t484*t557+t582+t521*t554;
- t487 = RATIONAL(-1.0,288.0);
- t559 = t487*t521;
- t590 = y*t526;
- t562 = t526*t524;
- t489 = RATIONAL(1.0,288.0);
- t579 = t521*t489;
- t610 = t487*t562+y*t559+t524*t579+t489*t590;
- t583 = t523*t484;
- t565 = t524*t484;
- t553 = t523*t488;
- t609 = t488*t592+x*t565+y*t583+t524*t553;
- t581 = t522*t523;
- t608 = t519*t581+t523*t554+t616*x;
- t558 = t524*t496;
- t517 = RATIONAL(-4.0,9.0);
- t518 = RATIONAL(1.0,9.0);
- t576 = t518*t523;
- t580 = t518*t524;
- t607 = y*t576+t523*t558+x*t580+t517*t592;
- t584 = t521*t522;
- t606 = (t574+t557)*RATIONAL(-1.0,576.0)+(t584+t563)*RATIONAL(1.0,576.0);
- t586 = x*t524;
- t605 = t490*t586+y*t570+t614;
- t561 = t507*t521;
- t564 = t521*t490;
- t604 = t526*t549+t524*t564+y*t561+t493*t590;
- t578 = t524*t488;
- t603 = t484*t592+x*t578+y*t553+t523*t565;
- t577 = t523*t490;
- t602 = t493*t586+y*t577+t614;
- t560 = t525*t493;
- t568 = t507*t523;
- t551 = t525*t492;
- t601 = x*t551+t523*t560+(t568+t588)*t522;
- t585 = y*t523;
- t550 = t519*t524;
- t600 = x*t550+t520*t592+t519*t585+t508*t523*t524;
- t599 = t487*t590+t489*t562+y*t579+t524*t559;
- t589 = x*t522;
- t555 = t525*t496;
- t598 = t522*t576+t517*t589+t523*t555+t518*x*t525;
- t597 = x*t560+t523*t551+t522*t577+t507*t589;
- t596 = x*t549+y*t568+t613;
- t595 = t492*t590+t524*t561+y*t564+t493*t562;
- t594 = t492*t585+t507*t586+t613;
- t593 = t525*t552+t519*t557+t616*t521;
- t591 = y*t521;
- t587 = x*t487;
- t506 = RATIONAL(1.0,24.0);
- t479 = t506*t525;
- t477 = t506*t526;
- t514 = RATIONAL(-1.0,6.0);
- t480 = t514*t526;
- t483 = t514*t525;
- t510 = RATIONAL(-5.0,6.0);
- t572 = t510*t522;
- t571 = RATIONAL(1.0,96.0)*t563+RATIONAL(5.0,96.0)*t584;
- t497 = RATIONAL(-1.0,24.0);
- t569 = t497*t563+t521*t572;
- t567 = t525*t489;
- t548 = t487*t523;
- t547 = t522*t489;
- t495 = RATIONAL(5.0,24.0);
- t546 = t495*t557+t480+t569+t615*t521;
- t499 = RATIONAL(-5.0,96.0);
- t502 = RATIONAL(-1.0,96.0);
- t545 = t499*t557+t497*t521+t477+t502*t574+t571;
- t544 = t499*t574+t497*t522+t479+t502*t557+t571;
- t543 = t513*t557+t512*t522+t495*t574+t483+t569;
- t505 = RATIONAL(-1.0,12.0);
- t500 = RATIONAL(-1.0,48.0);
- t503 = RATIONAL(5.0,48.0);
- t542 = t505+t503*t521+t500*t526;
- t501 = RATIONAL(1.0,48.0);
- t504 = RATIONAL(1.0,12.0);
- t498 = RATIONAL(-5.0,48.0);
- t541 = t504+t501*t526+t498*t521;
- t540 = t500*t525+t503*t522+t505;
- t539 = t501*t525+t498*t522+t504;
- t538 = t612+t521*t565+t488*t562+(t552+t556)*y;
- t537 = t496*t590+t508*t591+t484*t562+t521*t578+t612;
- t536 = t525*t583+t508*t589+t522*t553+x*t555+t611;
- t535 = t522*t587+t523*t547+x*t567+t525*t548+t606;
- t534 = t525*t587+t522*t548+t523*t567+x*t547+t606;
- t533 = t496*t589+x*t554+t525*t553+t484*t581+t611;
- t532 = t517*t591+t526*t558+t518*t590+t521*t580+t593;
- t531 = t521*t550+t520*t591+t519*t590+t524*t552+t593;
- t516 = RATIONAL(-2.0,3.0);
- t515 = RATIONAL(5.0,6.0);
- t511 = RATIONAL(1.0,4.0);
- t509 = RATIONAL(-5.0,4.0);
- t494 = RATIONAL(-5.0,24.0);
- t491 = RATIONAL(-5.0,16.0);
- coeffs_I->coeff_m2_m2 = t535+t603+t610;
- coeffs_I->coeff_m1_m2 = t536+t595+t605;
- coeffs_I->coeff_0_m2 = t542*t524+t541*y+t544;
- coeffs_I->coeff_p1_m2 = t533+t594+t595;
- coeffs_I->coeff_p2_m2 = t534+t609+t610;
- coeffs_I->coeff_m2_m1 = t537+t601+t602;
- coeffs_I->coeff_m1_m1 = t532+t598+t600;
- coeffs_I->coeff_0_m1 = (t513+t477+t494*t521)*t524+(t516+t515*t521+t480)*y
-+t543;
- coeffs_I->coeff_p1_m1 = t532+t607+t608;
- coeffs_I->coeff_p2_m1 = t537+t596+t597;
- coeffs_I->coeff_m2_0 = t540*t523+t539*x+t545;
- coeffs_I->coeff_m1_0 = (t479+t494*t522+t513)*t523+(t483+t516+t515*t522)*x
-+t546;
- coeffs_I->coeff_0_0 = t511*t525+RATIONAL(1.0,1.0)+t509*t522+(t491*t522+
-t511+RATIONAL(1.0,16.0)*t525)*t526+(t491*t525+RATIONAL(25.0,16.0)*t522+t509)*
-t521;
- coeffs_I->coeff_p1_0 = (t497*t525+t495*t522+t514)*t523+(t572+t615)*x+t546
-;
- coeffs_I->coeff_p2_0 = t539*t523+t540*x+t545;
- coeffs_I->coeff_m2_p1 = t538+t596+t601;
- coeffs_I->coeff_m1_p1 = t531+t598+t607;
- coeffs_I->coeff_0_p1 = (t497*t526+t495*t521+t514)*t524+(t512+t510*t521+
-t513*t526)*y+t543;
- coeffs_I->coeff_p1_p1 = t531+t600+t608;
- coeffs_I->coeff_p2_p1 = t538+t597+t602;
- coeffs_I->coeff_m2_p2 = t535+t599+t609;
- coeffs_I->coeff_m1_p2 = t536+t594+t604;
- coeffs_I->coeff_0_p2 = t541*t524+t542*y+t544;
- coeffs_I->coeff_p1_p2 = t533+t604+t605;
- coeffs_I->coeff_p2_p2 = t534+t599+t603;
+fp t652;
+fp t654;
+fp t674;
+fp t672;
+fp t771;
+fp t666;
+fp t667;
+fp t770;
+fp t675;
+fp t677;
+fp t655;
+fp t676;
+fp t708;
+fp t747;
+fp t769;
+fp t699;
+fp t651;
+fp t742;
+fp t768;
+fp t636;
+fp t707;
+fp t653;
+fp t720;
+fp t638;
+fp t732;
+fp t673;
+fp t719;
+fp t656;
+fp t735;
+fp t767;
+fp t724;
+fp t700;
+fp t766;
+fp t662;
+fp t713;
+fp t726;
+fp t744;
+fp t663;
+fp t746;
+fp t765;
+fp t722;
+fp t743;
+fp t739;
+fp t764;
+fp t741;
+fp t763;
+fp t737;
+fp t635;
+fp t640;
+fp t745;
+fp t701;
+fp t702;
+fp t762;
+fp t730;
+fp t761;
+fp t740;
+fp t733;
+fp t721;
+fp t728;
+fp t760;
+fp t705;
+fp t717;
+fp t759;
+fp t664;
+fp t738;
+fp t706;
+fp t665;
+fp t758;
+fp t757;
+fp t756;
+fp t714;
+fp t755;
+fp t754;
+fp t753;
+fp t729;
+fp t752;
+fp t703;
+fp t712;
+fp t710;
+fp t751;
+fp t750;
+fp t749;
+fp t748;
+fp t671;
+fp t711;
+fp t646;
+fp t715;
+fp t736;
+fp t734;
+fp t731;
+fp t727;
+fp t723;
+fp t658;
+fp t718;
+fp t709;
+fp t668;
+fp t632;
+fp t641;
+fp t629;
+fp t633;
+fp t627;
+fp t698;
+fp t697;
+fp t644;
+fp t643;
+fp t696;
+fp t695;
+fp t694;
+fp t642;
+fp t647;
+fp t648;
+fp t693;
+fp t645;
+fp t649;
+fp t650;
+fp t692;
+fp t691;
+fp t690;
+fp t689;
+fp t688;
+fp t687;
+fp t686;
+fp t685;
+fp t684;
+fp t683;
+fp t682;
+fp t670;
+fp t669;
+fp t661;
+fp t660;
+fp t659;
+fp t657;
+ t652 = RATIONAL(1.0,72.0);
+ t654 = RATIONAL(-1.0,18.0);
+ t674 = x*x;
+ t672 = t674*t674;
+ t771 = t654*t674+t652*t672;
+ t666 = RATIONAL(1.0,6.0);
+ t667 = RATIONAL(2.0,3.0);
+ t770 = t666*t672+t667;
+ t675 = y*y;
+ t677 = t675*y;
+ t655 = RATIONAL(-1.0,72.0);
+ t676 = t674*x;
+ t708 = t655*t676;
+ t747 = x*y;
+ t769 = t677*t708+t654*t747;
+ t699 = t652*t676;
+ t651 = RATIONAL(1.0,18.0);
+ t742 = x*t651;
+ t768 = t677*t699+y*t742;
+ t636 = RATIONAL(1.0,144.0);
+ t707 = t636*t674;
+ t653 = RATIONAL(1.0,36.0);
+ t720 = t653*t672;
+ t638 = RATIONAL(-1.0,144.0);
+ t732 = t674*t675;
+ t673 = t675*t675;
+ t719 = t673*t672;
+ t656 = RATIONAL(-1.0,36.0);
+ t735 = t638*t719+t656*t732;
+ t767 = t673*t707+t735+t675*t720;
+ t724 = t672*t675;
+ t700 = t653*t674;
+ t766 = t636*t724+t673*t700+t735;
+ t662 = RATIONAL(1.0,9.0);
+ t713 = t677*t662;
+ t726 = t656*t672;
+ t744 = y*t674;
+ t663 = RATIONAL(-4.0,9.0);
+ t746 = y*t672;
+ t765 = t674*t713+t662*t746+t663*t744+t677*t726;
+ t722 = t677*t676;
+ t743 = y*t676;
+ t739 = x*t663;
+ t764 = x*t713+t662*t743+t656*t722+y*t739;
+ t741 = t677*x;
+ t763 = t651*t743+t652*t741+t769;
+ t737 = x*t673;
+ t635 = RATIONAL(-1.0,288.0);
+ t640 = RATIONAL(1.0,288.0);
+ t745 = x*t675;
+ t701 = t635*t673;
+ t702 = t640*t676;
+ t762 = t675*t702+t635*t745+t676*t701+t640*t737;
+ t730 = t673*t674;
+ t761 = (t730+t724)*RATIONAL(-1.0,576.0)+(t719+t732)*RATIONAL(1.0,576.0);
+ t740 = x*t636;
+ t733 = t677*t638;
+ t721 = t636*t676;
+ t728 = t638*t676;
+ t760 = x*t733+t677*t721+(t728+t740)*y;
+ t705 = t655*t673;
+ t717 = t675*t676;
+ t759 = t673*t699+x*t705+t675*t742+t654*t717;
+ t664 = RATIONAL(4.0,9.0);
+ t738 = y*t664;
+ t706 = t653*t676;
+ t665 = RATIONAL(-1.0,9.0);
+ t758 = x*t738+t677*t706+(t743+t741)*t665;
+ t757 = t673*t702+t635*t717+x*t701+t640*t745;
+ t756 = t677*t740+y*t721+(t722+t747)*t638;
+ t714 = t665*t674;
+ t755 = t665*t746+t674*t738+(t714+t720)*t677;
+ t754 = t654*t745+t676*t705+t652*t737+t651*t717;
+ t753 = t655*t741+t654*t743+t768;
+ t729 = t656*t673;
+ t752 = t673*t721+t653*t745+x*t729+t638*t717;
+ t703 = t664*t675;
+ t712 = t653*t673;
+ t710 = t665*t675;
+ t751 = t674*t703+t673*t714+(t712+t710)*t672;
+ t750 = y*t699+t651*t741+t769;
+ t749 = t636*t717+t673*t728+x*t712+t656*t745;
+ t748 = y*t708+t654*t741+t768;
+ t671 = RATIONAL(-5.0,6.0);
+ t711 = t671*t674;
+ t646 = RATIONAL(-1.0,24.0);
+ t715 = t646*t672;
+ t736 = t673*t715+t675*t711;
+ t734 = t677*t672;
+ t731 = t651*t674;
+ t727 = t640*t672;
+ t723 = RATIONAL(1.0,96.0)*t719+RATIONAL(5.0,96.0)*t732;
+ t658 = RATIONAL(5.0,24.0);
+ t718 = t658*t674;
+ t709 = t640*t674;
+ t668 = RATIONAL(-1.0,6.0);
+ t632 = t668*t672;
+ t641 = RATIONAL(1.0,24.0);
+ t629 = t641*t672;
+ t633 = t668*t673;
+ t627 = t641*t673;
+ t698 = t677*t635;
+ t697 = t736+t673*t718+t633+t770*t675;
+ t644 = RATIONAL(-1.0,96.0);
+ t643 = RATIONAL(-5.0,96.0);
+ t696 = t627+t646*t675+t643*t730+t644*t724+t723;
+ t695 = t646*t674+t644*t730+t643*t724+t629+t723;
+ t694 = t658*t724+t667*t674+t666*t730+t632+t736;
+ t642 = RATIONAL(1.0,48.0);
+ t647 = RATIONAL(-5.0,48.0);
+ t648 = RATIONAL(1.0,12.0);
+ t693 = t647*t675+t642*t673+t648;
+ t645 = RATIONAL(5.0,48.0);
+ t649 = RATIONAL(-1.0,48.0);
+ t650 = RATIONAL(-1.0,12.0);
+ t692 = t650+t649*t672+t645*t674;
+ t691 = t642*t672+t648+t647*t674;
+ t690 = t650+t649*t673+t645*t675;
+ t689 = t677*t731+t655*t734+t766+t771*y;
+ t688 = t767+t636*t734+t674*t733+(t700+t726)*y;
+ t687 = t676*t710+x*t703+t673*t706+t665*t737+t751;
+ t686 = t751+t675*t739+t676*t729+(t737+t717)*t662;
+ t685 = t677*t727+t674*t698+y*t709+t635*t746+t761;
+ t684 = t677*t707+t672*t733+t656*t744+y*t720+t767;
+ t683 = t677*t709+y*t727+t672*t698+t635*t744+t761;
+ t682 = t766+y*t731+t655*t746+t771*t677;
+ t670 = RATIONAL(-2.0,3.0);
+ t669 = RATIONAL(5.0,6.0);
+ t661 = RATIONAL(-5.0,4.0);
+ t660 = RATIONAL(1.0,4.0);
+ t659 = RATIONAL(-5.0,24.0);
+ t657 = RATIONAL(-5.0,16.0);
+ coeffs_I->coeff_m2_m2 = t683+t760+t762;
+ coeffs_I->coeff_m1_m2 = t682+t750+t752;
+ coeffs_I->coeff_0_m2 = t691*y+t692*t677+t696;
+ coeffs_I->coeff_p1_m2 = t682+t748+t749;
+ coeffs_I->coeff_p2_m2 = t683+t756+t757;
+ coeffs_I->coeff_m2_m1 = t688+t759+t763;
+ coeffs_I->coeff_m1_m1 = t686+t758+t765;
+ coeffs_I->coeff_0_m1 = (t669*t674+t632+t670)*y+(t659*t674+t666+t629)*t677
++t697;
+ coeffs_I->coeff_p1_m1 = t687+t764+t765;
+ coeffs_I->coeff_p2_m1 = t688+t753+t754;
+ coeffs_I->coeff_m2_0 = t690*t676+t693*x+t695;
+ coeffs_I->coeff_m1_0 = (t666+t659*t675+t627)*t676+(t669*t675+t670+t633)*x
++t694;
+ coeffs_I->coeff_0_0 = RATIONAL(1.0,1.0)+t660*t672+t661*t674+(RATIONAL(
+25.0,16.0)*t674+t661+t657*t672)*t675+(t660+t657*t674+RATIONAL(1.0,16.0)*t672)*
+t673;
+ coeffs_I->coeff_p1_0 = (t658*t675+t646*t673+t668)*t676+(t667+t671*t675+
+t666*t673)*x+t694;
+ coeffs_I->coeff_p2_0 = t693*t676+t690*x+t695;
+ coeffs_I->coeff_m2_p1 = t684+t753+t759;
+ coeffs_I->coeff_m1_p1 = t686+t755+t764;
+ coeffs_I->coeff_0_p1 = (t711+t770)*y+(t715+t718+t668)*t677+t697;
+ coeffs_I->coeff_p1_p1 = t687+t755+t758;
+ coeffs_I->coeff_p2_p1 = t684+t754+t763;
+ coeffs_I->coeff_m2_p2 = t685+t756+t762;
+ coeffs_I->coeff_m1_p2 = t689+t748+t752;
+ coeffs_I->coeff_0_p2 = t692*y+t691*t677+t696;
+ coeffs_I->coeff_p1_p2 = t689+t749+t750;
+ coeffs_I->coeff_p2_p2 = t685+t757+t760;
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-19 21:48:07 +0000