From 80df27de110b6e739fda09ba4d9eae9fdbefc342 Mon Sep 17 00:00:00 2001
From: jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416>
Date: Mon, 2 Sep 2002 10:36:32 +0000
Subject: This is a major reworking of this interpolator. The interpolation
 code in template.c is now split up into subfunctions (in common/, Lagrange/,
 and Hermite/), so it should now compile in finite cpu/memory even with
 optimization.

git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@115 df1f8a13-aa1d-4dd4-9681-27ded5b42416
---
 .../2d.cube.order4.smooth0/coeffs-dyy.compute.c    | 107 +++++++++++++++++++++
 1 file changed, 107 insertions(+)
 create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c

(limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c')

diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c
new file mode 100644
index 0000000..29eef2a
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c
@@ -0,0 +1,107 @@
+fp t831;
+fp t817;
+fp t830;
+fp t848;
+fp t808;
+fp t826;
+fp t811;
+fp t847;
+fp t846;
+fp t813;
+fp t845;
+fp t844;
+fp t814;
+fp t816;
+fp t800;
+fp t812;
+fp t807;
+fp t843;
+fp t815;
+fp t842;
+fp t801;
+fp t841;
+fp t840;
+fp t839;
+fp t838;
+fp t806;
+fp t837;
+fp t809;
+fp t836;
+fp t835;
+fp t834;
+fp t833;
+fp t832;
+fp t810;
+fp t805;
+fp t799;
+fp t798;
+fp t797;
+fp t796;
+fp t795;
+fp t794;
+      t831 = y*y;
+      t817 = RATIONAL(-2.0,5.0)*t831;
+      t830 = x*x;
+      t848 = RATIONAL(1.0,98.0)*t830+t817+RATIONAL(181.0,735.0);
+      t808 = RATIONAL(-1.0,49.0)*t830;
+      t826 = RATIONAL(1.0,10.0);
+      t811 = t826*t831;
+      t847 = t808+t811+RATIONAL(71.0,2940.0);
+      t846 = t808+t817+RATIONAL(226.0,735.0);
+      t813 = RATIONAL(2.0,49.0)*t830;
+      t845 = t811+t813+RATIONAL(-289.0,2940.0);
+      t844 = x*y;
+      t814 = RATIONAL(-1.0,5.0)*y;
+      t816 = RATIONAL(1.0,5.0)*y;
+      t800 = t826*y;
+      t812 = RATIONAL(-1.0,10.0)*y;
+      t807 = RATIONAL(-2.0,49.0)*t830;
+      t843 = t807+t811+RATIONAL(191.0,2940.0);
+      t815 = RATIONAL(3.0,5.0)*t831;
+      t842 = t807+t815+RATIONAL(-41.0,98.0);
+      t801 = RATIONAL(1.0,49.0)*t830;
+      t841 = t801+t815+RATIONAL(-53.0,98.0);
+      t840 = t801+RATIONAL(166.0,735.0)+t817;
+      t839 = t800+t845;
+      t838 = t812+t845;
+      t806 = RATIONAL(1.0,35.0)*x;
+      t837 = t806+t847;
+      t809 = RATIONAL(-1.0,35.0)*x;
+      t836 = t809+t847;
+      t835 = t814+t846;
+      t834 = t816+t846;
+      t833 = RATIONAL(1.0,70.0)*x+t848;
+      t832 = RATIONAL(-1.0,70.0)*x+t848;
+      t810 = RATIONAL(-2.0,35.0)*x;
+      t805 = RATIONAL(2.0,35.0)*x;
+      t799 = x*t814;
+      t798 = x*t816;
+      t797 = RATIONAL(1.0,20.0)*t844;
+      t796 = x*t812;
+      t795 = RATIONAL(-1.0,20.0)*t844;
+      t794 = x*t800;
+      coeffs_dyy->coeff_m2_m2 = t810+t794+t838;
+      coeffs_dyy->coeff_m1_m2 = t797+t812+t836;
+      coeffs_dyy->coeff_0_m2 = t812+t843;
+      coeffs_dyy->coeff_p1_m2 = t812+t795+t837;
+      coeffs_dyy->coeff_p2_m2 = t796+t805+t838;
+      coeffs_dyy->coeff_m2_m1 = t799+t806+t834;
+      coeffs_dyy->coeff_m1_m1 = t796+t816+t833;
+      coeffs_dyy->coeff_0_m1 = t816+t840;
+      coeffs_dyy->coeff_p1_m1 = t794+t816+t832;
+      coeffs_dyy->coeff_p2_m1 = t809+t798+t834;
+      coeffs_dyy->coeff_m2_0 = t805+t842;
+      coeffs_dyy->coeff_m1_0 = t806+t841;
+      coeffs_dyy->coeff_0_0 = t813+t815+RATIONAL(-57.0,98.0);
+      coeffs_dyy->coeff_p1_0 = t809+t841;
+      coeffs_dyy->coeff_p2_0 = t810+t842;
+      coeffs_dyy->coeff_m2_p1 = t806+t798+t835;
+      coeffs_dyy->coeff_m1_p1 = t794+t814+t833;
+      coeffs_dyy->coeff_0_p1 = t814+t840;
+      coeffs_dyy->coeff_p1_p1 = t814+t796+t832;
+      coeffs_dyy->coeff_p2_p1 = t799+t809+t835;
+      coeffs_dyy->coeff_m2_p2 = t796+t810+t839;
+      coeffs_dyy->coeff_m1_p2 = t795+t800+t836;
+      coeffs_dyy->coeff_0_p2 = t800+t843;
+      coeffs_dyy->coeff_p1_p2 = t800+t797+t837;
+      coeffs_dyy->coeff_p2_p2 = t794+t805+t839;
-- 
cgit v1.2.3