aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c')
-rw-r--r--src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c510
1 files changed, 253 insertions, 257 deletions
diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c
index 46af5f2..74e206a 100644
--- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c
+++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c
@@ -1,261 +1,257 @@
-fp t779;
-fp t782;
-fp t742;
-fp t740;
-fp t851;
-fp t781;
-fp t741;
-fp t739;
-fp t850;
-fp t849;
-fp t774;
-fp t778;
-fp t769;
-fp t848;
-fp t758;
-fp t752;
-fp t847;
-fp t759;
-fp t800;
-fp t743;
-fp t846;
-fp t760;
-fp t720;
-fp t750;
-fp t802;
-fp t845;
-fp t821;
-fp t761;
-fp t748;
-fp t717;
-fp t780;
-fp t715;
-fp t812;
-fp t807;
-fp t844;
-fp t745;
-fp t714;
-fp t770;
-fp t843;
-fp t811;
-fp t754;
-fp t804;
-fp t842;
-fp t819;
-fp t798;
-fp t771;
-fp t735;
-fp t824;
-fp t841;
-fp t815;
-fp t840;
-fp t809;
-fp t839;
-fp t731;
-fp t801;
-fp t838;
-fp t777;
-fp t810;
-fp t822;
-fp t823;
-fp t837;
-fp t818;
-fp t836;
-fp t757;
-fp t776;
-fp t820;
-fp t772;
-fp t803;
-fp t835;
-fp t834;
-fp t833;
-fp t744;
-fp t832;
-fp t831;
-fp t830;
-fp t766;
-fp t799;
-fp t716;
-fp t829;
-fp t765;
-fp t828;
-fp t827;
-fp t718;
-fp t826;
-fp t825;
-fp t816;
-fp t817;
-fp t813;
-fp t773;
-fp t806;
-fp t805;
-fp t751;
-fp t797;
-fp t764;
-fp t749;
-fp t796;
-fp t746;
-fp t795;
-fp t767;
-fp t768;
-fp t794;
-fp t753;
-fp t793;
-fp t792;
-fp t791;
-fp t790;
-fp t789;
-fp t788;
-fp t787;
-fp t786;
-fp t763;
-fp t762;
-fp t747;
- t779 = y*y;
- t782 = t779*t779;
- t742 = RATIONAL(1.0,288.0);
- t740 = RATIONAL(-1.0,288.0);
- t851 = t740*t782+t742*t779;
- t781 = t779*y;
- t741 = RATIONAL(-1.0,144.0);
- t739 = RATIONAL(1.0,144.0);
- t850 = t741*t781+t739*y;
- t849 = t741*y+t739*t781;
- t774 = RATIONAL(1.0,3.0);
- t778 = x*x;
- t769 = RATIONAL(-4.0,9.0);
- t848 = t769+t774*t778;
- t758 = RATIONAL(-1.0,12.0);
- t752 = RATIONAL(-1.0,48.0);
- t847 = t752*t782+t758;
- t759 = RATIONAL(1.0,24.0);
- t800 = t759*t778;
- t743 = RATIONAL(1.0,18.0);
- t846 = t743*y+t781*t800;
- t760 = RATIONAL(-1.0,18.0);
- t720 = t760*y;
- t750 = RATIONAL(-1.0,24.0);
- t802 = t778*t750;
- t845 = t720+t781*t802;
- t821 = x*t743;
- t761 = RATIONAL(1.0,36.0);
- t748 = RATIONAL(-1.0,36.0);
- t717 = t748*t782;
- t780 = t778*x;
- t715 = t760*t779;
- t812 = t780*t717+x*t715;
- t807 = t779*t780;
- t844 = t761*t807+t812+t782*t821;
- t745 = RATIONAL(1.0,72.0);
- t714 = t745*t782;
- t770 = RATIONAL(1.0,9.0);
- t843 = x*t714+t812+t770*t807;
- t811 = t778*t782;
- t754 = RATIONAL(1.0,48.0);
- t804 = t778*t754;
- t842 = t761*t782+t752*t811+(t804+t748)*t779;
- t819 = y*t780;
- t798 = t780*t743;
- t771 = RATIONAL(-1.0,9.0);
- t735 = t771*t781;
- t824 = x*y;
- t841 = t760*t819+t770*t824+t781*t798+x*t735;
- t815 = t780*t782;
- t840 = t741*t807+t739*t815+t851*x;
- t809 = t778*t779;
- t839 = t717+t761*t779+t782*t804+t752*t809;
- t731 = t770*t781;
- t801 = t781*t780;
- t838 = y*t798+t760*t801+t771*t824+x*t731;
- t777 = RATIONAL(-2.0,9.0);
- t810 = t780*t770;
- t822 = x*RATIONAL(8.0,9.0);
- t823 = x*t782;
- t837 = t779*t822+t769*t807+t777*t823+t782*t810;
- t818 = x*t781;
- t836 = y*t822+t769*t819+t777*t818+t770*t801;
- t757 = RATIONAL(1.0,12.0);
- t776 = RATIONAL(-1.0,3.0);
- t820 = y*t778;
- t772 = RATIONAL(4.0,9.0);
- t803 = t778*t781;
- t835 = t757*t803+t776*t820+t735+t772*y;
- t834 = t731+t758*t803+t848*y;
- t833 = t772*t819+RATIONAL(2.0,9.0)*t818+t771*t801+RATIONAL(-8.0,9.0)*t824
+fp t935;
+fp t933;
+fp t894;
+fp t897;
+fp t1007;
+fp t920;
+fp t934;
+fp t936;
+fp t1006;
+fp t937;
+fp t895;
+fp t896;
+fp t1005;
+fp t1004;
+fp t932;
+fp t1003;
+fp t906;
+fp t907;
+fp t1002;
+fp t899;
+fp t956;
+fp t913;
+fp t1001;
+fp t911;
+fp t873;
+fp t955;
+fp t903;
+fp t1000;
+fp t908;
+fp t879;
+fp t898;
+fp t999;
+fp t919;
+fp t964;
+fp t909;
+fp t870;
+fp t965;
+fp t871;
+fp t971;
+fp t998;
+fp t910;
+fp t874;
+fp t972;
+fp t997;
+fp t927;
+fp t974;
+fp t954;
+fp t889;
+fp t996;
+fp t961;
+fp t995;
+fp t922;
+fp t994;
+fp t923;
+fp t973;
+fp t924;
+fp t892;
+fp t905;
+fp t993;
+fp t975;
+fp t992;
+fp t976;
+fp t978;
+fp t991;
+fp t925;
+fp t990;
+fp t968;
+fp t989;
+fp t988;
+fp t928;
+fp t960;
+fp t958;
+fp t987;
+fp t912;
+fp t930;
+fp t986;
+fp t985;
+fp t984;
+fp t876;
+fp t983;
+fp t869;
+fp t982;
+fp t981;
+fp t980;
+fp t979;
+fp t966;
+fp t970;
+fp t963;
+fp t921;
+fp t959;
+fp t957;
+fp t953;
+fp t929;
+fp t914;
+fp t952;
+fp t904;
+fp t951;
+fp t916;
+fp t950;
+fp t949;
+fp t948;
+fp t947;
+fp t946;
+fp t945;
+fp t944;
+fp t943;
+fp t942;
+fp t941;
+fp t918;
+fp t917;
+fp t915;
+ t935 = y*y;
+ t933 = t935*t935;
+ t894 = RATIONAL(-1.0,288.0);
+ t897 = RATIONAL(1.0,288.0);
+ t1007 = t897*t935+t894*t933;
+ t920 = RATIONAL(-4.0,9.0);
+ t934 = x*x;
+ t936 = t934*x;
+ t1006 = t920*t936+x*RATIONAL(8.0,9.0);
+ t937 = y*t935;
+ t895 = RATIONAL(1.0,144.0);
+ t896 = RATIONAL(-1.0,144.0);
+ t1005 = t896*t937+t895*y;
+ t1004 = t895*t937+t896*y;
+ t932 = RATIONAL(-2.0,3.0);
+ t1003 = t932+RATIONAL(5.0,6.0)*t935;
+ t906 = RATIONAL(-1.0,48.0);
+ t907 = RATIONAL(-1.0,12.0);
+ t1002 = t906*t933+t907;
+ t899 = RATIONAL(1.0,48.0);
+ t956 = t934*t899;
+ t913 = RATIONAL(-1.0,36.0);
+ t1001 = t913+t956;
+ t911 = RATIONAL(-1.0,18.0);
+ t873 = t911*y;
+ t955 = t934*t937;
+ t903 = RATIONAL(-1.0,24.0);
+ t1000 = t873+t903*t955;
+ t908 = RATIONAL(1.0,18.0);
+ t879 = t908*y;
+ t898 = RATIONAL(1.0,24.0);
+ t999 = t898*t955+t879;
+ t919 = RATIONAL(1.0,9.0);
+ t964 = t919*t936;
+ t909 = RATIONAL(1.0,72.0);
+ t870 = t909*t933;
+ t965 = t933*t936;
+ t871 = t911*t935;
+ t971 = x*t871+t913*t965;
+ t998 = t935*t964+x*t870+t971;
+ t910 = RATIONAL(1.0,36.0);
+ t874 = t910*t935;
+ t972 = x*t933;
+ t997 = t936*t874+t971+t908*t972;
+ t927 = RATIONAL(-1.0,9.0);
+ t974 = y*t927;
+ t954 = t936*t937;
+ t889 = t919*t937;
+ t996 = t911*t954+t936*t879+(t974+t889)*x;
+ t961 = t934*t898;
+ t995 = t908*t937+y*t961+t1000;
+ t922 = RATIONAL(-2.0,9.0);
+ t994 = t922*t972+t933*t964+t1006*t935;
+ t923 = RATIONAL(4.0,9.0);
+ t973 = y*t934;
+ t924 = RATIONAL(-1.0,3.0);
+ t892 = t927*t937;
+ t905 = RATIONAL(1.0,12.0);
+ t993 = t905*t955+t924*t973+t923*y+t892;
+ t975 = x*t937;
+ t992 = t922*t975+t919*t954+t1006*y;
+ t976 = y*t936;
+ t978 = x*y;
+ t991 = t927*t954+t923*t976+RATIONAL(2.0,9.0)*t975+RATIONAL(-8.0,9.0)*t978
;
- t744 = RATIONAL(-1.0,72.0);
- t832 = t744*t819+t745*t801+t850*x;
- t831 = t754*t803+t752*t820+t850;
- t830 = t744*t801+t745*t819+t849*x;
- t766 = RATIONAL(-1.0,6.0);
- t799 = t778*t766;
- t716 = t744*t781;
- t829 = y*t799+t716+t846;
- t765 = RATIONAL(1.0,6.0);
- t828 = t765*t809+t714+t715+t782*t802;
- t827 = (t800+t744)*t782+(t743+t799)*t779;
- t718 = t745*t781;
- t826 = t718+t765*t820+t845;
- t825 = t752*t803+y*t804+t849;
- t816 = RATIONAL(5.0,48.0)*t779;
- t817 = t759*t815+x*t816;
- t813 = t766*t782;
- t773 = RATIONAL(-5.0,3.0);
- t806 = t780*t813+t773*x*t779;
- t805 = t778*RATIONAL(-1.0,96.0);
- t751 = RATIONAL(-5.0,48.0);
- t797 = t751*t823+t750*t807+t817;
- t764 = RATIONAL(2.0,3.0);
- t749 = RATIONAL(5.0,12.0);
- t796 = t749*t823+t764*t807+t806;
- t746 = RATIONAL(-5.0,24.0);
- t795 = t817+t765*t780+t746*t807+t847*x;
- t767 = RATIONAL(5.0,6.0);
- t768 = RATIONAL(-2.0,3.0);
- t794 = t774*t823+RATIONAL(4.0,3.0)*x+t768*t780+t767*t807+t806;
- t753 = RATIONAL(1.0,96.0);
- t793 = t753*t811+t779*t805+t840+t851;
- t792 = t760*t781+y*t802+t844+t846;
- t791 = y*t810+t748*t801+t843+(t718+t720)*x;
- t790 = t776*t809+t771*t782+t772*t779+t757*t811+t837;
- t789 = t761*t801+y*t821+x*t716+t771*t819+t843;
- t788 = t753*t809+t740*t779+t840+(t742+t805)*t782;
- t787 = t837+t758*t811+t770*t782+t848*t779;
- t786 = t743*t781+y*t800+t844+t845;
- t763 = RATIONAL(1.0,4.0);
- t762 = RATIONAL(-5.0,8.0);
- t747 = RATIONAL(5.0,24.0);
- coeffs_dx->coeff_m2_m2 = t788+t830+t831;
- coeffs_dx->coeff_m1_m2 = t786+t839+t841;
- coeffs_dx->coeff_0_m2 = (t757*y+t758*t781)*t780+(t746*y+t747*t781)*x+t797
+ t925 = RATIONAL(1.0,3.0);
+ t990 = t907*t955+t925*t973+t920*y+t889;
+ t968 = t935*t936;
+ t989 = t895*t965+t896*t968+t1007*x;
+ t988 = x*t892+t908*t954+t919*t978+t911*t976;
+ t928 = RATIONAL(1.0,6.0);
+ t960 = t934*t903;
+ t958 = t934*t935;
+ t987 = t928*t958+t870+t933*t960+t871;
+ t912 = RATIONAL(-1.0,72.0);
+ t930 = RATIONAL(-1.0,6.0);
+ t986 = t908*t935+t930*t958+(t961+t912)*t933;
+ t985 = t912*t954+t909*t976+t1004*x;
+ t984 = t899*t955+t906*t973+t1005;
+ t876 = t909*t937;
+ t983 = y*t964+t913*t954+(t873+t876)*x;
+ t869 = t912*t937;
+ t982 = t910*t954+t936*t974+(t879+t869)*x;
+ t981 = y*t960+t911*t937+t999;
+ t980 = y*t956+t906*t955+t1004;
+ t979 = t909*t954+t912*t976+t1005*x;
+ t966 = RATIONAL(5.0,48.0)*t935;
+ t970 = t898*t965+x*t966;
+ t963 = t933*t934;
+ t921 = RATIONAL(-5.0,3.0);
+ t959 = t930*t965+t921*x*t935;
+ t957 = t934*RATIONAL(1.0,96.0);
+ t953 = t934*RATIONAL(-1.0,96.0);
+ t929 = RATIONAL(2.0,3.0);
+ t914 = RATIONAL(5.0,12.0);
+ t952 = t929*t968+t914*t972+t959;
+ t904 = RATIONAL(-5.0,48.0);
+ t951 = t904*t972+t903*t968+t970;
+ t916 = RATIONAL(-5.0,24.0);
+ t950 = t970+t916*t968+t928*t936+t1002*x;
+ t949 = t925*t972+RATIONAL(4.0,3.0)*x+t959+t1003*t936;
+ t948 = t997+t874+t906*t958+t1001*t933;
+ t947 = t930*t973+t869+t998+t999;
+ t946 = t989+(t957+t894)*t935+(t953+t897)*t933;
+ t945 = t997+t910*t933+t906*t963+t1001*t935;
+ t944 = t905*t963+t923*t935+t927*t933+t924*t958+t994;
+ t943 = t933*t957+t935*t953+t989+t1007;
+ t942 = t925*t958+t919*t933+t920*t935+t907*t963+t994;
+ t941 = t876+t928*t973+t998+t1000;
+ t918 = RATIONAL(1.0,4.0);
+ t917 = RATIONAL(-5.0,8.0);
+ t915 = RATIONAL(5.0,24.0);
+ coeffs_dx->coeff_m2_m2 = t946+t984+t985;
+ coeffs_dx->coeff_m1_m2 = t948+t988+t995;
+ coeffs_dx->coeff_0_m2 = (t907*t937+t905*y)*t936+(t916*y+t915*t937)*x+t951
;
- coeffs_dx->coeff_p1_m2 = t792+t841+t842;
- coeffs_dx->coeff_p2_m2 = t793+t825+t830;
- coeffs_dx->coeff_m2_m1 = t789+t826+t827;
- coeffs_dx->coeff_m1_m1 = t787+t833+t835;
- coeffs_dx->coeff_0_m1 = (t765*t781+t768*y)*t780+(RATIONAL(5.0,3.0)*y+
-RATIONAL(-5.0,12.0)*t781)*x+t796;
- coeffs_dx->coeff_p1_m1 = t790+t833+t834;
- coeffs_dx->coeff_p2_m1 = t789+t828+t829;
- coeffs_dx->coeff_m2_0 = t757+t754*t782+t751*t779+(RATIONAL(5.0,16.0)*t779
-+RATIONAL(-1.0,4.0)+RATIONAL(-1.0,16.0)*t782)*t778+t795;
- coeffs_dx->coeff_m1_0 = t768+t813+t767*t779+(RATIONAL(1.0,2.0)+t762*t779+
-RATIONAL(1.0,8.0)*t782)*t778+t794;
- coeffs_dx->coeff_0_0 = (1.0+t763*t782+RATIONAL(-5.0,4.0)*t779)*t780+(
-RATIONAL(25.0,8.0)*t779+RATIONAL(-5.0,2.0)+t762*t782)*x;
- coeffs_dx->coeff_p1_0 = RATIONAL(-5.0,6.0)*t779+t764+t765*t782+(RATIONAL(
-5.0,8.0)*t779+RATIONAL(-1.0,2.0)+RATIONAL(-1.0,8.0)*t782)*t778+t794;
- coeffs_dx->coeff_p2_0 = t816+(t763+RATIONAL(-5.0,16.0)*t779+RATIONAL(1.0,
-16.0)*t782)*t778+t795+t847;
- coeffs_dx->coeff_m2_p1 = t791+t827+t829;
- coeffs_dx->coeff_m1_p1 = t787+t834+t836;
- coeffs_dx->coeff_0_p1 = (t766*t781+t764*y)*t780+(t773*y+t749*t781)*x+t796
+ coeffs_dx->coeff_p1_m2 = t945+t981+t988;
+ coeffs_dx->coeff_p2_m2 = t943+t980+t985;
+ coeffs_dx->coeff_m2_m1 = t941+t982+t986;
+ coeffs_dx->coeff_m1_m1 = t942+t991+t993;
+ coeffs_dx->coeff_0_m1 = (t928*t937+t932*y)*t936+(RATIONAL(-5.0,12.0)*t937
++RATIONAL(5.0,3.0)*y)*x+t952;
+ coeffs_dx->coeff_p1_m1 = t944+t990+t991;
+ coeffs_dx->coeff_p2_m1 = t947+t982+t987;
+ coeffs_dx->coeff_m2_0 = t899*t933+t905+t904*t935+(RATIONAL(-1.0,4.0)+
+RATIONAL(-1.0,16.0)*t933+RATIONAL(5.0,16.0)*t935)*t934+t950;
+ coeffs_dx->coeff_m1_0 = t930*t933+(t917*t935+RATIONAL(1.0,2.0)+RATIONAL(
+1.0,8.0)*t933)*t934+t949+t1003;
+ coeffs_dx->coeff_0_0 = (t918*t933+1.0+RATIONAL(-5.0,4.0)*t935)*t936+(
+RATIONAL(-5.0,2.0)+t917*t933+RATIONAL(25.0,8.0)*t935)*x;
+ coeffs_dx->coeff_p1_0 = t929+t928*t933+RATIONAL(-5.0,6.0)*t935+(RATIONAL(
+-1.0,2.0)+RATIONAL(5.0,8.0)*t935+RATIONAL(-1.0,8.0)*t933)*t934+t949;
+ coeffs_dx->coeff_p2_0 = t966+(RATIONAL(1.0,16.0)*t933+t918+RATIONAL(-5.0,
+16.0)*t935)*t934+t950+t1002;
+ coeffs_dx->coeff_m2_p1 = t947+t983+t986;
+ coeffs_dx->coeff_m1_p1 = t942+t990+t992;
+ coeffs_dx->coeff_0_p1 = (t929*y+t930*t937)*t936+(t921*y+t914*t937)*x+t952
;
- coeffs_dx->coeff_p1_p1 = t790+t835+t836;
- coeffs_dx->coeff_p2_p1 = t791+t826+t828;
- coeffs_dx->coeff_m2_p2 = t788+t825+t832;
- coeffs_dx->coeff_m1_p2 = t792+t838+t839;
- coeffs_dx->coeff_0_p2 = (t757*t781+t758*y)*t780+(t747*y+t746*t781)*x+t797
+ coeffs_dx->coeff_p1_p1 = t944+t992+t993;
+ coeffs_dx->coeff_p2_p1 = t941+t983+t987;
+ coeffs_dx->coeff_m2_p2 = t946+t979+t980;
+ coeffs_dx->coeff_m1_p2 = t948+t981+t996;
+ coeffs_dx->coeff_0_p2 = (t907*y+t905*t937)*t936+(t916*t937+t915*y)*x+t951
;
- coeffs_dx->coeff_p1_p2 = t786+t838+t842;
- coeffs_dx->coeff_p2_p2 = t793+t831+t832;
+ coeffs_dx->coeff_p1_p2 = t945+t995+t996;
+ coeffs_dx->coeff_p2_p2 = t943+t979+t984;