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
---
 .../3d.cube.order3.smooth0/coeffs-dx.compute.c     | 518 +++++++++++++++++++++
 1 file changed, 518 insertions(+)
 create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeffs-dx.compute.c

(limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeffs-dx.compute.c')

diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeffs-dx.compute.c
new file mode 100644
index 0000000..1753d9b
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order3.smooth0/coeffs-dx.compute.c
@@ -0,0 +1,518 @@
+fp t911;
+fp t898;
+fp t808;
+fp t913;
+fp t809;
+fp t1006;
+fp t912;
+fp t826;
+fp t906;
+fp t804;
+fp t962;
+fp t1005;
+fp t840;
+fp t1004;
+fp t888;
+fp t797;
+fp t755;
+fp t1003;
+fp t829;
+fp t1002;
+fp t903;
+fp t812;
+fp t803;
+fp t837;
+fp t1001;
+fp t834;
+fp t805;
+fp t963;
+fp t1000;
+fp t999;
+fp t998;
+fp t997;
+fp t833;
+fp t996;
+fp t835;
+fp t995;
+fp t832;
+fp t759;
+fp t994;
+fp t993;
+fp t902;
+fp t816;
+fp t958;
+fp t813;
+fp t992;
+fp t838;
+fp t961;
+fp t991;
+fp t827;
+fp t990;
+fp t989;
+fp t822;
+fp t988;
+fp t987;
+fp t986;
+fp t985;
+fp t984;
+fp t983;
+fp t982;
+fp t886;
+fp t783;
+fp t889;
+fp t799;
+fp t894;
+fp t793;
+fp t758;
+fp t760;
+fp t981;
+fp t757;
+fp t762;
+fp t980;
+fp t979;
+fp t978;
+fp t761;
+fp t764;
+fp t976;
+fp t763;
+fp t975;
+fp t974;
+fp t973;
+fp t971;
+fp t970;
+fp t969;
+fp t967;
+fp t966;
+fp t965;
+fp t964;
+fp t960;
+fp t957;
+fp t955;
+fp t954;
+fp t751;
+fp t953;
+fp t756;
+fp t952;
+fp t951;
+fp t950;
+fp t830;
+fp t949;
+fp t752;
+fp t948;
+fp t947;
+fp t946;
+fp t836;
+fp t945;
+fp t944;
+fp t943;
+fp t753;
+fp t942;
+fp t941;
+fp t940;
+fp t939;
+fp t938;
+fp t749;
+fp t937;
+fp t754;
+fp t936;
+fp t935;
+fp t934;
+fp t933;
+fp t932;
+fp t750;
+fp t931;
+fp t930;
+fp t929;
+fp t928;
+fp t927;
+fp t926;
+fp t925;
+fp t924;
+fp t923;
+fp t922;
+fp t921;
+fp t920;
+fp t841;
+fp t919;
+fp t918;
+fp t824;
+fp t917;
+fp t916;
+fp t915;
+fp t914;
+fp t905;
+fp t904;
+fp t901;
+fp t900;
+fp t899;
+fp t896;
+fp t895;
+fp t893;
+fp t891;
+fp t890;
+fp t887;
+fp t885;
+fp t884;
+fp t880;
+fp t879;
+fp t878;
+fp t877;
+fp t875;
+fp t874;
+fp t871;
+fp t870;
+fp t867;
+fp t866;
+fp t865;
+fp t863;
+fp t861;
+fp t860;
+fp t859;
+fp t855;
+fp t854;
+fp t852;
+fp t849;
+fp t848;
+fp t846;
+fp t823;
+fp t821;
+fp t820;
+fp t819;
+fp t818;
+fp t817;
+fp t815;
+fp t814;
+fp t811;
+fp t810;
+fp t807;
+fp t806;
+fp t802;
+fp t801;
+fp t800;
+fp t798;
+fp t796;
+fp t795;
+fp t794;
+fp t792;
+fp t791;
+fp t790;
+fp t789;
+fp t788;
+fp t787;
+fp t786;
+fp t785;
+fp t784;
+fp t782;
+fp t781;
+fp t780;
+fp t779;
+fp t778;
+fp t777;
+fp t776;
+fp t775;
+fp t774;
+fp t773;
+fp t772;
+fp t771;
+fp t770;
+fp t769;
+fp t768;
+fp t767;
+fp t766;
+fp t765;
+      t911 = z*z;
+      t898 = RATIONAL(3.0,160.0);
+      t808 = t898*t911;
+      t913 = y*y;
+      t809 = t898*t913;
+      t1006 = t808+t809;
+      t912 = x*x;
+      t826 = RATIONAL(-1.0,32.0)*t912;
+      t906 = RATIONAL(-3.0,160.0);
+      t804 = t906*t913;
+      t962 = t808+t804;
+      t1005 = t826+t962;
+      t840 = RATIONAL(3.0,40.0)*x;
+      t1004 = t840+RATIONAL(-97.0,3000.0);
+      t888 = RATIONAL(-3.0,1000.0);
+      t797 = t888*y;
+      t755 = z*t797;
+      t1003 = RATIONAL(53.0,1000.0)+t840+t755;
+      t829 = RATIONAL(-1.0,80.0)*x;
+      t1002 = RATIONAL(-167.0,12000.0)+t829;
+      t903 = RATIONAL(-1.0,160.0);
+      t812 = t903*t913;
+      t803 = t903*t911;
+      t837 = RATIONAL(-3.0,32.0)*t912;
+      t1001 = t812+t837+t803;
+      t834 = RATIONAL(1.0,32.0)*t912;
+      t805 = t906*t911;
+      t963 = t804+t805;
+      t1000 = t834+t963;
+      t999 = t834+t1006;
+      t998 = RATIONAL(-79.0,2000.0)+RATIONAL(-1.0,10.0)*x;
+      t997 = t834+t962;
+      t833 = RATIONAL(-9.0,80.0)*x;
+      t996 = RATIONAL(-137.0,4000.0)+t833;
+      t835 = RATIONAL(1.0,40.0)*x;
+      t995 = RATIONAL(-17.0,3000.0)+t835;
+      t832 = RATIONAL(3.0,80.0)*x;
+      t759 = y*t832;
+      t994 = t759+t832;
+      t993 = RATIONAL(-49.0,2000.0)+RATIONAL(-3.0,20.0)*x;
+      t902 = RATIONAL(1.0,160.0);
+      t816 = t902*t913;
+      t958 = t816+t837;
+      t813 = t902*t911;
+      t992 = t958+t813;
+      t838 = RATIONAL(3.0,32.0)*t912;
+      t961 = t803+t838;
+      t991 = t812+t961;
+      t827 = RATIONAL(1.0,80.0)*x;
+      t990 = RATIONAL(-137.0,12000.0)+t827;
+      t989 = t826+t1006;
+      t822 = RATIONAL(-11.0,80.0)*x;
+      t988 = RATIONAL(-107.0,4000.0)+t822;
+      t987 = RATIONAL(-4.0,375.0)+RATIONAL(-1.0,40.0)*x;
+      t986 = t816+t813+t838;
+      t985 = t826+t963;
+      t984 = y+z;
+      t983 = y*z;
+      t982 = RATIONAL(-3.0,80.0)*x;
+      t886 = RATIONAL(-27.0,1000.0);
+      t783 = t886*y;
+      t889 = RATIONAL(-9.0,1000.0);
+      t799 = t889*y;
+      t894 = RATIONAL(-1.0,1000.0);
+      t793 = t894*z;
+      t758 = z*t832;
+      t760 = y*t829;
+      t981 = t758+t760;
+      t757 = y*t982;
+      t762 = z*t827;
+      t980 = t757+t762;
+      t979 = t816+t803;
+      t978 = t809+t805;
+      t761 = y*t827;
+      t764 = z*t982;
+      t976 = t761+t764;
+      t763 = z*t829;
+      t975 = t760+t763;
+      t974 = t758+t761;
+      t973 = t757+t764;
+      t971 = t757+t763;
+      t970 = t759+t763;
+      t969 = t759+t762;
+      t967 = t760+t764;
+      t966 = t761+t763;
+      t965 = t758+t759;
+      t964 = t760+t762;
+      t960 = t759+t764;
+      t957 = t761+t762;
+      t955 = t757+t758;
+      t954 = t812+t813;
+      t751 = z*t799;
+      t953 = t751+t978;
+      t756 = RATIONAL(3.0,1000.0)*t983;
+      t952 = t756+t954;
+      t951 = RATIONAL(1.0,10.0)*x+t973;
+      t950 = t832+t974;
+      t830 = RATIONAL(7.0,80.0)*x;
+      t949 = t830+t971;
+      t752 = RATIONAL(9.0,1000.0)*t983;
+      t948 = t752+t978;
+      t947 = t840+t975;
+      t946 = t756+t979;
+      t836 = RATIONAL(1.0,20.0)*x;
+      t945 = t836+t957;
+      t944 = t835+t965;
+      t943 = t751+t991;
+      t753 = y*t793;
+      t942 = t753+t986;
+      t941 = t838+t952;
+      t940 = t755+t989;
+      t939 = t762+t803+t958;
+      t938 = t834+t953;
+      t749 = z*t783;
+      t937 = t749+t985;
+      t754 = RATIONAL(1.0,1000.0)*t983;
+      t936 = t754+t1001;
+      t935 = t764+t997;
+      t934 = t762+t1005;
+      t933 = t758+t997;
+      t932 = t755+t838+t954;
+      t750 = RATIONAL(27.0,1000.0)*t983;
+      t931 = t750+t999;
+      t930 = t755+t816+t961;
+      t929 = t756+t1000;
+      t928 = t752+t826+t836+RATIONAL(-29.0,6000.0);
+      t927 = t755+RATIONAL(141.0,2000.0)+t836+t837;
+      t926 = t764+t837+t954;
+      t925 = t826+t953;
+      t924 = t834+t948;
+      t923 = t752+t992;
+      t922 = t838+t946;
+      t921 = t755+RATIONAL(79.0,3000.0)+t1000;
+      t920 = RATIONAL(-41.0,750.0)+t749+t999;
+      t841 = RATIONAL(-1.0,8.0)*x;
+      t919 = t752+RATIONAL(-33.0,2000.0)+t841+t991;
+      t918 = t841+RATIONAL(-87.0,2000.0)+t754+t986;
+      t824 = RATIONAL(1.0,16.0)*x;
+      t917 = t756+RATIONAL(-691.0,12000.0)+t824+t989;
+      t916 = t753+RATIONAL(299.0,4000.0)+t824+t1001;
+      t915 = RATIONAL(281.0,12000.0)+t750+t824+t985;
+      t914 = t751+t824+RATIONAL(191.0,4000.0)+t992;
+      t905 = RATIONAL(-7.0,250.0);
+      t904 = RATIONAL(-9.0,500.0);
+      t901 = RATIONAL(-9.0,250.0);
+      t900 = RATIONAL(-1.0,500.0);
+      t899 = RATIONAL(2.0,125.0);
+      t896 = RATIONAL(-1.0,250.0);
+      t895 = RATIONAL(-17.0,500.0);
+      t893 = RATIONAL(-3.0,2000.0);
+      t891 = RATIONAL(11.0,500.0);
+      t890 = RATIONAL(-1.0,2000.0);
+      t887 = RATIONAL(-7.0,2000.0);
+      t885 = RATIONAL(37.0,2000.0);
+      t884 = RATIONAL(13.0,1000.0);
+      t880 = RATIONAL(49.0,2000.0);
+      t879 = RATIONAL(27.0,2000.0);
+      t878 = RATIONAL(-37.0,1000.0);
+      t877 = RATIONAL(43.0,2000.0);
+      t875 = RATIONAL(31.0,2000.0);
+      t874 = RATIONAL(-19.0,2000.0);
+      t871 = RATIONAL(93.0,2000.0);
+      t870 = RATIONAL(-21.0,2000.0);
+      t867 = RATIONAL(19.0,1000.0);
+      t866 = RATIONAL(-23.0,2000.0);
+      t865 = RATIONAL(-39.0,2000.0);
+      t863 = RATIONAL(-17.0,2000.0);
+      t861 = RATIONAL(-31.0,1000.0);
+      t860 = RATIONAL(33.0,2000.0);
+      t859 = RATIONAL(-13.0,2000.0);
+      t855 = RATIONAL(183.0,4000.0);
+      t854 = RATIONAL(111.0,2000.0);
+      t852 = RATIONAL(129.0,2000.0);
+      t849 = RATIONAL(53.0,12000.0);
+      t848 = RATIONAL(323.0,4000.0);
+      t846 = RATIONAL(-607.0,12000.0);
+      t823 = t900*y;
+      t821 = t896*z;
+      t820 = t905*y;
+      t819 = t904*y;
+      t818 = t899*z;
+      t817 = t899*y;
+      t815 = t905*z;
+      t814 = t900*z;
+      t811 = t901*y;
+      t810 = t901*z;
+      t807 = t904*z;
+      t806 = t896*y;
+      t802 = t887*y;
+      t801 = t895*y;
+      t800 = t890*z;
+      t798 = t888*z;
+      t796 = t894*y;
+      t795 = t891*y;
+      t794 = t889*z;
+      t792 = t891*z;
+      t791 = t887*z;
+      t790 = t895*z;
+      t789 = t890*y;
+      t788 = t875*z;
+      t787 = t885*z;
+      t786 = t859*z;
+      t785 = t885*y;
+      t784 = t877*z;
+      t782 = t884*z;
+      t781 = t861*y;
+      t780 = t874*y;
+      t779 = t874*z;
+      t778 = t867*y;
+      t777 = t861*z;
+      t776 = t880*z;
+      t775 = t870*z;
+      t774 = t878*y;
+      t773 = t875*y;
+      t772 = t877*y;
+      t771 = t878*z;
+      t770 = t867*z;
+      t769 = t870*y;
+      t768 = t880*y;
+      t767 = t884*y;
+      t766 = t886*z;
+      t765 = t859*y;
+      coeffs_dx->coeff_m1_m1_m1 = RATIONAL(-47.0,1500.0)+t984*RATIONAL(147.0,
+2000.0)+t937+t951;
+      coeffs_dx->coeff_0_m1_m1 = RATIONAL(-39.0,4000.0)+t789+t800+RATIONAL(
+-13.0,80.0)*x+t943+t965;
+      coeffs_dx->coeff_p1_m1_m1 = RATIONAL(157.0,2000.0)+t774+t771+t923+t944;
+      coeffs_dx->coeff_p2_m1_m1 = t811+t832+t810+RATIONAL(-449.0,12000.0)+t931+
+t973;
+      coeffs_dx->coeff_m1_0_m1 = t830+t846+t852*z+t789+t925+t967;
+      coeffs_dx->coeff_0_0_m1 = t863*y+t791+t930+t974+t993;
+      coeffs_dx->coeff_p1_0_m1 = t837+t790+t806+t848+t950+t952;
+      coeffs_dx->coeff_p2_0_m1 = t766+t760+t767+t752+t935+t995;
+      coeffs_dx->coeff_m1_p1_m1 = t854*z+t774+t826+t948+t976+t1004;
+      coeffs_dx->coeff_0_p1_m1 = t786+t806+t922+t981+t988;
+      coeffs_dx->coeff_p1_p1_m1 = t777+t860*y+t927+t954+t981;
+      coeffs_dx->coeff_p2_p1_m1 = t761+t807+t751+t768+t935+t990;
+      coeffs_dx->coeff_m1_p2_m1 = t811+t871*z+t915+t960;
+      coeffs_dx->coeff_0_p2_m1 = t767+t779+t919+t955;
+      coeffs_dx->coeff_p1_p2_m1 = t815+t768+t914+t955;
+      coeffs_dx->coeff_p2_p2_m1 = t794+t893*y+t920+t960;
+      coeffs_dx->coeff_m1_m1_0 = t751+t852*y+t846+t800+t949+t1005;
+      coeffs_dx->coeff_0_m1_0 = t802+t863*z+t932+t969+t993;
+      coeffs_dx->coeff_p1_m1_0 = t821+t756+t801+t848+t939+t994;
+      coeffs_dx->coeff_p2_m1_0 = t783+t782+t924+t971+t995;
+      coeffs_dx->coeff_m1_0_0 = t791+RATIONAL(-437.0,6000.0)+t802+t940+t947;
+      coeffs_dx->coeff_0_0_0 = RATIONAL(-161.0,4000.0)+t822+t779+t780+t942+t957
+;
+      coeffs_dx->coeff_p1_0_0 = RATIONAL(21.0,250.0)+t797+t798+t936+t945;
+      coeffs_dx->coeff_p2_0_0 = t827+t817+RATIONAL(349.0,12000.0)+t818+t929+
+t975;
+      coeffs_dx->coeff_m1_p1_0 = t786+t801+t917+t966;
+      coeffs_dx->coeff_0_p1_0 = t775+t797+t918+t964;
+      coeffs_dx->coeff_p1_p1_0 = t773+t814+t916+t964;
+      coeffs_dx->coeff_p2_p1_0 = t772+t770+t921+t966;
+      coeffs_dx->coeff_m1_p2_0 = t783+t779+t928+t962+t970;
+      coeffs_dx->coeff_0_p2_0 = t817+t866*z+t941+t980+t996;
+      coeffs_dx->coeff_p1_p2_0 = t772+t757+t793+t939+t1003;
+      coeffs_dx->coeff_p2_p2_0 = t792+t769+t938+t970+t1002;
+      coeffs_dx->coeff_m1_m1_p1 = t752+t757+t854*y+t771+t934+t1004;
+      coeffs_dx->coeff_0_m1_p1 = t765+t821+t941+t970+t988;
+      coeffs_dx->coeff_p1_m1_p1 = t860*z+t781+t927+t970+t979;
+      coeffs_dx->coeff_p2_m1_p1 = t819+t776+t938+t980+t990;
+      coeffs_dx->coeff_m1_0_p1 = t790+t765+t917+t964;
+      coeffs_dx->coeff_0_0_p1 = t769+t798+t918+t966;
+      coeffs_dx->coeff_p1_0_p1 = t823+t788+t916+t966;
+      coeffs_dx->coeff_p2_0_p1 = t778+t784+t921+t964;
+      coeffs_dx->coeff_m1_p1_p1 = RATIONAL(-17.0,375.0)+t777+t781+t940+t945;
+      coeffs_dx->coeff_0_p1_p1 = RATIONAL(-191.0,4000.0)+t833+t814+t823+t942+
+t975;
+      coeffs_dx->coeff_p1_p1_p1 = RATIONAL(133.0,2000.0)+t984*RATIONAL(29.0,
+2000.0)+t936+t947;
+      coeffs_dx->coeff_p2_p1_p1 = RATIONAL(319.0,12000.0)+t829+t787+t785+t929+
+t957;
+      coeffs_dx->coeff_m1_p2_p1 = t849+t819+t751+t815+t934+t994;
+      coeffs_dx->coeff_0_p2_p1 = t793+t778+t932+t971+t998;
+      coeffs_dx->coeff_p1_p2_p1 = t855+t837+t879*z+t785+t946+t949;
+      coeffs_dx->coeff_p2_p2_p1 = t865*y+t788+t924+t969+t987;
+      coeffs_dx->coeff_m1_m1_p2 = t810+t871*y+t915+t955;
+      coeffs_dx->coeff_0_m1_p2 = t780+t782+t919+t960;
+      coeffs_dx->coeff_p1_m1_p2 = t776+t820+t914+t960;
+      coeffs_dx->coeff_p2_m1_p2 = t799+t893*z+t920+t955;
+      coeffs_dx->coeff_m1_0_p2 = t780+t766+t928+t978+t981;
+      coeffs_dx->coeff_0_0_p2 = t866*y+t818+t922+t976+t996;
+      coeffs_dx->coeff_p1_0_p2 = t796+t761+t784+t926+t1003;
+      coeffs_dx->coeff_p2_0_p2 = t775+t795+t760+t751+t933+t1002;
+      coeffs_dx->coeff_m1_p1_p2 = t849+t820+t807+t925+t950;
+      coeffs_dx->coeff_0_p1_p2 = t770+t796+t930+t967+t998;
+      coeffs_dx->coeff_p1_p1_p2 = t879*y+t760+t855+t830+t756+t787+t926;
+      coeffs_dx->coeff_p2_p1_p2 = t761+t752+t865*z+t773+t933+t987;
+      coeffs_dx->coeff_m1_p2_p2 = t794+RATIONAL(307.0,6000.0)+t799+t937+t944;
+      coeffs_dx->coeff_0_p2_p2 = t792+t795+RATIONAL(-7.0,80.0)*x+RATIONAL(
+-129.0,4000.0)+t943+t973;
+      coeffs_dx->coeff_p1_p2_p2 = RATIONAL(13.0,500.0)+t788+t773+t923+t951;
+      coeffs_dx->coeff_p2_p2_p2 = RATIONAL(-539.0,12000.0)+t982+t984*RATIONAL(
+-57.0,2000.0)+t931+t965;
-- 
cgit v1.2.3