diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0')
12 files changed, 516 insertions, 518 deletions
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 deleted file mode 100644 index a022b8a..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,196 +0,0 @@ -fp t160, - t131, - t159, - t133, - t158, - t117, - t66, - t132, - t128, - t87, - t157, - t129, - t84, - t156, - t124, - t68, - t130, - t116, - t77, - t155, - t126, - t70, - t125, - t83, - t154, - t71, - t153, - t108, - t122, - t69, - t152, - t151, - t63, - t150, - t67, - t86, - t149, - t65, - t148, - t85, - t147, - t79, - t146, - t82, - t145, - t144, - t143, - t142, - t123, - t141, - t101, - t140, - t139, - t138, - t137, - t136, - t127, - t121, - t119, - t118, - t115, - t113, - t112, - t111, - t109, - t106, - t105, - t103, - t102, - t100, - t98, - t97, - t96, - t95, - t94, - t93, - t92, - t91, - t90, - t89, - t88, - t81, - t80, - t78, - t76, - t75, - t74, - t72; - t160 = x*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; - 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; - 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); - 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 = t133+t131; - t127 = RATIONAL(1.0,80.0); - t121 = RATIONAL(9.0,80.0); - t119 = RATIONAL(-1.0,80.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(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(-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+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 = 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 deleted file mode 100644 index 5bfa6f4..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,95 +0,0 @@ -fp t201, - t161, - t200, - t184, - t188, - t174, - t187, - t180, - t199, - t164, - t198, - t183, - t175, - t197, - t162, - t196, - t181, - t195, - t194, - t163, - t193, - t186, - t176, - t182, - t192, - t185, - t172, - t191, - t179, - t190, - t189, - t178, - t177, - t169, - t167, - t166, - t165; - t201 = x*y; - t161 = RATIONAL(1.0,20.0)*t201; - t200 = t161+RATIONAL(9.0,40.0)*x; - t184 = RATIONAL(-1.0,40.0); - t188 = y*y; - t174 = t184*t188; - t187 = x*x; - 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 = t162+RATIONAL(11.0,40.0)*x; - t181 = RATIONAL(3.0,8.0)*t187; - 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); - t176 = t186*t188; - t182 = RATIONAL(-1.0,8.0)*t187; - 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; - 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+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)+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 deleted file mode 100644 index 29b5494..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c +++ /dev/null @@ -1,40 +0,0 @@ -fp t244, - t258, - t246, - t257, - t245, - t256, - t243, - t255, - t250, - t249, - t248, - t247; - 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 deleted file mode 100644 index 20e2106..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c +++ /dev/null @@ -1,52 +0,0 @@ -fp t276, - t275, - t274, - t273, - t272, - t271, - t270, - t269, - t268, - t267, - t266, - t265, - t264, - t263, - t262, - t261, - t260, - t259; - t276 = RATIONAL(-1.0,50.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(-1.0,100.0); - t268 = RATIONAL(7.0,100.0); - t267 = RATIONAL(-13.0,100.0); - t266 = t272*y; - t265 = t271*y; - t264 = t273*x; - t263 = t274*x; - t262 = t271*x; - t261 = t272*x; - 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 deleted file mode 100644 index c766157..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c +++ /dev/null @@ -1,95 +0,0 @@ -fp t242, - t225, - t228, - t212, - t229, - t220, - t241, - t203, - t240, - t204, - t239, - t205, - t238, - t224, - t213, - t237, - t202, - t236, - t226, - t215, - t222, - t235, - t227, - t218, - t234, - t221, - t233, - t232, - t223, - t231, - t230, - t219, - t217, - t209, - t208, - t207, - t206; - t242 = x*y; - t225 = RATIONAL(-1.0,40.0); - t228 = x*x; - t212 = t225*t228; - t229 = y*y; - 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; - t237 = t213+t220; - t202 = RATIONAL(3.0,20.0)*t242; - 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); - 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; - 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(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 deleted file mode 100644 index b449df9..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c +++ /dev/null @@ -1,40 +0,0 @@ -fp t277, - t292, - t279, - t291, - t280, - t290, - 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; - t291 = t279+RATIONAL(11.0,40.0); - t280 = RATIONAL(-3.0,20.0)*x; - t290 = t280+RATIONAL(13.0,40.0); - t278 = RATIONAL(1.0,20.0)*x; - t289 = t278+RATIONAL(9.0,40.0); - t284 = RATIONAL(-1.0,4.0)*y; - 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.order3.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..ac11344 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c @@ -0,0 +1,190 @@ +fp t158; +fp t131; +fp t157; +fp t132; +fp t156; +fp t117; +fp t63; +fp t129; +fp t130; +fp t84; +fp t155; +fp t154; +fp t128; +fp t87; +fp t153; +fp t126; +fp t64; +fp t133; +fp t116; +fp t76; +fp t152; +fp t67; +fp t86; +fp t151; +fp t123; +fp t71; +fp t150; +fp t124; +fp t65; +fp t149; +fp t125; +fp t81; +fp t148; +fp t85; +fp t147; +fp t146; +fp t112; +fp t69; +fp t145; +fp t110; +fp t66; +fp t144; +fp t143; +fp t70; +fp t120; +fp t142; +fp t122; +fp t141; +fp t140; +fp t101; +fp t139; +fp t138; +fp t137; +fp t136; +fp t127; +fp t121; +fp t118; +fp t115; +fp t113; +fp t109; +fp t107; +fp t106; +fp t104; +fp t103; +fp t102; +fp t100; +fp t98; +fp t97; +fp t96; +fp t95; +fp t94; +fp t93; +fp t92; +fp t91; +fp t90; +fp t89; +fp t88; +fp t83; +fp t82; +fp t80; +fp t78; +fp t77; +fp t74; +fp t73; + t158 = x*y; + t131 = x*x; + t157 = t131*y; + t132 = y*y; + t156 = x*t132; + t117 = RATIONAL(1.0,40.0); + t63 = t117*t157; + t129 = RATIONAL(1.0,8.0); + t130 = t132*y; + t84 = t129*t130; + t155 = t63+t84; + t154 = RATIONAL(-1.0,100.0)*t158+RATIONAL(37.0,200.0); + t128 = RATIONAL(-1.0,8.0); + t87 = t128*t130; + t153 = t63+t87; + t126 = RATIONAL(-3.0,40.0); + t64 = t126*t156; + t133 = t131*x; + t116 = RATIONAL(1.0,24.0); + t76 = t116*t133; + t152 = t64+t76; + t67 = t117*t156; + t86 = t129*t133; + t151 = t67+t86; + t123 = RATIONAL(-1.0,40.0); + t71 = t123*t156; + t150 = t71+t86; + t124 = RATIONAL(3.0,40.0); + t65 = t124*t156; + t149 = t65+t76; + t125 = RATIONAL(-1.0,24.0); + t81 = t125*t133; + t148 = t64+t81; + t85 = t128*t133; + t147 = t71+t85; + t146 = t65+t81; + t112 = RATIONAL(13.0,80.0); + t69 = t126*t157; + t145 = t69+t112*t131; + t110 = RATIONAL(11.0,80.0); + t66 = t123*t157; + t144 = t66+t110*t131; + t143 = t67+t85; + t70 = t124*t157; + t120 = RATIONAL(7.0,80.0); + t142 = t70+t120*t131; + t122 = RATIONAL(-1.0,50.0); + t141 = t122+RATIONAL(-13.0,100.0)*t158; + t140 = RATIONAL(2.0,25.0)*t158+RATIONAL(-2.0,25.0); + t101 = RATIONAL(-17.0,200.0); + t139 = RATIONAL(-9.0,100.0)*t158+t101; + t138 = RATIONAL(11.0,50.0)+t122*t158; + t137 = RATIONAL(7.0,100.0)*t158+RATIONAL(13.0,100.0); + t136 = t131+t132; + t127 = RATIONAL(1.0,80.0); + t121 = RATIONAL(9.0,80.0); + t118 = RATIONAL(-1.0,80.0); + t115 = RATIONAL(3.0,80.0); + t113 = RATIONAL(-17.0,80.0); + t109 = RATIONAL(-23.0,80.0); + t107 = RATIONAL(-19.0,80.0); + t106 = RATIONAL(-63.0,400.0); + t104 = 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(111.0,400.0); + t94 = RATIONAL(37.0,1200.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); + t83 = t121*t131; + t82 = t125*t130; + t80 = t121*t132; + t78 = t116*t130; + t77 = t120*t132; + t74 = t112*t132; + t73 = t110*t132; + coeffs_I->coeff_m1_m1 = t74+t82+t92*y+RATIONAL(-23.0,200.0)+(t92+RATIONAL +(6.0,25.0)*y)*x+t145+t148; + coeffs_I->coeff_0_m1 = t70+t91*y+t73+t103*x+t82+t109*t131+t138+t150; + coeffs_I->coeff_p1_m1 = t80+t82+t95*x+t90*y+t141+t142+t143; + coeffs_I->coeff_p2_m1 = t93*y+t89*x+t77+t82+t69+t115*t131+t139+t149; + coeffs_I->coeff_m1_0 = t109*t132+t84+t91*x+t103*y+t138+t144+t146; + coeffs_I->coeff_0_0 = RATIONAL(93.0,200.0)+t100*y+t136*RATIONAL(-21.0, +80.0)+(t100+RATIONAL(-1.0,25.0)*y)*x+t151+t155; + coeffs_I->coeff_p1_0 = t107*t132+t106*y+t83+t98*x+t147+t154+t155; + coeffs_I->coeff_p2_0 = t66+t113*t132+t127*t131+t84+t104*y+t97*x+t137+t152 +; + coeffs_I->coeff_m1_p1 = t95*y+t90*x+t77+t83+t141+t146+t153; + coeffs_I->coeff_0_p1 = t87+t106*x+t66+t80+t98*y+t107*t131+t151+t154; + coeffs_I->coeff_p1_p1 = t87+t73+t101+t96*y+(t96+RATIONAL(3.0,50.0)*y)*x+ +t144+t147; + coeffs_I->coeff_p2_p1 = t118*t131+t102*y+t74+t94*x+t140+t152+t153; + coeffs_I->coeff_m1_p2 = t93*x+t89*y+t115*t132+t78+t139+t142+t148; + coeffs_I->coeff_0_p2 = t97*y+t69+t127*t132+t113*t131+t78+t104*x+t137+t150 +; + coeffs_I->coeff_p1_p2 = t102*x+t94*y+t78+t118*t132+t140+t143+t145; + coeffs_I->coeff_p2_p2 = t70+t88*y+t78+RATIONAL(7.0,200.0)+t136*RATIONAL( +-3.0,80.0)+(t88+RATIONAL(-3.0,50.0)*y)*x+t149; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..0ea2d65 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c @@ -0,0 +1,97 @@ +fp t199; +fp t185; +fp t178; +fp t181; +fp t186; +fp t173; +fp t198; +fp t159; +fp t197; +fp t182; +fp t175; +fp t196; +fp t161; +fp t195; +fp t162; +fp t194; +fp t160; +fp t193; +fp t179; +fp t184; +fp t176; +fp t192; +fp t183; +fp t171; +fp t191; +fp t180; +fp t190; +fp t189; +fp t177; +fp t188; +fp t187; +fp t174; +fp t169; +fp t167; +fp t165; +fp t164; +fp t163; + t199 = x*y; + t185 = x*x; + t178 = RATIONAL(-3.0,8.0)*t185; + t181 = RATIONAL(1.0,40.0); + t186 = y*y; + t173 = t181*t186; + t198 = t178+t173; + t159 = RATIONAL(-3.0,20.0)*t199; + t197 = RATIONAL(13.0,40.0)*x+t159; + t182 = RATIONAL(-1.0,40.0); + t175 = t182*t186; + t196 = t178+t175; + t161 = RATIONAL(-1.0,20.0)*t199; + t195 = t161+RATIONAL(11.0,40.0)*x; + t162 = RATIONAL(3.0,20.0)*t199; + t194 = RATIONAL(7.0,40.0)*x+t162; + t160 = RATIONAL(1.0,20.0)*t199; + t193 = RATIONAL(9.0,40.0)*x+t160; + t179 = RATIONAL(-1.0,8.0)*t185; + t184 = RATIONAL(-3.0,40.0); + t176 = t184*t186; + t192 = t179+t176; + t183 = RATIONAL(3.0,40.0); + t171 = t183*t186; + t191 = t179+t171; + t180 = RATIONAL(3.0,8.0)*t185; + t190 = t180+t175; + t189 = t180+t173; + t177 = RATIONAL(1.0,8.0)*t185; + t188 = t171+t177; + t187 = t176+t177; + t174 = RATIONAL(-1.0,50.0)*y; + t169 = RATIONAL(2.0,25.0)*y; + t167 = RATIONAL(-9.0,100.0)*y; + t165 = RATIONAL(-1.0,100.0)*y; + t164 = RATIONAL(7.0,100.0)*y; + t163 = RATIONAL(-13.0,100.0)*y; + coeffs_dx->coeff_m1_m1 = RATIONAL(6.0,25.0)*y+RATIONAL(-109.0,1200.0)+ +t192+t197; + coeffs_dx->coeff_0_m1 = t162+t174+RATIONAL(-31.0,400.0)+RATIONAL(-23.0, +40.0)*x+t190; + coeffs_dx->coeff_p1_m1 = RATIONAL(111.0,400.0)+t163+t194+t198; + coeffs_dx->coeff_p2_m1 = t183*x+RATIONAL(-131.0,1200.0)+t167+t159+t188; + coeffs_dx->coeff_m1_0 = RATIONAL(-223.0,1200.0)+t174+t191+t195; + coeffs_dx->coeff_0_0 = RATIONAL(-1.0,25.0)*y+t160+RATIONAL(-21.0,40.0)*x+ +RATIONAL(-57.0,400.0)+t189; + coeffs_dx->coeff_p1_0 = t165+RATIONAL(117.0,400.0)+t193+t196; + coeffs_dx->coeff_p2_0 = t161+t181*x+RATIONAL(43.0,1200.0)+t164+t187; + coeffs_dx->coeff_m1_p1 = RATIONAL(-157.0,1200.0)+t163+t191+t193; + coeffs_dx->coeff_0_p1 = RATIONAL(-63.0,400.0)+t165+t161+RATIONAL(-19.0, +40.0)*x+t189; + coeffs_dx->coeff_p1_p1 = RATIONAL(103.0,400.0)+RATIONAL(3.0,50.0)*y+t195+ +t196; + coeffs_dx->coeff_p2_p1 = t160+t169+RATIONAL(37.0,1200.0)+t182*x+t187; + coeffs_dx->coeff_m1_p2 = RATIONAL(89.0,1200.0)+t167+t192+t194; + coeffs_dx->coeff_0_p2 = t159+t164+RATIONAL(-17.0,40.0)*x+RATIONAL(-49.0, +400.0)+t190; + coeffs_dx->coeff_p1_p2 = RATIONAL(69.0,400.0)+t169+t197+t198; + coeffs_dx->coeff_p2_p2 = RATIONAL(-149.0,1200.0)+t162+RATIONAL(-3.0,50.0) +*y+t184*x+t188; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..ba139bf --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,40 @@ +fp t244; +fp t256; +fp t242; +fp t255; +fp t241; +fp t254; +fp t243; +fp t253; +fp t248; +fp t247; +fp t246; +fp t245; + t244 = RATIONAL(-3.0,20.0)*y; + t256 = t244+RATIONAL(13.0,40.0); + t242 = RATIONAL(-1.0,20.0)*y; + t255 = t242+RATIONAL(11.0,40.0); + t241 = RATIONAL(1.0,20.0)*y; + t254 = t241+RATIONAL(9.0,40.0); + t243 = RATIONAL(3.0,20.0)*y; + t253 = t243+RATIONAL(7.0,40.0); + t248 = RATIONAL(-1.0,4.0)*x; + t247 = RATIONAL(1.0,4.0)*x; + t246 = RATIONAL(-3.0,4.0)*x; + t245 = RATIONAL(3.0,4.0)*x; + coeffs_dxx->coeff_m1_m1 = t248+t256; + coeffs_dxx->coeff_0_m1 = t245+t243+RATIONAL(-23.0,40.0); + coeffs_dxx->coeff_p1_m1 = t246+t253; + coeffs_dxx->coeff_p2_m1 = t244+RATIONAL(3.0,40.0)+t247; + coeffs_dxx->coeff_m1_0 = t248+t255; + coeffs_dxx->coeff_0_0 = RATIONAL(-21.0,40.0)+t241+t245; + coeffs_dxx->coeff_p1_0 = t246+t254; + coeffs_dxx->coeff_p2_0 = t247+t242+RATIONAL(1.0,40.0); + coeffs_dxx->coeff_m1_p1 = t248+t254; + coeffs_dxx->coeff_0_p1 = RATIONAL(-19.0,40.0)+t242+t245; + coeffs_dxx->coeff_p1_p1 = t246+t255; + coeffs_dxx->coeff_p2_p1 = RATIONAL(-1.0,40.0)+t241+t247; + coeffs_dxx->coeff_m1_p2 = t248+t253; + coeffs_dxx->coeff_0_p2 = t244+t245+RATIONAL(-17.0,40.0); + coeffs_dxx->coeff_p1_p2 = t246+t256; + coeffs_dxx->coeff_p2_p2 = RATIONAL(-3.0,40.0)+t247+t243; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c new file mode 100644 index 0000000..e9208a2 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c @@ -0,0 +1,52 @@ +fp t274; +fp t273; +fp t272; +fp t271; +fp t270; +fp t269; +fp t268; +fp t267; +fp t266; +fp t265; +fp t264; +fp t263; +fp t262; +fp t261; +fp t260; +fp t259; +fp t258; +fp t257; + t274 = RATIONAL(-3.0,20.0); + t273 = RATIONAL(3.0,20.0); + t272 = RATIONAL(2.0,25.0); + t271 = RATIONAL(-1.0,20.0); + t270 = RATIONAL(1.0,20.0); + t269 = RATIONAL(-1.0,50.0); + t268 = RATIONAL(-1.0,100.0); + t267 = RATIONAL(7.0,100.0); + t266 = RATIONAL(-9.0,100.0); + t265 = RATIONAL(-13.0,100.0); + t264 = t274*x; + t263 = t273*x; + t262 = t271*x; + t261 = t274*y; + t260 = t270*x; + t259 = t273*y; + t258 = t271*y; + t257 = t270*y; + coeffs_dxy->coeff_m1_m1 = t264+RATIONAL(6.0,25.0)+t261; + coeffs_dxy->coeff_0_m1 = t263+t269+t258; + coeffs_dxy->coeff_p1_m1 = t257+t265+t263; + coeffs_dxy->coeff_p2_m1 = t259+t264+t266; + coeffs_dxy->coeff_m1_0 = t259+t262+t269; + coeffs_dxy->coeff_0_0 = t260+RATIONAL(-1.0,25.0)+t257; + coeffs_dxy->coeff_p1_0 = t258+t268+t260; + coeffs_dxy->coeff_p2_0 = t261+t262+t267; + coeffs_dxy->coeff_m1_p1 = t259+t265+t260; + coeffs_dxy->coeff_0_p1 = t262+t257+t268; + coeffs_dxy->coeff_p1_p1 = t262+RATIONAL(3.0,50.0)+t258; + coeffs_dxy->coeff_p2_p1 = t261+t260+t272; + coeffs_dxy->coeff_m1_p2 = t266+t261+t263; + coeffs_dxy->coeff_0_p2 = t267+t258+t264; + coeffs_dxy->coeff_p1_p2 = t264+t257+t272; + coeffs_dxy->coeff_p2_p2 = t263+t259+RATIONAL(-3.0,50.0); diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c new file mode 100644 index 0000000..c12d9d1 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c @@ -0,0 +1,97 @@ +fp t240; +fp t227; +fp t219; +fp t226; +fp t222; +fp t211; +fp t239; +fp t201; +fp t238; +fp t223; +fp t213; +fp t237; +fp t200; +fp t236; +fp t203; +fp t235; +fp t202; +fp t234; +fp t224; +fp t216; +fp t221; +fp t233; +fp t225; +fp t217; +fp t232; +fp t218; +fp t231; +fp t230; +fp t220; +fp t229; +fp t228; +fp t215; +fp t210; +fp t207; +fp t206; +fp t205; +fp t204; + t240 = x*y; + t227 = y*y; + t219 = RATIONAL(-3.0,8.0)*t227; + t226 = x*x; + t222 = RATIONAL(1.0,40.0); + t211 = t222*t226; + t239 = t219+t211; + t201 = RATIONAL(1.0,20.0)*t240; + t238 = t201+RATIONAL(9.0,40.0)*y; + t223 = RATIONAL(-1.0,40.0); + t213 = t223*t226; + t237 = t219+t213; + t200 = RATIONAL(3.0,20.0)*t240; + t236 = t200+RATIONAL(7.0,40.0)*y; + t203 = RATIONAL(-3.0,20.0)*t240; + t235 = t203+RATIONAL(13.0,40.0)*y; + t202 = RATIONAL(-1.0,20.0)*t240; + t234 = t202+RATIONAL(11.0,40.0)*y; + t224 = RATIONAL(3.0,40.0); + t216 = t224*t226; + t221 = RATIONAL(-1.0,8.0)*t227; + t233 = t216+t221; + t225 = RATIONAL(-3.0,40.0); + t217 = t225*t226; + t232 = t217+t221; + t218 = RATIONAL(1.0,8.0)*t227; + t231 = t218+t217; + t230 = t216+t218; + t220 = RATIONAL(3.0,8.0)*t227; + t229 = t220+t213; + t228 = t220+t211; + t215 = RATIONAL(-1.0,50.0)*x; + t210 = RATIONAL(2.0,25.0)*x; + t207 = RATIONAL(-9.0,100.0)*x; + t206 = RATIONAL(-1.0,100.0)*x; + t205 = RATIONAL(7.0,100.0)*x; + t204 = RATIONAL(-13.0,100.0)*x; + coeffs_dy->coeff_m1_m1 = RATIONAL(6.0,25.0)*x+RATIONAL(-109.0,1200.0)+ +t232+t235; + coeffs_dy->coeff_0_m1 = t215+RATIONAL(-223.0,1200.0)+t233+t234; + coeffs_dy->coeff_p1_m1 = RATIONAL(-157.0,1200.0)+t204+t233+t238; + coeffs_dy->coeff_p2_m1 = t207+RATIONAL(89.0,1200.0)+t232+t236; + coeffs_dy->coeff_m1_0 = RATIONAL(-23.0,40.0)*y+RATIONAL(-31.0,400.0)+t200 ++t215+t229; + coeffs_dy->coeff_0_0 = RATIONAL(-21.0,40.0)*y+RATIONAL(-57.0,400.0)+ +RATIONAL(-1.0,25.0)*x+t201+t228; + coeffs_dy->coeff_p1_0 = t202+RATIONAL(-19.0,40.0)*y+RATIONAL(-63.0,400.0) ++t206+t228; + coeffs_dy->coeff_p2_0 = t203+RATIONAL(-49.0,400.0)+RATIONAL(-17.0,40.0)*y ++t205+t229; + coeffs_dy->coeff_m1_p1 = t204+RATIONAL(111.0,400.0)+t236+t239; + coeffs_dy->coeff_0_p1 = t206+RATIONAL(117.0,400.0)+t237+t238; + coeffs_dy->coeff_p1_p1 = RATIONAL(3.0,50.0)*x+RATIONAL(103.0,400.0)+t234+ +t237; + coeffs_dy->coeff_p2_p1 = t210+RATIONAL(69.0,400.0)+t235+t239; + coeffs_dy->coeff_m1_p2 = t207+t203+RATIONAL(-131.0,1200.0)+t224*y+t230; + coeffs_dy->coeff_0_p2 = t205+t202+t222*y+RATIONAL(43.0,1200.0)+t231; + coeffs_dy->coeff_p1_p2 = t201+RATIONAL(37.0,1200.0)+t223*y+t210+t231; + coeffs_dy->coeff_p2_p2 = t225*y+t200+RATIONAL(-3.0,50.0)*x+RATIONAL( +-149.0,1200.0)+t230; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c new file mode 100644 index 0000000..6bbdde7 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c @@ -0,0 +1,40 @@ +fp t276; +fp t290; +fp t277; +fp t289; +fp t278; +fp t288; +fp t275; +fp t287; +fp t282; +fp t281; +fp t280; +fp t279; + t276 = RATIONAL(-1.0,20.0)*x; + t290 = t276+RATIONAL(11.0,40.0); + t277 = RATIONAL(3.0,20.0)*x; + t289 = RATIONAL(7.0,40.0)+t277; + t278 = RATIONAL(-3.0,20.0)*x; + t288 = t278+RATIONAL(13.0,40.0); + t275 = RATIONAL(1.0,20.0)*x; + t287 = t275+RATIONAL(9.0,40.0); + t282 = RATIONAL(1.0,4.0)*y; + t281 = RATIONAL(-1.0,4.0)*y; + t280 = RATIONAL(3.0,4.0)*y; + t279 = RATIONAL(-3.0,4.0)*y; + coeffs_dyy->coeff_m1_m1 = t281+t288; + coeffs_dyy->coeff_0_m1 = t281+t290; + coeffs_dyy->coeff_p1_m1 = t281+t287; + coeffs_dyy->coeff_p2_m1 = t281+t289; + coeffs_dyy->coeff_m1_0 = t277+t280+RATIONAL(-23.0,40.0); + coeffs_dyy->coeff_0_0 = t280+t275+RATIONAL(-21.0,40.0); + coeffs_dyy->coeff_p1_0 = t280+RATIONAL(-19.0,40.0)+t276; + coeffs_dyy->coeff_p2_0 = RATIONAL(-17.0,40.0)+t280+t278; + coeffs_dyy->coeff_m1_p1 = t279+t289; + coeffs_dyy->coeff_0_p1 = t279+t287; + coeffs_dyy->coeff_p1_p1 = t279+t290; + coeffs_dyy->coeff_p2_p1 = t279+t288; + coeffs_dyy->coeff_m1_p2 = t278+t282+RATIONAL(3.0,40.0); + coeffs_dyy->coeff_0_p2 = t282+RATIONAL(1.0,40.0)+t276; + coeffs_dyy->coeff_p1_p2 = t275+t282+RATIONAL(-1.0,40.0); + coeffs_dyy->coeff_p2_p2 = t277+RATIONAL(-3.0,40.0)+t282; |