diff options
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.c | 510 |
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; |