aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c')
-rw-r--r--src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c530
1 files changed, 275 insertions, 255 deletions
diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c
index 59006b6..763a422 100644
--- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c
+++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c
@@ -1,255 +1,275 @@
-fp t1007;
-fp t1006;
-fp t1013;
-fp t1015;
-fp t964;
-fp t1082;
-fp t974;
-fp t976;
-fp t1081;
-fp t1080;
-fp t975;
-fp t977;
-fp t1017;
-fp t1079;
-fp t1004;
-fp t1005;
-fp t1078;
-fp t987;
-fp t993;
-fp t1077;
-fp t985;
-fp t1014;
-fp t1040;
-fp t995;
-fp t957;
-fp t1076;
-fp t994;
-fp t978;
-fp t952;
-fp t1075;
-fp t1052;
-fp t996;
-fp t1016;
-fp t1039;
-fp t983;
-fp t1042;
-fp t958;
-fp t1048;
-fp t1074;
-fp t980;
-fp t959;
-fp t1038;
-fp t1073;
-fp t972;
-fp t1054;
-fp t1050;
-fp t1043;
-fp t1072;
-fp t1071;
-fp t1049;
-fp t1010;
-fp t1070;
-fp t1053;
-fp t1009;
-fp t1069;
-fp t992;
-fp t1011;
-fp t1068;
-fp t1067;
-fp t979;
-fp t1066;
-fp t989;
-fp t1065;
-fp t955;
-fp t1064;
-fp t1001;
-fp t951;
-fp t1063;
-fp t1062;
-fp t1061;
-fp t1060;
-fp t1059;
-fp t1058;
-fp t1000;
-fp t1057;
-fp t1056;
-fp t1051;
-fp t1055;
-fp t1035;
-fp t1008;
-fp t1047;
-fp t1045;
-fp t1046;
-fp t1041;
-fp t1037;
-fp t1036;
-fp t1003;
-fp t973;
-fp t968;
-fp t984;
-fp t999;
-fp t1034;
-fp t986;
-fp t1033;
-fp t1032;
-fp t981;
-fp t1031;
-fp t1030;
-fp t1029;
-fp t1028;
-fp t988;
-fp t990;
-fp t1027;
-fp t1026;
-fp t1025;
-fp t1024;
-fp t1023;
-fp t1022;
-fp t1021;
-fp t998;
-fp t997;
- t1007 = RATIONAL(4.0,9.0);
- t1006 = RATIONAL(-1.0,9.0);
- t1013 = x*x;
- t1015 = t1013*x;
- t964 = t1006*t1015;
- t1082 = t964+t1007*x;
- t974 = RATIONAL(1.0,144.0);
- t976 = RATIONAL(-1.0,144.0);
- t1081 = t974*x+t976*t1015;
- t1080 = t974*t1015+t976*x;
- t975 = RATIONAL(-1.0,288.0);
- t977 = RATIONAL(1.0,288.0);
- t1017 = t1013*t1013;
- t1079 = t975*t1017+t977*t1013;
- t1004 = RATIONAL(-4.0,9.0);
- t1005 = RATIONAL(1.0,9.0);
- t1078 = t1004*t1013+t1005*t1017;
- t987 = RATIONAL(-1.0,48.0);
- t993 = RATIONAL(-1.0,12.0);
- t1077 = t987*t1017+t993;
- t985 = RATIONAL(-1.0,24.0);
- t1014 = y*y;
- t1040 = t1014*t1015;
- t995 = RATIONAL(-1.0,18.0);
- t957 = t995*x;
- t1076 = t985*t1040+t957;
- t994 = RATIONAL(1.0,24.0);
- t978 = RATIONAL(1.0,18.0);
- t952 = t978*x;
- t1075 = t994*t1040+t952;
- t1052 = y*t1017;
- t996 = RATIONAL(1.0,36.0);
- t1016 = t1014*y;
- t1039 = t1013*t1016;
- t983 = RATIONAL(-1.0,36.0);
- t1042 = t1017*t1016;
- t958 = t995*t1013;
- t1048 = t983*t1042+y*t958;
- t1074 = t978*t1052+t1048+t996*t1039;
- t980 = RATIONAL(1.0,72.0);
- t959 = t980*t1017;
- t1038 = t1005*t1016;
- t1073 = y*t959+t1048+t1013*t1038;
- t972 = t1005*t1015;
- t1054 = x*y;
- t1050 = x*t1016;
- t1043 = t1016*t1015;
- t1072 = t1006*t1054+y*t972+t978*t1050+t995*t1043;
- t1071 = RATIONAL(2.0,9.0)*y*t1015+RATIONAL(-8.0,9.0)*t1054+t1082*t1016;
- t1049 = y*RATIONAL(-2.0,9.0);
- t1010 = RATIONAL(8.0,9.0);
- t1070 = t1004*t1050+t1016*t972+t1010*t1054+t1015*t1049;
- t1053 = x*t1014;
- t1009 = RATIONAL(1.0,3.0);
- t1069 = t972+t1004*x+t993*t1040+t1009*t1053;
- t992 = RATIONAL(1.0,12.0);
- t1011 = RATIONAL(-1.0,3.0);
- t1068 = t1011*t1053+t992*t1040+t1082;
- t1067 = t995*t1050+y*t964+t978*t1043+t1005*t1054;
- t979 = RATIONAL(-1.0,72.0);
- t1066 = t980*t1043+t979*t1050+t1081*y;
- t989 = RATIONAL(1.0,48.0);
- t1065 = t987*t1053+t989*t1040+t1081;
- t955 = t980*t1015;
- t1064 = t983*t1043+x*t1038+(t957+t955)*y;
- t1001 = RATIONAL(-1.0,6.0);
- t951 = t979*t1015;
- t1063 = t1001*t1053+t951+t1075;
- t1062 = t996*t1043+t1006*t1050+(t952+t951)*y;
- t1061 = t979*t1043+t980*t1050+t1080*y;
- t1060 = t989*t1053+t987*t1040+t1080;
- t1059 = t978*t1015+t994*t1053+t1076;
- t1058 = t995*t1015+t985*t1053+t1075;
- t1000 = RATIONAL(1.0,6.0);
- t1057 = t955+t1000*t1053+t1076;
- t1056 = t974*t1042+t976*t1039+t1079*y;
- t1051 = y*t1013;
- t1055 = t1010*t1051+t1017*t1049+t1078*t1016;
- t1035 = t1001*t1017;
- t1008 = RATIONAL(-5.0,3.0);
- t1047 = t1008*t1051+t1016*t1035;
- t1045 = RATIONAL(5.0,48.0)*t1013;
- t1046 = t994*t1042+y*t1045;
- t1041 = t1017*t1014;
- t1037 = t1013*t1014;
- t1036 = RATIONAL(5.0,6.0)*t1013;
- t1003 = RATIONAL(-2.0,3.0);
- t973 = t1003*t1016;
- t968 = t1000*t1016;
- t984 = RATIONAL(5.0,12.0);
- t999 = RATIONAL(2.0,3.0);
- t1034 = t984*t1052+t999*t1039+t1047;
- t986 = RATIONAL(-5.0,48.0);
- t1033 = t986*t1052+t985*t1039+t1046;
- t1032 = RATIONAL(5.0,24.0)*y+t993*t1016;
- t981 = RATIONAL(-5.0,24.0);
- t1031 = t992*t1016+t981*y;
- t1030 = t973+t1009*t1052+t1016*t1036+RATIONAL(4.0,3.0)*y+t1047;
- t1029 = t968+t981*t1039+t1046+t1077*y;
- t1028 = t994*t1041+t1001*t1037+t978*t1013+t979*t1017+t1073;
- t988 = RATIONAL(1.0,96.0);
- t990 = RATIONAL(-1.0,96.0);
- t1027 = t988*t1037+t977*t1017+t975*t1013+t990*t1041+t1056;
- t1026 = t1007*t1013+t992*t1041+t1006*t1017+t1011*t1037+t1055;
- t1025 = t987*t1041+t983*t1013+t996*t1017+t989*t1037+t1074;
- t1024 = t993*t1041+t1009*t1037+t1055+t1078;
- t1023 = t987*t1037+t983*t1017+t996*t1013+t989*t1041+t1074;
- t1022 = t990*t1037+t988*t1041+t1056+t1079;
- t1021 = t959+t985*t1041+t1000*t1037+t958+t1073;
- t998 = RATIONAL(1.0,4.0);
- t997 = RATIONAL(-5.0,8.0);
- coeffs_dy->coeff_m2_m2 = t1027+t1061+t1065;
- coeffs_dy->coeff_m1_m2 = t1028+t1057+t1062;
- coeffs_dy->coeff_0_m2 = t986*t1013+t989*t1017+t992+(RATIONAL(5.0,16.0)*
-t1013+RATIONAL(-1.0,4.0)+RATIONAL(-1.0,16.0)*t1017)*t1014+t1029;
- coeffs_dy->coeff_p1_m2 = t1028+t1063+t1064;
- coeffs_dy->coeff_p2_m2 = t1027+t1060+t1066;
- coeffs_dy->coeff_m2_m1 = t1023+t1059+t1067;
- coeffs_dy->coeff_m1_m1 = t1024+t1068+t1071;
- coeffs_dy->coeff_0_m1 = t1036+t1003+t1035+(t997*t1013+RATIONAL(1.0,2.0)+
-RATIONAL(1.0,8.0)*t1017)*t1014+t1030;
- coeffs_dy->coeff_p1_m1 = t1024+t1069+t1070;
- coeffs_dy->coeff_p2_m1 = t1023+t1058+t1072;
- coeffs_dy->coeff_m2_0 = t1032*t1015+t1031*x+t1033;
- coeffs_dy->coeff_m1_0 = (RATIONAL(-5.0,12.0)*y+t968)*t1015+(RATIONAL(5.0,
-3.0)*y+t973)*x+t1034;
- coeffs_dy->coeff_0_0 = (RATIONAL(-5.0,4.0)*t1013+t998*t1017+1.0)*t1016+(
-t997*t1017+RATIONAL(25.0,8.0)*t1013+RATIONAL(-5.0,2.0))*y;
- coeffs_dy->coeff_p1_0 = (t1001*t1016+t984*y)*t1015+(t1008*y+t999*t1016)*x
-+t1034;
- coeffs_dy->coeff_p2_0 = t1031*t1015+t1032*x+t1033;
- coeffs_dy->coeff_m2_p1 = t1025+t1058+t1067;
- coeffs_dy->coeff_m1_p1 = t1026+t1069+t1071;
- coeffs_dy->coeff_0_p1 = t1000*t1017+t999+RATIONAL(-5.0,6.0)*t1013+(
-RATIONAL(5.0,8.0)*t1013+RATIONAL(-1.0,8.0)*t1017+RATIONAL(-1.0,2.0))*t1014+
-t1030;
- coeffs_dy->coeff_p1_p1 = t1026+t1068+t1070;
- coeffs_dy->coeff_p2_p1 = t1025+t1059+t1072;
- coeffs_dy->coeff_m2_p2 = t1022+t1060+t1061;
- coeffs_dy->coeff_m1_p2 = t1021+t1062+t1063;
- coeffs_dy->coeff_0_p2 = t1045+(RATIONAL(1.0,16.0)*t1017+RATIONAL(-5.0,
-16.0)*t1013+t998)*t1014+t1029+t1077;
- coeffs_dy->coeff_p1_p2 = t1021+t1057+t1064;
- coeffs_dy->coeff_p2_p2 = t1022+t1065+t1066;
+fp t1155;
+fp t1170;
+fp t1172;
+fp t1120;
+fp t1156;
+fp t1245;
+fp t1132;
+fp t1131;
+fp t1244;
+fp t1243;
+fp t1163;
+fp t1122;
+fp t1159;
+fp t1242;
+fp t1169;
+fp t1241;
+fp t1130;
+fp t1133;
+fp t1240;
+fp t1143;
+fp t1142;
+fp t1196;
+fp t1239;
+fp t1171;
+fp t1173;
+fp t1213;
+fp t1144;
+fp t1146;
+fp t1105;
+fp t1238;
+fp t1200;
+fp t1145;
+fp t1117;
+fp t1237;
+fp t1166;
+fp t1205;
+fp t1134;
+fp t1148;
+fp t1201;
+fp t1236;
+fp t1215;
+fp t1135;
+fp t1198;
+fp t1235;
+fp t1112;
+fp t1217;
+fp t1234;
+fp t1214;
+fp t1212;
+fp t1233;
+fp t1161;
+fp t1232;
+fp t1231;
+fp t1207;
+fp t1199;
+fp t1230;
+fp t1195;
+fp t1210;
+fp t1149;
+fp t1106;
+fp t1229;
+fp t1203;
+fp t1228;
+fp t1218;
+fp t1147;
+fp t1111;
+fp t1104;
+fp t1227;
+fp t1226;
+fp t1141;
+fp t1160;
+fp t1225;
+fp t1114;
+fp t1108;
+fp t1224;
+fp t1204;
+fp t1223;
+fp t1118;
+fp t1139;
+fp t1164;
+fp t1208;
+fp t1222;
+fp t1221;
+fp t1107;
+fp t1220;
+fp t1219;
+fp t1216;
+fp t1168;
+fp t1128;
+fp t1197;
+fp t1211;
+fp t1119;
+fp t1193;
+fp t1209;
+fp t1206;
+fp t1165;
+fp t1202;
+fp t1194;
+fp t1192;
+fp t1140;
+fp t1191;
+fp t1190;
+fp t1150;
+fp t1189;
+fp t1152;
+fp t1188;
+fp t1187;
+fp t1186;
+fp t1185;
+fp t1136;
+fp t1137;
+fp t1184;
+fp t1183;
+fp t1182;
+fp t1181;
+fp t1180;
+fp t1179;
+fp t1178;
+fp t1177;
+fp t1154;
+fp t1153;
+ t1155 = RATIONAL(1.0,9.0);
+ t1170 = x*x;
+ t1172 = x*t1170;
+ t1120 = t1155*t1172;
+ t1156 = RATIONAL(-4.0,9.0);
+ t1245 = t1156*x+t1120;
+ t1132 = RATIONAL(-1.0,144.0);
+ t1131 = RATIONAL(1.0,144.0);
+ t1244 = t1132*t1172+t1131*x;
+ t1243 = t1132*x+t1131*t1172;
+ t1163 = RATIONAL(-1.0,9.0);
+ t1122 = t1163*t1172;
+ t1159 = RATIONAL(4.0,9.0);
+ t1242 = t1122+t1159*x;
+ t1169 = t1170*t1170;
+ t1241 = t1156*t1170+t1155*t1169;
+ t1130 = RATIONAL(-1.0,288.0);
+ t1133 = RATIONAL(1.0,288.0);
+ t1240 = t1130*t1169+t1133*t1170;
+ t1143 = RATIONAL(-1.0,12.0);
+ t1142 = RATIONAL(-1.0,48.0);
+ t1196 = t1142*t1169;
+ t1239 = t1196+t1143;
+ t1171 = y*y;
+ t1173 = y*t1171;
+ t1213 = y*t1169;
+ t1144 = RATIONAL(1.0,18.0);
+ t1146 = RATIONAL(1.0,36.0);
+ t1105 = t1146*t1170;
+ t1238 = t1173*t1105+t1144*t1213;
+ t1200 = t1173*t1155;
+ t1145 = RATIONAL(1.0,72.0);
+ t1117 = t1145*t1169;
+ t1237 = t1170*t1200+y*t1117;
+ t1166 = RATIONAL(-1.0,6.0);
+ t1205 = t1166*t1171;
+ t1134 = RATIONAL(1.0,24.0);
+ t1148 = RATIONAL(-1.0,72.0);
+ t1201 = t1169*t1171;
+ t1236 = t1148*t1169+t1134*t1201+(t1144+t1205)*t1170;
+ t1215 = x*t1171;
+ t1135 = RATIONAL(1.0,48.0);
+ t1198 = t1171*t1172;
+ t1235 = t1142*t1215+t1135*t1198+t1244;
+ t1112 = t1145*t1172;
+ t1217 = x*t1173;
+ t1234 = t1173*t1112+t1148*t1217+t1244*y;
+ t1214 = y*RATIONAL(8.0,9.0);
+ t1212 = y*RATIONAL(-2.0,9.0);
+ t1233 = t1170*t1214+t1169*t1212+t1241*t1173;
+ t1161 = RATIONAL(1.0,3.0);
+ t1232 = t1143*t1198+t1161*t1215+t1245;
+ t1231 = x*t1214+t1172*t1212+t1245*t1173;
+ t1207 = t1173*t1170;
+ t1199 = t1173*t1169;
+ t1230 = t1131*t1199+t1132*t1207+t1240*y;
+ t1195 = t1135*t1169;
+ t1210 = t1170*t1171;
+ t1149 = RATIONAL(-1.0,36.0);
+ t1106 = t1149*t1169;
+ t1229 = t1142*t1210+t1171*t1195+t1105+t1106;
+ t1203 = t1135*t1171;
+ t1228 = t1171*t1196+t1146*t1169+(t1203+t1149)*t1170;
+ t1218 = x*y;
+ t1147 = RATIONAL(-1.0,18.0);
+ t1111 = t1147*x;
+ t1104 = t1144*t1172;
+ t1227 = y*t1122+t1155*t1218+(t1111+t1104)*t1173;
+ t1226 = RATIONAL(2.0,9.0)*y*t1172+RATIONAL(-8.0,9.0)*t1218+t1242*t1173;
+ t1141 = RATIONAL(1.0,12.0);
+ t1160 = RATIONAL(-1.0,3.0);
+ t1225 = t1160*t1215+t1141*t1198+t1242;
+ t1114 = t1147*t1172;
+ t1108 = t1144*x;
+ t1224 = y*t1120+t1163*t1218+(t1108+t1114)*t1173;
+ t1204 = t1173*t1172;
+ t1223 = t1149*t1204+x*t1200+(t1112+t1111)*y;
+ t1118 = t1147*t1170;
+ t1139 = RATIONAL(-1.0,24.0);
+ t1164 = RATIONAL(1.0,6.0);
+ t1208 = t1164*t1171;
+ t1222 = t1118+t1139*t1201+t1170*t1208+t1117;
+ t1221 = t1142*t1198+x*t1203+t1243;
+ t1107 = t1148*t1172;
+ t1220 = t1146*t1204+t1163*t1217+(t1108+t1107)*y;
+ t1219 = t1173*t1107+t1145*t1217+t1243*y;
+ t1216 = RATIONAL(-5.0,3.0)*y;
+ t1168 = RATIONAL(-2.0,3.0);
+ t1128 = t1168*t1173;
+ t1197 = t1166*t1169;
+ t1211 = t1170*t1216+t1173*t1197;
+ t1119 = t1164*t1173;
+ t1193 = RATIONAL(5.0,48.0)*t1170;
+ t1209 = t1134*t1199+y*t1193;
+ t1206 = y*t1118+t1173*t1106;
+ t1165 = RATIONAL(2.0,3.0);
+ t1202 = t1165*t1173;
+ t1194 = RATIONAL(5.0,6.0)*t1170;
+ t1192 = t1111+t1139*t1198+t1206;
+ t1140 = RATIONAL(-5.0,48.0);
+ t1191 = t1139*t1207+t1140*t1213+t1209;
+ t1190 = t1134*t1198+t1108+t1206;
+ t1150 = RATIONAL(5.0,12.0);
+ t1189 = t1170*t1202+t1150*t1213+t1211;
+ t1152 = RATIONAL(-5.0,24.0);
+ t1188 = t1141*t1173+t1152*y;
+ t1187 = RATIONAL(5.0,24.0)*y+t1143*t1173;
+ t1186 = t1128+t1173*t1194+RATIONAL(4.0,3.0)*y+t1161*t1213+t1211;
+ t1185 = t1119+t1152*t1207+t1209+t1239*y;
+ t1136 = RATIONAL(1.0,96.0);
+ t1137 = RATIONAL(-1.0,96.0);
+ t1184 = t1137*t1201+t1130*t1170+t1133*t1169+t1136*t1210+t1230;
+ t1183 = t1134*t1215+t1104+t1192+t1238;
+ t1182 = t1136*t1201+t1137*t1210+t1230+t1240;
+ t1181 = t1139*t1215+t1114+t1190+t1238;
+ t1180 = t1112+x*t1208+t1192+t1237;
+ t1179 = t1159*t1170+t1163*t1169+t1160*t1210+t1141*t1201+t1233;
+ t1178 = x*t1205+t1107+t1190+t1237;
+ t1177 = t1161*t1210+t1143*t1201+t1233+t1241;
+ t1154 = RATIONAL(1.0,4.0);
+ t1153 = RATIONAL(-5.0,8.0);
+ coeffs_dy->coeff_m2_m2 = t1184+t1219+t1235;
+ coeffs_dy->coeff_m1_m2 = t1180+t1220+t1236;
+ coeffs_dy->coeff_0_m2 = t1140*t1170+t1195+t1141+(RATIONAL(-1.0,16.0)*
+t1169+RATIONAL(5.0,16.0)*t1170+RATIONAL(-1.0,4.0))*t1171+t1185;
+ coeffs_dy->coeff_p1_m2 = t1178+t1223+t1236;
+ coeffs_dy->coeff_p2_m2 = t1184+t1221+t1234;
+ coeffs_dy->coeff_m2_m1 = t1183+t1227+t1229;
+ coeffs_dy->coeff_m1_m1 = t1177+t1225+t1226;
+ coeffs_dy->coeff_0_m1 = t1194+t1197+t1168+(RATIONAL(1.0,2.0)+t1153*t1170+
+RATIONAL(1.0,8.0)*t1169)*t1171+t1186;
+ coeffs_dy->coeff_p1_m1 = t1177+t1231+t1232;
+ coeffs_dy->coeff_p2_m1 = t1181+t1224+t1229;
+ coeffs_dy->coeff_m2_0 = t1187*t1172+t1188*x+t1191;
+ coeffs_dy->coeff_m1_0 = (RATIONAL(-5.0,12.0)*y+t1119)*t1172+(t1128+
+RATIONAL(5.0,3.0)*y)*x+t1189;
+ coeffs_dy->coeff_0_0 = (RATIONAL(-5.0,4.0)*t1170+t1154*t1169+1.0)*t1173+(
+RATIONAL(-5.0,2.0)+RATIONAL(25.0,8.0)*t1170+t1153*t1169)*y;
+ coeffs_dy->coeff_p1_0 = (t1166*t1173+t1150*y)*t1172+(t1202+t1216)*x+t1189
+;
+ coeffs_dy->coeff_p2_0 = t1188*t1172+t1187*x+t1191;
+ coeffs_dy->coeff_m2_p1 = t1181+t1227+t1228;
+ coeffs_dy->coeff_m1_p1 = t1179+t1226+t1232;
+ coeffs_dy->coeff_0_p1 = t1164*t1169+t1165+RATIONAL(-5.0,6.0)*t1170+(
+RATIONAL(5.0,8.0)*t1170+RATIONAL(-1.0,8.0)*t1169+RATIONAL(-1.0,2.0))*t1171+
+t1186;
+ coeffs_dy->coeff_p1_p1 = t1179+t1225+t1231;
+ coeffs_dy->coeff_p2_p1 = t1183+t1224+t1228;
+ coeffs_dy->coeff_m2_p2 = t1182+t1219+t1221;
+ coeffs_dy->coeff_m1_p2 = t1178+t1220+t1222;
+ coeffs_dy->coeff_0_p2 = t1193+(RATIONAL(1.0,16.0)*t1169+t1154+RATIONAL(
+-5.0,16.0)*t1170)*t1171+t1185+t1239;
+ coeffs_dy->coeff_p1_p2 = t1180+t1222+t1223;
+ coeffs_dy->coeff_p2_p2 = t1182+t1234+t1235;
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-20 00:46:36 +0000