diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs')
21 files changed, 0 insertions, 2355 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order1.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order1.smooth0/coeffs-I.compute.c deleted file mode 100644 index 1680b1f..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order1.smooth0/coeffs-I.compute.c +++ /dev/null @@ -1,8 +0,0 @@ -fp t1; -fp t2; - t1 = x*y; - t2 = y-t1; - coeffs_I->coeff_0_0 = RATIONAL(1.0,1.0)-x-t2; - coeffs_I->coeff_p1_0 = x-t1; - coeffs_I->coeff_0_p1 = t2; - coeffs_I->coeff_p1_p1 = t1; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order1.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order1.smooth0/coeffs-dx.compute.c deleted file mode 100644 index 13d4254..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order1.smooth0/coeffs-dx.compute.c +++ /dev/null @@ -1,4 +0,0 @@ - coeffs_dx->coeff_0_0 = y+RATIONAL(-1.0,1.0); - coeffs_dx->coeff_p1_0 = -y+RATIONAL(1.0,1.0); - coeffs_dx->coeff_0_p1 = -y; - coeffs_dx->coeff_p1_p1 = y; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order1.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order1.smooth0/coeffs-dy.compute.c deleted file mode 100644 index ca7e73e..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order1.smooth0/coeffs-dy.compute.c +++ /dev/null @@ -1,4 +0,0 @@ - coeffs_dy->coeff_0_0 = x+RATIONAL(-1.0,1.0); - coeffs_dy->coeff_p1_0 = -x; - coeffs_dy->coeff_0_p1 = RATIONAL(1.0,1.0)-x; - coeffs_dy->coeff_p1_p1 = x; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c deleted file mode 100644 index a9bdd03..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c +++ /dev/null @@ -1,53 +0,0 @@ -fp t17; -fp t15; -fp t18; -fp t19; -fp t4; -fp t16; -fp t24; -fp t26; -fp t25; -fp t14; -fp t12; -fp t13; -fp t20; -fp t9; -fp t23; -fp t11; -fp t22; -fp t21; -fp t8; -fp t6; -fp t5; -fp t3; - t17 = x*x; - t15 = RATIONAL(1.0,4.0); - t18 = y*y; - t19 = t15*t18; - t4 = t17*t19; - t16 = RATIONAL(-1.0,4.0); - t24 = y*t16; - t26 = t4+t17*t24; - t25 = y*t15; - t14 = RATIONAL(1.0,2.0); - t12 = t14*t18; - t13 = RATIONAL(-1.0,2.0); - t20 = t13*t18; - t9 = t17*t20; - t23 = t12+t9; - t11 = t14*t17; - t22 = t11+t9; - t21 = t4+t17*t25; - t8 = x*t19; - t6 = x*t24; - t5 = x*t25; - t3 = t16*x*t18; - coeffs_I->coeff_m1_m1 = t3+t5+t26; - coeffs_I->coeff_0_m1 = (t13+t11)*y+t23; - coeffs_I->coeff_p1_m1 = t8+t6+t26; - coeffs_I->coeff_m1_0 = (t12+t13)*x+t22; - coeffs_I->coeff_0_0 = -t18+t17*t18+RATIONAL(1.0,1.0)-t17; - coeffs_I->coeff_p1_0 = (t20+t14)*x+t22; - coeffs_I->coeff_m1_p1 = t6+t3+t21; - coeffs_I->coeff_0_p1 = (t13*t17+t14)*y+t23; - coeffs_I->coeff_p1_p1 = t5+t8+t21; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dx.compute.c deleted file mode 100644 index bc709b8..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dx.compute.c +++ /dev/null @@ -1,35 +0,0 @@ -fp t35; -fp t41; -fp t38; -fp t29; -fp t36; -fp t40; -fp t37; -fp t39; -fp t34; -fp t33; -fp t30; -fp t28; -fp t27; - t35 = RATIONAL(1.0,2.0); - t41 = x*t35; - t38 = y*y; - t29 = t38*t41; - t36 = RATIONAL(1.0,4.0); - t40 = t29+t36*t38; - t37 = RATIONAL(-1.0,4.0); - t39 = t29+t37*t38; - t34 = RATIONAL(-1.0,2.0); - t33 = t36*y; - t30 = t37*y; - t28 = t34*x*y; - t27 = y*t41; - coeffs_dx->coeff_m1_m1 = t33+t28+t39; - coeffs_dx->coeff_0_m1 = (y-t38)*x; - coeffs_dx->coeff_p1_m1 = t28+t30+t40; - coeffs_dx->coeff_m1_0 = x+t34+(t35-x)*t38; - coeffs_dx->coeff_0_0 = (RATIONAL(-2.0,1.0)+RATIONAL(2.0,1.0)*t38)*x; - coeffs_dx->coeff_p1_0 = x+t35+(t34-x)*t38; - coeffs_dx->coeff_m1_p1 = t30+t27+t39; - coeffs_dx->coeff_0_p1 = (-y-t38)*x; - coeffs_dx->coeff_p1_p1 = t27+t33+t40; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c deleted file mode 100644 index d67f268..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c +++ /dev/null @@ -1,21 +0,0 @@ -fp t62; -fp t61; -fp t60; -fp t59; -fp t58; -fp t57; - t62 = y*y; - t61 = RATIONAL(1.0,2.0); - t60 = t61*t62; - t59 = -t62+RATIONAL(1.0,1.0); - t58 = RATIONAL(-1.0,2.0)*y+t60; - t57 = t60+t61*y; - coeffs_dxx->coeff_m1_m1 = t58; - coeffs_dxx->coeff_0_m1 = y-t62; - coeffs_dxx->coeff_p1_m1 = t58; - coeffs_dxx->coeff_m1_0 = t59; - coeffs_dxx->coeff_0_0 = RATIONAL(2.0,1.0)*t62+RATIONAL(-2.0,1.0); - coeffs_dxx->coeff_p1_0 = t59; - coeffs_dxx->coeff_m1_p1 = t57; - coeffs_dxx->coeff_0_p1 = -t62-y; - coeffs_dxx->coeff_p1_p1 = t57; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c deleted file mode 100644 index 0fda87b..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c +++ /dev/null @@ -1,29 +0,0 @@ -fp t72; -fp t68; -fp t74; -fp t69; -fp t73; -fp t71; -fp t70; -fp t65; -fp t64; -fp t63; - t72 = x*y; - t68 = RATIONAL(-1.0,2.0); - t74 = t68*x+t72; - t69 = RATIONAL(1.0,2.0); - t73 = t69*x+t72; - t71 = RATIONAL(-1.0,4.0); - t70 = RATIONAL(1.0,4.0); - t65 = t68*y; - t64 = t69*y; - t63 = RATIONAL(-2.0,1.0)*t72; - coeffs_dxy->coeff_m1_m1 = t65+t70+t74; - coeffs_dxy->coeff_0_m1 = t63+x; - coeffs_dxy->coeff_p1_m1 = t71+t64+t74; - coeffs_dxy->coeff_m1_0 = t63+y; - coeffs_dxy->coeff_0_0 = RATIONAL(4.0,1.0)*t72; - coeffs_dxy->coeff_p1_0 = -y+t63; - coeffs_dxy->coeff_m1_p1 = t65+t71+t73; - coeffs_dxy->coeff_0_p1 = t63-x; - coeffs_dxy->coeff_p1_p1 = t70+t64+t73; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c deleted file mode 100644 index e373f6b..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c +++ /dev/null @@ -1,35 +0,0 @@ -fp t50; -fp t56; -fp t53; -fp t43; -fp t51; -fp t55; -fp t52; -fp t54; -fp t49; -fp t47; -fp t45; -fp t44; -fp t42; - t50 = RATIONAL(1.0,2.0); - t56 = y*t50; - t53 = x*x; - t43 = t53*t56; - t51 = RATIONAL(1.0,4.0); - t55 = t43+t51*x; - t52 = RATIONAL(-1.0,4.0); - t54 = t43+t52*x; - t49 = RATIONAL(-1.0,2.0); - t47 = t52*t53; - t45 = t51*t53; - t44 = x*t56; - t42 = t49*x*y; - coeffs_dy->coeff_m1_m1 = t42+t47+t55; - coeffs_dy->coeff_0_m1 = t49+y+(-y+t50)*t53; - coeffs_dy->coeff_p1_m1 = t44+t47+t54; - coeffs_dy->coeff_m1_0 = (x-t53)*y; - coeffs_dy->coeff_0_0 = (RATIONAL(-2.0,1.0)+RATIONAL(2.0,1.0)*t53)*y; - coeffs_dy->coeff_p1_0 = (-x-t53)*y; - coeffs_dy->coeff_m1_p1 = t45+t42+t54; - coeffs_dy->coeff_0_p1 = t50+y+(t49-y)*t53; - coeffs_dy->coeff_p1_p1 = t44+t45+t55; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c deleted file mode 100644 index 9a941cd..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c +++ /dev/null @@ -1,21 +0,0 @@ -fp t80; -fp t79; -fp t78; -fp t77; -fp t76; -fp t75; - t80 = x*x; - t79 = RATIONAL(1.0,2.0); - t78 = t79*t80; - t77 = RATIONAL(1.0,1.0)-t80; - t76 = t78+RATIONAL(-1.0,2.0)*x; - t75 = t79*x+t78; - coeffs_dyy->coeff_m1_m1 = t76; - coeffs_dyy->coeff_0_m1 = t77; - coeffs_dyy->coeff_p1_m1 = t75; - coeffs_dyy->coeff_m1_0 = -t80+x; - coeffs_dyy->coeff_0_0 = RATIONAL(2.0,1.0)*t80+RATIONAL(-2.0,1.0); - coeffs_dyy->coeff_p1_0 = -x-t80; - coeffs_dyy->coeff_m1_p1 = t76; - coeffs_dyy->coeff_0_p1 = t77; - coeffs_dyy->coeff_p1_p1 = t75; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c deleted file mode 100644 index cb6836d..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c +++ /dev/null @@ -1,162 +0,0 @@ -fp t128; -fp t116; -fp t105; -fp t120; -fp t152; -fp t127; -fp t151; -fp t150; -fp t118; -fp t149; -fp t148; -fp t121; -fp t147; -fp t146; -fp t137; -fp t138; -fp t145; -fp t139; -fp t141; -fp t144; -fp t129; -fp t130; -fp t142; -fp t143; -fp t115; -fp t104; -fp t96; -fp t119; -fp t136; -fp t110; -fp t135; -fp t101; -fp t114; -fp t134; -fp t133; -fp t113; -fp t112; -fp t111; -fp t109; -fp t108; -fp t100; -fp t98; -fp t91; -fp t88; -fp t87; -fp t86; -fp t85; -fp t84; -fp t83; -fp t82; - t128 = y*y; - t116 = RATIONAL(-3.0,40.0); - t105 = t116*y; - t120 = RATIONAL(-3.0,20.0); - t152 = t105+t120*t128; - t127 = x*x; - t151 = y*t127; - t150 = t128*x; - t118 = RATIONAL(-1.0,20.0); - t149 = t118*y; - t148 = RATIONAL(1.0,20.0)*y; - t121 = RATIONAL(3.0,20.0); - t147 = t121*y; - t146 = RATIONAL(-1.0,40.0)*y; - t137 = RATIONAL(1.0,8.0)*t127; - t138 = RATIONAL(1.0,4.0)*t127; - t145 = y*t137+t128*t138; - t139 = RATIONAL(-1.0,8.0)*t127; - t141 = RATIONAL(-1.0,4.0)*t127; - t144 = y*t139+t128*t141; - t129 = t127*x; - t130 = t128*y; - t142 = t129*t130; - t143 = RATIONAL(1.0,50.0)*t142+RATIONAL(-1.0,400.0)*t127; - t115 = RATIONAL(3.0,40.0); - t104 = t115*t128; - t96 = t128*t137; - t119 = RATIONAL(1.0,40.0); - t136 = t116*t128+t119*t150+t96; - t110 = RATIONAL(-11.0,40.0); - t135 = t110*t150+t104+t96; - t101 = t128*t139; - t114 = RATIONAL(11.0,40.0); - t134 = t114*t128+t101+x*t104; - t133 = RATIONAL(7.0,40.0)*t150+t101+RATIONAL(9.0,40.0)*t128; - t113 = RATIONAL(-1.0,300.0); - t112 = RATIONAL(11.0,20.0); - t111 = RATIONAL(-11.0,20.0); - t109 = RATIONAL(29.0,600.0); - t108 = t120*y; - t100 = y*t138; - t98 = y*t141; - t91 = x*t146; - t88 = RATIONAL(1.0,12.0)*t151; - t87 = RATIONAL(-1.0,12.0)*t151; - t86 = RATIONAL(-1.0,50.0)*t142; - t85 = RATIONAL(-1.0,24.0)*t151; - t84 = RATIONAL(1.0,24.0)*t151; - t83 = RATIONAL(1.0,100.0)*t142; - t82 = RATIONAL(-1.0,100.0)*t142; - coeffs_I->coeff_m1_m1 = RATIONAL(-111.0,2000.0)+RATIONAL(37.0,400.0)*t127 -+t87+RATIONAL(-37.0,600.0)*t129+t149+t83+(RATIONAL(-17.0,300.0)*t127+RATIONAL( --2.0,125.0))*t130+(RATIONAL(481.0,6000.0)+RATIONAL(11.0,60.0)*y+RATIONAL(59.0, -750.0)*t130)*x+t135; - coeffs_I->coeff_0_m1 = t88+RATIONAL(-111.0,400.0)*t127+RATIONAL(333.0, -2000.0)+RATIONAL(-11.0,60.0)*y+RATIONAL(37.0,200.0)*t129+(RATIONAL(-481.0, -2000.0)+t149)*x+(RATIONAL(7.0,500.0)*x+RATIONAL(-3.0,100.0)*t129+RATIONAL(13.0, -150.0)*t127+RATIONAL(-89.0,750.0))*t130+t134; - coeffs_I->coeff_p1_m1 = t88+RATIONAL(-333.0,2000.0)+RATIONAL(-37.0,200.0) -*t129+RATIONAL(111.0,400.0)*t127+t108+(RATIONAL(-7.0,60.0)*y+RATIONAL(481.0, -2000.0))*x+(RATIONAL(3.0,100.0)*t129+t113*t127+RATIONAL(-73.0,750.0)*x+RATIONAL -(-6.0,125.0))*t130+t133; - coeffs_I->coeff_p2_m1 = t82+t87+t148+RATIONAL(-37.0,400.0)*t127+RATIONAL( -111.0,2000.0)+RATIONAL(37.0,600.0)*t129+(RATIONAL(2.0,125.0)+RATIONAL(-2.0,75.0 -)*t127)*t130+(RATIONAL(7.0,1500.0)*t130+RATIONAL(-481.0,6000.0)+RATIONAL(-1.0, -60.0)*y)*x+t136; - coeffs_I->coeff_m1_0 = RATIONAL(201.0,2000.0)+RATIONAL(133.0,400.0)*t127+ -RATIONAL(-11.0,200.0)*t129+(RATIONAL(23.0,200.0)*t127+RATIONAL(1.0,150.0)*t129+ -RATIONAL(81.0,1000.0))*t130+(t114*y+RATIONAL(-957.0,2000.0)+t112*t128+RATIONAL( --851.0,3000.0)*t130)*x+t144+t152; - coeffs_I->coeff_0_0 = t86+RATIONAL(-199.0,400.0)*t127+t111*t128+t110*y+ -RATIONAL(33.0,200.0)*t129+RATIONAL(1397.0,2000.0)+(RATIONAL(-19.0,200.0)*t127+ -RATIONAL(257.0,1000.0))*t130+(RATIONAL(-129.0,2000.0)+RATIONAL(101.0,1000.0)* -t130+t152)*x+t145; - coeffs_I->coeff_p1_0 = RATIONAL(-9.0,40.0)*y+RATIONAL(-33.0,200.0)*t129+ -RATIONAL(603.0,2000.0)+RATIONAL(-9.0,20.0)*t128+(RATIONAL(243.0,1000.0)+ -RATIONAL(-31.0,200.0)*t127)*t130+(RATIONAL(149.0,1000.0)*t130+RATIONAL(-7.0, -40.0)*y+RATIONAL(-7.0,20.0)*t128+RATIONAL(1129.0,2000.0))*x+t143+t145; - coeffs_I->coeff_p2_0 = t115*y+RATIONAL(11.0,200.0)*t129+RATIONAL(67.0, -400.0)*t127+RATIONAL(-43.0,2000.0)*x+t91+RATIONAL(-201.0,2000.0)+(t118*x+t121)* -t128+(RATIONAL(-1.0,150.0)*t129+RATIONAL(101.0,3000.0)*x+RATIONAL(-81.0,1000.0) -+RATIONAL(27.0,200.0)*t127)*t130+t144; - coeffs_I->coeff_m1_p1 = RATIONAL(-29.0,600.0)*t129+RATIONAL(29.0,400.0)* -t127+RATIONAL(-87.0,2000.0)+t147+t100+(t111*y+RATIONAL(377.0,6000.0))*x+( -RATIONAL(-13.0,100.0)*t127+RATIONAL(1.0,300.0)*t129+RATIONAL(-9.0,125.0)+ -RATIONAL(203.0,750.0)*x)*t130+t135; - coeffs_I->coeff_0_p1 = RATIONAL(261.0,2000.0)+t98+RATIONAL(-87.0,400.0)* -t127+t112*y+t82+RATIONAL(29.0,200.0)*t129+(RATIONAL(-71.0,250.0)+RATIONAL(7.0, -50.0)*t127)*t130+(t147+RATIONAL(-31.0,500.0)*t130+RATIONAL(-377.0,2000.0))*x+ -t134; - coeffs_I->coeff_p1_p1 = RATIONAL(87.0,400.0)*t127+RATIONAL(9.0,20.0)*y+ -RATIONAL(-29.0,200.0)*t129+t98+t83+RATIONAL(-261.0,2000.0)+(RATIONAL(-27.0, -125.0)+RATIONAL(11.0,100.0)*t127)*t130+(RATIONAL(-47.0,250.0)*t130+RATIONAL(7.0 -,20.0)*y+RATIONAL(377.0,2000.0))*x+t133; - coeffs_I->coeff_p2_p1 = RATIONAL(87.0,2000.0)+RATIONAL(-29.0,400.0)*t127+ -t100+t109*t129+t108+(RATIONAL(-377.0,6000.0)+t148)*x+(RATIONAL(9.0,125.0)+ -RATIONAL(-31.0,1500.0)*x+t113*t129+RATIONAL(-3.0,25.0)*t127)*t130+t136; - coeffs_I->coeff_m1_p2 = RATIONAL(1.0,400.0)*t127+t146+t86+RATIONAL(-1.0, -600.0)*t129+t85+RATIONAL(-3.0,2000.0)+(RATIONAL(7.0,1000.0)+RATIONAL(43.0,600.0 -)*t127)*t130+(RATIONAL(11.0,120.0)*y+RATIONAL(-197.0,3000.0)*t130+RATIONAL(13.0 -,6000.0))*x; - coeffs_I->coeff_0_p2 = RATIONAL(-3.0,400.0)*t127+RATIONAL(-11.0,120.0)*y+ -t91+t84+RATIONAL(-13.0,2000.0)*x+RATIONAL(9.0,2000.0)+RATIONAL(1.0,200.0)*t129+ -(RATIONAL(-53.0,1000.0)*x+RATIONAL(3.0,50.0)*t129+RATIONAL(437.0,3000.0)+ -RATIONAL(-79.0,600.0)*t127)*t130; - coeffs_I->coeff_p1_p2 = RATIONAL(3.0,400.0)*t127+t105+t84+RATIONAL(-1.0, -200.0)*t129+RATIONAL(-9.0,2000.0)+(RATIONAL(13.0,2000.0)+RATIONAL(-7.0,120.0)*y -)*x+(RATIONAL(409.0,3000.0)*x+t109*t127+RATIONAL(21.0,1000.0)+RATIONAL(-3.0, -50.0)*t129)*t130; - coeffs_I->coeff_p2_p2 = RATIONAL(3.0,2000.0)+t119*y+t85+RATIONAL(1.0, -600.0)*t129+(RATIONAL(-7.0,1000.0)+RATIONAL(7.0,600.0)*t127)*t130+(RATIONAL( --53.0,3000.0)*t130+RATIONAL(-1.0,120.0)*y+RATIONAL(-13.0,6000.0))*x+t143; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c deleted file mode 100644 index a593111..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c +++ /dev/null @@ -1,131 +0,0 @@ -fp t207; -fp t193; -fp t206; -fp t205; -fp t204; -fp t188; -fp t168; -fp t203; -fp t167; -fp t202; -fp t201; -fp t191; -fp t200; -fp t199; -fp t192; -fp t194; -fp t198; -fp t197; -fp t183; -fp t196; -fp t187; -fp t186; -fp t185; -fp t184; -fp t182; -fp t181; -fp t180; -fp t179; -fp t178; -fp t177; -fp t176; -fp t175; -fp t174; -fp t173; -fp t172; -fp t171; -fp t169; -fp t158; -fp t156; -fp t155; -fp t154; -fp t153; - t207 = x*y; - t193 = y*y; - t206 = x*t193; - t205 = x*RATIONAL(-1.0,4.0); - t204 = x*RATIONAL(1.0,2.0); - t188 = RATIONAL(1.0,4.0); - t168 = t188*t206; - t203 = RATIONAL(-1.0,6.0)*t207+t168; - t167 = t193*t205; - t202 = RATIONAL(1.0,6.0)*t207+t167; - t201 = t188*t207+t193*t204; - t191 = RATIONAL(-1.0,2.0); - t200 = y*t205+t191*t206; - t199 = y*t204+t168; - t192 = x*x; - t194 = y*t193; - t198 = t192*t194; - t197 = t191*t207+t167; - t183 = RATIONAL(-1.0,200.0); - t196 = RATIONAL(3.0,50.0)*t198+t183*x; - t187 = RATIONAL(-1.0,20.0); - t186 = RATIONAL(3.0,200.0); - t185 = RATIONAL(-3.0,200.0); - t184 = RATIONAL(1.0,200.0); - t182 = RATIONAL(37.0,200.0); - t181 = RATIONAL(87.0,200.0); - t180 = RATIONAL(-37.0,200.0); - t179 = RATIONAL(-87.0,200.0); - t178 = RATIONAL(29.0,200.0); - t177 = RATIONAL(-29.0,200.0); - t176 = RATIONAL(111.0,200.0); - t175 = RATIONAL(-111.0,200.0); - t174 = RATIONAL(3.0,40.0)*t193; - t173 = RATIONAL(7.0,40.0)*t193; - t172 = RATIONAL(-1.0,40.0)*y; - t171 = RATIONAL(1.0,40.0)*t193; - t169 = RATIONAL(-11.0,40.0)*t193; - t158 = RATIONAL(-3.0,50.0)*t198; - t156 = RATIONAL(-1.0,12.0)*t207; - t155 = RATIONAL(1.0,12.0)*t207; - t154 = RATIONAL(3.0,100.0)*t198; - t153 = RATIONAL(-3.0,100.0)*t198; - coeffs_dx->coeff_m1_m1 = t180*t192+RATIONAL(481.0,6000.0)+t154+t169+ -RATIONAL(11.0,60.0)*y+RATIONAL(59.0,750.0)*t194+(RATIONAL(-17.0,150.0)*t194+ -t182)*x+t203; - coeffs_dx->coeff_0_m1 = t176*t192+t174+RATIONAL(-481.0,2000.0)+t175*x+ -t187*y+(RATIONAL(13.0,75.0)*x+RATIONAL(7.0,500.0)+RATIONAL(-9.0,100.0)*t192)* -t194+t202; - coeffs_dx->coeff_p1_m1 = t173+t176*x+RATIONAL(-7.0,60.0)*y+t175*t192+ -RATIONAL(481.0,2000.0)+(RATIONAL(-1.0,150.0)*x+RATIONAL(-73.0,750.0)+RATIONAL( -9.0,100.0)*t192)*t194+t202; - coeffs_dx->coeff_p2_m1 = RATIONAL(7.0,1500.0)*t194+t153+RATIONAL(-1.0, -60.0)*y+t171+t182*t192+RATIONAL(-481.0,6000.0)+(RATIONAL(-4.0,75.0)*t194+t180)* -x+t203; - coeffs_dx->coeff_m1_0 = RATIONAL(-33.0,200.0)*t192+RATIONAL(11.0,40.0)*y+ -RATIONAL(-957.0,2000.0)+RATIONAL(133.0,200.0)*x+RATIONAL(11.0,20.0)*t193+( -RATIONAL(-851.0,3000.0)+RATIONAL(23.0,100.0)*x+RATIONAL(1.0,50.0)*t192)*t194+ -t200; - coeffs_dx->coeff_0_0 = RATIONAL(-3.0,40.0)*y+RATIONAL(-129.0,2000.0)+t158 -+RATIONAL(101.0,1000.0)*t194+RATIONAL(-3.0,20.0)*t193+RATIONAL(99.0,200.0)*t192 -+(RATIONAL(-19.0,100.0)*t194+RATIONAL(-199.0,200.0))*x+t201; - coeffs_dx->coeff_p1_0 = RATIONAL(-7.0,40.0)*y+RATIONAL(1129.0,2000.0)+ -RATIONAL(-7.0,20.0)*t193+RATIONAL(-99.0,200.0)*t192+(RATIONAL(149.0,1000.0)+ -RATIONAL(-31.0,100.0)*x)*t194+t196+t201; - coeffs_dx->coeff_p2_0 = t172+t187*t193+RATIONAL(33.0,200.0)*t192+RATIONAL -(-43.0,2000.0)+RATIONAL(67.0,200.0)*x+(RATIONAL(-1.0,50.0)*t192+RATIONAL(27.0, -100.0)*x+RATIONAL(101.0,3000.0))*t194+t200; - coeffs_dx->coeff_m1_p1 = RATIONAL(-11.0,20.0)*y+t169+RATIONAL(377.0, -6000.0)+t178*x+t177*t192+(RATIONAL(-13.0,50.0)*x+RATIONAL(1.0,100.0)*t192+ -RATIONAL(203.0,750.0))*t194+t199; - coeffs_dx->coeff_0_p1 = t153+t181*t192+t174+RATIONAL(-31.0,500.0)*t194+ -RATIONAL(3.0,20.0)*y+RATIONAL(-377.0,2000.0)+(t179+RATIONAL(7.0,25.0)*t194)*x+ -t197; - coeffs_dx->coeff_p1_p1 = RATIONAL(7.0,20.0)*y+RATIONAL(-47.0,250.0)*t194+ -t179*t192+RATIONAL(377.0,2000.0)+t173+t154+(RATIONAL(11.0,50.0)*t194+t181)*x+ -t197; - coeffs_dx->coeff_p2_p1 = RATIONAL(1.0,20.0)*y+RATIONAL(-377.0,6000.0)+ -t178*t192+t171+t177*x+(RATIONAL(-31.0,1500.0)+RATIONAL(-6.0,25.0)*x+RATIONAL( --1.0,100.0)*t192)*t194+t199; - coeffs_dx->coeff_m1_p2 = RATIONAL(11.0,120.0)*y+t183*t192+RATIONAL(13.0, -6000.0)+RATIONAL(-197.0,3000.0)*t194+t156+t158+(RATIONAL(43.0,300.0)*t194+t184) -*x; - coeffs_dx->coeff_0_p2 = t186*t192+t172+t155+RATIONAL(-13.0,2000.0)+t185*x -+(RATIONAL(-53.0,1000.0)+RATIONAL(-79.0,300.0)*x+RATIONAL(9.0,50.0)*t192)*t194; - coeffs_dx->coeff_p1_p2 = t155+t185*t192+RATIONAL(-7.0,120.0)*y+t186*x+ -RATIONAL(13.0,2000.0)+(RATIONAL(409.0,3000.0)+RATIONAL(-9.0,50.0)*t192+RATIONAL -(29.0,300.0)*x)*t194; - coeffs_dx->coeff_p2_p2 = RATIONAL(-1.0,120.0)*y+t184*t192+RATIONAL(-13.0, -6000.0)+t156+(RATIONAL(-53.0,3000.0)+RATIONAL(7.0,300.0)*x)*t194+t196; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c deleted file mode 100644 index 8c53891..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c +++ /dev/null @@ -1,74 +0,0 @@ -fp t286; -fp t287; -fp t296; -fp t282; -fp t276; -fp t284; -fp t295; -fp t285; -fp t283; -fp t294; -fp t274; -fp t293; -fp t292; -fp t291; -fp t290; -fp t289; -fp t270; -fp t269; -fp t268; -fp t266; -fp t265; - t286 = y*y; - t287 = t286*y; - t296 = x*t287; - t282 = RATIONAL(1.0,4.0); - t276 = t282*t286; - t284 = RATIONAL(1.0,2.0); - t295 = t284*y+t276; - t285 = RATIONAL(-1.0,2.0); - t283 = RATIONAL(-1.0,4.0); - t294 = t283*y+t285*t286; - t274 = t283*t286; - t293 = t285*y+t274; - t292 = RATIONAL(3.0,25.0)*t296+RATIONAL(-1.0,200.0); - t291 = t284*t286+t282*y; - t290 = RATIONAL(-1.0,6.0)*y+t276; - t289 = RATIONAL(1.0,6.0)*y+t274; - t270 = RATIONAL(1.0,12.0)*y; - t269 = RATIONAL(-1.0,12.0)*y; - t268 = RATIONAL(-3.0,25.0)*t296; - t266 = RATIONAL(3.0,50.0)*t296; - t265 = RATIONAL(-3.0,50.0)*t296; - coeffs_dxx->coeff_m1_m1 = RATIONAL(-37.0,100.0)*x+RATIONAL(-17.0,150.0)* -t287+t266+RATIONAL(37.0,200.0)+t290; - coeffs_dxx->coeff_0_m1 = RATIONAL(-111.0,200.0)+RATIONAL(13.0,75.0)*t287+ -(RATIONAL(-9.0,50.0)*t287+RATIONAL(111.0,100.0))*x+t289; - coeffs_dxx->coeff_p1_m1 = RATIONAL(-1.0,150.0)*t287+RATIONAL(111.0,200.0) -+(RATIONAL(-111.0,100.0)+RATIONAL(9.0,50.0)*t287)*x+t289; - coeffs_dxx->coeff_p2_m1 = RATIONAL(-37.0,200.0)+RATIONAL(-4.0,75.0)*t287+ -t265+RATIONAL(37.0,100.0)*x+t290; - coeffs_dxx->coeff_m1_0 = RATIONAL(23.0,100.0)*t287+RATIONAL(133.0,200.0)+ -(RATIONAL(-33.0,100.0)+RATIONAL(1.0,25.0)*t287)*x+t294; - coeffs_dxx->coeff_0_0 = RATIONAL(-199.0,200.0)+RATIONAL(-19.0,100.0)*t287 -+RATIONAL(99.0,100.0)*x+t268+t291; - coeffs_dxx->coeff_p1_0 = RATIONAL(-99.0,100.0)*x+RATIONAL(-31.0,100.0)* -t287+t291+t292; - coeffs_dxx->coeff_p2_0 = RATIONAL(67.0,200.0)+RATIONAL(27.0,100.0)*t287+( -RATIONAL(-1.0,25.0)*t287+RATIONAL(33.0,100.0))*x+t294; - coeffs_dxx->coeff_m1_p1 = RATIONAL(-13.0,50.0)*t287+RATIONAL(29.0,200.0)+ -(RATIONAL(1.0,50.0)*t287+RATIONAL(-29.0,100.0))*x+t295; - coeffs_dxx->coeff_0_p1 = RATIONAL(-87.0,200.0)+t265+RATIONAL(87.0,100.0)* -x+RATIONAL(7.0,25.0)*t287+t293; - coeffs_dxx->coeff_p1_p1 = RATIONAL(-87.0,100.0)*x+t266+RATIONAL(87.0, -200.0)+RATIONAL(11.0,50.0)*t287+t293; - coeffs_dxx->coeff_p2_p1 = RATIONAL(-6.0,25.0)*t287+RATIONAL(-29.0,200.0)+ -(RATIONAL(-1.0,50.0)*t287+RATIONAL(29.0,100.0))*x+t295; - coeffs_dxx->coeff_m1_p2 = RATIONAL(43.0,300.0)*t287+RATIONAL(-1.0,100.0)* -x+t268+t269+RATIONAL(1.0,200.0); - coeffs_dxx->coeff_0_p2 = RATIONAL(-79.0,300.0)*t287+RATIONAL(-3.0,200.0)+ -t270+(RATIONAL(3.0,100.0)+RATIONAL(9.0,25.0)*t287)*x; - coeffs_dxx->coeff_p1_p2 = t270+RATIONAL(29.0,300.0)*t287+RATIONAL(3.0, -200.0)+(RATIONAL(-9.0,25.0)*t287+RATIONAL(-3.0,100.0))*x; - coeffs_dxx->coeff_p2_p2 = RATIONAL(1.0,100.0)*x+RATIONAL(7.0,300.0)*t287+ -t269+t292; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c deleted file mode 100644 index ed83ce7..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c +++ /dev/null @@ -1,87 +0,0 @@ -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; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c deleted file mode 100644 index 6c1ce7f..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c +++ /dev/null @@ -1,114 +0,0 @@ -fp t264; -fp t240; -fp t242; -fp t263; -fp t247; -fp t262; -fp t261; -fp t243; -fp t226; -fp t249; -fp t248; -fp t258; -fp t232; -fp t229; -fp t257; -fp t256; -fp t234; -fp t218; -fp t255; -fp t219; -fp t254; -fp t253; -fp t252; -fp t239; -fp t251; -fp t244; -fp t241; -fp t235; -fp t228; -fp t227; -fp t225; -fp t224; -fp t223; -fp t215; -fp t211; -fp t210; -fp t208; - t264 = RATIONAL(7.0,20.0)*x+RATIONAL(9.0,20.0); - t240 = RATIONAL(1.0,20.0); - t242 = RATIONAL(-3.0,20.0); - t263 = t240*x+t242; - t247 = x*x; - t262 = y*t247; - t261 = RATIONAL(-11.0,20.0)*x; - t243 = RATIONAL(3.0,20.0); - t226 = t243*y; - t249 = t247*x; - t248 = y*y; - t258 = t249*t248; - t232 = RATIONAL(1.0,4.0)*t247; - t229 = RATIONAL(-1.0,4.0)*t247; - t257 = RATIONAL(-1.0,8.0)*t247+RATIONAL(-1.0,2.0)*t262; - t256 = RATIONAL(1.0,8.0)*t247+RATIONAL(1.0,2.0)*t262; - t234 = RATIONAL(11.0,20.0); - t218 = y*t229; - t255 = x*t226+t234*y+t218; - t219 = y*t232; - t254 = t226+t219+y*t261; - t253 = t219+t263*y; - t252 = t218+t264*y; - t239 = RATIONAL(-3.0,40.0); - t251 = RATIONAL(-3.0,10.0)*y+t239; - t244 = RATIONAL(-1.0,40.0); - t241 = RATIONAL(-1.0,20.0); - t235 = RATIONAL(-1.0,100.0); - t228 = RATIONAL(-1.0,24.0)*t247; - t227 = t244*x; - t225 = RATIONAL(1.0,24.0)*t247; - t224 = RATIONAL(-1.0,12.0)*t247; - t223 = RATIONAL(1.0,12.0)*t247; - t215 = RATIONAL(3.0,50.0)*t258; - t211 = RATIONAL(-3.0,50.0)*t258; - t210 = RATIONAL(3.0,100.0)*t258; - t208 = RATIONAL(-3.0,100.0)*t258; - coeffs_dy->coeff_m1_m1 = t241+t210+RATIONAL(11.0,60.0)*x+t224+(RATIONAL( --17.0,100.0)*t247+RATIONAL(-6.0,125.0)+RATIONAL(59.0,250.0)*x)*t248+t254; - coeffs_dy->coeff_0_m1 = t223+t241*x+RATIONAL(-11.0,60.0)+(RATIONAL(21.0, -500.0)*x+RATIONAL(13.0,50.0)*t247+RATIONAL(-9.0,100.0)*t249+RATIONAL(-89.0, -250.0))*t248+t255; - coeffs_dy->coeff_p1_m1 = t242+t223+RATIONAL(-7.0,60.0)*x+(RATIONAL(-18.0, -125.0)+t235*t247+RATIONAL(9.0,100.0)*t249+RATIONAL(-73.0,250.0)*x)*t248+t252; - coeffs_dy->coeff_p2_m1 = RATIONAL(-1.0,60.0)*x+t240+t208+t224+(RATIONAL( --2.0,25.0)*t247+RATIONAL(6.0,125.0)+RATIONAL(7.0,500.0)*x)*t248+t253; - coeffs_dy->coeff_m1_0 = (RATIONAL(11.0,40.0)+RATIONAL(11.0,10.0)*y)*x+( -RATIONAL(-851.0,1000.0)*x+RATIONAL(69.0,200.0)*t247+RATIONAL(1.0,50.0)*t249+ -RATIONAL(243.0,1000.0))*t248+t251+t257; - coeffs_dy->coeff_0_0 = RATIONAL(-11.0,10.0)*y+t211+RATIONAL(-11.0,40.0)+( -RATIONAL(771.0,1000.0)+RATIONAL(-57.0,200.0)*t247)*t248+(RATIONAL(303.0,1000.0) -*t248+t251)*x+t256; - coeffs_dy->coeff_p1_0 = t215+RATIONAL(-9.0,10.0)*y+RATIONAL(-9.0,40.0)+( -RATIONAL(-93.0,200.0)*t247+RATIONAL(729.0,1000.0))*t248+(RATIONAL(-7.0,10.0)*y+ -RATIONAL(-7.0,40.0)+RATIONAL(447.0,1000.0)*t248)*x+t256; - coeffs_dy->coeff_p2_0 = RATIONAL(3.0,40.0)+t227+(RATIONAL(-1.0,10.0)*x+ -RATIONAL(3.0,10.0))*y+(RATIONAL(81.0,200.0)*t247+RATIONAL(101.0,1000.0)*x+ -RATIONAL(-243.0,1000.0)+RATIONAL(-1.0,50.0)*t249)*t248+t257; - coeffs_dy->coeff_m1_p1 = t243+t232+t261+(RATIONAL(-27.0,125.0)+RATIONAL( --39.0,100.0)*t247+RATIONAL(203.0,250.0)*x+RATIONAL(1.0,100.0)*t249)*t248+t254; - coeffs_dy->coeff_0_p1 = t229+t234+t208+t243*x+(RATIONAL(-213.0,250.0)+ -RATIONAL(21.0,50.0)*t247+RATIONAL(-93.0,500.0)*x)*t248+t255; - coeffs_dy->coeff_p1_p1 = t229+t210+(RATIONAL(33.0,100.0)*t247+RATIONAL( --141.0,250.0)*x+RATIONAL(-81.0,125.0))*t248+t252+t264; - coeffs_dy->coeff_p2_p1 = t232+(RATIONAL(-9.0,25.0)*t247+RATIONAL(-31.0, -500.0)*x+t235*t249+RATIONAL(27.0,125.0))*t248+t253+t263; - coeffs_dy->coeff_m1_p2 = t244+t228+RATIONAL(11.0,120.0)*x+t211+(RATIONAL( --197.0,1000.0)*x+RATIONAL(43.0,200.0)*t247+RATIONAL(21.0,1000.0))*t248; - coeffs_dy->coeff_0_p2 = t225+t227+RATIONAL(-11.0,120.0)+(RATIONAL(-159.0, -1000.0)*x+RATIONAL(9.0,50.0)*t249+RATIONAL(-79.0,200.0)*t247+RATIONAL(437.0, -1000.0))*t248; - coeffs_dy->coeff_p1_p2 = t239+RATIONAL(-7.0,120.0)*x+t225+(RATIONAL(29.0, -200.0)*t247+RATIONAL(-9.0,50.0)*t249+RATIONAL(63.0,1000.0)+RATIONAL(409.0, -1000.0)*x)*t248; - coeffs_dy->coeff_p2_p2 = RATIONAL(1.0,40.0)+t215+t228+RATIONAL(-1.0,120.0 -)*x+(RATIONAL(7.0,200.0)*t247+RATIONAL(-21.0,1000.0)+RATIONAL(-53.0,1000.0)*x)* -t248; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c deleted file mode 100644 index 1b690c6..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c +++ /dev/null @@ -1,71 +0,0 @@ -fp t350; -fp t351; -fp t357; -fp t342; -fp t356; -fp t343; -fp t355; -fp t348; -fp t354; -fp t353; -fp t349; -fp t345; -fp t341; -fp t340; -fp t335; -fp t334; -fp t333; -fp t332; - t350 = x*x; - t351 = t350*x; - t357 = y*t351; - t342 = RATIONAL(-1.0,4.0)*t350; - t356 = RATIONAL(7.0,20.0)*x+t342+RATIONAL(9.0,20.0); - t343 = RATIONAL(1.0,4.0)*t350; - t355 = RATIONAL(1.0,20.0)*x+t343+RATIONAL(-3.0,20.0); - t348 = RATIONAL(3.0,20.0); - t354 = t348*x+t342+RATIONAL(11.0,20.0); - t353 = t343+t348+RATIONAL(-11.0,20.0)*x; - t349 = RATIONAL(-3.0,10.0); - t345 = RATIONAL(-1.0,50.0); - t341 = RATIONAL(-1.0,2.0)*t350; - t340 = RATIONAL(1.0,2.0)*t350; - t335 = RATIONAL(3.0,50.0)*t357; - t334 = RATIONAL(-3.0,50.0)*t357; - t333 = RATIONAL(-3.0,25.0)*t357; - t332 = RATIONAL(3.0,25.0)*t357; - coeffs_dyy->coeff_m1_m1 = t335+(RATIONAL(59.0,125.0)*x+RATIONAL(-12.0, -125.0)+RATIONAL(-17.0,50.0)*t350)*y+t353; - coeffs_dyy->coeff_0_m1 = (RATIONAL(-89.0,125.0)+RATIONAL(21.0,250.0)*x+ -RATIONAL(13.0,25.0)*t350+RATIONAL(-9.0,50.0)*t351)*y+t354; - coeffs_dyy->coeff_p1_m1 = (RATIONAL(-73.0,125.0)*x+RATIONAL(-36.0,125.0)+ -RATIONAL(9.0,50.0)*t351+t345*t350)*y+t356; - coeffs_dyy->coeff_p2_m1 = t334+(RATIONAL(12.0,125.0)+RATIONAL(-4.0,25.0)* -t350+RATIONAL(7.0,250.0)*x)*y+t355; - coeffs_dyy->coeff_m1_0 = RATIONAL(11.0,10.0)*x+t349+t341+(RATIONAL(243.0, -500.0)+RATIONAL(1.0,25.0)*t351+RATIONAL(69.0,100.0)*t350+RATIONAL(-851.0,500.0) -*x)*y; - coeffs_dyy->coeff_0_0 = t340+t349*x+t333+RATIONAL(-11.0,10.0)+(RATIONAL( --57.0,100.0)*t350+RATIONAL(771.0,500.0)+RATIONAL(303.0,500.0)*x)*y; - coeffs_dyy->coeff_p1_0 = t332+t340+RATIONAL(-9.0,10.0)+RATIONAL(-7.0,10.0 -)*x+(RATIONAL(729.0,500.0)+RATIONAL(-93.0,100.0)*t350+RATIONAL(447.0,500.0)*x)* -y; - coeffs_dyy->coeff_p2_0 = t341+RATIONAL(3.0,10.0)+RATIONAL(-1.0,10.0)*x+( -RATIONAL(-243.0,500.0)+RATIONAL(101.0,500.0)*x+RATIONAL(81.0,100.0)*t350+ -RATIONAL(-1.0,25.0)*t351)*y; - coeffs_dyy->coeff_m1_p1 = (RATIONAL(1.0,50.0)*t351+RATIONAL(-54.0,125.0)+ -RATIONAL(203.0,125.0)*x+RATIONAL(-39.0,50.0)*t350)*y+t353; - coeffs_dyy->coeff_0_p1 = t334+(RATIONAL(-213.0,125.0)+RATIONAL(21.0,25.0) -*t350+RATIONAL(-93.0,250.0)*x)*y+t354; - coeffs_dyy->coeff_p1_p1 = t335+(RATIONAL(-141.0,125.0)*x+RATIONAL(33.0, -50.0)*t350+RATIONAL(-162.0,125.0))*y+t356; - coeffs_dyy->coeff_p2_p1 = (RATIONAL(54.0,125.0)+RATIONAL(-18.0,25.0)*t350 -+RATIONAL(-31.0,250.0)*x+t345*t351)*y+t355; - coeffs_dyy->coeff_m1_p2 = t333+(RATIONAL(21.0,500.0)+RATIONAL(-197.0, -500.0)*x+RATIONAL(43.0,100.0)*t350)*y; - coeffs_dyy->coeff_0_p2 = (RATIONAL(9.0,25.0)*t351+RATIONAL(-159.0,500.0)* -x+RATIONAL(-79.0,100.0)*t350+RATIONAL(437.0,500.0))*y; - coeffs_dyy->coeff_p1_p2 = (RATIONAL(29.0,100.0)*t350+RATIONAL(-9.0,25.0)* -t351+RATIONAL(63.0,500.0)+RATIONAL(409.0,500.0)*x)*y; - coeffs_dyy->coeff_p2_p2 = t332+(RATIONAL(7.0,100.0)*t350+RATIONAL(-53.0, -500.0)*x+RATIONAL(-21.0,500.0))*y; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c deleted file mode 100644 index 1c237bc..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c +++ /dev/null @@ -1,287 +0,0 @@ -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; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c deleted file mode 100644 index 46af5f2..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c +++ /dev/null @@ -1,261 +0,0 @@ -fp t779; -fp t782; -fp t742; -fp t740; -fp t851; -fp t781; -fp t741; -fp t739; -fp t850; -fp t849; -fp t774; -fp t778; -fp t769; -fp t848; -fp t758; -fp t752; -fp t847; -fp t759; -fp t800; -fp t743; -fp t846; -fp t760; -fp t720; -fp t750; -fp t802; -fp t845; -fp t821; -fp t761; -fp t748; -fp t717; -fp t780; -fp t715; -fp t812; -fp t807; -fp t844; -fp t745; -fp t714; -fp t770; -fp t843; -fp t811; -fp t754; -fp t804; -fp t842; -fp t819; -fp t798; -fp t771; -fp t735; -fp t824; -fp t841; -fp t815; -fp t840; -fp t809; -fp t839; -fp t731; -fp t801; -fp t838; -fp t777; -fp t810; -fp t822; -fp t823; -fp t837; -fp t818; -fp t836; -fp t757; -fp t776; -fp t820; -fp t772; -fp t803; -fp t835; -fp t834; -fp t833; -fp t744; -fp t832; -fp t831; -fp t830; -fp t766; -fp t799; -fp t716; -fp t829; -fp t765; -fp t828; -fp t827; -fp t718; -fp t826; -fp t825; -fp t816; -fp t817; -fp t813; -fp t773; -fp t806; -fp t805; -fp t751; -fp t797; -fp t764; -fp t749; -fp t796; -fp t746; -fp t795; -fp t767; -fp t768; -fp t794; -fp t753; -fp t793; -fp t792; -fp t791; -fp t790; -fp t789; -fp t788; -fp t787; -fp t786; -fp t763; -fp t762; -fp t747; - t779 = y*y; - t782 = t779*t779; - t742 = RATIONAL(1.0,288.0); - t740 = RATIONAL(-1.0,288.0); - t851 = t740*t782+t742*t779; - t781 = t779*y; - t741 = RATIONAL(-1.0,144.0); - t739 = RATIONAL(1.0,144.0); - t850 = t741*t781+t739*y; - t849 = t741*y+t739*t781; - t774 = RATIONAL(1.0,3.0); - t778 = x*x; - t769 = RATIONAL(-4.0,9.0); - t848 = t769+t774*t778; - t758 = RATIONAL(-1.0,12.0); - t752 = RATIONAL(-1.0,48.0); - t847 = t752*t782+t758; - t759 = RATIONAL(1.0,24.0); - t800 = t759*t778; - t743 = RATIONAL(1.0,18.0); - t846 = t743*y+t781*t800; - t760 = RATIONAL(-1.0,18.0); - t720 = t760*y; - t750 = RATIONAL(-1.0,24.0); - t802 = t778*t750; - t845 = t720+t781*t802; - t821 = x*t743; - t761 = RATIONAL(1.0,36.0); - t748 = RATIONAL(-1.0,36.0); - t717 = t748*t782; - t780 = t778*x; - t715 = t760*t779; - t812 = t780*t717+x*t715; - t807 = t779*t780; - t844 = t761*t807+t812+t782*t821; - t745 = RATIONAL(1.0,72.0); - t714 = t745*t782; - t770 = RATIONAL(1.0,9.0); - t843 = x*t714+t812+t770*t807; - t811 = t778*t782; - t754 = RATIONAL(1.0,48.0); - t804 = t778*t754; - t842 = t761*t782+t752*t811+(t804+t748)*t779; - t819 = y*t780; - t798 = t780*t743; - t771 = RATIONAL(-1.0,9.0); - t735 = t771*t781; - t824 = x*y; - t841 = t760*t819+t770*t824+t781*t798+x*t735; - t815 = t780*t782; - t840 = t741*t807+t739*t815+t851*x; - t809 = t778*t779; - t839 = t717+t761*t779+t782*t804+t752*t809; - t731 = t770*t781; - t801 = t781*t780; - t838 = y*t798+t760*t801+t771*t824+x*t731; - t777 = RATIONAL(-2.0,9.0); - t810 = t780*t770; - t822 = x*RATIONAL(8.0,9.0); - t823 = x*t782; - t837 = t779*t822+t769*t807+t777*t823+t782*t810; - t818 = x*t781; - t836 = y*t822+t769*t819+t777*t818+t770*t801; - t757 = RATIONAL(1.0,12.0); - t776 = RATIONAL(-1.0,3.0); - t820 = y*t778; - t772 = RATIONAL(4.0,9.0); - t803 = t778*t781; - t835 = t757*t803+t776*t820+t735+t772*y; - t834 = t731+t758*t803+t848*y; - t833 = t772*t819+RATIONAL(2.0,9.0)*t818+t771*t801+RATIONAL(-8.0,9.0)*t824 -; - t744 = RATIONAL(-1.0,72.0); - t832 = t744*t819+t745*t801+t850*x; - t831 = t754*t803+t752*t820+t850; - t830 = t744*t801+t745*t819+t849*x; - t766 = RATIONAL(-1.0,6.0); - t799 = t778*t766; - t716 = t744*t781; - t829 = y*t799+t716+t846; - t765 = RATIONAL(1.0,6.0); - t828 = t765*t809+t714+t715+t782*t802; - t827 = (t800+t744)*t782+(t743+t799)*t779; - t718 = t745*t781; - t826 = t718+t765*t820+t845; - t825 = t752*t803+y*t804+t849; - t816 = RATIONAL(5.0,48.0)*t779; - t817 = t759*t815+x*t816; - t813 = t766*t782; - t773 = RATIONAL(-5.0,3.0); - t806 = t780*t813+t773*x*t779; - t805 = t778*RATIONAL(-1.0,96.0); - t751 = RATIONAL(-5.0,48.0); - t797 = t751*t823+t750*t807+t817; - t764 = RATIONAL(2.0,3.0); - t749 = RATIONAL(5.0,12.0); - t796 = t749*t823+t764*t807+t806; - t746 = RATIONAL(-5.0,24.0); - t795 = t817+t765*t780+t746*t807+t847*x; - t767 = RATIONAL(5.0,6.0); - t768 = RATIONAL(-2.0,3.0); - t794 = t774*t823+RATIONAL(4.0,3.0)*x+t768*t780+t767*t807+t806; - t753 = RATIONAL(1.0,96.0); - t793 = t753*t811+t779*t805+t840+t851; - t792 = t760*t781+y*t802+t844+t846; - t791 = y*t810+t748*t801+t843+(t718+t720)*x; - t790 = t776*t809+t771*t782+t772*t779+t757*t811+t837; - t789 = t761*t801+y*t821+x*t716+t771*t819+t843; - t788 = t753*t809+t740*t779+t840+(t742+t805)*t782; - t787 = t837+t758*t811+t770*t782+t848*t779; - t786 = t743*t781+y*t800+t844+t845; - t763 = RATIONAL(1.0,4.0); - t762 = RATIONAL(-5.0,8.0); - t747 = RATIONAL(5.0,24.0); - coeffs_dx->coeff_m2_m2 = t788+t830+t831; - coeffs_dx->coeff_m1_m2 = t786+t839+t841; - coeffs_dx->coeff_0_m2 = (t757*y+t758*t781)*t780+(t746*y+t747*t781)*x+t797 -; - coeffs_dx->coeff_p1_m2 = t792+t841+t842; - coeffs_dx->coeff_p2_m2 = t793+t825+t830; - coeffs_dx->coeff_m2_m1 = t789+t826+t827; - coeffs_dx->coeff_m1_m1 = t787+t833+t835; - coeffs_dx->coeff_0_m1 = (t765*t781+t768*y)*t780+(RATIONAL(5.0,3.0)*y+ -RATIONAL(-5.0,12.0)*t781)*x+t796; - coeffs_dx->coeff_p1_m1 = t790+t833+t834; - coeffs_dx->coeff_p2_m1 = t789+t828+t829; - coeffs_dx->coeff_m2_0 = t757+t754*t782+t751*t779+(RATIONAL(5.0,16.0)*t779 -+RATIONAL(-1.0,4.0)+RATIONAL(-1.0,16.0)*t782)*t778+t795; - coeffs_dx->coeff_m1_0 = t768+t813+t767*t779+(RATIONAL(1.0,2.0)+t762*t779+ -RATIONAL(1.0,8.0)*t782)*t778+t794; - coeffs_dx->coeff_0_0 = (1.0+t763*t782+RATIONAL(-5.0,4.0)*t779)*t780+( -RATIONAL(25.0,8.0)*t779+RATIONAL(-5.0,2.0)+t762*t782)*x; - coeffs_dx->coeff_p1_0 = RATIONAL(-5.0,6.0)*t779+t764+t765*t782+(RATIONAL( -5.0,8.0)*t779+RATIONAL(-1.0,2.0)+RATIONAL(-1.0,8.0)*t782)*t778+t794; - coeffs_dx->coeff_p2_0 = t816+(t763+RATIONAL(-5.0,16.0)*t779+RATIONAL(1.0, -16.0)*t782)*t778+t795+t847; - coeffs_dx->coeff_m2_p1 = t791+t827+t829; - coeffs_dx->coeff_m1_p1 = t787+t834+t836; - coeffs_dx->coeff_0_p1 = (t766*t781+t764*y)*t780+(t773*y+t749*t781)*x+t796 -; - coeffs_dx->coeff_p1_p1 = t790+t835+t836; - coeffs_dx->coeff_p2_p1 = t791+t826+t828; - coeffs_dx->coeff_m2_p2 = t788+t825+t832; - coeffs_dx->coeff_m1_p2 = t792+t838+t839; - coeffs_dx->coeff_0_p2 = (t757*t781+t758*y)*t780+(t747*y+t746*t781)*x+t797 -; - coeffs_dx->coeff_p1_p2 = t786+t838+t842; - coeffs_dx->coeff_p2_p2 = t793+t831+t832; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c deleted file mode 100644 index 15d553b..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c +++ /dev/null @@ -1,224 +0,0 @@ -fp t1210; -fp t1181; -fp t1212; -fp t1182; -fp t1274; -fp t1185; -fp t1249; -fp t1180; -fp t1245; -fp t1273; -fp t1184; -fp t1244; -fp t1211; -fp t1248; -fp t1183; -fp t1106; -fp t1272; -fp t1271; -fp t1270; -fp t1250; -fp t1197; -fp t1252; -fp t1269; -fp t1268; -fp t1155; -fp t1267; -fp t1098; -fp t1253; -fp t1203; -fp t1266; -fp t1265; -fp t1264; -fp t1205; -fp t1246; -fp t1263; -fp t1262; -fp t1236; -fp t1240; -fp t1261; -fp t1247; -fp t1198; -fp t1260; -fp t1259; -fp t1258; -fp t1176; -fp t1186; -fp t1209; -fp t1242; -fp t1257; -fp t1175; -fp t1229; -fp t1256; -fp t1255; -fp t1237; -fp t1238; -fp t1193; -fp t1204; -fp t1206; -fp t1254; -fp t1251; -fp t1233; -fp t1243; -fp t1187; -fp t1232; -fp t1202; -fp t1239; -fp t1189; -fp t1235; -fp t1231; -fp t1230; -fp t1191; -fp t1171; -fp t1169; -fp t1194; -fp t1228; -fp t1196; -fp t1227; -fp t1188; -fp t1226; -fp t1179; -fp t1225; -fp t1190; -fp t1224; -fp t1223; -fp t1222; -fp t1200; -fp t1201; -fp t1221; -fp t1220; -fp t1219; -fp t1218; -fp t1217; -fp t1174; -fp t1173; -fp t1216; -fp t1215; - t1210 = y*y; - t1181 = RATIONAL(-1.0,48.0); - t1212 = t1210*t1210; - t1182 = RATIONAL(1.0,48.0); - t1274 = t1212*t1182+t1210*t1181; - t1185 = RATIONAL(1.0,24.0); - t1249 = x*t1185; - t1180 = RATIONAL(-1.0,24.0); - t1245 = x*t1180; - t1273 = t1212*t1249+t1210*t1245; - t1184 = RATIONAL(-1.0,12.0); - t1244 = x*t1184; - t1211 = y*t1210; - t1248 = x*t1211; - t1183 = RATIONAL(1.0,12.0); - t1106 = t1183*t1248; - t1272 = t1106+y*t1244; - t1271 = t1274*x; - t1270 = y*t1249+t1211*t1245; - t1250 = x*RATIONAL(-2.0,3.0); - t1197 = RATIONAL(1.0,6.0); - t1252 = x*t1212; - t1269 = t1210*t1250+t1197*t1252; - t1268 = y*t1250+t1197*t1248; - t1155 = t1181*t1212; - t1267 = (t1182*t1210+t1155)*x; - t1098 = t1211*t1244; - t1253 = x*y; - t1203 = RATIONAL(1.0,3.0); - t1266 = t1203*t1253+t1098; - t1265 = t1185*t1248+y*t1245; - t1264 = t1183*t1253+t1098; - t1205 = RATIONAL(-1.0,3.0); - t1246 = x*t1205; - t1263 = t1210*t1246+t1183*t1252; - t1262 = t1212*t1245+t1210*t1249; - t1236 = t1210*t1203; - t1240 = t1212*t1184; - t1261 = (t1236+t1240)*x; - t1247 = x*RATIONAL(2.0,3.0); - t1198 = RATIONAL(-1.0,6.0); - t1260 = t1210*t1247+t1198*t1252; - t1259 = y*t1247+t1198*t1248; - t1258 = y*t1246+t1106; - t1176 = RATIONAL(1.0,72.0); - t1186 = RATIONAL(-1.0,18.0); - t1209 = x*x; - t1242 = t1209*t1240+t1186*t1210; - t1257 = t1242+t1209*t1236+t1176*t1212; - t1175 = RATIONAL(1.0,18.0); - t1229 = t1183*t1209; - t1256 = t1242+t1175*t1212+t1210*t1229; - t1255 = RATIONAL(-1.0,288.0)*t1212+RATIONAL(1.0,288.0)*t1210+t1274*t1209; - t1237 = t1203*t1209; - t1238 = t1210*t1209; - t1193 = RATIONAL(-4.0,3.0); - t1204 = RATIONAL(8.0,9.0); - t1206 = RATIONAL(-2.0,9.0); - t1254 = t1193*t1238+t1204*t1210+(t1237+t1206)*t1212; - t1251 = y*t1209; - t1233 = RATIONAL(1.0,8.0)*t1212; - t1243 = RATIONAL(5.0,48.0)*t1210+t1209*t1233; - t1187 = RATIONAL(-1.0,2.0); - t1232 = t1187*t1209; - t1202 = RATIONAL(-5.0,3.0); - t1239 = t1202*t1210+t1212*t1232; - t1189 = RATIONAL(-5.0,8.0); - t1235 = t1189*t1210; - t1231 = t1211*t1209; - t1230 = RATIONAL(2.0,1.0)*t1209; - t1191 = RATIONAL(1.0,2.0); - t1171 = t1191*t1209; - t1169 = RATIONAL(-2.0,1.0)*t1209; - t1194 = RATIONAL(1.0,4.0); - t1228 = RATIONAL(-5.0,24.0)+t1194*t1209; - t1196 = RATIONAL(-1.0,4.0); - t1227 = t1196*t1209+RATIONAL(5.0,24.0); - t1188 = RATIONAL(-1.0,8.0); - t1226 = RATIONAL(-5.0,48.0)*t1212+t1188*t1238+t1243; - t1179 = RATIONAL(5.0,12.0); - t1225 = t1179*t1212+t1210*t1230+t1239; - t1190 = RATIONAL(4.0,3.0); - t1224 = t1169+t1203*t1212+RATIONAL(5.0,2.0)*t1238+t1190+t1239; - t1223 = t1171+t1209*t1235+t1184+t1155+t1243; - t1222 = t1204*y+t1206*t1211+t1193*t1251+t1203*t1231+t1254; - t1200 = RATIONAL(1.0,9.0); - t1201 = RATIONAL(-1.0,9.0); - t1221 = t1197*t1251+t1201*y+t1200*t1211+t1198*t1231+t1256; - t1220 = RATIONAL(2.0,9.0)*t1211+t1205*t1231+t1190*t1251+RATIONAL(-8.0,9.0 -)*y+t1254; - t1219 = t1201*t1211+t1198*t1251+t1197*t1231+t1200*y+t1256; - t1218 = t1205*t1251+t1175*y+t1257+(RATIONAL(-1.0,72.0)+t1229)*t1211; - t1217 = t1257+t1184*t1231+t1176*t1211+(t1186+t1237)*y; - t1174 = RATIONAL(-1.0,144.0); - t1173 = RATIONAL(1.0,144.0); - t1216 = t1174*y+t1185*t1251+t1180*t1231+t1173*t1211+t1255; - t1215 = t1185*t1231+t1180*t1251+t1173*y+t1174*t1211+t1255; - coeffs_dxx->coeff_m2_m2 = t1216+t1265+t1267; - coeffs_dxx->coeff_m1_m2 = t1219+t1264+t1273; - coeffs_dxx->coeff_0_m2 = t1227*t1211+t1228*y+t1226; - coeffs_dxx->coeff_p1_m2 = t1219+t1262+t1272; - coeffs_dxx->coeff_p2_m2 = t1216+t1270+t1271; - coeffs_dxx->coeff_m2_m1 = t1218+t1263+t1266; - coeffs_dxx->coeff_m1_m1 = t1220+t1260+t1268; - coeffs_dxx->coeff_0_m1 = (t1171+RATIONAL(-5.0,12.0))*t1211+(t1169+ -RATIONAL(5.0,3.0))*y+t1225; - coeffs_dxx->coeff_p1_m1 = t1220+t1259+t1269; - coeffs_dxx->coeff_p2_m1 = t1218+t1258+t1261; - coeffs_dxx->coeff_m2_0 = (RATIONAL(5.0,8.0)*t1210+t1187+t1188*t1212)*x+ -t1223; - coeffs_dxx->coeff_m1_0 = (t1194*t1212+1.0+RATIONAL(-5.0,4.0)*t1210)*x+ -t1224; - coeffs_dxx->coeff_0_0 = t1189*t1212+RATIONAL(-5.0,2.0)+RATIONAL(25.0,8.0) -*t1210+(RATIONAL(3.0,1.0)+RATIONAL(3.0,4.0)*t1212+RATIONAL(-15.0,4.0)*t1210)* -t1209; - coeffs_dxx->coeff_p1_0 = (t1196*t1212+RATIONAL(5.0,4.0)*t1210-1.0)*x+ -t1224; - coeffs_dxx->coeff_p2_0 = (t1191+t1233+t1235)*x+t1223; - coeffs_dxx->coeff_m2_p1 = t1217+t1258+t1263; - coeffs_dxx->coeff_m1_p1 = t1222+t1259+t1260; - coeffs_dxx->coeff_0_p1 = (t1179+t1232)*t1211+(t1202+t1230)*y+t1225; - coeffs_dxx->coeff_p1_p1 = t1222+t1268+t1269; - coeffs_dxx->coeff_p2_p1 = t1217+t1261+t1266; - coeffs_dxx->coeff_m2_p2 = t1215+t1267+t1270; - coeffs_dxx->coeff_m1_p2 = t1221+t1272+t1273; - coeffs_dxx->coeff_0_p2 = t1228*t1211+t1227*y+t1226; - coeffs_dxx->coeff_p1_p2 = t1221+t1262+t1264; - coeffs_dxx->coeff_p2_p2 = t1215+t1265+t1271; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c deleted file mode 100644 index c0d0630..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c +++ /dev/null @@ -1,256 +0,0 @@ -fp t1387; -fp t1428; -fp t1430; -fp t1390; -fp t1485; -fp t1396; -fp t1484; -fp t1397; -fp t1388; -fp t1483; -fp t1426; -fp t1403; -fp t1482; -fp t1458; -fp t1427; -fp t1429; -fp t1389; -fp t1361; -fp t1425; -fp t1481; -fp t1343; -fp t1460; -fp t1480; -fp t1450; -fp t1418; -fp t1398; -fp t1395; -fp t1344; -fp t1479; -fp t1448; -fp t1414; -fp t1417; -fp t1423; -fp t1369; -fp t1401; -fp t1478; -fp t1459; -fp t1416; -fp t1422; -fp t1421; -fp t1402; -fp t1477; -fp t1366; -fp t1367; -fp t1413; -fp t1452; -fp t1476; -fp t1391; -fp t1475; -fp t1350; -fp t1474; -fp t1356; -fp t1400; -fp t1473; -fp t1373; -fp t1472; -fp t1362; -fp t1471; -fp t1351; -fp t1380; -fp t1470; -fp t1455; -fp t1469; -fp t1399; -fp t1355; -fp t1468; -fp t1386; -fp t1415; -fp t1467; -fp t1466; -fp t1449; -fp t1347; -fp t1465; -fp t1461; -fp t1464; -fp t1412; -fp t1463; -fp t1411; -fp t1447; -fp t1462; -fp t1457; -fp t1456; -fp t1454; -fp t1453; -fp t1393; -fp t1451; -fp t1446; -fp t1445; -fp t1409; -fp t1420; -fp t1444; -fp t1394; -fp t1443; -fp t1442; -fp t1441; -fp t1440; -fp t1439; -fp t1438; -fp t1437; -fp t1436; -fp t1435; -fp t1434; -fp t1433; -fp t1424; -fp t1419; -fp t1410; -fp t1408; -fp t1407; -fp t1406; -fp t1405; -fp t1404; -fp t1392; - t1387 = RATIONAL(1.0,144.0); - t1428 = y*y; - t1430 = y*t1428; - t1390 = RATIONAL(-1.0,72.0); - t1485 = t1387*y+t1390*t1430; - t1396 = RATIONAL(-1.0,48.0); - t1484 = t1396*t1428+t1387; - t1397 = RATIONAL(1.0,48.0); - t1388 = RATIONAL(-1.0,144.0); - t1483 = t1397*t1428+t1388; - t1426 = RATIONAL(-8.0,9.0); - t1403 = RATIONAL(4.0,9.0); - t1482 = t1426*y+t1403*t1430; - t1458 = x*t1430; - t1427 = x*x; - t1429 = t1427*x; - t1389 = RATIONAL(1.0,18.0); - t1361 = t1389*t1429; - t1425 = RATIONAL(2.0,9.0); - t1481 = y*t1361+t1425*t1458; - t1343 = t1389*t1430; - t1460 = y*t1429; - t1480 = x*t1343+t1425*t1460; - t1450 = t1428*t1429; - t1418 = RATIONAL(-1.0,9.0); - t1398 = RATIONAL(1.0,12.0); - t1395 = RATIONAL(-1.0,24.0); - t1344 = t1395*t1428; - t1479 = t1398*t1450+t1418*t1429+(t1389+t1344)*x; - t1448 = t1427*t1430; - t1414 = RATIONAL(1.0,6.0); - t1417 = RATIONAL(1.0,9.0); - t1423 = RATIONAL(-1.0,3.0); - t1369 = t1423*t1427; - t1401 = RATIONAL(-1.0,18.0); - t1478 = t1401*t1430+t1414*t1448+(t1369+t1417)*y; - t1459 = x*t1428; - t1416 = RATIONAL(-2.0,3.0); - t1422 = RATIONAL(8.0,9.0); - t1421 = RATIONAL(1.0,3.0); - t1402 = RATIONAL(-4.0,9.0); - t1477 = t1416*t1459+t1421*t1450+t1422*x+t1402*t1429; - t1366 = t1421*t1427; - t1367 = t1421*t1428; - t1413 = RATIONAL(-1.0,4.0); - t1452 = t1413*t1427; - t1476 = t1366+t1367+t1428*t1452+t1402; - t1391 = RATIONAL(1.0,72.0); - t1475 = t1391*t1429+t1395*t1450+t1483*x; - t1350 = t1397*t1427; - t1474 = t1391*t1430+t1395*t1448+(t1350+t1388)*y; - t1356 = t1396*t1427; - t1400 = RATIONAL(1.0,24.0); - t1473 = t1400*t1448+y*t1356+t1485; - t1373 = t1423*t1428; - t1472 = t1401*t1429+t1414*t1450+(t1373+t1417)*x; - t1362 = t1390*t1429; - t1471 = t1362+t1400*t1450+t1484*x; - t1351 = t1395*t1427; - t1380 = t1418*t1430; - t1470 = t1398*t1448+t1380+(t1351+t1389)*y; - t1455 = t1429*t1430; - t1469 = y*t1362+RATIONAL(1.0,36.0)*t1455+t1485*x; - t1399 = RATIONAL(-1.0,12.0); - t1355 = t1400*t1427; - t1468 = t1417*t1430+t1399*t1448+(t1355+t1401)*y; - t1386 = t1418*x; - t1415 = RATIONAL(-1.0,6.0); - t1467 = t1386+t1361+t1415*t1450+x*t1367; - t1466 = t1343+t1415*t1448+(t1366+t1418)*y; - t1449 = t1399*t1429; - t1347 = t1400*t1428; - t1465 = t1417*t1429+t1428*t1449+(t1401+t1347)*x; - t1461 = x*y; - t1464 = RATIONAL(16.0,9.0)*t1461+t1426*t1458+t1482*t1429; - t1412 = RATIONAL(2.0,3.0); - t1463 = t1412*t1459+t1423*t1450+t1426*x+t1403*t1429; - t1411 = RATIONAL(1.0,4.0); - t1447 = t1411*t1428; - t1462 = t1427*t1447+t1369+t1373+t1403; - t1457 = y*t1427; - t1456 = t1429*t1380+y*t1386; - t1454 = t1428*t1427; - t1453 = t1416*t1455+RATIONAL(-10.0,3.0)*t1461; - t1393 = RATIONAL(5.0,24.0); - t1451 = t1393*t1461+t1414*t1455; - t1446 = RATIONAL(1.0,8.0)*t1454+t1389+t1456; - t1445 = t1401+RATIONAL(-1.0,8.0)*t1454+t1456; - t1409 = RATIONAL(4.0,3.0); - t1420 = RATIONAL(5.0,3.0); - t1444 = t1420*t1458+t1409*t1460+t1453; - t1394 = RATIONAL(-5.0,12.0); - t1443 = t1394*t1458+y*t1449+t1451; - t1442 = t1420*t1460+t1409*t1458+t1453; - t1441 = t1399*t1458+t1394*t1460+t1451; - t1440 = t1350+RATIONAL(-1.0,16.0)*t1454+t1469+t1483; - t1439 = RATIONAL(1.0,16.0)*t1454+t1356+t1469+t1484; - t1438 = t1351+t1415*t1428+t1446+t1481; - t1437 = t1344+t1415*t1427+t1446+t1480; - t1436 = t1347+t1414*t1427+t1445+t1480; - t1435 = t1402*t1430+t1422*y+t1416*t1457+t1421*t1448+t1464; - t1434 = t1414*t1428+t1355+t1445+t1481; - t1433 = t1423*t1448+t1412*t1457+t1464+t1482; - t1424 = RATIONAL(5.0,4.0); - t1419 = RATIONAL(-5.0,3.0); - t1410 = RATIONAL(1.0,2.0); - t1408 = RATIONAL(-5.0,8.0); - t1407 = RATIONAL(-5.0,4.0); - t1406 = RATIONAL(-5.0,2.0); - t1405 = RATIONAL(-1.0,2.0); - t1404 = RATIONAL(5.0,8.0); - t1392 = RATIONAL(-5.0,24.0); - coeffs_dxy->coeff_m2_m2 = t1439+t1474+t1475; - coeffs_dxy->coeff_m1_m2 = t1434+t1470+t1472; - coeffs_dxy->coeff_0_m2 = (t1413*t1428+t1398)*t1429+(t1404*t1428+t1392)*x+ -t1443; - coeffs_dxy->coeff_p1_m2 = t1438+t1468+t1472; - coeffs_dxy->coeff_p2_m2 = t1440+t1473+t1475; - coeffs_dxy->coeff_m2_m1 = t1436+t1478+t1479; - coeffs_dxy->coeff_m1_m1 = t1433+t1462+t1463; - coeffs_dxy->coeff_0_m1 = (t1410*t1428+t1416)*t1429+(t1420+t1407*t1428)*x+ -t1444; - coeffs_dxy->coeff_p1_m1 = t1435+t1463+t1476; - coeffs_dxy->coeff_p2_m1 = t1437+t1466+t1479; - coeffs_dxy->coeff_m2_0 = (t1452+t1398)*t1430+(t1392+t1404*t1427)*y+t1441; - coeffs_dxy->coeff_m1_0 = (t1410*t1427+t1416)*t1430+(t1407*t1427+t1420)*y+ -t1442; - coeffs_dxy->coeff_0_0 = (t1406*y+t1430)*t1429+(RATIONAL(25.0,4.0)*y+t1406 -*t1430)*x; - coeffs_dxy->coeff_p1_0 = (t1405*t1427+t1412)*t1430+(t1424*t1427+t1419)*y+ -t1442; - coeffs_dxy->coeff_p2_0 = (t1399+t1411*t1427)*t1430+(t1393+t1408*t1427)*y+ -t1441; - coeffs_dxy->coeff_m2_p1 = t1437+t1465+t1478; - coeffs_dxy->coeff_m1_p1 = t1433+t1476+t1477; - coeffs_dxy->coeff_0_p1 = (t1405*t1428+t1412)*t1429+(t1419+t1424*t1428)*x+ -t1444; - coeffs_dxy->coeff_p1_p1 = t1435+t1462+t1477; - coeffs_dxy->coeff_p2_p1 = t1436+t1465+t1466; - coeffs_dxy->coeff_m2_p2 = t1440+t1471+t1474; - coeffs_dxy->coeff_m1_p2 = t1438+t1467+t1470; - coeffs_dxy->coeff_0_p2 = (t1447+t1399)*t1429+(t1393+t1408*t1428)*x+t1443; - coeffs_dxy->coeff_p1_p2 = t1434+t1467+t1468; - coeffs_dxy->coeff_p2_p2 = t1439+t1471+t1473; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c deleted file mode 100644 index 59006b6..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c +++ /dev/null @@ -1,255 +0,0 @@ -fp t1007; -fp t1006; -fp t1013; -fp t1015; -fp t964; -fp t1082; -fp t974; -fp t976; -fp t1081; -fp t1080; -fp t975; -fp t977; -fp t1017; -fp t1079; -fp t1004; -fp t1005; -fp t1078; -fp t987; -fp t993; -fp t1077; -fp t985; -fp t1014; -fp t1040; -fp t995; -fp t957; -fp t1076; -fp t994; -fp t978; -fp t952; -fp t1075; -fp t1052; -fp t996; -fp t1016; -fp t1039; -fp t983; -fp t1042; -fp t958; -fp t1048; -fp t1074; -fp t980; -fp t959; -fp t1038; -fp t1073; -fp t972; -fp t1054; -fp t1050; -fp t1043; -fp t1072; -fp t1071; -fp t1049; -fp t1010; -fp t1070; -fp t1053; -fp t1009; -fp t1069; -fp t992; -fp t1011; -fp t1068; -fp t1067; -fp t979; -fp t1066; -fp t989; -fp t1065; -fp t955; -fp t1064; -fp t1001; -fp t951; -fp t1063; -fp t1062; -fp t1061; -fp t1060; -fp t1059; -fp t1058; -fp t1000; -fp t1057; -fp t1056; -fp t1051; -fp t1055; -fp t1035; -fp t1008; -fp t1047; -fp t1045; -fp t1046; -fp t1041; -fp t1037; -fp t1036; -fp t1003; -fp t973; -fp t968; -fp t984; -fp t999; -fp t1034; -fp t986; -fp t1033; -fp t1032; -fp t981; -fp t1031; -fp t1030; -fp t1029; -fp t1028; -fp t988; -fp t990; -fp t1027; -fp t1026; -fp t1025; -fp t1024; -fp t1023; -fp t1022; -fp t1021; -fp t998; -fp t997; - t1007 = RATIONAL(4.0,9.0); - t1006 = RATIONAL(-1.0,9.0); - t1013 = x*x; - t1015 = t1013*x; - t964 = t1006*t1015; - t1082 = t964+t1007*x; - t974 = RATIONAL(1.0,144.0); - t976 = RATIONAL(-1.0,144.0); - t1081 = t974*x+t976*t1015; - t1080 = t974*t1015+t976*x; - t975 = RATIONAL(-1.0,288.0); - t977 = RATIONAL(1.0,288.0); - t1017 = t1013*t1013; - t1079 = t975*t1017+t977*t1013; - t1004 = RATIONAL(-4.0,9.0); - t1005 = RATIONAL(1.0,9.0); - t1078 = t1004*t1013+t1005*t1017; - t987 = RATIONAL(-1.0,48.0); - t993 = RATIONAL(-1.0,12.0); - t1077 = t987*t1017+t993; - t985 = RATIONAL(-1.0,24.0); - t1014 = y*y; - t1040 = t1014*t1015; - t995 = RATIONAL(-1.0,18.0); - t957 = t995*x; - t1076 = t985*t1040+t957; - t994 = RATIONAL(1.0,24.0); - t978 = RATIONAL(1.0,18.0); - t952 = t978*x; - t1075 = t994*t1040+t952; - t1052 = y*t1017; - t996 = RATIONAL(1.0,36.0); - t1016 = t1014*y; - t1039 = t1013*t1016; - t983 = RATIONAL(-1.0,36.0); - t1042 = t1017*t1016; - t958 = t995*t1013; - t1048 = t983*t1042+y*t958; - t1074 = t978*t1052+t1048+t996*t1039; - t980 = RATIONAL(1.0,72.0); - t959 = t980*t1017; - t1038 = t1005*t1016; - t1073 = y*t959+t1048+t1013*t1038; - t972 = t1005*t1015; - t1054 = x*y; - t1050 = x*t1016; - t1043 = t1016*t1015; - t1072 = t1006*t1054+y*t972+t978*t1050+t995*t1043; - t1071 = RATIONAL(2.0,9.0)*y*t1015+RATIONAL(-8.0,9.0)*t1054+t1082*t1016; - t1049 = y*RATIONAL(-2.0,9.0); - t1010 = RATIONAL(8.0,9.0); - t1070 = t1004*t1050+t1016*t972+t1010*t1054+t1015*t1049; - t1053 = x*t1014; - t1009 = RATIONAL(1.0,3.0); - t1069 = t972+t1004*x+t993*t1040+t1009*t1053; - t992 = RATIONAL(1.0,12.0); - t1011 = RATIONAL(-1.0,3.0); - t1068 = t1011*t1053+t992*t1040+t1082; - t1067 = t995*t1050+y*t964+t978*t1043+t1005*t1054; - t979 = RATIONAL(-1.0,72.0); - t1066 = t980*t1043+t979*t1050+t1081*y; - t989 = RATIONAL(1.0,48.0); - t1065 = t987*t1053+t989*t1040+t1081; - t955 = t980*t1015; - t1064 = t983*t1043+x*t1038+(t957+t955)*y; - t1001 = RATIONAL(-1.0,6.0); - t951 = t979*t1015; - t1063 = t1001*t1053+t951+t1075; - t1062 = t996*t1043+t1006*t1050+(t952+t951)*y; - t1061 = t979*t1043+t980*t1050+t1080*y; - t1060 = t989*t1053+t987*t1040+t1080; - t1059 = t978*t1015+t994*t1053+t1076; - t1058 = t995*t1015+t985*t1053+t1075; - t1000 = RATIONAL(1.0,6.0); - t1057 = t955+t1000*t1053+t1076; - t1056 = t974*t1042+t976*t1039+t1079*y; - t1051 = y*t1013; - t1055 = t1010*t1051+t1017*t1049+t1078*t1016; - t1035 = t1001*t1017; - t1008 = RATIONAL(-5.0,3.0); - t1047 = t1008*t1051+t1016*t1035; - t1045 = RATIONAL(5.0,48.0)*t1013; - t1046 = t994*t1042+y*t1045; - t1041 = t1017*t1014; - t1037 = t1013*t1014; - t1036 = RATIONAL(5.0,6.0)*t1013; - t1003 = RATIONAL(-2.0,3.0); - t973 = t1003*t1016; - t968 = t1000*t1016; - t984 = RATIONAL(5.0,12.0); - t999 = RATIONAL(2.0,3.0); - t1034 = t984*t1052+t999*t1039+t1047; - t986 = RATIONAL(-5.0,48.0); - t1033 = t986*t1052+t985*t1039+t1046; - t1032 = RATIONAL(5.0,24.0)*y+t993*t1016; - t981 = RATIONAL(-5.0,24.0); - t1031 = t992*t1016+t981*y; - t1030 = t973+t1009*t1052+t1016*t1036+RATIONAL(4.0,3.0)*y+t1047; - t1029 = t968+t981*t1039+t1046+t1077*y; - t1028 = t994*t1041+t1001*t1037+t978*t1013+t979*t1017+t1073; - t988 = RATIONAL(1.0,96.0); - t990 = RATIONAL(-1.0,96.0); - t1027 = t988*t1037+t977*t1017+t975*t1013+t990*t1041+t1056; - t1026 = t1007*t1013+t992*t1041+t1006*t1017+t1011*t1037+t1055; - t1025 = t987*t1041+t983*t1013+t996*t1017+t989*t1037+t1074; - t1024 = t993*t1041+t1009*t1037+t1055+t1078; - t1023 = t987*t1037+t983*t1017+t996*t1013+t989*t1041+t1074; - t1022 = t990*t1037+t988*t1041+t1056+t1079; - t1021 = t959+t985*t1041+t1000*t1037+t958+t1073; - t998 = RATIONAL(1.0,4.0); - t997 = RATIONAL(-5.0,8.0); - coeffs_dy->coeff_m2_m2 = t1027+t1061+t1065; - coeffs_dy->coeff_m1_m2 = t1028+t1057+t1062; - coeffs_dy->coeff_0_m2 = t986*t1013+t989*t1017+t992+(RATIONAL(5.0,16.0)* -t1013+RATIONAL(-1.0,4.0)+RATIONAL(-1.0,16.0)*t1017)*t1014+t1029; - coeffs_dy->coeff_p1_m2 = t1028+t1063+t1064; - coeffs_dy->coeff_p2_m2 = t1027+t1060+t1066; - coeffs_dy->coeff_m2_m1 = t1023+t1059+t1067; - coeffs_dy->coeff_m1_m1 = t1024+t1068+t1071; - coeffs_dy->coeff_0_m1 = t1036+t1003+t1035+(t997*t1013+RATIONAL(1.0,2.0)+ -RATIONAL(1.0,8.0)*t1017)*t1014+t1030; - coeffs_dy->coeff_p1_m1 = t1024+t1069+t1070; - coeffs_dy->coeff_p2_m1 = t1023+t1058+t1072; - coeffs_dy->coeff_m2_0 = t1032*t1015+t1031*x+t1033; - coeffs_dy->coeff_m1_0 = (RATIONAL(-5.0,12.0)*y+t968)*t1015+(RATIONAL(5.0, -3.0)*y+t973)*x+t1034; - coeffs_dy->coeff_0_0 = (RATIONAL(-5.0,4.0)*t1013+t998*t1017+1.0)*t1016+( -t997*t1017+RATIONAL(25.0,8.0)*t1013+RATIONAL(-5.0,2.0))*y; - coeffs_dy->coeff_p1_0 = (t1001*t1016+t984*y)*t1015+(t1008*y+t999*t1016)*x -+t1034; - coeffs_dy->coeff_p2_0 = t1031*t1015+t1032*x+t1033; - coeffs_dy->coeff_m2_p1 = t1025+t1058+t1067; - coeffs_dy->coeff_m1_p1 = t1026+t1069+t1071; - coeffs_dy->coeff_0_p1 = t1000*t1017+t999+RATIONAL(-5.0,6.0)*t1013+( -RATIONAL(5.0,8.0)*t1013+RATIONAL(-1.0,8.0)*t1017+RATIONAL(-1.0,2.0))*t1014+ -t1030; - coeffs_dy->coeff_p1_p1 = t1026+t1068+t1070; - coeffs_dy->coeff_p2_p1 = t1025+t1059+t1072; - coeffs_dy->coeff_m2_p2 = t1022+t1060+t1061; - coeffs_dy->coeff_m1_p2 = t1021+t1062+t1063; - coeffs_dy->coeff_0_p2 = t1045+(RATIONAL(1.0,16.0)*t1017+RATIONAL(-5.0, -16.0)*t1013+t998)*t1014+t1029+t1077; - coeffs_dy->coeff_p1_p2 = t1021+t1057+t1064; - coeffs_dy->coeff_p2_p2 = t1022+t1065+t1066; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c deleted file mode 100644 index b1a740e..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c +++ /dev/null @@ -1,223 +0,0 @@ -fp t1613; -fp t1675; -fp t1611; -fp t1645; -fp t1612; -fp t1614; -fp t1586; -fp t1644; -fp t1504; -fp t1674; -fp t1585; -fp t1615; -fp t1649; -fp t1584; -fp t1635; -fp t1673; -fp t1583; -fp t1646; -fp t1648; -fp t1588; -fp t1672; -fp t1650; -fp t1671; -fp t1555; -fp t1654; -fp t1670; -fp t1651; -fp t1604; -fp t1647; -fp t1669; -fp t1603; -fp t1605; -fp t1668; -fp t1587; -fp t1643; -fp t1494; -fp t1667; -fp t1666; -fp t1633; -fp t1665; -fp t1664; -fp t1663; -fp t1609; -fp t1652; -fp t1662; -fp t1661; -fp t1660; -fp t1659; -fp t1589; -fp t1632; -fp t1636; -fp t1642; -fp t1579; -fp t1658; -fp t1580; -fp t1657; -fp t1601; -fp t1656; -fp t1655; -fp t1653; -fp t1641; -fp t1598; -fp t1640; -fp t1608; -fp t1591; -fp t1639; -fp t1568; -fp t1597; -fp t1573; -fp t1595; -fp t1637; -fp t1634; -fp t1599; -fp t1631; -fp t1602; -fp t1630; -fp t1592; -fp t1629; -fp t1578; -fp t1628; -fp t1596; -fp t1627; -fp t1626; -fp t1576; -fp t1577; -fp t1625; -fp t1624; -fp t1623; -fp t1622; -fp t1621; -fp t1620; -fp t1606; -fp t1607; -fp t1619; -fp t1618; - t1613 = y*y; - t1675 = RATIONAL(8.0,9.0)+t1613*RATIONAL(-4.0,3.0); - t1611 = RATIONAL(-1.0,3.0); - t1645 = t1611*x; - t1612 = x*x; - t1614 = t1612*x; - t1586 = RATIONAL(1.0,12.0); - t1644 = y*t1586; - t1504 = t1614*t1644; - t1674 = y*t1645+t1504; - t1585 = RATIONAL(1.0,48.0); - t1615 = t1612*t1612; - t1649 = y*t1615; - t1584 = RATIONAL(-1.0,48.0); - t1635 = t1612*t1584; - t1673 = t1585*t1649+y*t1635; - t1583 = RATIONAL(-1.0,24.0); - t1646 = x*t1583; - t1648 = y*t1614; - t1588 = RATIONAL(1.0,24.0); - t1672 = y*t1646+t1588*t1648; - t1650 = x*t1588; - t1671 = y*t1650+t1583*t1648; - t1555 = t1584*t1615; - t1654 = y*t1612; - t1670 = y*t1555+t1585*t1654; - t1651 = RATIONAL(2.0,3.0)*y; - t1604 = RATIONAL(-1.0,6.0); - t1647 = y*t1604; - t1669 = t1614*t1647+x*t1651; - t1603 = RATIONAL(1.0,6.0); - t1605 = RATIONAL(-2.0,3.0); - t1668 = t1605*t1654+t1603*t1649; - t1587 = RATIONAL(-1.0,12.0); - t1643 = y*t1587; - t1494 = t1614*t1643; - t1667 = x*t1644+t1494; - t1666 = t1583*t1654+t1588*t1649; - t1633 = t1603*t1614; - t1665 = (t1605*x+t1633)*y; - t1664 = t1615*t1647+t1612*t1651; - t1663 = t1615*t1644+t1611*t1654; - t1609 = RATIONAL(1.0,3.0); - t1652 = y*t1609; - t1662 = t1494+x*t1652; - t1661 = t1583*t1649+t1588*t1654; - t1660 = x*t1643+t1504; - t1659 = t1612*t1652+t1615*t1643; - t1589 = RATIONAL(-1.0,18.0); - t1632 = t1613*t1615; - t1636 = t1589*t1612+t1587*t1632; - t1642 = t1613*t1612; - t1579 = RATIONAL(1.0,18.0); - t1658 = t1636+t1579*t1615+t1586*t1642; - t1580 = RATIONAL(1.0,72.0); - t1657 = t1636+t1580*t1615+t1609*t1642; - t1601 = RATIONAL(-2.0,9.0); - t1656 = t1601*t1615+t1609*t1632+t1675*t1612; - t1655 = RATIONAL(-1.0,288.0)*t1615+t1613*t1635+RATIONAL(1.0,288.0)*t1612+ -t1585*t1632; - t1653 = x*t1613; - t1641 = RATIONAL(2.0,1.0)*t1613; - t1598 = RATIONAL(1.0,8.0); - t1640 = RATIONAL(5.0,48.0)*t1612+t1598*t1632; - t1608 = RATIONAL(-5.0,3.0); - t1591 = RATIONAL(-1.0,2.0); - t1639 = t1608*t1612+t1591*t1632; - t1568 = RATIONAL(-2.0,1.0)*t1613; - t1597 = RATIONAL(1.0,2.0); - t1573 = t1597*t1613; - t1595 = RATIONAL(-5.0,8.0); - t1637 = t1595*t1612; - t1634 = t1613*t1614; - t1599 = RATIONAL(1.0,4.0); - t1631 = t1599*t1613+RATIONAL(-5.0,24.0); - t1602 = RATIONAL(-1.0,4.0); - t1630 = t1602*t1613+RATIONAL(5.0,24.0); - t1592 = RATIONAL(-1.0,8.0); - t1629 = t1592*t1642+RATIONAL(-5.0,48.0)*t1615+t1640; - t1578 = RATIONAL(5.0,12.0); - t1628 = t1612*t1641+t1578*t1615+t1639; - t1596 = RATIONAL(4.0,3.0); - t1627 = t1609*t1615+RATIONAL(5.0,2.0)*t1642+t1568+t1596+t1639; - t1626 = t1555+t1587+t1573+t1613*t1637+t1640; - t1576 = RATIONAL(1.0,144.0); - t1577 = RATIONAL(-1.0,144.0); - t1625 = t1613*t1646+t1577*t1614+t1576*x+t1588*t1634+t1655; - t1624 = t1583*t1634+t1577*x+t1613*t1650+t1576*t1614+t1655; - t1623 = t1609*t1653+t1589*x+t1580*t1614+t1587*t1634+t1657; - t1622 = t1579*x+RATIONAL(-1.0,72.0)*t1614+t1586*t1634+t1613*t1645+t1657; - t1621 = t1601*t1614+t1609*t1634+t1656+t1675*x; - t1620 = RATIONAL(-8.0,9.0)*x+RATIONAL(2.0,9.0)*t1614+t1596*t1653+t1611* -t1634+t1656; - t1606 = RATIONAL(1.0,9.0); - t1607 = RATIONAL(-1.0,9.0); - t1619 = t1607*x+t1604*t1634+t1606*t1614+t1603*t1653+t1658; - t1618 = t1606*x+t1607*t1614+t1613*t1633+t1604*t1653+t1658; - coeffs_dyy->coeff_m2_m2 = t1624+t1670+t1672; - coeffs_dyy->coeff_m1_m2 = t1622+t1662+t1663; - coeffs_dyy->coeff_0_m2 = (RATIONAL(5.0,8.0)*t1612+t1592*t1615+t1591)*y+ -t1626; - coeffs_dyy->coeff_p1_m2 = t1623+t1663+t1674; - coeffs_dyy->coeff_p2_m2 = t1625+t1670+t1671; - coeffs_dyy->coeff_m2_m1 = t1618+t1666+t1667; - coeffs_dyy->coeff_m1_m1 = t1620+t1664+t1665; - coeffs_dyy->coeff_0_m1 = (1.0+t1599*t1615+RATIONAL(-5.0,4.0)*t1612)*y+ -t1627; - coeffs_dyy->coeff_p1_m1 = t1621+t1664+t1669; - coeffs_dyy->coeff_p2_m1 = t1619+t1660+t1666; - coeffs_dyy->coeff_m2_0 = t1630*t1614+t1631*x+t1629; - coeffs_dyy->coeff_m1_0 = (t1573+RATIONAL(-5.0,12.0))*t1614+(RATIONAL(5.0, -3.0)+t1568)*x+t1628; - coeffs_dyy->coeff_0_0 = RATIONAL(25.0,8.0)*t1612+RATIONAL(-5.0,2.0)+t1595 -*t1615+(RATIONAL(3.0,1.0)+RATIONAL(3.0,4.0)*t1615+RATIONAL(-15.0,4.0)*t1612)* -t1613; - coeffs_dyy->coeff_p1_0 = (t1591*t1613+t1578)*t1614+(t1641+t1608)*x+t1628; - coeffs_dyy->coeff_p2_0 = t1631*t1614+t1630*x+t1629; - coeffs_dyy->coeff_m2_p1 = t1618+t1660+t1661; - coeffs_dyy->coeff_m1_p1 = t1620+t1668+t1669; - coeffs_dyy->coeff_0_p1 = (t1602*t1615-1.0+RATIONAL(5.0,4.0)*t1612)*y+ -t1627; - coeffs_dyy->coeff_p1_p1 = t1621+t1665+t1668; - coeffs_dyy->coeff_p2_p1 = t1619+t1661+t1667; - coeffs_dyy->coeff_m2_p2 = t1624+t1671+t1673; - coeffs_dyy->coeff_m1_p2 = t1622+t1659+t1674; - coeffs_dyy->coeff_0_p2 = (t1637+t1597+t1598*t1615)*y+t1626; - coeffs_dyy->coeff_p1_p2 = t1623+t1659+t1662; - coeffs_dyy->coeff_p2_p2 = t1625+t1672+t1673; |