diff options
author | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-08-20 16:31:28 +0000 |
---|---|---|
committer | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-08-20 16:31:28 +0000 |
commit | d5b49cab5a99a22dd26bd55fcf65659aa7f7e020 (patch) | |
tree | e3df72179dc33241852470ee9f86b6d9e67b43e4 /src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs | |
parent | d1aaafaad12880241613d3465da77245ab51bbff (diff) |
add $Header:$ lines, switch a few $Id:$ --> $Header:$
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@92 df1f8a13-aa1d-4dd4-9681-27ded5b42416
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs')
17 files changed, 1152 insertions, 1144 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c index 34f402b..52b6d48 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c @@ -1,49 +1,53 @@ -fp t25, - t28, - t20, - t27, +fp t27, + t26, t18, + t28, + t21, t34, t33, t23, + t15, t22, - t17, t32, - t15, + t17, t31, + t14, t30, + t25, + t20, t29, - t26, - t21, t19, t16, - t14; - t25 = RATIONAL(1.0,6.0); - t28 = y*y; - t20 = t25*t28; + t13, + t12; t27 = x*x; - t18 = t25*t27; - t34 = t20+RATIONAL(-1.0,9.0)+t18; + t26 = RATIONAL(1.0,6.0); + t18 = t26*t27; + t28 = y*y; + t21 = t26*t28; + t34 = t18+t21+RATIONAL(-1.0,9.0); t33 = x*y; - t23 = RATIONAL(2.0,9.0); - t22 = RATIONAL(-1.0,3.0); - t17 = t22*t27; - t32 = t20+t23+t17; - t15 = t22*t28; - t31 = t15+t23+t18; - t30 = RATIONAL(-1.0,4.0)*t33+t34; - t29 = RATIONAL(1.0,4.0)*t33+t34; - t26 = RATIONAL(-1.0,6.0); - t21 = t26*y; - t19 = t26*x; - t16 = t25*x; - t14 = t25*y; - coeff_I_m1_m1 = t21+t19+t29; - coeff_I_0_m1 = t21+t32; - coeff_I_p1_m1 = t21+t16+t30; - coeff_I_m1_0 = t19+t31; - coeff_I_0_0 = t15+t17+RATIONAL(5.0,9.0); - coeff_I_p1_0 = t16+t31; - coeff_I_m1_p1 = t14+t19+t30; - coeff_I_0_p1 = t14+t32; - coeff_I_p1_p1 = t14+t16+t29; + t23 = RATIONAL(-1.0,3.0); + t15 = t23*t28; + t22 = RATIONAL(2.0,9.0); + t32 = t15+t18+t22; + t17 = t23*t27; + t31 = t17+t21+t22; + t14 = t26*y; + t30 = t14+t34; + t25 = RATIONAL(-1.0,6.0); + t20 = t25*y; + t29 = t20+t34; + t19 = t25*x; + t16 = t26*x; + t13 = RATIONAL(1.0,4.0)*t33; + t12 = RATIONAL(-1.0,4.0)*t33; + coeff_I_m1_m1 = t13+t19+t29; + coeff_I_0_m1 = t20+t31; + coeff_I_p1_m1 = t12+t16+t29; + coeff_I_m1_0 = t19+t32; + coeff_I_0_0 = t15+RATIONAL(5.0,9.0)+t17; + coeff_I_p1_0 = t16+t32; + coeff_I_m1_p1 = t19+t12+t30; + coeff_I_0_p1 = t14+t31; + coeff_I_p1_p1 = t16+t13+t30; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c index 22a3439..14b1a1c 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c @@ -3,19 +3,19 @@ fp t37, t41, t40, t39, - t35; + t36; t37 = RATIONAL(1.0,3.0)*x; - t42 = t37+RATIONAL(1.0,4.0)*y; - t41 = RATIONAL(-1.0,4.0)*y+t37; - t40 = RATIONAL(-1.0,6.0); - t39 = RATIONAL(1.0,6.0); - t35 = RATIONAL(-2.0,3.0)*x; - coeff_dx_m1_m1 = t40+t42; - coeff_dx_0_m1 = t35; - coeff_dx_p1_m1 = t39+t41; - coeff_dx_m1_0 = t40+t37; - coeff_dx_0_0 = t35; - coeff_dx_p1_0 = t39+t37; - coeff_dx_m1_p1 = t40+t41; - coeff_dx_0_p1 = t35; - coeff_dx_p1_p1 = t39+t42; + t42 = RATIONAL(1.0,4.0)*y+t37; + t41 = t37+RATIONAL(-1.0,4.0)*y; + t40 = RATIONAL(1.0,6.0); + t39 = RATIONAL(-1.0,6.0); + t36 = RATIONAL(-2.0,3.0)*x; + coeff_dx_m1_m1 = t39+t42; + coeff_dx_0_m1 = t36; + coeff_dx_p1_m1 = t40+t41; + coeff_dx_m1_0 = t39+t37; + coeff_dx_0_0 = t36; + coeff_dx_p1_0 = t40+t37; + coeff_dx_m1_p1 = t39+t41; + coeff_dx_0_p1 = t36; + coeff_dx_p1_p1 = t40+t42; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c index 5f6427a..3381569 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c @@ -1,13 +1,13 @@ fp t52, t51; - t52 = RATIONAL(1.0,3.0); - t51 = RATIONAL(-2.0,3.0); - coeff_dxx_m1_m1 = t52; - coeff_dxx_0_m1 = t51; - coeff_dxx_p1_m1 = t52; - coeff_dxx_m1_0 = t52; - coeff_dxx_0_0 = t51; - coeff_dxx_p1_0 = t52; - coeff_dxx_m1_p1 = t52; - coeff_dxx_0_p1 = t51; - coeff_dxx_p1_p1 = t52; + t52 = RATIONAL(-2.0,3.0); + t51 = RATIONAL(1.0,3.0); + coeff_dxx_m1_m1 = t51; + coeff_dxx_0_m1 = t52; + coeff_dxx_p1_m1 = t51; + coeff_dxx_m1_0 = t51; + coeff_dxx_0_0 = t52; + coeff_dxx_p1_0 = t51; + coeff_dxx_m1_p1 = t51; + coeff_dxx_0_p1 = t52; + coeff_dxx_p1_p1 = t51; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c index 4afa754..65e74bd 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c @@ -5,17 +5,17 @@ fp t46, t47, t43; t46 = RATIONAL(1.0,3.0)*y; - t50 = RATIONAL(-1.0,4.0)*x+t46; - t49 = RATIONAL(1.0,4.0)*x+t46; - t48 = RATIONAL(-1.0,6.0); - t47 = RATIONAL(1.0,6.0); + t50 = RATIONAL(1.0,4.0)*x+t46; + t49 = RATIONAL(-1.0,4.0)*x+t46; + t48 = RATIONAL(1.0,6.0); + t47 = RATIONAL(-1.0,6.0); t43 = RATIONAL(-2.0,3.0)*y; - coeff_dy_m1_m1 = t48+t49; - coeff_dy_0_m1 = t46+t48; - coeff_dy_p1_m1 = t48+t50; + coeff_dy_m1_m1 = t47+t50; + coeff_dy_0_m1 = t47+t46; + coeff_dy_p1_m1 = t47+t49; coeff_dy_m1_0 = t43; coeff_dy_0_0 = t43; coeff_dy_p1_0 = t43; - coeff_dy_m1_p1 = t47+t50; - coeff_dy_0_p1 = t47+t46; - coeff_dy_p1_p1 = t47+t49; + coeff_dy_m1_p1 = t48+t49; + coeff_dy_0_p1 = t46+t48; + coeff_dy_p1_p1 = t48+t50; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c index a643108..5d0b896 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c @@ -1,13 +1,13 @@ fp t57, t56; - t57 = RATIONAL(1.0,3.0); - t56 = RATIONAL(-2.0,3.0); - coeff_dyy_m1_m1 = t57; - coeff_dyy_0_m1 = t57; - coeff_dyy_p1_m1 = t57; - coeff_dyy_m1_0 = t56; - coeff_dyy_0_0 = t56; - coeff_dyy_p1_0 = t56; - coeff_dyy_m1_p1 = t57; - coeff_dyy_0_p1 = t57; - coeff_dyy_p1_p1 = t57; + t57 = RATIONAL(-2.0,3.0); + t56 = RATIONAL(1.0,3.0); + coeff_dyy_m1_m1 = t56; + coeff_dyy_0_m1 = t56; + coeff_dyy_p1_m1 = t56; + coeff_dyy_m1_0 = t57; + coeff_dyy_0_0 = t57; + coeff_dyy_p1_0 = t57; + coeff_dyy_m1_p1 = t56; + coeff_dyy_0_p1 = t56; + coeff_dyy_p1_p1 = t56; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c index 55cc52e..a022b8a 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c @@ -1,73 +1,73 @@ fp t160, t131, t159, - t130, - t158, - t124, - t68, t133, - t126, - t82, - t157, - t71, + t158, + t117, + t66, t132, - t83, - t156, - t121, - t70, t128, t87, - t155, - t67, + t157, t129, - t86, + t84, + t156, + t124, + t68, + t130, + t116, + t77, + t155, + t126, + t70, + t125, + t83, t154, - t85, + t71, t153, - t117, - t77, + t108, + t122, + t69, t152, - t127, - t62, t151, t63, t150, - t125, - t64, + t67, + t86, t149, - t79, + t65, t148, + t85, t147, + t79, t146, - t69, - t110, + t82, t145, - t84, t144, t143, - t123, t142, - t104, + t123, t141, + t101, t140, t139, t138, t137, t136, - t122, - t120, + t127, + t121, t119, - t116, - t114, + t118, + t115, t113, + t112, t111, - t108, - t107, + t109, + t106, t105, t103, t102, - t101, - t99, + t100, t98, t97, t96, @@ -85,112 +85,112 @@ fp t160, t76, t75, t74, - t73; + t72; t160 = x*y; - t131 = y*y; - t159 = x*t131; - t130 = x*x; - t158 = t130*y; + t131 = x*x; + t159 = t131*y; + t133 = y*y; + t158 = x*t133; + t117 = RATIONAL(1.0,40.0); + t66 = t117*t159; + t132 = t133*y; + t128 = RATIONAL(-1.0,8.0); + t87 = t128*t132; + t157 = t66+t87; + t129 = RATIONAL(1.0,8.0); + t84 = t129*t132; + t156 = t66+t84; t124 = RATIONAL(3.0,40.0); t68 = t124*t158; - t133 = t131*y; - t126 = RATIONAL(-1.0,24.0); - t82 = t126*t133; - t157 = t68+t82; + t130 = t131*x; + t116 = RATIONAL(1.0,24.0); + t77 = t116*t130; + t155 = t68+t77; + t126 = RATIONAL(-3.0,40.0); + t70 = t126*t159; + t125 = RATIONAL(-1.0,24.0); + t83 = t125*t132; + t154 = t70+t83; t71 = t124*t159; - t132 = t130*x; - t83 = t126*t132; - t156 = t71+t83; - t121 = RATIONAL(1.0,40.0); - t70 = t121*t158; - t128 = RATIONAL(-1.0,8.0); - t87 = t128*t133; - t155 = t70+t87; - t67 = t121*t159; - t129 = RATIONAL(1.0,8.0); - t86 = t129*t132; - t154 = t67+t86; - t85 = t128*t132; - t153 = t67+t85; - t117 = RATIONAL(1.0,24.0); - t77 = t117*t132; - t152 = t71+t77; - t127 = RATIONAL(-3.0,40.0); - t62 = t127*t158; - t151 = t62+t82; - t63 = t127*t159; - t150 = t63+t77; - t125 = RATIONAL(-1.0,40.0); - t64 = t125*t159; - t149 = t64+t86; - t79 = t117*t133; - t148 = t68+t79; - t147 = t63+t83; - t146 = t62+t79; - t69 = t125*t158; - t110 = RATIONAL(11.0,80.0); - t145 = t69+t110*t130; - t84 = t129*t133; - t144 = t70+t84; - t143 = t64+t85; + t153 = t71+t83; + t108 = RATIONAL(11.0,80.0); + t122 = RATIONAL(-1.0,40.0); + t69 = t122*t158; + t152 = t69+t108*t133; + t151 = RATIONAL(-1.0,100.0)*t160+RATIONAL(37.0,200.0); + t63 = t122*t159; + t150 = t63+t87; + t67 = t117*t158; + t86 = t129*t130; + t149 = t67+t86; + t65 = t126*t158; + t148 = t65+t77; + t85 = t128*t130; + t147 = t67+t85; + t79 = t116*t132; + t146 = t70+t79; + t82 = t125*t130; + t145 = t68+t82; + t144 = t63+t84; + t143 = t65+t82; + t142 = t71+t79; t123 = RATIONAL(-1.0,50.0); - t142 = RATIONAL(-13.0,100.0)*t160+t123; - t104 = RATIONAL(-17.0,200.0); - t141 = RATIONAL(-9.0,100.0)*t160+t104; - t140 = t123*t160+RATIONAL(11.0,50.0); - t139 = RATIONAL(-1.0,100.0)*t160+RATIONAL(37.0,200.0); - t138 = RATIONAL(-2.0,25.0)+RATIONAL(2.0,25.0)*t160; + t141 = t123+RATIONAL(-13.0,100.0)*t160; + t101 = RATIONAL(-17.0,200.0); + t140 = RATIONAL(-9.0,100.0)*t160+t101; + t139 = RATIONAL(11.0,50.0)+t123*t160; + t138 = RATIONAL(2.0,25.0)*t160+RATIONAL(-2.0,25.0); t137 = RATIONAL(7.0,100.0)*t160+RATIONAL(13.0,100.0); - t136 = t130+t131; - t122 = RATIONAL(9.0,80.0); - t120 = RATIONAL(7.0,80.0); + t136 = t133+t131; + t127 = RATIONAL(1.0,80.0); + t121 = RATIONAL(9.0,80.0); t119 = RATIONAL(-1.0,80.0); - t116 = RATIONAL(3.0,80.0); - t114 = RATIONAL(1.0,80.0); - t113 = RATIONAL(13.0,80.0); - t111 = RATIONAL(-23.0,80.0); - t108 = RATIONAL(-19.0,80.0); - t107 = RATIONAL(-17.0,80.0); - t105 = RATIONAL(-31.0,400.0); - t103 = RATIONAL(69.0,400.0); - t102 = RATIONAL(-57.0,400.0); - t101 = RATIONAL(-63.0,400.0); - t99 = RATIONAL(-49.0,400.0); - t98 = RATIONAL(43.0,1200.0); - t97 = RATIONAL(103.0,400.0); - t96 = RATIONAL(111.0,400.0); + t118 = RATIONAL(7.0,80.0); + t115 = RATIONAL(3.0,80.0); + t113 = RATIONAL(-19.0,80.0); + t112 = RATIONAL(-17.0,80.0); + t111 = RATIONAL(13.0,80.0); + t109 = RATIONAL(-23.0,80.0); + t106 = RATIONAL(-63.0,400.0); + t105 = RATIONAL(-49.0,400.0); + t103 = RATIONAL(-31.0,400.0); + t102 = RATIONAL(69.0,400.0); + t100 = RATIONAL(-57.0,400.0); + t98 = RATIONAL(117.0,400.0); + t97 = RATIONAL(43.0,1200.0); + t96 = RATIONAL(103.0,400.0); t95 = RATIONAL(37.0,1200.0); - t94 = RATIONAL(89.0,1200.0); - t93 = RATIONAL(117.0,400.0); + t94 = RATIONAL(111.0,400.0); + t93 = RATIONAL(89.0,1200.0); t92 = RATIONAL(-109.0,1200.0); t91 = RATIONAL(-223.0,1200.0); t90 = RATIONAL(-157.0,1200.0); - t89 = RATIONAL(-131.0,1200.0); - t88 = RATIONAL(-149.0,1200.0); - t81 = t122*t131; - t80 = t120*t130; - t78 = t120*t131; - t76 = t122*t130; - t75 = t113*t131; - t74 = t113*t130; - t73 = t110*t131; + t89 = RATIONAL(-149.0,1200.0); + t88 = RATIONAL(-131.0,1200.0); + t81 = t121*t133; + t80 = t118*t131; + t78 = t118*t133; + t76 = t121*t131; + t75 = t111*t131; + t74 = t111*t133; + t72 = t108*t131; coeff_I_m1_m1 = t75+RATIONAL(-23.0,200.0)+t74+t92*y+(t92+RATIONAL(6.0, -25.0)*y)*x+t147+t151; - coeff_I_0_m1 = t111*t130+t105*x+t73+t91*y+t140+t149+t157; - coeff_I_p1_m1 = t81+t96*x+t90*y+t80+t142+t153+t157; - coeff_I_p2_m1 = t78+t94*y+t89*x+t116*t130+t141+t151+t152; - coeff_I_m1_0 = t105*y+t91*x+t111*t131+t84+t140+t145+t156; - coeff_I_0_0 = RATIONAL(93.0,200.0)+t102*y+t136*RATIONAL(-21.0,80.0)+( -RATIONAL(-1.0,25.0)*y+t102)*x+t144+t154; - coeff_I_p1_0 = t108*t131+t101*y+t76+t93*x+t139+t143+t144; - coeff_I_p2_0 = t107*t131+t114*t130+t99*y+t84+t69+t98*x+t137+t150; - coeff_I_m1_p1 = t90*x+t96*y+t78+t76+t142+t155+t156; - coeff_I_0_p1 = t69+t101*x+t108*t130+t87+t93*y+t81+t139+t154; - coeff_I_p1_p1 = t97*y+t73+t104+t87+(RATIONAL(3.0,50.0)*y+t97)*x+t143+t145 +25.0)*y)*x+t143+t154; + coeff_I_0_m1 = t109*t131+t91*y+t86+t103*x+t139+t152+t153; + coeff_I_p1_m1 = t81+t94*x+t90*y+t80+t141+t147+t153; + coeff_I_p2_m1 = t78+t93*y+t88*x+t115*t131+t140+t154+t155; + coeff_I_m1_0 = t72+t109*t133+t91*x+t103*y+t139+t144+t145; + coeff_I_0_0 = RATIONAL(93.0,200.0)+t100*y+t136*RATIONAL(-21.0,80.0)+( +RATIONAL(-1.0,25.0)*y+t100)*x+t149+t156; + coeff_I_p1_0 = t69+t76+t98*x+t85+t113*t133+t106*y+t151+t156; + coeff_I_p2_0 = t97*x+t127*t131+t105*y+t112*t133+t137+t144+t148; + coeff_I_m1_p1 = t76+t90*x+t78+t94*y+t141+t145+t157; + coeff_I_0_p1 = t113*t131+t81+t106*x+t98*y+t149+t150+t151; + coeff_I_p1_p1 = t72+t101+t96*y+t85+(t96+RATIONAL(3.0,50.0)*y)*x+t150+t152 ; - coeff_I_p2_p1 = t95*x+t119*t130+t75+t103*y+t138+t150+t155; - coeff_I_m1_p2 = t116*t131+t94*x+t89*y+t80+t141+t147+t148; - coeff_I_0_p2 = t99*x+t98*y+t114*t131+t107*t130+t137+t146+t149; - coeff_I_p1_p2 = t103*x+t119*t131+t95*y+t74+t138+t146+t153; - coeff_I_p2_p2 = RATIONAL(7.0,200.0)+t88*y+t136*RATIONAL(-3.0,80.0)+(t88+ -RATIONAL(-3.0,50.0)*y)*x+t148+t152; + coeff_I_p2_p1 = t102*y+t74+t119*t131+t95*x+t138+t148+t157; + coeff_I_m1_p2 = t80+t88*y+t93*x+t115*t133+t140+t142+t143; + coeff_I_0_p2 = t127*t133+t86+t97*y+t105*x+t112*t131+t69+t137+t146; + coeff_I_p1_p2 = t95*y+t75+t119*t133+t102*x+t138+t146+t147; + coeff_I_p2_p2 = t89*y+RATIONAL(7.0,200.0)+t136*RATIONAL(-3.0,80.0)+(t89+ +RATIONAL(-3.0,50.0)*y)*x+t142+t155; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c index 4650764..5bfa6f4 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c @@ -1,95 +1,95 @@ fp t201, - t164, + t161, t200, - t163, - t199, - t183, + t184, t188, - t173, + t174, t187, - t179, + t180, + t199, + t164, t198, - t161, + t183, + t175, t197, t162, t196, - t185, - t176, - t195, t181, + t195, t194, + t163, t193, t186, - t177, + t176, t182, t192, - t180, - t191, - t184, + t185, t172, + t191, + t179, t190, t189, t178, - t175, - t170, - t168, + t177, + t169, + t167, t166, t165; t201 = x*y; - t164 = RATIONAL(-3.0,20.0)*t201; - t200 = RATIONAL(13.0,40.0)*x+t164; - t163 = RATIONAL(3.0,20.0)*t201; - t199 = RATIONAL(7.0,40.0)*x+t163; - t183 = RATIONAL(1.0,40.0); + t161 = RATIONAL(1.0,20.0)*t201; + t200 = t161+RATIONAL(9.0,40.0)*x; + t184 = RATIONAL(-1.0,40.0); t188 = y*y; - t173 = t183*t188; + t174 = t184*t188; t187 = x*x; - t179 = RATIONAL(-3.0,8.0)*t187; - t198 = t173+t179; - t161 = RATIONAL(1.0,20.0)*t201; - t197 = RATIONAL(9.0,40.0)*x+t161; + t180 = RATIONAL(-3.0,8.0)*t187; + t199 = t174+t180; + t164 = RATIONAL(-3.0,20.0)*t201; + t198 = t164+RATIONAL(13.0,40.0)*x; + t183 = RATIONAL(1.0,40.0); + t175 = t183*t188; + t197 = t175+t180; t162 = RATIONAL(-1.0,20.0)*t201; - t196 = RATIONAL(11.0,40.0)*x+t162; - t185 = RATIONAL(-1.0,40.0); - t176 = t185*t188; - t195 = t176+t179; + t196 = t162+RATIONAL(11.0,40.0)*x; t181 = RATIONAL(3.0,8.0)*t187; - t194 = t173+t181; - t193 = t176+t181; + t195 = t174+t181; + t194 = t175+t181; + t163 = RATIONAL(3.0,20.0)*t201; + t193 = RATIONAL(7.0,40.0)*x+t163; t186 = RATIONAL(-3.0,40.0); - t177 = t186*t188; + t176 = t186*t188; t182 = RATIONAL(-1.0,8.0)*t187; - t192 = t177+t182; - t180 = RATIONAL(1.0,8.0)*t187; - t191 = t177+t180; - t184 = RATIONAL(3.0,40.0); - t172 = t184*t188; - t190 = t172+t180; - t189 = t172+t182; + t192 = t176+t182; + t185 = RATIONAL(3.0,40.0); + t172 = t185*t188; + t191 = t172+t182; + t179 = RATIONAL(1.0,8.0)*t187; + t190 = t172+t179; + t189 = t176+t179; t178 = RATIONAL(2.0,25.0)*y; - t175 = RATIONAL(-1.0,50.0)*y; - t170 = RATIONAL(7.0,100.0)*y; - t168 = RATIONAL(-9.0,100.0)*y; - t166 = RATIONAL(-1.0,100.0)*y; + t177 = RATIONAL(-1.0,50.0)*y; + t169 = RATIONAL(-9.0,100.0)*y; + t167 = RATIONAL(-1.0,100.0)*y; + t166 = RATIONAL(7.0,100.0)*y; t165 = RATIONAL(-13.0,100.0)*y; - coeff_dx_m1_m1 = RATIONAL(-109.0,1200.0)+RATIONAL(6.0,25.0)*y+t192+t200; - coeff_dx_0_m1 = RATIONAL(-31.0,400.0)+RATIONAL(-23.0,40.0)*x+t175+t163+ -t193; - coeff_dx_p1_m1 = RATIONAL(111.0,400.0)+t165+t198+t199; - coeff_dx_p2_m1 = t168+t164+RATIONAL(-131.0,1200.0)+t184*x+t190; - coeff_dx_m1_0 = RATIONAL(-223.0,1200.0)+t175+t189+t196; - coeff_dx_0_0 = RATIONAL(-57.0,400.0)+RATIONAL(-21.0,40.0)*x+RATIONAL(-1.0 -,25.0)*y+t161+t194; - coeff_dx_p1_0 = RATIONAL(117.0,400.0)+t166+t195+t197; - coeff_dx_p2_0 = RATIONAL(43.0,1200.0)+t183*x+t162+t170+t191; - coeff_dx_m1_p1 = RATIONAL(-157.0,1200.0)+t165+t189+t197; - coeff_dx_0_p1 = RATIONAL(-63.0,400.0)+t166+t162+RATIONAL(-19.0,40.0)*x+ + coeff_dx_m1_m1 = RATIONAL(-109.0,1200.0)+RATIONAL(6.0,25.0)*y+t192+t198; + coeff_dx_0_m1 = RATIONAL(-31.0,400.0)+t177+RATIONAL(-23.0,40.0)*x+t163+ +t195; + coeff_dx_p1_m1 = RATIONAL(111.0,400.0)+t165+t193+t197; + coeff_dx_p2_m1 = t169+RATIONAL(-131.0,1200.0)+t164+t185*x+t190; + coeff_dx_m1_0 = RATIONAL(-223.0,1200.0)+t177+t191+t196; + coeff_dx_0_0 = RATIONAL(-1.0,25.0)*y+t161+RATIONAL(-21.0,40.0)*x+RATIONAL +(-57.0,400.0)+t194; + coeff_dx_p1_0 = t167+RATIONAL(117.0,400.0)+t199+t200; + coeff_dx_p2_0 = t183*x+t162+RATIONAL(43.0,1200.0)+t166+t189; + coeff_dx_m1_p1 = RATIONAL(-157.0,1200.0)+t165+t191+t200; + coeff_dx_0_p1 = RATIONAL(-19.0,40.0)*x+RATIONAL(-63.0,400.0)+t162+t167+ t194; - coeff_dx_p1_p1 = RATIONAL(3.0,50.0)*y+RATIONAL(103.0,400.0)+t195+t196; - coeff_dx_p2_p1 = t185*x+t161+RATIONAL(37.0,1200.0)+t178+t191; - coeff_dx_m1_p2 = t168+RATIONAL(89.0,1200.0)+t192+t199; - coeff_dx_0_p2 = t164+t170+RATIONAL(-17.0,40.0)*x+RATIONAL(-49.0,400.0)+ -t193; - coeff_dx_p1_p2 = RATIONAL(69.0,400.0)+t178+t198+t200; + coeff_dx_p1_p1 = RATIONAL(3.0,50.0)*y+RATIONAL(103.0,400.0)+t196+t199; + coeff_dx_p2_p1 = t161+t184*x+RATIONAL(37.0,1200.0)+t178+t189; + coeff_dx_m1_p2 = RATIONAL(89.0,1200.0)+t169+t192+t193; + coeff_dx_0_p2 = RATIONAL(-17.0,40.0)*x+t166+RATIONAL(-49.0,400.0)+t164+ +t195; + coeff_dx_p1_p2 = t178+RATIONAL(69.0,400.0)+t197+t198; coeff_dx_p2_p2 = t163+RATIONAL(-149.0,1200.0)+t186*x+RATIONAL(-3.0,50.0)* y+t190; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c index 395eceb..29b5494 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c @@ -1,40 +1,40 @@ -fp t243, +fp t244, t258, - t245, - t257, t246, + t257, + t245, t256, - t244, + t243, t255, t250, t249, t248, t247; - t243 = RATIONAL(1.0,20.0)*y; - t258 = t243+RATIONAL(9.0,40.0); - t245 = RATIONAL(3.0,20.0)*y; - t257 = t245+RATIONAL(7.0,40.0); - t246 = RATIONAL(-3.0,20.0)*y; - t256 = t246+RATIONAL(13.0,40.0); - t244 = RATIONAL(-1.0,20.0)*y; - t255 = t244+RATIONAL(11.0,40.0); - t250 = RATIONAL(-1.0,4.0)*x; - t249 = RATIONAL(3.0,4.0)*x; - t248 = RATIONAL(-3.0,4.0)*x; - t247 = RATIONAL(1.0,4.0)*x; - coeff_dxx_m1_m1 = t250+t256; - coeff_dxx_0_m1 = RATIONAL(-23.0,40.0)+t245+t249; - coeff_dxx_p1_m1 = t248+t257; - coeff_dxx_p2_m1 = t246+t247+RATIONAL(3.0,40.0); - coeff_dxx_m1_0 = t250+t255; - coeff_dxx_0_0 = t243+RATIONAL(-21.0,40.0)+t249; - coeff_dxx_p1_0 = t248+t258; - coeff_dxx_p2_0 = RATIONAL(1.0,40.0)+t244+t247; - coeff_dxx_m1_p1 = t250+t258; - coeff_dxx_0_p1 = t244+RATIONAL(-19.0,40.0)+t249; - coeff_dxx_p1_p1 = t248+t255; - coeff_dxx_p2_p1 = t243+t247+RATIONAL(-1.0,40.0); - coeff_dxx_m1_p2 = t250+t257; - coeff_dxx_0_p2 = RATIONAL(-17.0,40.0)+t246+t249; - coeff_dxx_p1_p2 = t248+t256; - coeff_dxx_p2_p2 = RATIONAL(-3.0,40.0)+t247+t245; + t244 = RATIONAL(-3.0,20.0)*y; + t258 = t244+RATIONAL(13.0,40.0); + t246 = RATIONAL(-1.0,20.0)*y; + t257 = RATIONAL(11.0,40.0)+t246; + t245 = RATIONAL(1.0,20.0)*y; + t256 = RATIONAL(9.0,40.0)+t245; + t243 = RATIONAL(3.0,20.0)*y; + t255 = t243+RATIONAL(7.0,40.0); + t250 = RATIONAL(1.0,4.0)*x; + t249 = RATIONAL(-1.0,4.0)*x; + t248 = RATIONAL(3.0,4.0)*x; + t247 = RATIONAL(-3.0,4.0)*x; + coeff_dxx_m1_m1 = t249+t258; + coeff_dxx_0_m1 = t248+RATIONAL(-23.0,40.0)+t243; + coeff_dxx_p1_m1 = t247+t255; + coeff_dxx_p2_m1 = RATIONAL(3.0,40.0)+t250+t244; + coeff_dxx_m1_0 = t249+t257; + coeff_dxx_0_0 = RATIONAL(-21.0,40.0)+t245+t248; + coeff_dxx_p1_0 = t247+t256; + coeff_dxx_p2_0 = RATIONAL(1.0,40.0)+t246+t250; + coeff_dxx_m1_p1 = t249+t256; + coeff_dxx_0_p1 = t248+RATIONAL(-19.0,40.0)+t246; + coeff_dxx_p1_p1 = t247+t257; + coeff_dxx_p2_p1 = RATIONAL(-1.0,40.0)+t250+t245; + coeff_dxx_m1_p2 = t249+t255; + coeff_dxx_0_p2 = t244+RATIONAL(-17.0,40.0)+t248; + coeff_dxx_p1_p2 = t247+t258; + coeff_dxx_p2_p2 = RATIONAL(-3.0,40.0)+t250+t243; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c index b48e5bc..20e2106 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c @@ -17,36 +17,36 @@ fp t276, t260, t259; t276 = RATIONAL(-1.0,50.0); - t275 = RATIONAL(-3.0,20.0); - t274 = RATIONAL(3.0,20.0); - t273 = RATIONAL(2.0,25.0); + t275 = RATIONAL(2.0,25.0); + t274 = RATIONAL(-3.0,20.0); + t273 = RATIONAL(3.0,20.0); t272 = RATIONAL(-1.0,20.0); t271 = RATIONAL(1.0,20.0); t270 = RATIONAL(-9.0,100.0); - t269 = RATIONAL(7.0,100.0); - t268 = RATIONAL(-1.0,100.0); + t269 = RATIONAL(-1.0,100.0); + t268 = RATIONAL(7.0,100.0); t267 = RATIONAL(-13.0,100.0); - t266 = t275*y; - t265 = t274*y; - t264 = t275*x; + t266 = t272*y; + t265 = t271*y; + t264 = t273*x; t263 = t274*x; t262 = t271*x; t261 = t272*x; - t260 = t272*y; - t259 = t271*y; - coeff_dxy_m1_m1 = t266+t264+RATIONAL(6.0,25.0); - coeff_dxy_0_m1 = t263+t260+t276; - coeff_dxy_p1_m1 = t263+t259+t267; - coeff_dxy_p2_m1 = t270+t264+t265; - coeff_dxy_m1_0 = t276+t261+t265; - coeff_dxy_0_0 = RATIONAL(-1.0,25.0)+t259+t262; - coeff_dxy_p1_0 = t268+t262+t260; - coeff_dxy_p2_0 = t261+t269+t266; - coeff_dxy_m1_p1 = t267+t262+t265; - coeff_dxy_0_p1 = t261+t268+t259; - coeff_dxy_p1_p1 = RATIONAL(3.0,50.0)+t261+t260; - coeff_dxy_p2_p1 = t266+t262+t273; - coeff_dxy_m1_p2 = t266+t263+t270; - coeff_dxy_0_p2 = t269+t264+t260; - coeff_dxy_p1_p2 = t264+t259+t273; - coeff_dxy_p2_p2 = t263+t265+RATIONAL(-3.0,50.0); + t260 = t274*y; + t259 = t273*y; + coeff_dxy_m1_m1 = t260+RATIONAL(6.0,25.0)+t263; + coeff_dxy_0_m1 = t276+t264+t266; + coeff_dxy_p1_m1 = t264+t267+t265; + coeff_dxy_p2_m1 = t263+t270+t259; + coeff_dxy_m1_0 = t276+t259+t261; + coeff_dxy_0_0 = t262+RATIONAL(-1.0,25.0)+t265; + coeff_dxy_p1_0 = t269+t262+t266; + coeff_dxy_p2_0 = t260+t268+t261; + coeff_dxy_m1_p1 = t262+t267+t259; + coeff_dxy_0_p1 = t269+t265+t261; + coeff_dxy_p1_p1 = RATIONAL(3.0,50.0)+t266+t261; + coeff_dxy_p2_p1 = t260+t275+t262; + coeff_dxy_m1_p2 = t264+t260+t270; + coeff_dxy_0_p2 = t268+t263+t266; + coeff_dxy_p1_p2 = t263+t275+t265; + coeff_dxy_p2_p2 = t259+RATIONAL(-3.0,50.0)+t264; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c index 678f775..c766157 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c @@ -1,95 +1,95 @@ fp t242, + t225, t228, - t226, t212, t229, t220, t241, - t224, - t213, + t203, t240, - t202, + t204, t239, - t223, + t205, t238, - t204, + t224, + t213, t237, + t202, t236, - t227, - t217, - t221, + t226, + t215, + t222, t235, - t205, + t227, + t218, t234, - t203, + t221, t233, - t225, - t218, t232, - t222, + t223, t231, t230, t219, - t214, + t217, t209, t208, t207, t206; t242 = x*y; + t225 = RATIONAL(-1.0,40.0); t228 = x*x; - t226 = RATIONAL(-1.0,40.0); - t212 = t226*t228; + t212 = t225*t228; t229 = y*y; - t220 = RATIONAL(3.0,8.0)*t229; + t220 = RATIONAL(-3.0,8.0)*t229; t241 = t212+t220; + t203 = RATIONAL(1.0,20.0)*t242; + t240 = t203+RATIONAL(9.0,40.0)*y; + t204 = RATIONAL(-1.0,20.0)*t242; + t239 = RATIONAL(11.0,40.0)*y+t204; + t205 = RATIONAL(-3.0,20.0)*t242; + t238 = t205+RATIONAL(13.0,40.0)*y; t224 = RATIONAL(1.0,40.0); t213 = t224*t228; - t240 = t213+t220; + t237 = t213+t220; t202 = RATIONAL(3.0,20.0)*t242; - t239 = t202+RATIONAL(7.0,40.0)*y; - t223 = RATIONAL(-3.0,8.0)*t229; - t238 = t213+t223; - t204 = RATIONAL(1.0,20.0)*t242; - t237 = t204+RATIONAL(9.0,40.0)*y; - t236 = t212+t223; + t236 = RATIONAL(7.0,40.0)*y+t202; + t226 = RATIONAL(3.0,40.0); + t215 = t226*t228; + t222 = RATIONAL(-1.0,8.0)*t229; + t235 = t215+t222; t227 = RATIONAL(-3.0,40.0); - t217 = t227*t228; - t221 = RATIONAL(-1.0,8.0)*t229; - t235 = t217+t221; - t205 = RATIONAL(-3.0,20.0)*t242; - t234 = t205+RATIONAL(13.0,40.0)*y; - t203 = RATIONAL(-1.0,20.0)*t242; - t233 = t203+RATIONAL(11.0,40.0)*y; - t225 = RATIONAL(3.0,40.0); - t218 = t225*t228; - t232 = t218+t221; - t222 = RATIONAL(1.0,8.0)*t229; - t231 = t218+t222; - t230 = t217+t222; + t218 = t227*t228; + t234 = t218+t222; + t221 = RATIONAL(3.0,8.0)*t229; + t233 = t212+t221; + t232 = t213+t221; + t223 = RATIONAL(1.0,8.0)*t229; + t231 = t218+t223; + t230 = t215+t223; t219 = RATIONAL(2.0,25.0)*x; - t214 = RATIONAL(-1.0,50.0)*x; - t209 = RATIONAL(-1.0,100.0)*x; - t208 = RATIONAL(-9.0,100.0)*x; - t207 = RATIONAL(7.0,100.0)*x; + t217 = RATIONAL(-1.0,50.0)*x; + t209 = RATIONAL(-9.0,100.0)*x; + t208 = RATIONAL(7.0,100.0)*x; + t207 = RATIONAL(-1.0,100.0)*x; t206 = RATIONAL(-13.0,100.0)*x; - coeff_dy_m1_m1 = RATIONAL(-109.0,1200.0)+RATIONAL(6.0,25.0)*x+t234+t235; - coeff_dy_0_m1 = t214+RATIONAL(-223.0,1200.0)+t232+t233; - coeff_dy_p1_m1 = t206+RATIONAL(-157.0,1200.0)+t232+t237; - coeff_dy_p2_m1 = t208+RATIONAL(89.0,1200.0)+t235+t239; - coeff_dy_m1_0 = RATIONAL(-23.0,40.0)*y+t214+RATIONAL(-31.0,400.0)+t202+ -t241; - coeff_dy_0_0 = RATIONAL(-21.0,40.0)*y+RATIONAL(-57.0,400.0)+t204+RATIONAL -(-1.0,25.0)*x+t240; - coeff_dy_p1_0 = t209+RATIONAL(-63.0,400.0)+t203+RATIONAL(-19.0,40.0)*y+ -t240; - coeff_dy_p2_0 = t207+RATIONAL(-49.0,400.0)+RATIONAL(-17.0,40.0)*y+t205+ -t241; - coeff_dy_m1_p1 = RATIONAL(111.0,400.0)+t206+t238+t239; - coeff_dy_0_p1 = RATIONAL(117.0,400.0)+t209+t236+t237; - coeff_dy_p1_p1 = RATIONAL(103.0,400.0)+RATIONAL(3.0,50.0)*x+t233+t236; - coeff_dy_p2_p1 = t219+RATIONAL(69.0,400.0)+t234+t238; - coeff_dy_m1_p2 = t225*y+t208+RATIONAL(-131.0,1200.0)+t205+t231; - coeff_dy_0_p2 = t207+t203+t224*y+RATIONAL(43.0,1200.0)+t230; - coeff_dy_p1_p2 = t219+t226*y+RATIONAL(37.0,1200.0)+t204+t230; - coeff_dy_p2_p2 = t202+RATIONAL(-3.0,50.0)*x+RATIONAL(-149.0,1200.0)+t227* -y+t231; + coeff_dy_m1_m1 = RATIONAL(6.0,25.0)*x+RATIONAL(-109.0,1200.0)+t234+t238; + coeff_dy_0_m1 = RATIONAL(-223.0,1200.0)+t217+t235+t239; + coeff_dy_p1_m1 = t206+RATIONAL(-157.0,1200.0)+t235+t240; + coeff_dy_p2_m1 = t209+RATIONAL(89.0,1200.0)+t234+t236; + coeff_dy_m1_0 = RATIONAL(-23.0,40.0)*y+t202+RATIONAL(-31.0,400.0)+t217+ +t233; + coeff_dy_0_0 = RATIONAL(-1.0,25.0)*x+RATIONAL(-57.0,400.0)+t203+RATIONAL( +-21.0,40.0)*y+t232; + coeff_dy_p1_0 = RATIONAL(-63.0,400.0)+t207+t204+RATIONAL(-19.0,40.0)*y+ +t232; + coeff_dy_p2_0 = t205+RATIONAL(-17.0,40.0)*y+RATIONAL(-49.0,400.0)+t208+ +t233; + coeff_dy_m1_p1 = RATIONAL(111.0,400.0)+t206+t236+t237; + coeff_dy_0_p1 = RATIONAL(117.0,400.0)+t207+t240+t241; + coeff_dy_p1_p1 = RATIONAL(3.0,50.0)*x+RATIONAL(103.0,400.0)+t239+t241; + coeff_dy_p2_p1 = t219+RATIONAL(69.0,400.0)+t237+t238; + coeff_dy_m1_p2 = t209+RATIONAL(-131.0,1200.0)+t205+t226*y+t230; + coeff_dy_0_p2 = t204+t208+t224*y+RATIONAL(43.0,1200.0)+t231; + coeff_dy_p1_p2 = t225*y+t219+t203+RATIONAL(37.0,1200.0)+t231; + coeff_dy_p2_p2 = RATIONAL(-149.0,1200.0)+t202+RATIONAL(-3.0,50.0)*x+t227* +y+t230; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c index dea29f0..b449df9 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c @@ -1,40 +1,40 @@ -fp t279, +fp t277, t292, - t280, + t279, t291, - t278, + t280, t290, - t277, + t278, t289, t284, t283, t282, t281; + t277 = RATIONAL(3.0,20.0)*x; + t292 = t277+RATIONAL(7.0,40.0); t279 = RATIONAL(-1.0,20.0)*x; - t292 = t279+RATIONAL(11.0,40.0); + t291 = t279+RATIONAL(11.0,40.0); t280 = RATIONAL(-3.0,20.0)*x; - t291 = RATIONAL(13.0,40.0)+t280; + t290 = t280+RATIONAL(13.0,40.0); t278 = RATIONAL(1.0,20.0)*x; - t290 = t278+RATIONAL(9.0,40.0); - t277 = RATIONAL(3.0,20.0)*x; - t289 = t277+RATIONAL(7.0,40.0); + t289 = t278+RATIONAL(9.0,40.0); t284 = RATIONAL(-1.0,4.0)*y; - t283 = RATIONAL(1.0,4.0)*y; - t282 = RATIONAL(3.0,4.0)*y; - t281 = RATIONAL(-3.0,4.0)*y; - coeff_dyy_m1_m1 = t284+t291; - coeff_dyy_0_m1 = t284+t292; - coeff_dyy_p1_m1 = t284+t290; - coeff_dyy_p2_m1 = t284+t289; - coeff_dyy_m1_0 = t282+RATIONAL(-23.0,40.0)+t277; - coeff_dyy_0_0 = RATIONAL(-21.0,40.0)+t278+t282; - coeff_dyy_p1_0 = RATIONAL(-19.0,40.0)+t282+t279; - coeff_dyy_p2_0 = t280+t282+RATIONAL(-17.0,40.0); - coeff_dyy_m1_p1 = t281+t289; - coeff_dyy_0_p1 = t281+t290; - coeff_dyy_p1_p1 = t281+t292; - coeff_dyy_p2_p1 = t281+t291; - coeff_dyy_m1_p2 = t280+t283+RATIONAL(3.0,40.0); - coeff_dyy_0_p2 = t283+RATIONAL(1.0,40.0)+t279; - coeff_dyy_p1_p2 = t278+RATIONAL(-1.0,40.0)+t283; - coeff_dyy_p2_p2 = t277+t283+RATIONAL(-3.0,40.0); + t283 = RATIONAL(-3.0,4.0)*y; + t282 = RATIONAL(1.0,4.0)*y; + t281 = RATIONAL(3.0,4.0)*y; + coeff_dyy_m1_m1 = t284+t290; + coeff_dyy_0_m1 = t284+t291; + coeff_dyy_p1_m1 = t284+t289; + coeff_dyy_p2_m1 = t284+t292; + coeff_dyy_m1_0 = RATIONAL(-23.0,40.0)+t281+t277; + coeff_dyy_0_0 = t281+RATIONAL(-21.0,40.0)+t278; + coeff_dyy_p1_0 = RATIONAL(-19.0,40.0)+t279+t281; + coeff_dyy_p2_0 = t281+RATIONAL(-17.0,40.0)+t280; + coeff_dyy_m1_p1 = t283+t292; + coeff_dyy_0_p1 = t283+t289; + coeff_dyy_p1_p1 = t283+t291; + coeff_dyy_p2_p1 = t283+t290; + coeff_dyy_m1_p2 = t282+RATIONAL(3.0,40.0)+t280; + coeff_dyy_0_p2 = t282+RATIONAL(1.0,40.0)+t279; + coeff_dyy_p1_p2 = RATIONAL(-1.0,40.0)+t278+t282; + coeff_dyy_p2_p2 = t277+t282+RATIONAL(-3.0,40.0); diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c index 49655f0..cfe391a 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c @@ -1,239 +1,243 @@ -fp t413, +fp t417, t418, - t466, - t403, - t417, - t358, + t468, t414, + t411, + t361, + t416, t373, - t438, - t444, - t465, - t405, - t364, + t435, + t446, + t467, + t410, t367, - t439, + t368, + t437, + t466, + t465, + t441, t464, - t436, t463, - t410, - t359, - t442, - t409, - t318, - t392, - t462, - t411, - t360, - t441, + t395, + t445, t400, - t298, - t461, + t403, t366, - t443, - t395, + t462, + t397, + t404, + t362, t401, + t461, + t389, + t409, + t316, t460, - t440, - t296, - t399, + t402, + t312, + t384, t459, - t372, + t412, + t443, + t407, + t309, t458, - t333, + t295, + t399, + t442, t457, - t396, - t412, - t327, + t406, + t405, + t319, t456, - t407, - t331, - t391, + t444, + t370, t455, - t402, - t313, - t397, + t315, t454, t325, t453, - t324, + t358, + t398, + t294, t452, - t311, t451, + t413, + t364, t450, - t387, - t398, + t357, + t380, + t381, t449, + t332, t415, - t371, - t416, - t357, - t383, - t382, - t448, - t317, t356, + t448, t447, - t362, - t446, - t445, - t437, - t435, - t393, - t404, - t369, - t376, + t440, + t360, + t439, + t438, + t436, + t390, + t408, + t363, + t392, t434, - t365, - t384, t386, + t374, t433, + t372, t432, t431, + t317, t430, t429, + t326, t428, t427, - t406, + t385, t426, t425, - t408, + t388, t424, t423, t394, - t389, - t388, - t385; - t413 = x*x; + t393, + t391, + t387; + t417 = x*x; t418 = y*y; - t466 = t418+t413; - t403 = RATIONAL(-1.0,60.0); - t417 = y*t418; - t358 = t403*t417; - t414 = t413*x; - t373 = t403*t414; - t438 = y*t373+x*t358; - t444 = x*y; - t465 = t438+RATIONAL(37.0,300.0)*t444; - t405 = RATIONAL(1.0,60.0); - t364 = t405*t414; - t367 = t405*t417; - t439 = y*t364+x*t367; - t464 = t439+RATIONAL(-37.0,300.0)*t444; - t436 = t413*t418; - t463 = RATIONAL(-99.0,1225.0)+RATIONAL(-1.0,98.0)*t436; - t410 = RATIONAL(-1.0,30.0); - t359 = t410*t417; - t442 = y*t413; - t409 = RATIONAL(1.0,70.0); - t318 = t409*t442; - t392 = RATIONAL(11.0,105.0); - t462 = t359+t318+t392*y; - t411 = RATIONAL(1.0,30.0); - t360 = t411*t417; - t441 = y*t414; - t400 = RATIONAL(-1.0,120.0); - t298 = RATIONAL(-21.0,200.0)*t444; - t461 = x*t360+t298+t400*t441; - t366 = t410*t414; - t443 = x*t418; - t395 = RATIONAL(31.0,210.0); + t468 = t418+t417; + t414 = t417*x; + t411 = RATIONAL(1.0,60.0); + t361 = t411*t414; + t416 = y*t418; + t373 = t411*t416; + t435 = x*t373+y*t361; + t446 = x*y; + t467 = RATIONAL(-11.0,150.0)*t446+t435; + t410 = RATIONAL(-1.0,60.0); + t367 = t410*t416; + t368 = t410*t414; + t437 = x*t367+y*t368; + t466 = RATIONAL(37.0,300.0)*t446+t437; + t465 = t437+RATIONAL(11.0,150.0)*t446; + t441 = t417*t418; + t464 = RATIONAL(-1.0,98.0)*t441+RATIONAL(-99.0,1225.0); + t463 = RATIONAL(-37.0,300.0)*t446+t435; + t395 = RATIONAL(-31.0,210.0); + t445 = x*t418; + t400 = RATIONAL(1.0,140.0); + t403 = RATIONAL(1.0,30.0); + t366 = t403*t414; + t462 = t400*t445+t395*x+t366; + t397 = RATIONAL(31.0,210.0); + t404 = RATIONAL(-1.0,30.0); + t362 = t404*t414; t401 = RATIONAL(-1.0,140.0); - t460 = t395*x+t366+t401*t443; - t440 = x*t417; - t296 = RATIONAL(21.0,200.0)*t444; + t461 = t362+t401*t445+t397*x; + t389 = RATIONAL(31.0,420.0); + t409 = RATIONAL(-1.0,35.0); + t316 = t409*t445; + t460 = t389*x+t316+t368; + t402 = RATIONAL(1.0,35.0); + t312 = t402*t445; + t384 = RATIONAL(-31.0,420.0); + t459 = t312+t384*x+t361; + t412 = RATIONAL(1.0,84.0); + t443 = y*t417; + t407 = RATIONAL(-1.0,70.0); + t309 = t407*t443; + t458 = t309+t412*y+t373; + t295 = RATIONAL(21.0,200.0)*t446; t399 = RATIONAL(1.0,120.0); - t459 = t296+y*t366+t399*t440; - t372 = t411*t414; - t458 = t400*t440+y*t372+t298; - t333 = t409*t443; - t457 = t366+t333+t392*x; - t396 = RATIONAL(-11.0,105.0); - t412 = RATIONAL(-1.0,70.0); - t327 = t412*t443; - t456 = t396*x+t372+t327; - t407 = RATIONAL(1.0,35.0); - t331 = t407*t442; - t391 = RATIONAL(-31.0,420.0); - t455 = t331+t391*y+t367; - t402 = RATIONAL(-1.0,35.0); - t313 = t402*t442; - t397 = RATIONAL(31.0,420.0); - t454 = t358+t397*y+t313; - t325 = t402*t443; - t453 = t397*x+t325+t373; - t324 = t407*t443; - t452 = t364+t324+t391*x; - t311 = t412*t442; - t451 = t360+t311+t396*y; - t450 = x*t359+t296+t399*t441; - t387 = RATIONAL(-31.0,210.0); - t398 = RATIONAL(1.0,140.0); - t449 = t372+t387*x+t398*t443; - t415 = t413*t413; - t371 = t410*t415; - t416 = t418*t418; - t357 = t399*t416; - t383 = RATIONAL(71.0,5880.0); - t382 = RATIONAL(113.0,735.0); - t448 = t371+t382*t413+t383*t418+t357+t463; - t317 = RATIONAL(1.0,49.0)*t436; + t442 = x*t416; + t457 = t295+t399*t442+y*t362; + t406 = RATIONAL(-1.0,84.0); + t405 = RATIONAL(1.0,70.0); + t319 = t405*t443; + t456 = t406*y+t319+t367; + t444 = t414*y; + t370 = t404*t416; + t455 = t295+t399*t444+x*t370; + t315 = t405*t445; + t454 = t406*x+t315+t368; + t325 = t407*t445; + t453 = t361+t412*x+t325; + t358 = t403*t416; + t398 = RATIONAL(-1.0,120.0); + t294 = RATIONAL(-21.0,200.0)*t446; + t452 = t398*t444+t294+x*t358; + t451 = y*t366+t398*t442+t294; + t413 = t417*t417; + t364 = t404*t413; + t450 = RATIONAL(-24.0,1225.0)+t364+RATIONAL(1.0,196.0)*t441+t468*RATIONAL +(181.0,1470.0); + t357 = t399*t413; + t380 = RATIONAL(71.0,5880.0); + t381 = RATIONAL(113.0,735.0); + t449 = t357+t380*t417+t381*t418+t464; + t332 = RATIONAL(1.0,49.0)*t441; + t415 = t418*t418; t356 = t399*t415; - t447 = t317+RATIONAL(51.0,1225.0)+t356+t357+t466*RATIONAL(-289.0,5880.0); - t362 = t410*t416; - t446 = t382*t418+t356+t383*t413+t362+t463; - t445 = t362+RATIONAL(1.0,196.0)*t436+t371+RATIONAL(-24.0,1225.0)+t466* -RATIONAL(181.0,1470.0); - t437 = RATIONAL(-1.0,49.0)*t436+RATIONAL(96.0,1225.0); - t435 = RATIONAL(1.0,98.0)*t436+RATIONAL(246.0,1225.0); - t393 = RATIONAL(-41.0,196.0); - t404 = RATIONAL(1.0,20.0); - t369 = t404*t415; - t376 = RATIONAL(191.0,5880.0); - t434 = t376*t418+t369+t357+t393*t413+t437; - t365 = t404*t416; - t384 = RATIONAL(-53.0,196.0); - t386 = RATIONAL(83.0,735.0); - t433 = t365+t371+t384*t418+t386*t413+t435; - t432 = t365+t376*t413+t356+t393*t418+t437; - t431 = t369+t384*t413+t362+t386*t418+t435; - t430 = RATIONAL(-11.0,150.0)*t444+t439+t447; - t429 = t359+t395*y+t401*t442+t445; - t428 = t360+t387*y+t398*t442+t445; - t427 = RATIONAL(11.0,150.0)*t444+t438+t447; - t406 = RATIONAL(1.0,84.0); - t426 = t367+t311+t406*y+t448; - t425 = t364+t406*x+t327+t446; - t408 = RATIONAL(-1.0,84.0); - t424 = t358+t318+t408*y+t448; - t423 = t333+t408*x+t373+t446; - t394 = RATIONAL(-17.0,420.0); - t389 = RATIONAL(17.0,105.0); - t388 = RATIONAL(17.0,420.0); - t385 = RATIONAL(-17.0,105.0); - coeff_I_m2_m2 = t430+t453+t454; - coeff_I_m1_m2 = t424+t456+t459; - coeff_I_0_m2 = t331+t394*y+t358+t434; - coeff_I_p1_m2 = t424+t457+t458; - coeff_I_p2_m2 = t427+t452+t454; - coeff_I_m2_m1 = t423+t450+t451; - coeff_I_m1_m1 = t428+t449+t465; - coeff_I_0_m1 = t360+t385*y+t318+t431; - coeff_I_p1_m1 = t428+t460+t464; - coeff_I_p2_m1 = t425+t451+t461; - coeff_I_m2_0 = t373+t394*x+t324+t432; - coeff_I_m1_0 = t385*x+t333+t372+t433; - coeff_I_0_0 = t317+t365+t369+RATIONAL(541.0,1225.0)+t466*RATIONAL(-57.0, + t448 = RATIONAL(51.0,1225.0)+t357+t332+t356+t468*RATIONAL(-289.0,5880.0); + t447 = t380*t418+t364+t381*t417+t356+t464; + t440 = RATIONAL(1.0,98.0)*t441+RATIONAL(246.0,1225.0); + t360 = t404*t415; + t439 = t360+t370; + t438 = t358+t360; + t436 = RATIONAL(-1.0,49.0)*t441+RATIONAL(96.0,1225.0); + t390 = RATIONAL(83.0,735.0); + t408 = RATIONAL(1.0,20.0); + t363 = t408*t413; + t392 = RATIONAL(-53.0,196.0); + t434 = t360+t390*t418+t363+t392*t417+t440; + t386 = RATIONAL(-41.0,196.0); + t374 = RATIONAL(191.0,5880.0); + t433 = t386*t417+t374*t418+t363+t356+t436; + t372 = t408*t415; + t432 = t357+t374*t417+t386*t418+t372+t436; + t431 = t364+t392*t418+t390*t417+t372+t440; + t317 = t402*t443; + t430 = t384*y+t317+t373+t448; + t429 = t395*y+t400*t443+t438+t450; + t326 = t409*t443; + t428 = t389*y+t367+t326+t448; + t427 = t397*y+t401*t443+t439+t450; + t385 = RATIONAL(11.0,105.0); + t426 = t385*x+t315+t362+t447; + t425 = t385*y+t319+t439+t449; + t388 = RATIONAL(-11.0,105.0); + t424 = t388*y+t309+t438+t449; + t423 = t388*x+t366+t325+t447; + t394 = RATIONAL(17.0,420.0); + t393 = RATIONAL(17.0,105.0); + t391 = RATIONAL(-17.0,105.0); + t387 = RATIONAL(-17.0,420.0); + coeff_I_m2_m2 = t428+t460+t467; + coeff_I_m1_m2 = t423+t456+t457; + coeff_I_0_m2 = t317+t367+t387*y+t433; + coeff_I_p1_m2 = t426+t451+t456; + coeff_I_p2_m2 = t428+t459+t465; + coeff_I_m2_m1 = t424+t454+t455; + coeff_I_m1_m1 = t429+t462+t466; + coeff_I_0_m1 = t319+t391*y+t358+t434; + coeff_I_p1_m1 = t429+t461+t463; + coeff_I_p2_m1 = t424+t452+t453; + coeff_I_m2_0 = t368+t312+t387*x+t432; + coeff_I_m1_0 = t366+t315+t391*x+t431; + coeff_I_0_0 = RATIONAL(541.0,1225.0)+t363+t372+t332+t468*RATIONAL(-57.0, 196.0); - coeff_I_p1_0 = t389*x+t327+t366+t433; - coeff_I_p2_0 = t325+t364+t388*x+t432; - coeff_I_m2_p1 = t423+t461+t462; - coeff_I_m1_p1 = t429+t449+t464; - coeff_I_0_p1 = t359+t389*y+t311+t431; - coeff_I_p1_p1 = t429+t460+t465; - coeff_I_p2_p1 = t425+t450+t462; - coeff_I_m2_p2 = t427+t453+t455; - coeff_I_m1_p2 = t426+t456+t458; - coeff_I_0_p2 = t313+t388*y+t367+t434; - coeff_I_p1_p2 = t426+t457+t459; - coeff_I_p2_p2 = t430+t452+t455; + coeff_I_p1_0 = t325+t362+t393*x+t431; + coeff_I_p2_0 = t394*x+t316+t361+t432; + coeff_I_m2_p1 = t425+t452+t454; + coeff_I_m1_p1 = t427+t462+t463; + coeff_I_0_p1 = t309+t370+t393*y+t434; + coeff_I_p1_p1 = t427+t461+t466; + coeff_I_p2_p1 = t425+t453+t455; + coeff_I_m2_p2 = t430+t460+t465; + coeff_I_m1_p2 = t423+t451+t458; + coeff_I_0_p2 = t394*y+t326+t373+t433; + coeff_I_p1_p2 = t426+t457+t458; + coeff_I_p2_p2 = t430+t459+t467; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c index 83ba345..e56ee74 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c @@ -1,48 +1,54 @@ -fp t531, - t534, - t535, - t513, - t570, - t508, - t569, +fp t538, + t507, + t539, t537, + t559, + t573, + t560, + t572, + t514, + t554, + t571, + t526, t536, - t502, - t552, + t516, + t570, + t531, + t517, + t569, + t481, t568, - t512, - t555, + t561, + t494, t567, - t556, + t495, t566, - t509, - t533, - t510, - t565, - t530, + t525, t506, - t507, + t565, + t510, t564, + t511, t563, - t471, t562, - t561, - t523, - t511, - t560, - t528, - t504, - t559, + t475, + t527, + t509, t558, + t513, t557, + t529, + t512, + t556, + t555, + t518, t478, - t554, t553, - t516, - t476, + t479, + t552, t551, - t483, t550, + t490, t549, t548, t547, @@ -51,97 +57,99 @@ fp t531, t544, t543, t542, - t541, - t540, - t484, + t486, + t482, t477, - t475, - t473, - t469, - t467; - t531 = RATIONAL(1.0,30.0); - t534 = x*x; - t535 = t534*x; - t513 = t531*t535; - t570 = t513+RATIONAL(71.0,2940.0)*x; - t508 = RATIONAL(-2.0,15.0)*t535; - t569 = t508+RATIONAL(226.0,735.0)*x; - t537 = y*y; - t536 = y*t537; - t502 = RATIONAL(1.0,20.0)*t534; - t552 = RATIONAL(1.0,60.0)*t536+y*t502; - t568 = t552+RATIONAL(-37.0,300.0)*y; - t512 = RATIONAL(-1.0,20.0)*t534; - t555 = y*t512+RATIONAL(-1.0,60.0)*t536; - t567 = t555+RATIONAL(37.0,300.0)*y; - t556 = x*t537; - t566 = t508+RATIONAL(181.0,735.0)*x+RATIONAL(1.0,98.0)*t556; - t509 = RATIONAL(1.0,10.0)*t534; - t533 = RATIONAL(-1.0,70.0); - t510 = t533*t537; - t565 = t509+t510+RATIONAL(-11.0,105.0); - t530 = RATIONAL(1.0,70.0); - t506 = t530*t537; - t507 = RATIONAL(-1.0,10.0)*t534; - t564 = t506+t507+RATIONAL(11.0,105.0); - t563 = t510+RATIONAL(1.0,84.0)+t502; - t471 = RATIONAL(2.0,49.0)*t556; - t562 = t513+t471+RATIONAL(-289.0,2940.0)*x; - t561 = t506+t512+RATIONAL(-1.0,84.0); - t523 = RATIONAL(-1.0,35.0); - t511 = t523*t537; - t560 = t511+t512+RATIONAL(31.0,420.0); - t528 = RATIONAL(1.0,35.0); - t504 = t528*t537; - t559 = RATIONAL(-31.0,420.0)+t502+t504; - t558 = x*y; - t557 = y*t534; - t478 = RATIONAL(-1.0,49.0)*t556; - t554 = RATIONAL(-21.0,200.0)*y+t478; - t553 = RATIONAL(21.0,200.0)*y+t478; - t516 = RATIONAL(1.0,5.0)*t535; - t476 = RATIONAL(-2.0,49.0)*t556; - t551 = t516+t476+RATIONAL(-41.0,98.0)*x; - t483 = RATIONAL(1.0,49.0)*t556; - t550 = t516+t483+RATIONAL(-53.0,98.0)*x; - t549 = t508+RATIONAL(166.0,735.0)*x+t483; - t548 = t513+t476+RATIONAL(191.0,2940.0)*x; - t547 = t509+RATIONAL(-31.0,210.0)+RATIONAL(1.0,140.0)*t537+t566; - t546 = RATIONAL(11.0,150.0)*y+t555+t562; - t545 = y*t509+RATIONAL(-1.0,120.0)*t536+t554+t569; - t544 = t507+RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t537+t566; - t543 = y*t507+RATIONAL(1.0,120.0)*t536+t553+t569; - t542 = t531*t536+RATIONAL(-1.0,40.0)*t557+t554+t570; - t541 = RATIONAL(-11.0,150.0)*y+t552+t562; - t540 = RATIONAL(1.0,40.0)*t557+RATIONAL(-1.0,30.0)*t536+t553+t570; - t484 = t530*t558; - t477 = t528*t558; - t475 = RATIONAL(2.0,35.0)*t558; - t473 = t533*t558; - t469 = RATIONAL(-2.0,35.0)*t558; - t467 = t523*t558; - coeff_dx_m2_m2 = t469+t541+t560; - coeff_dx_m1_m2 = t477+t543+t565; - coeff_dx_0_m2 = t475+t551; - coeff_dx_p1_m2 = t477+t545+t564; - coeff_dx_p2_m2 = t469+t546+t559; - coeff_dx_m2_m1 = t467+t540+t561; - coeff_dx_m1_m1 = t484+t547+t567; - coeff_dx_0_m1 = t477+t550; - coeff_dx_p1_m1 = t484+t544+t568; - coeff_dx_p2_m1 = t467+t542+t563; - coeff_dx_m2_0 = t504+t512+RATIONAL(-17.0,420.0)+t548; - coeff_dx_m1_0 = RATIONAL(-17.0,105.0)+t509+t506+t549; - coeff_dx_0_0 = RATIONAL(-57.0,98.0)*x+t516+t471; - coeff_dx_p1_0 = t510+RATIONAL(17.0,105.0)+t507+t549; - coeff_dx_p2_0 = t511+RATIONAL(17.0,420.0)+t502+t548; - coeff_dx_m2_p1 = t477+t542+t561; - coeff_dx_m1_p1 = t473+t547+t568; - coeff_dx_0_p1 = t467+t550; - coeff_dx_p1_p1 = t473+t544+t567; - coeff_dx_p2_p1 = t477+t540+t563; - coeff_dx_m2_p2 = t475+t546+t560; - coeff_dx_m1_p2 = t467+t545+t565; - coeff_dx_0_p2 = t469+t551; - coeff_dx_p1_p2 = t467+t543+t564; - coeff_dx_p2_p2 = t475+t541+t559; + t474, + t471, + t469; + t538 = x*x; + t507 = RATIONAL(1.0,20.0)*t538; + t539 = y*y; + t537 = y*t539; + t559 = y*t507+RATIONAL(1.0,60.0)*t537; + t573 = RATIONAL(-37.0,300.0)*y+t559; + t560 = x*t539; + t572 = RATIONAL(1.0,98.0)*t560+RATIONAL(181.0,735.0)*x; + t514 = RATIONAL(-1.0,20.0)*t538; + t554 = y*t514+RATIONAL(-1.0,60.0)*t537; + t571 = RATIONAL(37.0,300.0)*y+t554; + t526 = RATIONAL(1.0,30.0); + t536 = x*t538; + t516 = t526*t536; + t570 = RATIONAL(71.0,2940.0)*x+t516; + t531 = RATIONAL(-1.0,35.0); + t517 = t531*t539; + t569 = RATIONAL(31.0,420.0)+t514+t517; + t481 = RATIONAL(2.0,49.0)*t560; + t568 = RATIONAL(-289.0,2940.0)*x+t481+t516; + t561 = y*t538; + t494 = RATIONAL(-21.0,200.0)*y; + t567 = RATIONAL(-1.0,40.0)*t561+t494+t526*t537; + t495 = RATIONAL(21.0,200.0)*y; + t566 = RATIONAL(1.0,40.0)*t561+t495+RATIONAL(-1.0,30.0)*t537; + t525 = RATIONAL(1.0,35.0); + t506 = t525*t539; + t565 = RATIONAL(-31.0,420.0)+t507+t506; + t510 = RATIONAL(-1.0,10.0)*t538; + t564 = y*t510+t495+RATIONAL(1.0,120.0)*t537; + t511 = RATIONAL(1.0,10.0)*t538; + t563 = y*t511+t494+RATIONAL(-1.0,120.0)*t537; + t562 = x*y; + t475 = RATIONAL(-1.0,49.0)*t560; + t527 = RATIONAL(1.0,70.0); + t509 = t527*t539; + t558 = t475+t509; + t513 = RATIONAL(-2.0,15.0)*t536; + t557 = t510+t513; + t529 = RATIONAL(-1.0,70.0); + t512 = t529*t539; + t556 = t475+t512; + t555 = t511+t513; + t518 = RATIONAL(1.0,5.0)*t536; + t478 = RATIONAL(-2.0,49.0)*t560; + t553 = t518+RATIONAL(-41.0,98.0)*x+t478; + t479 = RATIONAL(1.0,49.0)*t560; + t552 = t518+RATIONAL(-53.0,98.0)*x+t479; + t551 = RATIONAL(191.0,2940.0)*x+t478+t516; + t550 = RATIONAL(166.0,735.0)*x+t479+t513; + t490 = RATIONAL(226.0,735.0)*x; + t549 = RATIONAL(11.0,105.0)+t490+t557+t558; + t548 = t490+RATIONAL(-11.0,105.0)+t555+t556; + t547 = RATIONAL(-31.0,210.0)+RATIONAL(1.0,140.0)*t539+t555+t572; + t546 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t539+t557+t572; + t545 = RATIONAL(1.0,84.0)+t507+t556+t570; + t544 = RATIONAL(-11.0,150.0)*y+t559+t568; + t543 = RATIONAL(-1.0,84.0)+t514+t558+t570; + t542 = RATIONAL(11.0,150.0)*y+t554+t568; + t486 = t525*t562; + t482 = RATIONAL(-2.0,35.0)*t562; + t477 = RATIONAL(2.0,35.0)*t562; + t474 = t531*t562; + t471 = t527*t562; + t469 = t529*t562; + coeff_dx_m2_m2 = t482+t544+t569; + coeff_dx_m1_m2 = t486+t548+t564; + coeff_dx_0_m2 = t477+t553; + coeff_dx_p1_m2 = t486+t549+t563; + coeff_dx_p2_m2 = t482+t542+t565; + coeff_dx_m2_m1 = t474+t543+t566; + coeff_dx_m1_m1 = t471+t547+t571; + coeff_dx_0_m1 = t486+t552; + coeff_dx_p1_m1 = t471+t546+t573; + coeff_dx_p2_m1 = t474+t545+t567; + coeff_dx_m2_0 = t514+RATIONAL(-17.0,420.0)+t506+t551; + coeff_dx_m1_0 = t511+RATIONAL(-17.0,105.0)+t509+t550; + coeff_dx_0_0 = t481+t518+RATIONAL(-57.0,98.0)*x; + coeff_dx_p1_0 = RATIONAL(17.0,105.0)+t510+t512+t550; + coeff_dx_p2_0 = t507+RATIONAL(17.0,420.0)+t517+t551; + coeff_dx_m2_p1 = t486+t543+t567; + coeff_dx_m1_p1 = t469+t547+t573; + coeff_dx_0_p1 = t474+t552; + coeff_dx_p1_p1 = t469+t546+t571; + coeff_dx_p2_p1 = t486+t545+t566; + coeff_dx_m2_p2 = t477+t542+t569; + coeff_dx_m1_p2 = t474+t548+t563; + coeff_dx_0_p2 = t482+t553; + coeff_dx_p1_p2 = t474+t549+t564; + coeff_dx_p2_p2 = t477+t544+t565; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c index 385d8a7..3ac1e8c 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c @@ -1,111 +1,107 @@ -fp t712, - t692, - t707, - t711, - t693, - t729, +fp t713, + t683, + t708, + t712, t687, + t730, + t696, + t729, t728, - t697, + t690, t727, t726, - t725, + t697, t686, - t694, - t695, t698, - t689, + t691, + t693, + t725, + t699, t724, - t681, + t682, t723, - t696, t722, + t695, t721, - t688, t720, - t691, + t689, t719, - t685, t718, t717, t716, - t690, + t678, t715, + t676, t714, - t713, - t684, - t683, + t694, + t685, + t681, t680, t679, - t678, - t677, - t676, - t675; - t712 = y*y; - t692 = RATIONAL(2.0,49.0)*t712; - t707 = RATIONAL(1.0,10.0); - t711 = x*x; - t693 = t707*t711; - t729 = t692+t693+RATIONAL(-289.0,2940.0); - t687 = RATIONAL(-1.0,49.0)*t712; - t728 = t693+RATIONAL(71.0,2940.0)+t687; - t697 = RATIONAL(-2.0,5.0)*t711; - t727 = t697+RATIONAL(226.0,735.0)+t687; - t726 = t697+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t712; - t725 = x*y; - t686 = t707*x; - t694 = RATIONAL(-1.0,10.0)*x; - t695 = RATIONAL(-1.0,5.0)*x; - t698 = RATIONAL(1.0,5.0)*x; - t689 = RATIONAL(-2.0,49.0)*t712; - t724 = t693+RATIONAL(191.0,2940.0)+t689; - t681 = RATIONAL(1.0,49.0)*t712; - t723 = t697+RATIONAL(166.0,735.0)+t681; - t696 = RATIONAL(3.0,5.0)*t711; - t722 = t696+RATIONAL(-53.0,98.0)+t681; - t721 = t696+RATIONAL(-41.0,98.0)+t689; - t688 = RATIONAL(2.0,35.0)*y; - t720 = t688+t729; - t691 = RATIONAL(-2.0,35.0)*y; - t719 = t691+t729; - t685 = RATIONAL(-1.0,35.0)*y; - t718 = t685+t727; - t717 = t686+t728; - t716 = t695+t726; - t690 = RATIONAL(1.0,35.0)*y; - t715 = t690+t727; - t714 = t698+t726; - t713 = t694+t728; - t684 = RATIONAL(1.0,70.0)*y; - t683 = RATIONAL(-1.0,70.0)*y; - t680 = y*t695; - t679 = y*t698; - t678 = RATIONAL(-1.0,20.0)*t725; - t677 = y*t686; - t676 = y*t694; - t675 = RATIONAL(1.0,20.0)*t725; - coeff_dxx_m2_m2 = t677+t694+t719; - coeff_dxx_m1_m2 = t680+t698+t715; - coeff_dxx_0_m2 = t688+t721; - coeff_dxx_p1_m2 = t679+t695+t715; - coeff_dxx_p2_m2 = t686+t676+t719; - coeff_dxx_m2_m1 = t685+t675+t713; - coeff_dxx_m1_m1 = t676+t684+t714; - coeff_dxx_0_m1 = t690+t722; - coeff_dxx_p1_m1 = t684+t677+t716; - coeff_dxx_p2_m1 = t685+t678+t717; - coeff_dxx_m2_0 = t694+t724; - coeff_dxx_m1_0 = t698+t723; - coeff_dxx_0_0 = t692+RATIONAL(-57.0,98.0)+t696; - coeff_dxx_p1_0 = t695+t723; - coeff_dxx_p2_0 = t686+t724; - coeff_dxx_m2_p1 = t690+t678+t713; - coeff_dxx_m1_p1 = t677+t683+t714; - coeff_dxx_0_p1 = t685+t722; - coeff_dxx_p1_p1 = t683+t676+t716; - coeff_dxx_p2_p1 = t690+t675+t717; - coeff_dxx_m2_p2 = t676+t694+t720; - coeff_dxx_m1_p2 = t698+t679+t718; - coeff_dxx_0_p2 = t691+t721; - coeff_dxx_p1_p2 = t695+t680+t718; - coeff_dxx_p2_p2 = t686+t677+t720; + t677; + t713 = y*y; + t683 = RATIONAL(-1.0,49.0)*t713; + t708 = RATIONAL(1.0,10.0); + t712 = x*x; + t687 = t708*t712; + t730 = RATIONAL(71.0,2940.0)+t683+t687; + t696 = RATIONAL(-2.0,5.0)*t712; + t729 = t696+RATIONAL(226.0,735.0)+t683; + t728 = t696+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t713; + t690 = RATIONAL(2.0,49.0)*t713; + t727 = t690+RATIONAL(-289.0,2940.0)+t687; + t726 = x*y; + t697 = RATIONAL(1.0,5.0)*x; + t686 = RATIONAL(-1.0,10.0)*x; + t698 = RATIONAL(-1.0,5.0)*x; + t691 = t708*x; + t693 = RATIONAL(-2.0,49.0)*t713; + t725 = t693+RATIONAL(191.0,2940.0)+t687; + t699 = RATIONAL(3.0,5.0)*t712; + t724 = t693+t699+RATIONAL(-41.0,98.0); + t682 = RATIONAL(1.0,49.0)*t713; + t723 = t699+RATIONAL(-53.0,98.0)+t682; + t722 = t696+RATIONAL(166.0,735.0)+t682; + t695 = RATIONAL(1.0,35.0)*y; + t721 = t695+t730; + t720 = RATIONAL(-1.0,70.0)*y+t728; + t689 = RATIONAL(-1.0,35.0)*y; + t719 = t689+t730; + t718 = t698+t729; + t717 = RATIONAL(1.0,70.0)*y+t728; + t716 = t697+t729; + t678 = y*t686; + t715 = t678+t727; + t676 = y*t691; + t714 = t676+t727; + t694 = RATIONAL(2.0,35.0)*y; + t685 = RATIONAL(-2.0,35.0)*y; + t681 = y*t697; + t680 = y*t698; + t679 = RATIONAL(1.0,20.0)*t726; + t677 = RATIONAL(-1.0,20.0)*t726; + coeff_dxx_m2_m2 = t686+t685+t714; + coeff_dxx_m1_m2 = t680+t695+t716; + coeff_dxx_0_m2 = t694+t724; + coeff_dxx_p1_m2 = t695+t681+t718; + coeff_dxx_p2_m2 = t691+t685+t715; + coeff_dxx_m2_m1 = t679+t686+t719; + coeff_dxx_m1_m1 = t697+t678+t717; + coeff_dxx_0_m1 = t695+t723; + coeff_dxx_p1_m1 = t698+t676+t717; + coeff_dxx_p2_m1 = t691+t677+t719; + coeff_dxx_m2_0 = t686+t725; + coeff_dxx_m1_0 = t697+t722; + coeff_dxx_0_0 = RATIONAL(-57.0,98.0)+t690+t699; + coeff_dxx_p1_0 = t698+t722; + coeff_dxx_p2_0 = t691+t725; + coeff_dxx_m2_p1 = t677+t686+t721; + coeff_dxx_m1_p1 = t697+t676+t720; + coeff_dxx_0_p1 = t689+t723; + coeff_dxx_p1_p1 = t678+t698+t720; + coeff_dxx_p2_p1 = t691+t679+t721; + coeff_dxx_m2_p2 = t694+t686+t715; + coeff_dxx_m1_p2 = t689+t681+t716; + coeff_dxx_0_p2 = t685+t724; + coeff_dxx_p1_p2 = t680+t689+t718; + coeff_dxx_p2_p2 = t694+t691+t714; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c index 9388906..0505d86 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c @@ -1,119 +1,119 @@ -fp t778, - t777, - t796, - t766, - t735, +fp t779, + t778, + t797, t772, - t755, + t774, + t796, t795, - t751, - t748, + t770, + t739, + t766, + t752, t794, + t742, + t759, t793, t792, - t791, t733, + t791, t790, t789, t788, t787, - t765, - t770, + t735, t786, - t734, + t732, t785, - t767, - t771, + t773, + t777, t784, t783, t782, - t730, t781, t780, - t779, t776, - t775, - t773, - t769, - t756, - t754, + t771, + t768, + t767, + t758, + t757, + t755, t750, - t749, - t746, - t745, t743, t738, - t732, + t737, + t736, + t734, t731; - t778 = y*y; - t777 = x*x; - t796 = t778+t777; - t766 = RATIONAL(-1.0,35.0); - t735 = t766*x; - t772 = RATIONAL(1.0,35.0); - t755 = t772*y; - t795 = t735+t755; - t751 = t766*y; - t748 = t772*x; - t794 = t751+t748; - t793 = t735+t751; - t792 = t755+t748; - t791 = x*y; - t733 = RATIONAL(-2.0,49.0)*t791; - t790 = t733+RATIONAL(-21.0,200.0); - t789 = t733+RATIONAL(21.0,200.0); - t788 = t796*RATIONAL(-1.0,20.0); - t787 = t796*RATIONAL(1.0,20.0); - t765 = RATIONAL(-1.0,40.0); - t770 = RATIONAL(1.0,10.0); - t786 = t765*t777+t770*t778+t790; - t734 = RATIONAL(4.0,49.0)*t791; - t785 = t734+RATIONAL(11.0,150.0)+t788; - t767 = RATIONAL(1.0,40.0); - t771 = RATIONAL(-1.0,10.0); - t784 = t771*t778+t767*t777+t789; - t783 = t771*t777+t767*t778+t789; - t782 = t734+RATIONAL(-11.0,150.0)+t787; - t730 = RATIONAL(1.0,49.0)*t791; - t781 = RATIONAL(37.0,300.0)+t730+t788; - t780 = RATIONAL(-37.0,300.0)+t730+t787; - t779 = t765*t778+t770*t777+t790; + t779 = y*y; + t778 = x*x; + t797 = t778+t779; + t772 = RATIONAL(1.0,10.0); + t774 = RATIONAL(-1.0,40.0); + t796 = t772*t778+t774*t779; + t795 = t774*t778+t772*t779; + t770 = RATIONAL(-1.0,35.0); + t739 = t770*y; + t766 = RATIONAL(1.0,35.0); + t752 = t766*x; + t794 = t739+t752; + t742 = t770*x; + t759 = t766*y; + t793 = t742+t759; + t792 = x*y; + t733 = RATIONAL(-2.0,49.0)*t792; + t791 = t733+RATIONAL(-21.0,200.0); + t790 = t797*RATIONAL(-1.0,20.0); + t789 = t797*RATIONAL(1.0,20.0); + t788 = t733+RATIONAL(21.0,200.0); + t787 = t752+t759+t791; + t735 = RATIONAL(4.0,49.0)*t792; + t786 = t735+RATIONAL(11.0,150.0)+t790; + t732 = RATIONAL(1.0,49.0)*t792; + t785 = t732+RATIONAL(37.0,300.0)+t790; + t773 = RATIONAL(-1.0,10.0); + t777 = RATIONAL(1.0,40.0); + t784 = t773*t779+t777*t778+t788; + t783 = t735+RATIONAL(-11.0,150.0)+t789; + t782 = t732+RATIONAL(-37.0,300.0)+t789; + t781 = t773*t778+t777*t779+t788; + t780 = t739+t742+t791; t776 = RATIONAL(-2.0,35.0); - t775 = RATIONAL(-1.0,70.0); - t773 = RATIONAL(1.0,70.0); - t769 = RATIONAL(2.0,35.0); - t756 = t776*y; - t754 = t769*y; - t750 = t776*x; - t749 = t775*x; - t746 = t773*y; - t745 = t775*y; - t743 = t769*x; - t738 = t773*x; - t732 = RATIONAL(-4.0,49.0)*t791; - t731 = RATIONAL(2.0,49.0)*t791; - coeff_dxy_m2_m2 = t756+t750+t782; - coeff_dxy_m1_m2 = t783+t794; - coeff_dxy_0_m2 = t743+t732; - coeff_dxy_p1_m2 = t779+t792; - coeff_dxy_p2_m2 = t750+t754+t785; - coeff_dxy_m2_m1 = t784+t795; - coeff_dxy_m1_m1 = t738+t746+t781; - coeff_dxy_0_m1 = t748+t731; - coeff_dxy_p1_m1 = t745+t738+t780; - coeff_dxy_p2_m1 = t786+t793; - coeff_dxy_m2_0 = t754+t732; - coeff_dxy_m1_0 = t755+t731; - coeff_dxy_0_0 = t734; - coeff_dxy_p1_0 = t751+t731; - coeff_dxy_p2_0 = t732+t756; - coeff_dxy_m2_p1 = t786+t792; - coeff_dxy_m1_p1 = t749+t746+t780; - coeff_dxy_0_p1 = t735+t731; - coeff_dxy_p1_p1 = t749+t745+t781; + t771 = RATIONAL(2.0,35.0); + t768 = RATIONAL(-1.0,70.0); + t767 = RATIONAL(1.0,70.0); + t758 = t771*y; + t757 = t768*y; + t755 = t776*x; + t750 = t771*x; + t743 = t767*x; + t738 = t767*y; + t737 = t776*y; + t736 = t768*x; + t734 = RATIONAL(-4.0,49.0)*t792; + t731 = RATIONAL(2.0,49.0)*t792; + coeff_dxy_m2_m2 = t755+t737+t783; + coeff_dxy_m1_m2 = t781+t794; + coeff_dxy_0_m2 = t750+t734; + coeff_dxy_p1_m2 = t787+t796; + coeff_dxy_p2_m2 = t758+t755+t786; + coeff_dxy_m2_m1 = t784+t793; + coeff_dxy_m1_m1 = t738+t743+t785; + coeff_dxy_0_m1 = t752+t731; + coeff_dxy_p1_m1 = t757+t743+t782; + coeff_dxy_p2_m1 = t780+t795; + coeff_dxy_m2_0 = t734+t758; + coeff_dxy_m1_0 = t731+t759; + coeff_dxy_0_0 = t735; + coeff_dxy_p1_0 = t731+t739; + coeff_dxy_p2_0 = t734+t737; + coeff_dxy_m2_p1 = t787+t795; + coeff_dxy_m1_p1 = t738+t736+t782; + coeff_dxy_0_p1 = t731+t742; + coeff_dxy_p1_p1 = t757+t736+t785; coeff_dxy_p2_p1 = t784+t794; - coeff_dxy_m2_p2 = t743+t756+t785; - coeff_dxy_m1_p2 = t779+t793; - coeff_dxy_0_p2 = t750+t732; - coeff_dxy_p1_p2 = t783+t795; - coeff_dxy_p2_p2 = t754+t743+t782; + coeff_dxy_m2_p2 = t750+t737+t786; + coeff_dxy_m1_p2 = t780+t796; + coeff_dxy_0_p2 = t755+t734; + coeff_dxy_p1_p2 = t781+t793; + coeff_dxy_p2_p2 = t758+t750+t783; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c index 03c68c3..f4baa3a 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c @@ -1,50 +1,51 @@ -fp t641, - t611, - t638, - t639, - t657, +fp t644, + t609, + t630, + t643, + t616, + t675, + t636, + t621, + t622, t674, - t608, - t656, - t673, t661, - t589, + t601, + t641, + t673, + t631, + t600, t672, + t613, t671, - t635, - t640, - t609, + t662, + t642, + t611, t670, - t574, + t587, t669, - t614, - t598, + t619, + t590, t668, - t632, - t615, - t667, t634, - t619, + t612, + t667, t666, - t617, - t596, + t610, t665, - t627, - t606, + t632, + t618, t664, - t637, - t616, t663, - t662, t660, - t618, t659, + t574, t658, t586, + t657, + t623, + t656, t655, - t578, t654, - t620, t653, t652, t651, @@ -52,100 +53,95 @@ fp t641, t649, t648, t647, - t646, - t645, - t644, - t588, + t592, + t584, t583, t582, - t579, - t576, - t572; - t641 = y*y; - t611 = RATIONAL(1.0,20.0)*t641; - t638 = x*x; - t639 = t638*x; - t657 = RATIONAL(1.0,60.0)*t639+x*t611; - t674 = t657+RATIONAL(-37.0,300.0)*x; - t608 = RATIONAL(-1.0,20.0)*t641; - t656 = RATIONAL(-1.0,60.0)*t639+x*t608; - t673 = t656+RATIONAL(37.0,300.0)*x; - t661 = y*t638; - t589 = RATIONAL(-1.0,49.0)*t661; - t672 = t589+RATIONAL(226.0,735.0)*y; - t671 = RATIONAL(1.0,98.0)*t661+RATIONAL(181.0,735.0)*y; - t635 = RATIONAL(1.0,30.0); - t640 = y*t641; - t609 = t635*t640; - t670 = t589+t609+RATIONAL(71.0,2940.0)*y; - t574 = RATIONAL(2.0,49.0)*t661; - t669 = RATIONAL(-289.0,2940.0)*y+t574+t609; - t614 = RATIONAL(-1.0,10.0)*t641; - t598 = RATIONAL(21.0,200.0)*x; - t668 = x*t614+RATIONAL(1.0,120.0)*t639+t598; - t632 = RATIONAL(1.0,35.0); - t615 = t632*t638; - t667 = t611+RATIONAL(-31.0,420.0)+t615; - t634 = RATIONAL(1.0,70.0); - t619 = t634*t638; - t666 = t619+RATIONAL(-1.0,84.0)+t608; - t617 = RATIONAL(1.0,10.0)*t641; - t596 = RATIONAL(-21.0,200.0)*x; - t665 = RATIONAL(-1.0,120.0)*t639+x*t617+t596; - t627 = RATIONAL(-1.0,35.0); - t606 = t627*t638; - t664 = RATIONAL(31.0,420.0)+t606+t608; - t637 = RATIONAL(-1.0,70.0); - t616 = t637*t638; - t663 = t611+t616+RATIONAL(1.0,84.0); - t662 = x*y; - t660 = x*t641; - t618 = RATIONAL(-2.0,15.0)*t640; - t659 = t617+t618; - t658 = t618+t614; - t586 = RATIONAL(-2.0,49.0)*t661; - t655 = t586+t609+RATIONAL(191.0,2940.0)*y; - t578 = RATIONAL(1.0,49.0)*t661; - t654 = t618+RATIONAL(166.0,735.0)*y+t578; - t620 = RATIONAL(1.0,5.0)*t640; - t653 = RATIONAL(-53.0,98.0)*y+t620+t578; - t652 = t620+RATIONAL(-41.0,98.0)*y+t586; - t651 = RATIONAL(1.0,140.0)*t638+RATIONAL(-31.0,210.0)+t659+t671; - t650 = t619+RATIONAL(11.0,105.0)+t658+t672; - t649 = RATIONAL(-11.0,150.0)*x+t657+t669; - t648 = RATIONAL(11.0,150.0)*x+t656+t669; - t647 = t616+RATIONAL(-11.0,105.0)+t659+t672; - t646 = t635*t639+RATIONAL(-1.0,40.0)*t660+t596+t670; - t645 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t638+t658+t671; - t644 = RATIONAL(-1.0,30.0)*t639+RATIONAL(1.0,40.0)*t660+t598+t670; - t588 = t627*t662; - t583 = t632*t662; - t582 = RATIONAL(2.0,35.0)*t662; - t579 = t634*t662; - t576 = t637*t662; - t572 = RATIONAL(-2.0,35.0)*t662; - coeff_dy_m2_m2 = t572+t649+t664; - coeff_dy_m1_m2 = t588+t644+t666; - coeff_dy_0_m2 = RATIONAL(-17.0,420.0)+t615+t608+t655; - coeff_dy_p1_m2 = t583+t646+t666; - coeff_dy_p2_m2 = t582+t648+t664; - coeff_dy_m2_m1 = t583+t647+t668; - coeff_dy_m1_m1 = t579+t651+t673; - coeff_dy_0_m1 = RATIONAL(-17.0,105.0)+t617+t619+t654; - coeff_dy_p1_m1 = t576+t651+t674; - coeff_dy_p2_m1 = t588+t647+t665; - coeff_dy_m2_0 = t582+t652; - coeff_dy_m1_0 = t583+t653; - coeff_dy_0_0 = t574+RATIONAL(-57.0,98.0)*y+t620; - coeff_dy_p1_0 = t588+t653; - coeff_dy_p2_0 = t572+t652; - coeff_dy_m2_p1 = t583+t650+t665; - coeff_dy_m1_p1 = t579+t645+t674; - coeff_dy_0_p1 = RATIONAL(17.0,105.0)+t616+t614+t654; - coeff_dy_p1_p1 = t576+t645+t673; - coeff_dy_p2_p1 = t588+t650+t668; - coeff_dy_m2_p2 = t572+t648+t667; - coeff_dy_m1_p2 = t588+t646+t663; - coeff_dy_0_p2 = t611+RATIONAL(17.0,420.0)+t606+t655; - coeff_dy_p1_p2 = t583+t644+t663; - coeff_dy_p2_p2 = t582+t649+t667; + t578, + t576; + t644 = y*y; + t609 = RATIONAL(1.0,20.0)*t644; + t630 = RATIONAL(1.0,35.0); + t643 = x*x; + t616 = t630*t643; + t675 = t609+t616+RATIONAL(-31.0,420.0); + t636 = RATIONAL(-1.0,35.0); + t621 = t636*t643; + t622 = RATIONAL(-1.0,20.0)*t644; + t674 = t621+t622+RATIONAL(31.0,420.0); + t661 = x*t644; + t601 = RATIONAL(21.0,200.0)*x; + t641 = x*t643; + t673 = RATIONAL(1.0,40.0)*t661+t601+RATIONAL(-1.0,30.0)*t641; + t631 = RATIONAL(1.0,30.0); + t600 = RATIONAL(-21.0,200.0)*x; + t672 = t600+RATIONAL(-1.0,40.0)*t661+t631*t641; + t613 = RATIONAL(1.0,10.0)*t644; + t671 = RATIONAL(1.0,140.0)*t643+t613+RATIONAL(-31.0,210.0); + t662 = y*t643; + t642 = y*t644; + t611 = RATIONAL(-2.0,15.0)*t642; + t670 = RATIONAL(181.0,735.0)*y+t611+RATIONAL(1.0,98.0)*t662; + t587 = RATIONAL(-1.0,49.0)*t662; + t669 = RATIONAL(226.0,735.0)*y+t587+t611; + t619 = t631*t642; + t590 = RATIONAL(2.0,49.0)*t662; + t668 = RATIONAL(-289.0,2940.0)*y+t619+t590; + t634 = RATIONAL(-1.0,70.0); + t612 = t634*t643; + t667 = t613+t612+RATIONAL(-11.0,105.0); + t666 = RATIONAL(71.0,2940.0)*y+t587+t619; + t610 = RATIONAL(-1.0,10.0)*t644; + t665 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t643+t610; + t632 = RATIONAL(1.0,70.0); + t618 = t632*t643; + t664 = t610+t618+RATIONAL(11.0,105.0); + t663 = x*y; + t660 = x*t622+RATIONAL(-1.0,60.0)*t641; + t659 = x*t609+RATIONAL(1.0,60.0)*t641; + t574 = RATIONAL(1.0,49.0)*t662; + t658 = RATIONAL(166.0,735.0)*y+t611+t574; + t586 = RATIONAL(-2.0,49.0)*t662; + t657 = RATIONAL(191.0,2940.0)*y+t619+t586; + t623 = RATIONAL(1.0,5.0)*t642; + t656 = RATIONAL(-41.0,98.0)*y+t623+t586; + t655 = RATIONAL(-53.0,98.0)*y+t623+t574; + t654 = RATIONAL(37.0,300.0)*x+t660+t670; + t653 = RATIONAL(11.0,150.0)*x+t660+t668; + t652 = RATIONAL(-37.0,300.0)*x+t659+t670; + t651 = RATIONAL(-1.0,84.0)+t618+t622+t666; + t650 = x*t610+t601+RATIONAL(1.0,120.0)*t641+t669; + t649 = RATIONAL(1.0,84.0)+t609+t612+t666; + t648 = RATIONAL(-11.0,150.0)*x+t659+t668; + t647 = t600+RATIONAL(-1.0,120.0)*t641+x*t613+t669; + t592 = RATIONAL(-2.0,35.0)*t663; + t584 = t630*t663; + t583 = t636*t663; + t582 = RATIONAL(2.0,35.0)*t663; + t578 = t634*t663; + t576 = t632*t663; + coeff_dy_m2_m2 = t592+t648+t674; + coeff_dy_m1_m2 = t583+t651+t673; + coeff_dy_0_m2 = t616+RATIONAL(-17.0,420.0)+t622+t657; + coeff_dy_p1_m2 = t584+t651+t672; + coeff_dy_p2_m2 = t582+t653+t674; + coeff_dy_m2_m1 = t584+t650+t667; + coeff_dy_m1_m1 = t576+t654+t671; + coeff_dy_0_m1 = t618+t613+RATIONAL(-17.0,105.0)+t658; + coeff_dy_p1_m1 = t578+t652+t671; + coeff_dy_p2_m1 = t583+t647+t667; + coeff_dy_m2_0 = t582+t656; + coeff_dy_m1_0 = t584+t655; + coeff_dy_0_0 = RATIONAL(-57.0,98.0)*y+t590+t623; + coeff_dy_p1_0 = t583+t655; + coeff_dy_p2_0 = t592+t656; + coeff_dy_m2_p1 = t584+t647+t664; + coeff_dy_m1_p1 = t576+t652+t665; + coeff_dy_0_p1 = t612+t610+RATIONAL(17.0,105.0)+t658; + coeff_dy_p1_p1 = t578+t654+t665; + coeff_dy_p2_p1 = t583+t650+t664; + coeff_dy_m2_p2 = t592+t653+t675; + coeff_dy_m1_p2 = t583+t649+t672; + coeff_dy_0_p2 = t621+RATIONAL(17.0,420.0)+t609+t657; + coeff_dy_p1_p2 = t584+t649+t673; + coeff_dy_p2_p2 = t582+t648+t675; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c index cc0ccef..9e72985 100644 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c @@ -1,115 +1,115 @@ -fp t833, - t808, +fp t834, + t813, + t854, t853, - t807, + t830, + t835, + t814, t852, - t834, - t818, + t816, t851, t850, - t849, - t829, t815, - t817, - t806, + t808, + t818, t820, - t804, + t821, + t849, t848, + t810, t847, - t812, + t804, t846, t819, t845, - t803, t844, t843, - t798, t842, t841, + t801, t840, + t798, t839, t838, t837, - t800, t836, - t835, - t816, - t814, + t817, + t812, t811, - t810, t809, + t806, t805, + t803, t802, - t801, - t799, - t797; - t833 = x*x; - t808 = RATIONAL(-1.0,49.0)*t833; - t853 = t808+RATIONAL(71.0,2940.0); - t807 = RATIONAL(2.0,49.0)*t833; - t852 = t807+RATIONAL(-289.0,2940.0); - t834 = y*y; - t818 = RATIONAL(-2.0,5.0)*t834; - t851 = t818+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t833; - t850 = t818+t808+RATIONAL(226.0,735.0); - t849 = x*y; - t829 = RATIONAL(1.0,10.0); - t815 = t829*y; - t817 = RATIONAL(-1.0,5.0)*y; - t806 = RATIONAL(-1.0,10.0)*y; + t800, + t799; + t834 = x*x; + t813 = RATIONAL(-1.0,49.0)*t834; + t854 = t813+RATIONAL(226.0,735.0); + t853 = RATIONAL(1.0,98.0)*t834+RATIONAL(181.0,735.0); + t830 = RATIONAL(1.0,10.0); + t835 = y*y; + t814 = t830*t835; + t852 = t813+t814+RATIONAL(71.0,2940.0); + t816 = RATIONAL(2.0,49.0)*t834; + t851 = t814+t816+RATIONAL(-289.0,2940.0); + t850 = x*y; + t815 = RATIONAL(-1.0,10.0)*y; + t808 = t830*y; + t818 = RATIONAL(-1.0,5.0)*y; t820 = RATIONAL(1.0,5.0)*y; - t804 = t829*t834; - t848 = t815+t804; - t847 = t806+t804; - t812 = RATIONAL(1.0,49.0)*t833; - t846 = t818+RATIONAL(166.0,735.0)+t812; - t819 = RATIONAL(3.0,5.0)*t834; - t845 = t819+RATIONAL(-53.0,98.0)+t812; - t803 = RATIONAL(-2.0,49.0)*t833; - t844 = t803+RATIONAL(191.0,2940.0)+t804; - t843 = t803+t819+RATIONAL(-41.0,98.0); + t821 = RATIONAL(-2.0,5.0)*t835; + t849 = t818+t821; + t848 = t820+t821; + t810 = RATIONAL(-2.0,49.0)*t834; + t847 = t810+t814+RATIONAL(191.0,2940.0); + t804 = RATIONAL(1.0,49.0)*t834; + t846 = t804+t821+RATIONAL(166.0,735.0); + t819 = RATIONAL(3.0,5.0)*t835; + t845 = t810+t819+RATIONAL(-41.0,98.0); + t844 = t804+t819+RATIONAL(-53.0,98.0); + t843 = t849+t853; + t842 = t848+t853; + t841 = t849+t854; + t801 = x*t808; + t840 = t801+t851; t798 = x*t815; - t842 = t798+t851; - t841 = t848+t853; - t840 = t847+t853; - t839 = t847+t852; - t838 = t820+t850; - t837 = t817+t850; - t800 = x*t806; - t836 = t800+t851; - t835 = t848+t852; - t816 = RATIONAL(2.0,35.0)*x; - t814 = RATIONAL(1.0,70.0)*x; - t811 = RATIONAL(-1.0,70.0)*x; - t810 = RATIONAL(1.0,35.0)*x; - t809 = RATIONAL(-1.0,35.0)*x; - t805 = RATIONAL(-2.0,35.0)*x; - t802 = x*t817; - t801 = x*t820; - t799 = RATIONAL(-1.0,20.0)*t849; - t797 = RATIONAL(1.0,20.0)*t849; - coeff_dyy_m2_m2 = t798+t805+t839; - coeff_dyy_m1_m2 = t797+t809+t840; - coeff_dyy_0_m2 = t806+t844; - coeff_dyy_p1_m2 = t810+t799+t840; - coeff_dyy_p2_m2 = t800+t816+t839; - coeff_dyy_m2_m1 = t810+t802+t838; - coeff_dyy_m1_m1 = t820+t814+t836; + t839 = t798+t851; + t838 = t808+t852; + t837 = t848+t854; + t836 = t815+t852; + t817 = RATIONAL(-2.0,35.0)*x; + t812 = RATIONAL(-1.0,35.0)*x; + t811 = RATIONAL(1.0,35.0)*x; + t809 = RATIONAL(2.0,35.0)*x; + t806 = RATIONAL(1.0,70.0)*x; + t805 = RATIONAL(-1.0,70.0)*x; + t803 = x*t818; + t802 = x*t820; + t800 = RATIONAL(1.0,20.0)*t850; + t799 = RATIONAL(-1.0,20.0)*t850; + coeff_dyy_m2_m2 = t815+t817+t840; + coeff_dyy_m1_m2 = t800+t812+t836; + coeff_dyy_0_m2 = t815+t847; + coeff_dyy_p1_m2 = t799+t811+t836; + coeff_dyy_p2_m2 = t815+t809+t839; + coeff_dyy_m2_m1 = t803+t811+t837; + coeff_dyy_m1_m1 = t806+t798+t842; coeff_dyy_0_m1 = t820+t846; - coeff_dyy_p1_m1 = t811+t820+t842; - coeff_dyy_p2_m1 = t809+t801+t838; - coeff_dyy_m2_0 = t816+t843; - coeff_dyy_m1_0 = t810+t845; - coeff_dyy_0_0 = t819+RATIONAL(-57.0,98.0)+t807; - coeff_dyy_p1_0 = t809+t845; - coeff_dyy_p2_0 = t805+t843; - coeff_dyy_m2_p1 = t810+t801+t837; - coeff_dyy_m1_p1 = t817+t814+t842; - coeff_dyy_0_p1 = t817+t846; - coeff_dyy_p1_p1 = t811+t817+t836; - coeff_dyy_p2_p1 = t809+t802+t837; - coeff_dyy_m2_p2 = t805+t800+t835; - coeff_dyy_m1_p2 = t809+t799+t841; - coeff_dyy_0_p2 = t815+t844; - coeff_dyy_p1_p2 = t797+t810+t841; - coeff_dyy_p2_p2 = t816+t798+t835; + coeff_dyy_p1_m1 = t805+t801+t842; + coeff_dyy_p2_m1 = t802+t812+t837; + coeff_dyy_m2_0 = t809+t845; + coeff_dyy_m1_0 = t811+t844; + coeff_dyy_0_0 = RATIONAL(-57.0,98.0)+t819+t816; + coeff_dyy_p1_0 = t812+t844; + coeff_dyy_p2_0 = t817+t845; + coeff_dyy_m2_p1 = t802+t811+t841; + coeff_dyy_m1_p1 = t806+t801+t843; + coeff_dyy_0_p1 = t818+t846; + coeff_dyy_p1_p1 = t805+t798+t843; + coeff_dyy_p2_p1 = t812+t803+t841; + coeff_dyy_m2_p2 = t817+t808+t839; + coeff_dyy_m1_p2 = t812+t799+t838; + coeff_dyy_0_p2 = t808+t847; + coeff_dyy_p1_p2 = t811+t800+t838; + coeff_dyy_p2_p2 = t808+t809+t840; |