diff options
author | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-05-14 15:54:41 +0000 |
---|---|---|
committer | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-05-14 15:54:41 +0000 |
commit | 3f0ca6200461eaca761dcc73da676a9b373aca83 (patch) | |
tree | 5404cf6520498018bcc99679c31c1a80ae915111 /src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0 | |
parent | fc338426770f9069e68714142122c91d45fbb267 (diff) |
generate coeff-*.store.c files with non-I coefficients multiplied by "factor"
--> needed to handle 1/dx factors for derivatives
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@43 df1f8a13-aa1d-4dd4-9681-27ded5b42416
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0')
3 files changed, 70 insertions, 70 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c index d6cc566..1f7b51b 100644 --- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c +++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c @@ -1,38 +1,38 @@ -fp t104, - t108, +fp t96, + t100, + t101, + t104, t107, - t103, + t102, + t106, t113, - t101, - t96, - t97, t112, t111, - t102, - t100, + t103, t99, - t98; - t104 = x*x; - t108 = t104*t104; - t107 = x*t104; - t103 = t107*t107; - t113 = RATIONAL(-1.0,120.0)*t103+RATIONAL(1.0,12.0)*t108+RATIONAL(-3.0, -40.0)*t104; - t101 = RATIONAL(3.0,4.0); + t98, + t97; t96 = RATIONAL(-13.0,48.0); - t97 = RATIONAL(1.0,48.0); - t112 = t96*t108+t101*t104+t97*t103; - t111 = RATIONAL(-1.0,144.0)*t108+RATIONAL(1.0,180.0)*t104+RATIONAL(1.0, -720.0)*t103; - t102 = x*t108; - t100 = RATIONAL(1.0,60.0); - t99 = RATIONAL(-1.0,48.0); - t98 = RATIONAL(-1.0,60.0); - coeff_I_m3 = t98*x+t97*t107+RATIONAL(-1.0,240.0)*t102+t111; - coeff_I_m2 = RATIONAL(3.0,20.0)*x+RATIONAL(-1.0,6.0)*t107+t100*t102+t113; - coeff_I_m1 = RATIONAL(-3.0,4.0)*x+RATIONAL(13.0,48.0)*t107+t99*t102+t112; - coeff_I_0 = RATIONAL(7.0,18.0)*t108+RATIONAL(-49.0,36.0)*t104+RATIONAL( -1.0,1.0)+RATIONAL(-1.0,36.0)*t103; - coeff_I_p1 = t96*t107+t101*x+t97*t102+t112; - coeff_I_p2 = RATIONAL(1.0,6.0)*t107+RATIONAL(-3.0,20.0)*x+t98*t102+t113; - coeff_I_p3 = t99*t107+t100*x+RATIONAL(1.0,240.0)*t102+t111; + t100 = RATIONAL(1.0,48.0); + t101 = RATIONAL(3.0,4.0); + t104 = x*x; + t107 = t104*x; + t102 = t107*t107; + t106 = t104*t104; + t113 = t100*t102+t101*t104+t96*t106; + t112 = RATIONAL(1.0,720.0)*t102+RATIONAL(1.0,180.0)*t104+RATIONAL(-1.0, +144.0)*t106; + t111 = RATIONAL(1.0,12.0)*t106+RATIONAL(-3.0,40.0)*t104+RATIONAL(-1.0, +120.0)*t102; + t103 = t104*t107; + t99 = RATIONAL(-1.0,60.0); + t98 = RATIONAL(1.0,60.0); + t97 = RATIONAL(-1.0,48.0); + coeff_I_m3 = t99*x+t100*t107+RATIONAL(-1.0,240.0)*t103+t112; + coeff_I_m2 = RATIONAL(3.0,20.0)*x+RATIONAL(-1.0,6.0)*t107+t98*t103+t111; + coeff_I_m1 = RATIONAL(-3.0,4.0)*x+RATIONAL(13.0,48.0)*t107+t97*t103+t113; + coeff_I_0 = RATIONAL(7.0,18.0)*t106+RATIONAL(-49.0,36.0)*t104+RATIONAL( +1.0,1.0)+RATIONAL(-1.0,36.0)*t102; + coeff_I_p1 = t96*t107+t101*x+t100*t103+t113; + coeff_I_p2 = RATIONAL(1.0,6.0)*t107+RATIONAL(-3.0,20.0)*x+t99*t103+t111; + coeff_I_p3 = t97*t107+t98*x+RATIONAL(1.0,240.0)*t103+t112; diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c index 252c658..738bb32 100644 --- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c +++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c @@ -1,31 +1,31 @@ fp t125, - t126, - t129, + t128, t124, - t123, t133, + t123, t132, - t131; + t131, + t127; t125 = x*x; - t126 = x*t125; - t129 = t125*t125; - t124 = x*t129; + t128 = x*t125; + t124 = t128*t125; + t133 = RATIONAL(-1.0,36.0)*t128+RATIONAL(1.0,90.0)*x+RATIONAL(1.0,120.0)* +t124; t123 = RATIONAL(-3.0,20.0); - t133 = RATIONAL(-1.0,20.0)*t124+t123*x+RATIONAL(1.0,3.0)*t126; - t132 = RATIONAL(-13.0,12.0)*t126+RATIONAL(1.0,8.0)*t124+RATIONAL(3.0,2.0) -*x; - t131 = RATIONAL(1.0,90.0)*x+RATIONAL(-1.0,36.0)*t126+RATIONAL(1.0,120.0)* + t132 = RATIONAL(-1.0,20.0)*t124+t123*x+RATIONAL(1.0,3.0)*t128; + t131 = RATIONAL(-13.0,12.0)*t128+RATIONAL(3.0,2.0)*x+RATIONAL(1.0,8.0)* t124; - coeff_dx_m3 = RATIONAL(1.0,16.0)*t125+RATIONAL(-1.0,48.0)*t129+RATIONAL( --1.0,60.0)+t131; - coeff_dx_m2 = RATIONAL(1.0,12.0)*t129+RATIONAL(3.0,20.0)+RATIONAL(-1.0, -2.0)*t125+t133; - coeff_dx_m1 = RATIONAL(-3.0,4.0)+RATIONAL(-5.0,48.0)*t129+RATIONAL(13.0, -16.0)*t125+t132; - coeff_dx_0 = RATIONAL(-1.0,6.0)*t124+RATIONAL(-49.0,18.0)*x+RATIONAL(14.0 -,9.0)*t126; - coeff_dx_p1 = RATIONAL(-13.0,16.0)*t125+RATIONAL(5.0,48.0)*t129+RATIONAL( -3.0,4.0)+t132; - coeff_dx_p2 = t123+RATIONAL(-1.0,12.0)*t129+RATIONAL(1.0,2.0)*t125+t133; - coeff_dx_p3 = RATIONAL(1.0,60.0)+RATIONAL(-1.0,16.0)*t125+RATIONAL(1.0, -48.0)*t129+t131; + t127 = t125*t125; + coeff_dx_m3 = RATIONAL(1.0,16.0)*t125+RATIONAL(-1.0,48.0)*t127+RATIONAL( +-1.0,60.0)+t133; + coeff_dx_m2 = RATIONAL(-1.0,2.0)*t125+RATIONAL(3.0,20.0)+RATIONAL(1.0, +12.0)*t127+t132; + coeff_dx_m1 = RATIONAL(-3.0,4.0)+RATIONAL(13.0,16.0)*t125+RATIONAL(-5.0, +48.0)*t127+t131; + coeff_dx_0 = RATIONAL(-1.0,6.0)*t124+RATIONAL(14.0,9.0)*t128+RATIONAL( +-49.0,18.0)*x; + coeff_dx_p1 = RATIONAL(-13.0,16.0)*t125+RATIONAL(3.0,4.0)+RATIONAL(5.0, +48.0)*t127+t131; + coeff_dx_p2 = t123+RATIONAL(-1.0,12.0)*t127+RATIONAL(1.0,2.0)*t125+t132; + coeff_dx_p3 = RATIONAL(1.0,48.0)*t127+RATIONAL(-1.0,16.0)*t125+RATIONAL( +1.0,60.0)+t133; diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c index 55b1f9d..e63d0f2 100644 --- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c +++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c @@ -1,22 +1,22 @@ fp t143, t145, t150, - t141, + t140, t149, t148, t144; t143 = x*x; t145 = t143*t143; - t150 = RATIONAL(-1.0,4.0)*t145+t143+RATIONAL(-3.0,20.0); - t141 = RATIONAL(-1.0,12.0); - t149 = t141*t143+RATIONAL(1.0,24.0)*t145+RATIONAL(1.0,90.0); - t148 = RATIONAL(3.0,2.0)+RATIONAL(-13.0,4.0)*t143+RATIONAL(5.0,8.0)*t145; - t144 = x*t143; - coeff_dxx_m3 = t141*t144+RATIONAL(1.0,8.0)*x+t149; - coeff_dxx_m2 = RATIONAL(1.0,3.0)*t144-x+t150; - coeff_dxx_m1 = RATIONAL(-5.0,12.0)*t144+RATIONAL(13.0,8.0)*x+t148; - coeff_dxx_0 = RATIONAL(-5.0,6.0)*t145+RATIONAL(14.0,3.0)*t143+RATIONAL( --49.0,18.0); - coeff_dxx_p1 = RATIONAL(5.0,12.0)*t144+RATIONAL(-13.0,8.0)*x+t148; - coeff_dxx_p2 = RATIONAL(-1.0,3.0)*t144+x+t150; - coeff_dxx_p3 = RATIONAL(-1.0,8.0)*x+RATIONAL(1.0,12.0)*t144+t149; + t150 = RATIONAL(-13.0,4.0)*t143+RATIONAL(5.0,8.0)*t145+RATIONAL(3.0,2.0); + t140 = RATIONAL(-1.0,12.0); + t149 = t140*t143+RATIONAL(1.0,24.0)*t145+RATIONAL(1.0,90.0); + t148 = RATIONAL(-3.0,20.0)+RATIONAL(-1.0,4.0)*t145+t143; + t144 = t143*x; + coeff_dxx_m3 = t140*t144+RATIONAL(1.0,8.0)*x+t149; + coeff_dxx_m2 = RATIONAL(1.0,3.0)*t144-x+t148; + coeff_dxx_m1 = RATIONAL(-5.0,12.0)*t144+RATIONAL(13.0,8.0)*x+t150; + coeff_dxx_0 = RATIONAL(14.0,3.0)*t143+RATIONAL(-49.0,18.0)+RATIONAL(-5.0, +6.0)*t145; + coeff_dxx_p1 = RATIONAL(5.0,12.0)*t144+RATIONAL(-13.0,8.0)*x+t150; + coeff_dxx_p2 = x+RATIONAL(-1.0,3.0)*t144+t148; + coeff_dxx_p3 = RATIONAL(1.0,12.0)*t144+RATIONAL(-1.0,8.0)*x+t149; |