diff options
Diffstat (limited to 'src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c')
-rw-r--r-- | src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c | 205 |
1 files changed, 118 insertions, 87 deletions
diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c index ed83ce7..56cfd48 100644 --- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c +++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c @@ -1,87 +1,118 @@ -fp t311; -fp t314; -fp t318; -fp t301; -fp t331; -fp t322; -fp t330; -fp t315; -fp t302; -fp t329; -fp t328; -fp t316; -fp t327; -fp t319; -fp t326; -fp t323; -fp t324; -fp t325; -fp t317; -fp t310; -fp t309; -fp t307; -fp t306; -fp t300; -fp t299; -fp t298; -fp t297; - t311 = RATIONAL(-1.0,2.0)*x; - t314 = RATIONAL(1.0,2.0)*x; - t318 = RATIONAL(7.0,20.0); - t301 = y*t311; - t331 = t318*y+t301; - t322 = x*y; - t330 = t322+RATIONAL(1.0,4.0)*x; - t315 = RATIONAL(-11.0,20.0); - t302 = y*t314; - t329 = t302+t315*y; - t328 = -t322+RATIONAL(-1.0,4.0)*x; - t316 = RATIONAL(1.0,20.0); - t327 = t302+t316*y; - t319 = RATIONAL(3.0,20.0); - t326 = t319*y+t301; - t323 = x*x; - t324 = y*y; - t325 = t323*t324; - t317 = RATIONAL(-1.0,40.0); - t310 = RATIONAL(-1.0,6.0)*x; - t309 = RATIONAL(1.0,6.0)*x; - t307 = RATIONAL(-1.0,12.0)*x; - t306 = RATIONAL(1.0,12.0)*x; - t300 = RATIONAL(-9.0,50.0)*t325; - t299 = RATIONAL(9.0,50.0)*t325; - t298 = RATIONAL(-9.0,100.0)*t325; - t297 = RATIONAL(9.0,100.0)*t325; - coeffs_dxy->coeff_m1_m1 = t310+RATIONAL(11.0,60.0)+t297+(RATIONAL(59.0, -250.0)+RATIONAL(-17.0,50.0)*x)*t324+t329; - coeffs_dxy->coeff_0_m1 = RATIONAL(-1.0,20.0)+t309+(RATIONAL(21.0,500.0)+ -RATIONAL(13.0,25.0)*x+RATIONAL(-27.0,100.0)*t323)*t324+t326; - coeffs_dxy->coeff_p1_m1 = t309+RATIONAL(-7.0,60.0)+(RATIONAL(-1.0,50.0)*x -+RATIONAL(-73.0,250.0)+RATIONAL(27.0,100.0)*t323)*t324+t331; - coeffs_dxy->coeff_p2_m1 = t298+t310+RATIONAL(-1.0,60.0)+(RATIONAL(-4.0, -25.0)*x+RATIONAL(7.0,500.0))*t324+t327; - coeffs_dxy->coeff_m1_0 = RATIONAL(11.0,10.0)*y+RATIONAL(11.0,40.0)+( -RATIONAL(3.0,50.0)*t323+RATIONAL(-851.0,1000.0)+RATIONAL(69.0,100.0)*x)*t324+ -t328; - coeffs_dxy->coeff_0_0 = t300+RATIONAL(-3.0,10.0)*y+RATIONAL(-3.0,40.0)+( -RATIONAL(303.0,1000.0)+RATIONAL(-57.0,100.0)*x)*t324+t330; - coeffs_dxy->coeff_p1_0 = RATIONAL(-7.0,10.0)*y+RATIONAL(-7.0,40.0)+t299+( -RATIONAL(447.0,1000.0)+RATIONAL(-93.0,100.0)*x)*t324+t330; - coeffs_dxy->coeff_p2_0 = t317+RATIONAL(-1.0,10.0)*y+(RATIONAL(101.0, -1000.0)+RATIONAL(-3.0,50.0)*t323+RATIONAL(81.0,100.0)*x)*t324+t328; - coeffs_dxy->coeff_m1_p1 = t315+t314+(RATIONAL(203.0,250.0)+RATIONAL(-39.0 -,50.0)*x+RATIONAL(3.0,100.0)*t323)*t324+t329; - coeffs_dxy->coeff_0_p1 = t319+t311+t298+(RATIONAL(-93.0,500.0)+RATIONAL( -21.0,25.0)*x)*t324+t326; - coeffs_dxy->coeff_p1_p1 = t311+t297+t318+(RATIONAL(33.0,50.0)*x+RATIONAL( --141.0,250.0))*t324+t331; - coeffs_dxy->coeff_p2_p1 = t316+t314+(RATIONAL(-3.0,100.0)*t323+RATIONAL( --31.0,500.0)+RATIONAL(-18.0,25.0)*x)*t324+t327; - coeffs_dxy->coeff_m1_p2 = t307+RATIONAL(11.0,120.0)+t300+(RATIONAL(-197.0 -,1000.0)+RATIONAL(43.0,100.0)*x)*t324; - coeffs_dxy->coeff_0_p2 = t317+t306+(RATIONAL(-159.0,1000.0)+RATIONAL(27.0 -,50.0)*t323+RATIONAL(-79.0,100.0)*x)*t324; - coeffs_dxy->coeff_p1_p2 = t306+RATIONAL(-7.0,120.0)+(RATIONAL(-27.0,50.0) -*t323+RATIONAL(29.0,100.0)*x+RATIONAL(409.0,1000.0))*t324; - coeffs_dxy->coeff_p2_p2 = t307+RATIONAL(-1.0,120.0)+t299+(RATIONAL(-53.0, -1000.0)+RATIONAL(7.0,100.0)*x)*t324; +fp t435; +fp t424; +fp t439; +fp t409; +fp t460; +fp t438; +fp t416; +fp t459; +fp t451; +fp t429; +fp t445; +fp t436; +fp t403; +fp t458; +fp t423; +fp t408; +fp t457; +fp t456; +fp t437; +fp t455; +fp t454; +fp t427; +fp t453; +fp t405; +fp t452; +fp t421; +fp t422; +fp t390; +fp t450; +fp t449; +fp t448; +fp t434; +fp t447; +fp t446; +fp t425; +fp t444; +fp t443; +fp t426; +fp t391; +fp t442; +fp t441; +fp t420; +fp t440; +fp t433; +fp t432; +fp t431; +fp t430; +fp t428; +fp t394; +fp t389; + t435 = RATIONAL(-1.0,6.0); + t424 = RATIONAL(1.0,2.0); + t439 = y*y; + t409 = t424*t439; + t460 = t435+t409; + t438 = x*x; + t416 = t424*t438; + t459 = t416+t435; + t451 = t439+t438; + t429 = RATIONAL(-1.0,3.0); + t445 = t439*t438; + t436 = RATIONAL(3.0,4.0); + t403 = t436*t445; + t458 = t429+t403; + t423 = RATIONAL(-1.0,2.0); + t408 = t423*t438; + t457 = (t429+t408)*y; + t456 = t423+RATIONAL(-9.0,4.0)*t445; + t437 = x*y; + t455 = t438*y; + t454 = y+RATIONAL(9.0,4.0)*t445; + t427 = RATIONAL(3.0,2.0); + t453 = t427*t439; + t405 = t423*t439; + t452 = t427*t438; + t421 = RATIONAL(1.0,4.0); + t422 = RATIONAL(-3.0,4.0); + t390 = t422*t445; + t450 = t421*t439+t390; + t449 = t459*y; + t448 = t460*x; + t434 = RATIONAL(-1.0,4.0); + t447 = RATIONAL(1.0,18.0)+t434*t445; + t446 = t421*t438+t390; + t425 = RATIONAL(-3.0,2.0); + t444 = x+t437+t425*x*t439; + t443 = t437+(t429+t405)*x; + t426 = RATIONAL(-2.0,1.0); + t391 = t426*t437; + t442 = t391+(t423+t452)*y; + t441 = t391+(t453+t423)*x; + t420 = RATIONAL(1.0,12.0); + t440 = t420+t403+t451*t434; + t433 = RATIONAL(2.0,3.0); + t432 = RATIONAL(1.0,6.0); + t431 = RATIONAL(3.0,1.0); + t430 = RATIONAL(1.0,3.0); + t428 = RATIONAL(-3.0,1.0); + t394 = t421*t445; + t389 = t425*t455; + coeffs_dxy->coeff_m1_m1 = t394+RATIONAL(1.0,9.0)+t451*t432+t443+t457; + coeffs_dxy->coeff_0_m1 = t408+t432+(t433+t439)*x+t442+t450; + coeffs_dxy->coeff_p1_m1 = t405+t389+y+t416+t443+t458; + coeffs_dxy->coeff_p2_m1 = t420*t439+t435*t438+t447+t449; + coeffs_dxy->coeff_m1_0 = t405+t432+(t433+t438)*y+t441+t446; + coeffs_dxy->coeff_0_0 = t428*t455+t421+t451*t422+(t428*t439+1.0+RATIONAL( +4.0,1.0)*y)*x+t454; + coeffs_dxy->coeff_p1_0 = t436*t438+t453+(t431*t438+t426)*y+t441+t456; + coeffs_dxy->coeff_p2_0 = (-t438+t430)*y+t440; + coeffs_dxy->coeff_m1_p1 = t408+t409+t444+t457+t458; + coeffs_dxy->coeff_0_p1 = t452+t436*t439+(t426+t431*t439)*x+t442+t456; + coeffs_dxy->coeff_p1_p1 = t389+RATIONAL(1.0,1.0)+t451*t425+t444+t454; + coeffs_dxy->coeff_p2_p1 = t449+t450+t459; + coeffs_dxy->coeff_m1_p2 = t420*t438+t435*t439+t447+t448; + coeffs_dxy->coeff_0_p2 = (-t439+t430)*x+t440; + coeffs_dxy->coeff_p1_p2 = t446+t448+t460; + coeffs_dxy->coeff_p2_p2 = RATIONAL(1.0,36.0)+t394+t451*RATIONAL(-1.0,12.0 +); |