diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c')
-rw-r--r-- | src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c | 53 |
1 files changed, 53 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..a9bdd03 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c @@ -0,0 +1,53 @@ +fp t17; +fp t15; +fp t18; +fp t19; +fp t4; +fp t16; +fp t24; +fp t26; +fp t25; +fp t14; +fp t12; +fp t13; +fp t20; +fp t9; +fp t23; +fp t11; +fp t22; +fp t21; +fp t8; +fp t6; +fp t5; +fp t3; + t17 = x*x; + t15 = RATIONAL(1.0,4.0); + t18 = y*y; + t19 = t15*t18; + t4 = t17*t19; + t16 = RATIONAL(-1.0,4.0); + t24 = y*t16; + t26 = t4+t17*t24; + t25 = y*t15; + t14 = RATIONAL(1.0,2.0); + t12 = t14*t18; + t13 = RATIONAL(-1.0,2.0); + t20 = t13*t18; + t9 = t17*t20; + t23 = t12+t9; + t11 = t14*t17; + t22 = t11+t9; + t21 = t4+t17*t25; + t8 = x*t19; + t6 = x*t24; + t5 = x*t25; + t3 = t16*x*t18; + coeffs_I->coeff_m1_m1 = t3+t5+t26; + coeffs_I->coeff_0_m1 = (t13+t11)*y+t23; + coeffs_I->coeff_p1_m1 = t8+t6+t26; + coeffs_I->coeff_m1_0 = (t12+t13)*x+t22; + coeffs_I->coeff_0_0 = -t18+t17*t18+RATIONAL(1.0,1.0)-t17; + coeffs_I->coeff_p1_0 = (t20+t14)*x+t22; + coeffs_I->coeff_m1_p1 = t6+t3+t21; + coeffs_I->coeff_0_p1 = (t13*t17+t14)*y+t23; + coeffs_I->coeff_p1_p1 = t5+t8+t21; |