diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c')
-rw-r--r-- | src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c | 35 |
1 files changed, 35 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c new file mode 100644 index 0000000..e373f6b --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c @@ -0,0 +1,35 @@ +fp t50; +fp t56; +fp t53; +fp t43; +fp t51; +fp t55; +fp t52; +fp t54; +fp t49; +fp t47; +fp t45; +fp t44; +fp t42; + t50 = RATIONAL(1.0,2.0); + t56 = y*t50; + t53 = x*x; + t43 = t53*t56; + t51 = RATIONAL(1.0,4.0); + t55 = t43+t51*x; + t52 = RATIONAL(-1.0,4.0); + t54 = t43+t52*x; + t49 = RATIONAL(-1.0,2.0); + t47 = t52*t53; + t45 = t51*t53; + t44 = x*t56; + t42 = t49*x*y; + coeffs_dy->coeff_m1_m1 = t42+t47+t55; + coeffs_dy->coeff_0_m1 = t49+y+(-y+t50)*t53; + coeffs_dy->coeff_p1_m1 = t44+t47+t54; + coeffs_dy->coeff_m1_0 = (x-t53)*y; + coeffs_dy->coeff_0_0 = (RATIONAL(-2.0,1.0)+RATIONAL(2.0,1.0)*t53)*y; + coeffs_dy->coeff_p1_0 = (-x-t53)*y; + coeffs_dy->coeff_m1_p1 = t45+t42+t54; + coeffs_dy->coeff_0_p1 = t50+y+(t49-y)*t53; + coeffs_dy->coeff_p1_p1 = t44+t45+t55; |