← Back to branch summary

~pmdj/ubuntu/trusty/qemu/2.9+applesmc+fadtv3

« back to all changes in this revision

Viewing changes to pixman/demos/quad2quad.c

  • Committer: Phil Dennis-Jordan
  • Date: 2017-07-21 08:03:43 UTC
  • mfrom: (1.1.1)
  • Revision ID: phil@philjordan.eu-20170721080343-2yr2vdj7713czahv
New upstream release 2.9.0.

Show diffs side-by-side

added added

removed removed

Lines of Context:
pixman/demos/quad2quad.c
 
1
#include <math.h>
 
2
#include <stdio.h>
 
3
#include <pixman.h>
 
4
 
 
5
/* This code is basically the output of Maxima translated into C.
 
6
 *
 
7
 * See http://maxima.sourceforge.net/
 
8
 */
 
9
static void
 
10
quad_to_quad (double x0, double y0,
 
11
              double x1, double y1,
 
12
              double x2, double y2,
 
13
              double x3, double y3,
 
14
 
 
15
              double px0, double py0,
 
16
              double px1, double py1,
 
17
              double px2, double py2,
 
18
              double px3, double py3,
 
19
 
 
20
              struct pixman_f_transform *trans)
 
21
{
 
22
    double
 
23
        t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15, t16, t17, t18,
 
24
        t19, t20, t21, t22, t23, t24, t25, t26, t27, t28, t29, t30, t31, t32, t33, t34,
 
25
        t35, t36, t37, t38, t39, t40, t41, t42, t43, t44, t45, t46, t47, t48, t49, t50,
 
26
        t51, t52, t53, t54, t55, t56, t57, t58, t59, t60, t61, t62, t63, t64, t65, t66,
 
27
        t67, t68, t69, t70, t71, t72, t73, t74, t75, t76, t77, t78, t79, t80, t81, t82,
 
28
        t83, t84, t85, t86, t87, t88, t89, t90, t91, t92, t93, t94, t95, t96, t97, t98,
 
29
        t99, t100, t101, t102, t103, t104, t105, t106, t107, t108, t109, t110, t111,
 
30
        t112, t113, t114, t115, t116, t117, t118, t119, t120, t121, t122, t123,
 
31
        t124, t125, t126, t127, t128, t129, t130, t131, t132, t133, t134, t135,
 
32
        t136, t137, t138, t139, t140, t141, t142, t143, t144, t145, t146, t147,
 
33
        t148, t149, t150, t151, t152, t153, t154, t155, t156, t157, t158, t159,
 
34
        t160, t161, t162, t163, t164, t165, t166, t167, t168, t169, t170, t171,
 
35
        t172, t173, t174, t175, t176, t177, t178, t179, t180, t181, t182, t183,
 
36
        t184, t185, t186, t187, t188, t189, t190, t191, t192, t193, t194, t195,
 
37
        t196, t197, t198, t199, t200, t201, t202, t203, t204, t205, t206, t207,
 
38
        t208, t209, t210, t211, t212, t213, t214, t215, t216, t217, t218, t219,
 
39
        t220, t221, t222, t223, t224, t225, t226, t227, t228, t229, t230, t231,
 
40
        t232, t233, t234, t235, t236, t237, t238, t239, t240, t241, t242, t243,
 
41
        t244, t245, t246, t247, t248, t249, t250, t251, t252, t253, t254, t255,
 
42
        t256, t257, t258, t259, t260, t261, t262, t263, t264, t265, t266, t267,
 
43
        t268, t269, t270, t271, t272, t273, t274, t275, t276, t277, t278, t279,
 
44
        t280, t281, t282, t283, t284, t285, t286, t287, t288, t289, t290, t291,
 
45
        t292, t293, t294, t295, t296, t297, t298, t299, t300, t301, t302, t303,
 
46
        t304, t305, t306, t307, t308, t309, t310, t311, t312, t313, t314, t315,
 
47
        t316, t317, t318, t319, t320, t321, t322, t323, t324, t325, t326, t327,
 
48
        t328, t329, t330, t331, t332, t333, t334, t335, t336, t337, t338, t339,
 
49
        t340, t341, t342, t343, t344, t345, t346, t347, t348, t349, t350, t351,
 
50
        t352, t353, t354, t355, t356, t357, t358, t359, t360, t361, t362, t363,
 
51
        t364, t365, t366, t367, t368, t369, t370, t371, t372, t373, t374, t375,
 
52
        t376, t377, t378, t379, t380, t381, t382, t383, t384, t385, t386, t387,
 
53
        t388, t389, t390, t391, t392, t393, t394, t395, t396, t397, t398, t399,
 
54
        t400, t401, t402, t403, t404, t405, t406, t407, t408, t409, t410, t411,
 
55
        t412, t413, t414, t415, t416, t417, t418, t419, t420, t421, t422, t423,
 
56
        t424, t425, t426, t427, t428, t429, t430, t431, t432, t433, t434, t435,
 
57
        t436, t437, t438, t439, t440, t441, t442, t443, t444, t445, t446, t447,
 
58
        t448, t449, t450, t451, t452, t453, t454, t455, t456, t457, t458, t459,
 
59
        t460, t461, t462, t463, t464, t465, t466, t467, t468, t469, t470, t471,
 
60
        t472, t473, t474, t475, t476, t477, t478, t479, t480, t481, t482, t483,
 
61
        t484, t485, t486, t487, t488, t489, t490, t491, t492, t493, t494, t495,
 
62
        t496, t497, t498, t499, t500, t501, t502, t503, t504, t505, t506, t507,
 
63
        t508, t509, t510, t511, t512, t513, t514, t515, t516, t517, t518, t519,
 
64
        t520, t521, t522, t523, t524, t525, t526, t527, t528, t529, t530, t531,
 
65
        t532, t533, t534, t535, t536, t537, t538, t539, t540, t541, t542, t543,
 
66
        t544, t545, t546, t547, t548, t549, t550, t551, t552, t553, t554, t555,
 
67
        t556, t557, t558, t559, t560, t561, t562, t563, t564, t565, t566, t567,
 
68
        t568, t569, t570, t571, t572, t573, t574, t575, t576, t577, t578, t579,
 
69
        t580, t581, t582, t583, t584, t585, t586, t587, t588, t589, t590, t591,
 
70
        t592, t593, t594, t595, t596, t597, t598, t599, t600, t601, t602, t603,
 
71
        t604, t605, t606, t607, t608, t609, t610, t611, t612, t613, t614, t615,
 
72
        t616, t617, t618, t619, t620, t621, t622, t623, t624, t625, t626, t627,
 
73
        t628, t629, t630, t631, t632, t633, t634, t635, t636, t637, t638, t639,
 
74
        t640, t641, t642, t643, t644, t645, t646, t647, t648, t649, t650, t651,
 
75
        t652, t653, t654, t655, t656, t657, t658, t659, t660, t661, t662, t663,
 
76
        t664, t665, t666, t667, t668, t669, t670, t671, t672, t673, t674, t675,
 
77
        t676, t677, t678, t679, t680, t681, t682, t683, t684, t685, t686, t687,
 
78
        t688, t689, t690, t691, t692, t693, t694, t695, t696, t697, t698, t699,
 
79
        t700, t701, t702, t703, t704, t705, t706, t707, t708, t709, t710, t711,
 
80
        t712, t713, t714, t715, t716, t717, t718, t719, t720, t721, t722, t723,
 
81
        t724, t725, t726, t727, t728, t729, t730, t731, t732, t733, t734, t735,
 
82
        t736, t737, t738, t739, t740, t741, t742, t743, t744, t745, t746, t747,
 
83
        t748, t749, t750, t751, t752, t753, t754, t755, t756, t757, t758, t759,
 
84
        t760, t761, t762, t763, t764, t765, t766, t767, t768, t769, t770, t771,
 
85
        t772, t773, t774, t775, t776, t777, t778, t779, t780, t781, t782, t783,
 
86
        t784, t785, t786, t787, t788, t789, t790, t791, t792, t793, t794, t795,
 
87
        t796, t797, t798, t799, t800, t801, t802, t803, t804, t805, t806, t807,
 
88
        t808, t809, t810, t811, t812, t813, t814, t815, t816, t817, t818, t819,
 
89
        t820, t821, t822, t823, t824, t825, t826, t827, t828, t829, t830, t831,
 
90
        t832, t833, t834, t835, t836, t837, t838, t839, t840, t841, t842, t843,
 
91
        t844, t845, t846, t847, t848, t849, t850, t851, t852, t853, t854, t855,
 
92
        t856, t857, t858, t859, t860, t861, t862, t863, t864, t865, t866, t867,
 
93
        t868, t869, t870, t871, t872, t873, t874, t875, t876, t877, t878, t879,
 
94
        t880, t881, t882, t883, t884, t885, t886, t887, t888, t889, t890, t891,
 
95
        t892, t893, t894, t895, t896, t897, t898, t899, t900, t901, t902, t903,
 
96
        t904, t905, t906, t907, t908, t909, t910, t911, t912, t913, t914, t915,
 
97
        t916, t917, t918, t919, t920, t921, t922, t923, t924, t925, t926, t927,
 
98
        t928, t929, t930, t931, t932, t933, t934, t935, t936, t937, t938, t939,
 
99
        t940, t941, t942, t943, t944, t945, t946, t947, t948, t949, t950, t951,
 
100
        t952, t953, t954, t955, t956, t957, t958, t959, t960, t961, t962, t963,
 
101
        t964, t965, t966, t967, t968, t969, t970, t971, t972, t973, t974, t975,
 
102
        t976, t977, t978, t979, t980, t981, t982, t983, t984, t985, t986, t987,
 
103
        t988, t989, t990, t991, t992, t993, t994, t995, t996, t997, t998, t999,
 
104
        t1000, t1001, t1002, t1003, t1004, t1005, t1006, t1007, t1008, t1009,
 
105
        t1010, t1011, t1012, t1013, t1014, t1015, t1016, t1017, t1018, t1019,
 
106
        t1020, t1021, t1022, t1023, t1024, t1025, t1026, t1027, t1028, t1029,
 
107
        t1030, t1031, t1032, t1033, t1034, t1035, t1036, t1037, t1038, t1039,
 
108
        t1040, t1041, t1042, t1043, t1044, t1045, t1046, t1047, t1048, t1049,
 
109
        t1050, t1051, t1052, t1053, t1054, t1055, t1056, t1057, t1058, t1059,
 
110
        t1060, t1061, t1062, t1063, t1064, t1065, t1066, t1067, t1068, t1069,
 
111
        t1070, t1071, t1072, t1073;
 
112
 
 
113
    t1 = y1 * y1;
 
114
    t2 = x3 * x3;
 
115
    t3 = px2 * px3 * t2;
 
116
    t4 = (t3 - px2 * px3 * x2 * x3) * y2;
 
117
    t5 = x2 * x2;
 
118
    t6 = px2 * px3 * t5 * y3;
 
119
 
 
120
    t7 = - px2 * px3 * x2 * x3 * y3;
 
121
    t8 = py1 * (t7 + t6 + t4);
 
122
    t9 = px3 * py2 * x2 * x3;
 
123
 
 
124
    t10 = - px3 * py2 * t2;
 
125
    t11 = (t10 + t9) * y2;
 
126
    t12 = - px2 * py3 * t5 * y3;
 
127
 
 
128
    t13 = px2 * py3 * x2 * x3 * y3;
 
129
    t14 = y0 * y0;
 
130
    t15 = - px3 * py2;
 
131
    t16 = px2 * py3;
 
132
 
 
133
    t17 = t16 + t15;
 
134
    t18 = t17 * x2;
 
135
    t19 = px3 * py2 * x3;
 
136
    t20 = - px2 * py3 * x3;
 
137
 
 
138
    t21 = t20 + t19 + t18;
 
139
    t22 = px2 * px3 * t5;
 
140
    t23 = - 2 * px2 * px3 * x2 * x3;
 
141
 
 
142
    t24 = py1 * (t3 + t23 + t22);
 
143
    t25 = - px2 * py3 * t5;
 
144
    t26 = px2 * py3 * x3;
 
145
 
 
146
    t27 = x2 * (t26 + t19);
 
147
    t28 = t10 + t27 + t25;
 
148
    t29 = x1 * x1;
 
149
    t30 = px3 * py2;
 
150
 
 
151
    t31 = - px2 * py3;
 
152
    t32 = t31 + t30;
 
153
    t33 = t32 * y2;
 
154
    t34 = - px3 * py2 * y3;
 
155
 
 
156
    t35 = px2 * py3 * y3;
 
157
    t36 = t35 + t34 + t33;
 
158
    t37 = - px2 * px3 * t2;
 
159
 
 
160
    t38 = (t37 + px2 * px3 * x2 * x3) * y2;
 
161
    t39 = - px2 * px3 * t5 * y3;
 
162
 
 
163
    t40 = px2 * px3 * x2 * x3 * y3;
 
164
    t41 = py1 * (t40 + t39 + t38);
 
165
    t42 = - px2 * py3 * x2 * x3;
 
166
 
 
167
    t43 = px3 * py2 * t2;
 
168
    t44 = (t43 + t42) * y2;
 
169
    t45 = px2 * py3 * t5 * y3;
 
170
 
 
171
    t46 = - px3 * py2 * x2 * x3 * y3;
 
172
    t47 = (px2 * px3 * x3 - px2 * px3 * x2) * y2;
 
173
 
 
174
    t48 = px2 * px3 * x2 * y3;
 
175
    t49 = - px2 * px3 * x3 * y3;
 
176
    t50 = py1 * (t49 + t48 + t47);
 
177
 
 
178
    t51 = px2 * py3 * x2;
 
179
    t52 = - 2 * px3 * py2 * x3;
 
180
    t53 = (t26 + t52 + t51) * y2;
 
181
 
 
182
    t54 = px3 * py2 * x3 * y3;
 
183
    t55 = px3 * py2 * y3;
 
184
    t56 = - 2 * px2 * py3 * y3;
 
185
    t57 = t56 + t55;
 
186
 
 
187
    t58 = x2 * t57;
 
188
    t59 = - px2 * px3 * t5;
 
189
    t60 = 2 * px2 * px3 * x2 * x3;
 
190
    t61 = - px2;
 
191
 
 
192
    t62 = px3 + t61;
 
193
    t63 = t62 * x2;
 
194
    t64 = px2 * x3;
 
195
    t65 = - px3 * x3;
 
196
    t66 = t65 + t64 + t63;
 
197
 
 
198
    t67 = px2 * t5;
 
199
    t68 = - px2 * x3;
 
200
    t69 = x2 * (t65 + t68);
 
201
    t70 = px3 * t2;
 
202
 
 
203
    t71 = t70 + t69 + t67;
 
204
    t72 = - px3;
 
205
    t73 = t72 + px2;
 
206
    t74 = - px2 * y3;
 
207
    t75 = px3 * y3;
 
208
 
 
209
    t76 = t75 + t74 + t73 * y2;
 
210
    t77 = px2 * x2 * x3;
 
211
    t78 = - px3 * t2;
 
212
    t79 = - px2 * t5 * y3;
 
213
 
 
214
    t80 = px3 * x2 * x3 * y3;
 
215
    t81 = t80 + t79 + (t78 + t77) * y2;
 
216
 
 
217
    t82 = (px2 * px3 * x2 - px2 * px3 * x3) * y2;
 
218
    t83 = - px2 * px3 * x2 * y3;
 
219
 
 
220
    t84 = px2 * px3 * x3 * y3;
 
221
    t85 = - px2 * x2;
 
222
    t86 = 2 * px3 * x3;
 
223
    t87 = - px3 * x3 * y3;
 
224
 
 
225
    t88 = 2 * px2 * y3;
 
226
    t89 = - px3 * y3;
 
227
    t90 = t89 + t88;
 
228
    t91 = x2 * t90;
 
229
 
 
230
    t92 = t91 + t87 + (t86 + t68 + t85) * y2;
 
231
    t93 = px2 * py3 * t5;
 
232
    t94 = - px3 * py2 * x3;
 
233
 
 
234
    t95 = x2 * (t20 + t94);
 
235
    t96 = t32 * x2;
 
236
    t97 = t73 * x2;
 
237
    t98 = px3 * x3;
 
238
 
 
239
    t99 = t98 + t68 + t97;
 
240
    t100 = py1 * t99;
 
241
    t101 = - px2 * t5;
 
242
    t102 = x2 * (t98 + t64);
 
243
 
 
244
    t103 = t78 + t102 + t101;
 
245
    t104 = py1 * t103;
 
246
    t105 = - py2;
 
247
    t106 = py3 + t105;
 
248
 
 
249
    t107 = py2 * y3;
 
250
    t108 = - py3 * y3;
 
251
    t109 = t108 + t107 + t106 * y2;
 
252
    t110 = - px3 * x2 * x3;
 
253
 
 
254
    t111 = px2 * t5 * y3;
 
255
    t112 = - px2 * x2 * x3 * y3;
 
256
    t113 = t112 + t111 + (t70 + t110) * y2;
 
257
 
 
258
    t114 = - py2 * x3;
 
259
    t115 = py3 * x3;
 
260
    t116 = t115 + t114;
 
261
    t117 = py2 * x3 * y3;
 
262
 
 
263
    t118 = - py3 * x3 * y3;
 
264
    t119 = t118 + t117;
 
265
    t120 = x2 * t119;
 
266
 
 
267
    t121 = px1 * (t120 + x2 * t116 * y2);
 
268
    t122 = - px3 * py2 * x2;
 
269
    t123 = (t19 + t122) * y2;
 
270
 
 
271
    t124 = px2 * py3 * x2 * y3;
 
272
    t125 = - px2 * py3 * x3 * y3;
 
273
    t126 = px3 * x2;
 
274
 
 
275
    t127 = - px2 * x2 * y3;
 
276
    t128 = px2 * x3 * y3;
 
277
    t129 = t128 + t127 + (t65 + t126) * y2;
 
278
 
 
279
    t130 = - py3;
 
280
    t131 = t130 + py2;
 
281
    t132 = t131 * x2;
 
282
    t133 = py2 * x3;
 
283
    t134 = - py3 * x3;
 
284
 
 
285
    t135 = - py2 * x3 * y3;
 
286
    t136 = py3 * x3 * y3;
 
287
    t137 = - py2 * y3;
 
288
    t138 = py3 * y3;
 
289
 
 
290
    t139 = t138 + t137;
 
291
    t140 = x2 * t139;
 
292
 
 
293
    t141 = px1 * (t140 + t136 + t135 + (t134 + t133 + t132) * y2);
 
294
    t142 = y2 * y2;
 
295
 
 
296
    t143 = - px3 * py2 * x3 * y3;
 
297
    t144 = px2 * py3 * x3 * y3;
 
298
    t145 = t144 + t143;
 
299
 
 
300
    t146 = t142 * t145;
 
301
    t147 = y3 * y3;
 
302
    t148 = px3 * py2 * t147;
 
303
    t149 = - px2 * py3 * t147;
 
304
 
 
305
    t150 = t149 + t148;
 
306
    t151 = x2 * y2 * t150;
 
307
    t152 = t151 + t146;
 
308
    t153 = - px2 * py3 * y3;
 
309
 
 
310
    t154 = t153 + t55;
 
311
    t155 = t142 * t154;
 
312
    t156 = - px3 * py2 * t147;
 
313
 
 
314
    t157 = px2 * py3 * t147;
 
315
    t158 = t157 + t156;
 
316
    t159 = y2 * t158;
 
317
    t160 = t159 + t155;
 
318
 
 
319
    t161 = x0 * x0;
 
320
    t162 = py1 * t76;
 
321
    t163 = px1 * t109;
 
322
    t164 = px2 * y3;
 
323
    t165 = t89 + t164;
 
324
 
 
325
    t166 = - px2 * t147;
 
326
    t167 = px3 * t147;
 
327
    t168 = t167 + t166;
 
328
 
 
329
    t169 = y2 * t168 + t142 * t165;
 
330
    t170 = py1 * t169;
 
331
    t171 = py2 * t147;
 
332
 
 
333
    t172 = - py3 * t147;
 
334
    t173 = t172 + t171;
 
335
    t174 = y2 * t173 + t142 * t139;
 
336
 
 
337
    t175 = px1 * t174;
 
338
    t176 = t17 * t142;
 
339
    t177 = px2 * t147;
 
340
    t178 = - px3 * t147;
 
341
 
 
342
    t179 = t178 + t177 + t62 * t142;
 
343
    t180 = - py2 * t147;
 
344
    t181 = py3 * t147;
 
345
 
 
346
    t182 = t181 + t180 + t131 * t142;
 
347
 
 
348
    t183 = y1 * (px1 * t182 + py1 * t179 + t149 + t148 + t176)
 
349
        + t175 + t170 + t159 + t1 * (t163 + t162 + t35 + t34 + t33) + t155;
 
350
 
 
351
    t184 = - px2 * px3 * t2 * t142;
 
352
    t185 = 2 * px2 * px3 * x2 * x3 * y2 * y3;
 
353
 
 
354
    t186 = - px2 * px3 * t5 * t147;
 
355
    t187 = py1 * (t186 + t185 + t184);
 
356
 
 
357
    t188 = px3 * py2 * t2 * t142;
 
358
    t189 = x2 * y2 * (t125 + t143);
 
359
    t190 = px2 * py3 * t5 * t147;
 
360
 
 
361
    t191 = t190 + t189 + t188;
 
362
    t192 = px2 * px3 * x3 * t142;
 
363
    t193 = y2 * (t49 + t83);
 
364
 
 
365
    t194 = px2 * px3 * x2 * t147;
 
366
    t195 = py1 * (t194 + t193 + t192);
 
367
 
 
368
    t196 = - px3 * py2 * x3 * t142;
 
369
    t197 = 2 * px3 * py2 * x3 * y3;
 
370
    t198 = 2 * px2 * py3 * y3;
 
371
 
 
372
    t199 = t198 + t34;
 
373
    t200 = x2 * t199;
 
374
    t201 = y2 * (t200 + t125 + t197);
 
375
 
 
376
    t202 = - px2 * py3 * x2 * t147;
 
377
    t203 = - px2 * x3 * y3;
 
378
    t204 = px3 * x3 * y3;
 
379
 
 
380
    t205 = t204 + t203;
 
381
    t206 = t142 * t205;
 
382
    t207 = t178 + t177;
 
383
    t208 = x2 * y2 * t207;
 
384
 
 
385
    t209 = t208 + t206;
 
386
    t210 = px2 * px3 * t2 * t142;
 
387
    t211 = - 2 * px2 * px3 * x2 * x3 * y2 * y3;
 
388
 
 
389
    t212 = px2 * px3 * t5 * t147;
 
390
    t213 = - px3 * t2 * t142;
 
391
    t214 = x2 * y2 * (t204 + t128);
 
392
 
 
393
    t215 = - px2 * t5 * t147;
 
394
    t216 = t215 + t214 + t213;
 
395
    t217 = - px2 * px3 * x3 * t142;
 
396
 
 
397
    t218 = y2 * (t84 + t48);
 
398
    t219 = - px2 * px3 * x2 * t147;
 
399
    t220 = px3 * x3 * t142;
 
400
 
 
401
    t221 = - 2 * px3 * x3 * y3;
 
402
    t222 = - 2 * px2 * y3;
 
403
    t223 = t75 + t222;
 
404
    t224 = x2 * t223;
 
405
 
 
406
    t225 = y2 * (t224 + t221 + t128);
 
407
    t226 = px2 * x2 * t147;
 
408
    t227 = t226 + t225 + t220;
 
409
 
 
410
    t228 = t125 + t54;
 
411
    t229 = t142 * t228;
 
412
    t230 = x2 * y2 * t158;
 
413
    t231 = t87 + t128;
 
414
 
 
415
    t232 = t142 * t231;
 
416
    t233 = x2 * y2 * t168;
 
417
    t234 = t233 + t232;
 
418
    t235 = py1 * t234;
 
419
 
 
420
    t236 = - px3 * py2 * t2 * t142;
 
421
    t237 = x2 * y2 * (t144 + t54);
 
422
 
 
423
    t238 = - px2 * py3 * t5 * t147;
 
424
    t239 = px3 * t2 * t142;
 
425
    t240 = x2 * y2 * (t87 + t203);
 
426
 
 
427
    t241 = px2 * t5 * t147;
 
428
    t242 = t241 + t240 + t239;
 
429
    t243 = py1 * t242;
 
430
 
 
431
    t244 = px2 * py3 * x3 * t142;
 
432
    t245 = - px2 * py3 * x2 * y3;
 
433
    t246 = y2 * (t143 + t245);
 
434
 
 
435
    t247 = px3 * py2 * x2 * t147;
 
436
    t248 = - px2 * x3 * t142;
 
437
    t249 = px2 * x2 * y3;
 
438
 
 
439
    t250 = y2 * (t204 + t249);
 
440
    t251 = - px3 * x2 * t147;
 
441
    t252 = t251 + t250 + t248;
 
442
 
 
443
    t253 = t134 + t133;
 
444
    t254 = t253 * t142;
 
445
    t255 = t108 + t107;
 
446
    t256 = x2 * t255;
 
447
 
 
448
    t257 = t256 + t136 + t135;
 
449
    t258 = y2 * t257;
 
450
    t259 = t181 + t180;
 
451
    t260 = x2 * t259;
 
452
 
 
453
    t261 = px1 * (t260 + t258 + t254);
 
454
    t262 = py1 * (t37 + t60 + t59);
 
455
 
 
456
    t263 = t43 + t95 + t93;
 
457
    t264 = px1 * t263;
 
458
    t265 = t26 + t94;
 
459
    t266 = x2 * t265 * y2;
 
460
 
 
461
    t267 = x2 * t228;
 
462
    t268 = t267 + t266;
 
463
    t269 = py1 * (t84 + t83 + t82);
 
464
 
 
465
    t270 = - 2 * px2 * py3;
 
466
    t271 = (t26 + (t270 + t30) * x2) * y2;
 
467
    t272 = px3 * py2 * x2 * y3;
 
468
 
 
469
    t273 = - 2 * px3 * py2 * x3 * y3;
 
470
    t274 = t149 + t148 + t176;
 
471
 
 
472
    t275 = py1 * (t212 + t211 + t210);
 
473
    t276 = t238 + t237 + t236;
 
474
    t277 = px1 * t276;
 
475
 
 
476
    t278 = py1 * (t219 + t218 + t217);
 
477
    t279 = 2 * px3 * py2 * x3;
 
478
    t280 = t20 + t279;
 
479
 
 
480
    t281 = t280 * t142;
 
481
    t282 = - px3 * py2 * x2 * y3;
 
482
    t283 = y2 * (t125 + t282);
 
483
 
 
484
    t284 = 2 * px2 * py3 * t147;
 
485
    t285 = x2 * (t284 + t156);
 
486
    t286 = px1 * t103;
 
487
 
 
488
    t287 = t98 + t68;
 
489
    t288 = x2 * t287 * y2;
 
490
    t289 = x2 * t231;
 
491
    t290 = t289 + t288;
 
492
 
 
493
    t291 = 2 * px2;
 
494
    t292 = - px3 * x2 * y3;
 
495
    t293 = 2 * px3 * x3 * y3;
 
496
 
 
497
    t294 = t293 + t203 + t292 + (t68 + (t72 + t291) * x2) * y2;
 
498
    t295 = px1 * t242;
 
499
 
 
500
    t296 = - 2 * px3 * x3;
 
501
    t297 = t296 + t64;
 
502
    t298 = px3 * x2 * y3;
 
503
    t299 = y2 * (t128 + t298);
 
504
 
 
505
    t300 = - 2 * px2 * t147;
 
506
    t301 = x2 * (t167 + t300) + t299 + t297 * t142;
 
507
    t302 = py1 * t71;
 
508
 
 
509
    t303 = py1 * t290;
 
510
    t304 = 2 * py2 * x3;
 
511
    t305 = - 2 * py3 * x3;
 
512
    t306 = - 2 * py2 * x3 * y3;
 
513
 
 
514
    t307 = 2 * py3 * x3 * y3;
 
515
    t308 = t307 + t306;
 
516
    t309 = - 2 * px2 * py3 * x3;
 
517
 
 
518
    t310 = (t309 + t19 + t51) * y2;
 
519
    t311 = - 2 * px3 * py2 * y3;
 
520
    t312 = t35 + t311;
 
521
 
 
522
    t313 = x2 * t312;
 
523
    t314 = 2 * px2 * x3;
 
524
    t315 = 2 * px3 * y3;
 
525
    t316 = t315 + t74;
 
526
 
 
527
    t317 = x2 * t316;
 
528
    t318 = t317 + t87 + (t65 + t314 + t85) * y2;
 
529
    t319 = t106 * x2;
 
530
 
 
531
    t320 = px1 * (t256 + t118 + t117 + (t115 + t114 + t319) * y2);
 
532
    t321 = py1 * t216;
 
533
 
 
534
    t322 = 2 * px2 * py3 * x3 * y3;
 
535
    t323 = 2 * px3 * py2 * y3;
 
536
    t324 = t153 + t323;
 
537
 
 
538
    t325 = x2 * t324;
 
539
    t326 = y2 * (t325 + t322 + t143);
 
540
    t327 = - 2 * px2 * x3 * y3;
 
541
 
 
542
    t328 = - 2 * px3 * y3;
 
543
    t329 = t328 + t164;
 
544
    t330 = x2 * t329;
 
545
 
 
546
    t331 = y2 * (t330 + t204 + t327);
 
547
    t332 = t226 + t331 + t220;
 
548
    t333 = t116 * t142;
 
549
 
 
550
    t334 = t140 + t118 + t117;
 
551
    t335 = y2 * t334;
 
552
    t336 = x2 * t173;
 
553
 
 
554
    t337 = px1 * (t336 + t335 + t333);
 
555
    t338 = t26 + t94 + t96;
 
556
    t339 = t17 * y2;
 
557
 
 
558
    t340 = t153 + t55 + t339;
 
559
    t341 = px2 * px3 * t142;
 
560
    t342 = - 2 * px2 * px3 * y2 * y3;
 
561
 
 
562
    t343 = px2 * px3 * t147;
 
563
    t344 = py1 * (t343 + t342 + t341);
 
564
    t345 = - px2 * py3 * t142;
 
565
 
 
566
    t346 = y2 * (t35 + t55);
 
567
    t347 = t156 + t346 + t345;
 
568
    t348 = px1 * t347 + t344;
 
569
 
 
570
    t349 = t89 + t164 + t62 * y2;
 
571
    t350 = - px2 * px3 * t142;
 
572
    t351 = 2 * px2 * px3 * y2 * y3;
 
573
 
 
574
    t352 = - px2 * px3 * t147;
 
575
    t353 = px2 * t142;
 
576
    t354 = y2 * (t89 + t74);
 
577
 
 
578
    t355 = t167 + t354 + t353;
 
579
    t356 = px1 * t355 + t352 + t351 + t350;
 
580
    t357 = py1 * t66;
 
581
 
 
582
    t358 = py1 * t349;
 
583
    t359 = 2 * py2;
 
584
    t360 = - 2 * py3;
 
585
    t361 = - 2 * py2 * y3;
 
586
 
 
587
    t362 = 2 * py3 * y3;
 
588
    t363 = px3 * py2 * t142;
 
589
    t364 = y2 * (t153 + t34);
 
590
 
 
591
    t365 = - px3 * t142;
 
592
    t366 = y2 * (t75 + t164);
 
593
    t367 = t166 + t366 + t365;
 
594
 
 
595
    t368 = py1 * t367;
 
596
    t369 = px1 * (t172 + t171 + t106 * t142);
 
597
    t370 = t35 + t34;
 
598
 
 
599
    t371 = t142 * t370;
 
600
    t372 = y2 * t150;
 
601
    t373 = t372 + t371;
 
602
    t374 = t230 + t229;
 
603
 
 
604
    t375 = py1 * (t352 + t351 + t350);
 
605
    t376 = t157 + t364 + t363;
 
606
    t377 = px1 * t376 + t375;
 
607
 
 
608
    t378 = t75 + t74;
 
609
    t379 = y2 * t207 + t142 * t378;
 
610
    t380 = px1 * t367 + t343 + t342 + t341;
 
611
 
 
612
    t381 = py1 * t209;
 
613
    t382 = py1 * t355;
 
614
    t383 = py1 * t379;
 
615
    t384 = 2 * py2 * y3;
 
616
 
 
617
    t385 = - 2 * py3 * y3;
 
618
    t386 = t385 + t384;
 
619
    t387 = - 2 * py2 * t147;
 
620
    t388 = 2 * py3 * t147;
 
621
 
 
622
    t389 = px2 * py3 * t2;
 
623
    t390 = t389 + t10;
 
624
    t391 = x2 * t390 * y2;
 
625
    t392 = t5 * t228;
 
626
 
 
627
    t393 = - px2 * t2;
 
628
    t394 = t70 + t393;
 
629
    t395 = x2 * t394 * y2;
 
630
    t396 = t5 * t231;
 
631
 
 
632
    t397 = t396 + t395;
 
633
    t398 = py1 * t397;
 
634
    t399 = py2 * t2;
 
635
    t400 = - py3 * t2;
 
636
 
 
637
    t401 = t400 + t399;
 
638
    t402 = x2 * t401 * y2;
 
639
    t403 = t136 + t135;
 
640
    t404 = t5 * t403;
 
641
 
 
642
    t405 = t404 + t402;
 
643
    t406 = px1 * t405;
 
644
    t407 = t1 * (t406 + t398 + t392 + t391);
 
645
 
 
646
    t408 = t65 + t64;
 
647
    t409 = t5 * t408;
 
648
    t410 = x2 * t394;
 
649
    t411 = t410 + t409;
 
650
 
 
651
    t412 = py1 * t411;
 
652
    t413 = t5 * t116;
 
653
    t414 = x2 * t401;
 
654
    t415 = t414 + t413;
 
655
 
 
656
    t416 = px1 * t415;
 
657
    t417 = py2 * t5;
 
658
    t418 = x2 * (t134 + t114);
 
659
    t419 = py3 * t2;
 
660
 
 
661
    t420 = t419 + t418 + t417;
 
662
    t421 = px1 * t420;
 
663
    t422 = t265 * y2;
 
664
    t423 = x2 * t154;
 
665
 
 
666
    t424 = px2 * x2;
 
667
    t425 = (t68 + t424) * y2;
 
668
    t426 = - py2 * x2;
 
669
    t427 = (t133 + t426) * y2;
 
670
 
 
671
    t428 = py3 * x2 * y3;
 
672
    t429 = t20 + t19;
 
673
    t430 = x2 * t429;
 
674
    t431 = - px2 * py3 * t2;
 
675
 
 
676
    t432 = (t431 + t43 + t430) * y2;
 
677
    t433 = t5 * t370;
 
678
    t434 = x2 * t145;
 
679
 
 
680
    t435 = - px2 * x2 * x3;
 
681
    t436 = px2 * t2;
 
682
    t437 = (t436 + t435) * y2;
 
683
    t438 = px3 * t5 * y3;
 
684
 
 
685
    t439 = - px3 * x2 * x3 * y3;
 
686
    t440 = py2 * x2 * x3;
 
687
    t441 = - py2 * t2;
 
688
 
 
689
    t442 = (t441 + t440) * y2;
 
690
    t443 = - py3 * t5 * y3;
 
691
    t444 = py3 * x2 * x3 * y3;
 
692
 
 
693
    t445 = t5 * t287;
 
694
    t446 = t78 + t436;
 
695
    t447 = x2 * t446;
 
696
    t448 = - t2;
 
697
 
 
698
    t449 = t448 + 2 * x2 * x3 - t5;
 
699
    t450 = px1 * t449;
 
700
    t451 = (t98 + t85) * y2;
 
701
    t452 = - x2 * y3;
 
702
 
 
703
    t453 = x3 * y3;
 
704
    t454 = t453 + t452 + (x2 - x3) * y2;
 
705
    t455 = px1 * t454;
 
706
    t456 = t65 + t314;
 
707
 
 
708
    t457 = x2 * t456;
 
709
    t458 = (t78 + t457) * y2;
 
710
    t459 = x2 * (t293 + t203);
 
711
 
 
712
    t460 = - x2 * x3 * y3 + t5 * y3 + (t2 - x2 * x3) * y2;
 
713
    t461 = px1 * t460;
 
714
    t462 = t5 * t253;
 
715
 
 
716
    t463 = t419 + t441;
 
717
    t464 = x2 * t463;
 
718
    t465 = - py2 * t5;
 
719
    t466 = x2 * (t115 + t133);
 
720
 
 
721
    t467 = t2 - 2 * x2 * x3 + t5;
 
722
    t468 = py1 * t467;
 
723
    t469 = py2 * x2;
 
724
    t470 = (t134 + t469) * y2;
 
725
 
 
726
    t471 = - py2 * x2 * y3;
 
727
    t472 = x2 * y3;
 
728
    t473 = - x3 * y3;
 
729
    t474 = t473 + t472 + (x3 - x2) * y2;
 
730
 
 
731
    t475 = py1 * t474;
 
732
    t476 = - 2 * py2 * x3;
 
733
    t477 = t115 + t476;
 
734
    t478 = x2 * t477;
 
735
 
 
736
    t479 = (t419 + t478) * y2;
 
737
    t480 = py2 * t5 * y3;
 
738
    t481 = - 2 * py3 * x3 * y3;
 
739
 
 
740
    t482 = x2 * (t481 + t117);
 
741
    t483 = x2 * x3 * y3 - t5 * y3 + (t448 + x2 * x3) * y2;
 
742
 
 
743
    t484 = py1 * t483;
 
744
    t485 = t431 + t43;
 
745
    t486 = t485 * t142;
 
746
    t487 = t5 * t158;
 
747
 
 
748
    t488 = t446 * t142;
 
749
    t489 = t5 * t168;
 
750
    t490 = t489 + t488;
 
751
    t491 = py1 * t490;
 
752
 
 
753
    t492 = t463 * t142;
 
754
    t493 = t5 * t173;
 
755
    t494 = t493 + t492;
 
756
    t495 = px1 * t494;
 
757
 
 
758
    t496 = x1 * y1 * (t495 + t491 + t487 + t486);
 
759
    t497 = t142 * t119;
 
760
    t498 = x2 * y2 * t259;
 
761
 
 
762
    t499 = t498 + t497;
 
763
    t500 = px1 * t499;
 
764
    t501 = t29 * (t500 + t381 + t151 + t146);
 
765
 
 
766
    t502 = t429 * t142;
 
767
    t503 = x2 * t370;
 
768
    t504 = y2 * (t503 + t125 + t54);
 
769
    t505 = x2 * t158;
 
770
 
 
771
    t506 = - px3 * x3 * t142;
 
772
    t507 = - px2 * x2 * t147;
 
773
    t508 = py3 * x3 * t142;
 
774
 
 
775
    t509 = y2 * (t118 + t471);
 
776
    t510 = py2 * x2 * t147;
 
777
    t511 = - py2 * t142;
 
778
 
 
779
    t512 = y2 * (t138 + t107);
 
780
    t513 = t172 + t512 + t511;
 
781
    t514 = px1 * t513;
 
782
 
 
783
    t515 = y2 * t259 + t142 * t255;
 
784
    t516 = px1 * t515;
 
785
    t517 = py1 * t454;
 
786
 
 
787
    t518 = - py2 * x3 * t142;
 
788
    t519 = t108 + t384;
 
789
    t520 = x2 * t519;
 
790
 
 
791
    t521 = y2 * (t520 + t307 + t135);
 
792
    t522 = - py3 * x2 * t147;
 
793
    t523 = py2 * t142;
 
794
 
 
795
    t524 = y2 * (t108 + t137);
 
796
    t525 = - t147 + 2 * y2 * y3 - t142;
 
797
    t526 = py1 * t525;
 
798
 
 
799
    t527 = x2 * t147 + y2 * (t473 + t452) + x3 * t142;
 
800
    t528 = py1 * t527;
 
801
    t529 = px1 * t474;
 
802
 
 
803
    t530 = px2 * x3 * t142;
 
804
    t531 = px3 * x2 * t147;
 
805
 
 
806
    t532 = - x2 * t147 + y2 * (t453 + t472) - x3 * t142;
 
807
    t533 = px1 * t532;
 
808
 
 
809
    t534 = - px2 * t142;
 
810
    t535 = t147 - 2 * y2 * y3 + t142;
 
811
    t536 = px1 * t535;
 
812
 
 
813
    t537 = t447 + t445;
 
814
    t538 = py1 * t537;
 
815
    t539 = t464 + t462;
 
816
    t540 = px1 * t539;
 
817
 
 
818
    t541 = 2 * px3 * py2 * t2;
 
819
    t542 = - 2 * px2 * py3 * t2;
 
820
    t543 = x2 * t446 * y2;
 
821
 
 
822
    t544 = t5 * t205;
 
823
    t545 = t544 + t543;
 
824
    t546 = py1 * t545;
 
825
    t547 = x2 * t463 * y2;
 
826
 
 
827
    t548 = t5 * t119;
 
828
    t549 = t548 + t547;
 
829
    t550 = px1 * t549;
 
830
    t551 = x2 * t265;
 
831
 
 
832
    t552 = (t389 + t10 + t551) * y2;
 
833
    t553 = t5 * t154;
 
834
    t554 = 2 * px3 * t2;
 
835
 
 
836
    t555 = (t554 + t393 + t110) * y2;
 
837
    t556 = t5 * t90;
 
838
    t557 = py3 * x2 * x3;
 
839
 
 
840
    t558 = - 2 * py3 * t2;
 
841
    t559 = (t558 + t399 + t557) * y2;
 
842
    t560 = py2 * x2 * x3 * y3;
 
843
 
 
844
    t561 = t138 + t361;
 
845
    t562 = t5 * t561;
 
846
    t563 = t390 * t142;
 
847
    t564 = t5 * t150;
 
848
 
 
849
    t565 = - px2 * t2 * t142;
 
850
    t566 = - px3 * t5 * t147;
 
851
    t567 = t566 + t214 + t565;
 
852
 
 
853
    t568 = py1 * t567;
 
854
    t569 = py2 * t2 * t142;
 
855
    t570 = x2 * y2 * (t118 + t135);
 
856
 
 
857
    t571 = py3 * t5 * t147;
 
858
    t572 = t571 + t570 + t569;
 
859
    t573 = px1 * t572;
 
860
    t574 = t86 + t68;
 
861
 
 
862
    t575 = x2 * t574;
 
863
    t576 = (t78 + t575) * y2;
 
864
    t577 = 2 * px2 * x3 * y3;
 
865
 
 
866
    t578 = x2 * (t87 + t577);
 
867
    t579 = px1 * t527;
 
868
 
 
869
    t580 = - t5 * t147 + 2 * x2 * x3 * y2 * y3 - t2 * t142;
 
870
    t581 = px1 * t580;
 
871
    t582 = t305 + t133;
 
872
 
 
873
    t583 = x2 * t582;
 
874
    t584 = (t419 + t583) * y2;
 
875
    t585 = x2 * (t136 + t306);
 
876
 
 
877
    t586 = py1 * t532;
 
878
    t587 = - py3 * t2 * t142;
 
879
    t588 = x2 * y2 * (t136 + t117);
 
880
 
 
881
    t589 = - py2 * t5 * t147;
 
882
    t590 = t5 * t147 - 2 * x2 * x3 * y2 * y3 + t2 * t142;
 
883
 
 
884
    t591 = py1 * t590;
 
885
    t592 = t400 + t466 + t465;
 
886
    t593 = px1 * t592;
 
887
    t594 = t309 + t279;
 
888
 
 
889
    t595 = t198 + t311;
 
890
    t596 = x2 * t378;
 
891
    t597 = t596 + t408 * y2;
 
892
    t598 = py1 * t597;
 
893
 
 
894
    t599 = t256 + t116 * y2;
 
895
    t600 = px1 * t599;
 
896
    t601 = t178 + t366 + t534;
 
897
 
 
898
    t602 = py1 * t601;
 
899
    t603 = t181 + t524 + t523;
 
900
    t604 = px1 * t603;
 
901
    t605 = t265 * t142;
 
902
 
 
903
    t606 = t423 + t144 + t143;
 
904
    t607 = y2 * t606;
 
905
    t608 = x2 * t150;
 
906
    t609 = 2 * py2 * x3 * y3;
 
907
 
 
908
    t610 = t362 + t137;
 
909
    t611 = x2 * t610;
 
910
    t612 = y2 * (t611 + t118 + t609);
 
911
 
 
912
    t613 = py1 * t449;
 
913
    t614 = t419 + t613 + t418 + t417;
 
914
    t615 = py1 * t460;
 
915
 
 
916
    t616 = py1 * t535;
 
917
    t617 = t616 + t172 + t512 + t511;
 
918
    t618 = t134 + t304;
 
919
 
 
920
    t619 = t618 * t142;
 
921
    t620 = - py3 * x2 * y3;
 
922
    t621 = y2 * (t135 + t620);
 
923
 
 
924
    t622 = x2 * (t388 + t180);
 
925
    t623 = px1 * t467;
 
926
    t624 = t623 + t78 + t102 + t101;
 
927
 
 
928
    t625 = px1 * t483;
 
929
    t626 = px1 * t525;
 
930
    t627 = t167 + t626 + t354 + t353;
 
931
 
 
932
    t628 = - 2 * px2 * x3;
 
933
    t629 = t98 + t628;
 
934
    t630 = t629 * t142;
 
935
    t631 = - 2 * px3 * t147;
 
936
 
 
937
    t632 = x2 * (t631 + t177);
 
938
    t633 = - 2 * px2 * py3 * x3 * y3;
 
939
    t634 = t633 + t197;
 
940
 
 
941
    t635 = - 2 * px3 * py2 * t147;
 
942
    t636 = t142 * t403;
 
943
    t637 = x2 * y2 * t173;
 
944
 
 
945
    t638 = t637 + t636;
 
946
    t639 = px1 * t638;
 
947
    t640 = t589 + t588 + t587;
 
948
    t641 = px1 * t640;
 
949
 
 
950
    t642 = px1 * t590;
 
951
    t643 = py1 * t580;
 
952
 
 
953
    t644 = (x0 * (px0 * (y1 * (x1 * (t528 + t522 + t612 + t518)
 
954
                               + t643 + t571 + t570 + t569)
 
955
                         + t29 * t515 + x1 * t638 + t1 * (t615 + t444 + t443 + t442))
 
956
                  + py0 * (y1 * (x1 * (t533 + t531 + t331 + t530)
 
957
                                 + t642 + t566 + t214 + t565)
 
958
                           + x1 * t234 + t29 * t379 + t1 * (t625 + t439 + t438 + t437))
 
959
                  + y1 * (x1 * (px1 * (t622 + t621 + t619) + py1 * (t632 + t299 + t630)
 
960
                                + t608 + t607 + t605)
 
961
                          + t641 + t243 + t564 + t563)
 
962
                  + x1 * (t639 + t235 + x2 * y2 * (t284 + t635) + t142 * t634)
 
963
                  + t29 * (t175 + t170)
 
964
                  + t1 * (px1 * (t482 + t480 + t479) + py1 * (t459 + t79 + t458) + t434
 
965
                          + t433 + t432))
 
966
            + y0 * (x0 * (py0 * (x1 * (t579 + t632 + t299 + t630)
 
967
                                 + t489 + t29 * t627
 
968
                                 + y1 * (x1 * t597 + t625 + t556 + t112 + t555) + t488
 
969
                                 + t624 * t1)
 
970
                          + px0 * (x1 * (t586 + t622 + t621 + t619)
 
971
                                   + t29 * t617 + t493
 
972
                                   + y1 * (x1 * t599 + t615 + t562 + t560 + t559) + t492
 
973
                                   + t614 * t1)
 
974
                          + x1 * (px1 * (t522 + t612 + t518) + py1 * (t531 + t331 + t530)
 
975
                                  + t608 + t607 + t605)
 
976
                          + t29 * (t604 + t602) + t487
 
977
                          + y1 * (x1 * (t600 + t598 + x2 * t595 + t594 * y2)
 
978
                                  + px1 * (t585 + t480 + t584) + py1 * (t578 + t79 + t576) + t267
 
979
                                  + t553 + t552) + t486 + (t593 + t302) * t1)
 
980
                    + px0 * (x1 * (t591 + t589 + t588 + t587)
 
981
                             + t29 * (t586 + t510 + t509 + t508)
 
982
                             + y1 * (x1 * (t484 + t585 + t480 + t584) + t548 + t547) + t415 * t1)
 
983
                    + py0 * (x1 * (t581 + t241 + t240 + t239)
 
984
                             + t29 * (t579 + t507 + t250 + t506)
 
985
                             + y1 * (x1 * (t461 + t578 + t79 + t576) + t544 + t543) + t411 * t1)
 
986
                    + x1 * (t573 + t568 + t564 + t563)
 
987
                    + t29 * (px1 * (t522 + t521 + t518) + py1 * (t531 + t225 + t530) + t505
 
988
                             + t504 + t502)
 
989
                    + y1 * (x1 * (px1 * (t562 + t560 + t559) + py1 * (t556 + t112 + t555)
 
990
                                  + t267 + t553 + t552)
 
991
                            + t550 + t546 + t5 * (t322 + t273) + x2 * (t542 + t541) * y2)
 
992
                    + (t540 + t538) * t1)
 
993
            + t161 * (py0 * (y1 * (x1 * (t536 + t178 + t366 + t534)
 
994
                                   + t533 + t531 + t225 + t530)
 
995
                             + x1 * t169 + t208 + t1 * (t529 + t204 + t292 + t425) + t206)
 
996
                      + px0 * (y1 * (t528 + x1 * (t181 + t526 + t524 + t523) + t522 + t521
 
997
                                     + t518)
 
998
                               + x1 * t174 + t498 + t1 * (t517 + t118 + t428 + t427) + t497)
 
999
                      + x1 * (t516 + t383)
 
1000
                      + y1 * (x1 * (t514 + t382) + px1 * (t510 + t509 + t508)
 
1001
                              + py1 * (t507 + t250 + t506) + t505 + t504
 
1002
                              + t502) + t151
 
1003
                      + t1 * (px1 * (t136 + t471 + t470) + py1 * (t87 + t249 + t451) + t423
 
1004
                              + t422) + t146) + t501 + t496
 
1005
            + t14 * (px0 * (x1 * (t484 + t482 + t480 + t479)
 
1006
                            + t29 * (t475 + t136 + t471 + t470) + t404 + t402
 
1007
                            + (x1 * (t468 + t400 + t466 + t465) + t464 + t462) * y1)
 
1008
                     + py0 * (x1 * (t461 + t459 + t79 + t458)
 
1009
                              + t29 * (t455 + t87 + t249 + t451) + t396 + t395
 
1010
                              + (x1 * (t70 + t450 + t69 + t67) + t447 + t445) * y1)
 
1011
                     + x1 * (px1 * (t444 + t443 + t442) + py1 * (t439 + t438 + t437) + t434
 
1012
                             + t433 + t432)
 
1013
                     + t29 * (px1 * (t118 + t428 + t427) + py1 * (t204 + t292 + t425) + t423
 
1014
                              + t422) + t392 + t391
 
1015
                     + (x1 * (t421 + t104) + t416 + t412) * y1) + t407);
 
1016
    t645 = t5 * t265;
 
1017
 
 
1018
    t646 = t115 + t114 + t132;
 
1019
    t647 = px1 * t646;
 
1020
    t648 = x2 * t485;
 
1021
    t649 = t32 * t5;
 
1022
 
 
1023
    t650 = t70 + t393 + t73 * t5;
 
1024
    t651 = t400 + t399 + t106 * t5;
 
1025
 
 
1026
    t652 = t540 + x1 * (px1 * t651 + py1 * t650 + t389 + t10 + t649) + t538 + t648
 
1027
        + t29 * (t647 + t357 + t20 + t19 + t18) + t645;
 
1028
    t653 = t648 + t645;
 
1029
 
 
1030
    t654 = t392 + t391;
 
1031
    t655 = px1 * t654;
 
1032
    t656 = t309 + t19;
 
1033
    t657 = x2 * t656;
 
1034
 
 
1035
    t658 = (t389 + t657) * y2;
 
1036
    t659 = px3 * py2 * t5 * y3;
 
1037
    t660 = x2 * (t144 + t273);
 
1038
 
 
1039
    t661 = - px3 * py2 * t5;
 
1040
    t662 = t431 + t27 + t661;
 
1041
    t663 = px1 * t662 + t24;
 
1042
 
 
1043
    t664 = t5 * t429;
 
1044
    t665 = x2 * t390;
 
1045
    t666 = t665 + t664;
 
1046
    t667 = px3 * py2 * x2;
 
1047
 
 
1048
    t668 = (t20 + t667) * y2;
 
1049
    t669 = x2 * t485 * y2;
 
1050
    t670 = t5 * t145;
 
1051
    t671 = t670 + t669;
 
1052
 
 
1053
    t672 = px1 * t671;
 
1054
    t673 = t26 + t52;
 
1055
    t674 = x2 * t673;
 
1056
    t675 = (t389 + t674) * y2;
 
1057
 
 
1058
    t676 = x2 * (t633 + t54);
 
1059
    t677 = px3 * t5;
 
1060
    t678 = t436 + t69 + t677;
 
1061
 
 
1062
    t679 = px1 * t678 + t37 + t60 + t59;
 
1063
    t680 = - px3 * x2;
 
1064
 
 
1065
    t681 = t203 + t298 + (t64 + t680) * y2;
 
1066
    t682 = px1 * t545;
 
1067
    t683 = - px3 * t5 * y3;
 
1068
 
 
1069
    t684 = t578 + t683 + (t393 + t575) * y2;
 
1070
    t685 = 2 * py3 * x3;
 
1071
    t686 = t685 + t476;
 
1072
 
 
1073
    t687 = 2 * py2 * t2;
 
1074
    t688 = px1 * (t419 + t441 + t131 * t5);
 
1075
    t689 = - px2 * py3 * x2;
 
1076
 
 
1077
    t690 = 2 * px2 * py3 * x3;
 
1078
    t691 = (t690 + t94 + t689) * y2;
 
1079
 
 
1080
    t692 = t330 + t204 + (t98 + t628 + t424) * y2;
 
1081
    t693 = t134 + t133 + t319;
 
1082
 
 
1083
    t694 = px1 * (t140 + t118 + t117 + t693 * y2);
 
1084
    t695 = (t542 + t43 + t9) * y2;
 
1085
 
 
1086
    t696 = t5 * t312;
 
1087
    t697 = 2 * px2 * t2;
 
1088
    t698 = t5 * t316 + t112 + (t78 + t697 + t110) * y2;
 
1089
 
 
1090
    t699 = x2 * t253;
 
1091
    t700 = t5 * t255;
 
1092
    t701 = x2 * t403;
 
1093
 
 
1094
    t702 = px1 * (t701 + t700 + (t419 + t441 + t699) * y2);
 
1095
    t703 = px2 * py3 * x2 * x3;
 
1096
 
 
1097
    t704 = (t10 + t703) * y2;
 
1098
    t705 = px3 * py2 * x2 * x3 * y3;
 
1099
    t706 = (t20 + t279 + t689) * y2;
 
1100
 
 
1101
    t707 = t439 + t111 + (t70 + t435) * y2;
 
1102
    t708 = t224 + t204 + (t296 + t64 + t424) * y2;
 
1103
 
 
1104
    t709 = - 2 * py2;
 
1105
    t710 = 2 * py3;
 
1106
    t711 = py1 * t678;
 
1107
 
 
1108
    t712 = t459 + t683 + (t393 + t457) * y2;
 
1109
    t713 = x2 * t116;
 
1110
    t714 = t5 * t139;
 
1111
 
 
1112
    t715 = px1 * (t120 + t714 + (t400 + t399 + t713) * y2);
 
1113
    t716 = 2 * px2 * py3;
 
1114
 
 
1115
    t717 = (t94 + (t716 + t15) * x2) * y2;
 
1116
    t718 = - 2 * px2;
 
1117
 
 
1118
    t719 = t221 + t128 + t249 + (t98 + (px3 + t718) * x2) * y2;
 
1119
 
 
1120
    t720 = px1 * (t256 + t136 + t135 + t646 * y2);
 
1121
    t721 = - px2 * py3 * t2 * t142;
 
1122
 
 
1123
    t722 = - px3 * py2 * t5 * t147;
 
1124
    t723 = t722 + t237 + t721;
 
1125
    t724 = - px2 * py3 * x3 * t142;
 
1126
 
 
1127
    t725 = y2 * (t54 + t124);
 
1128
    t726 = px1 * y2 * t257;
 
1129
    t727 = - px3 * py2 * x2 * t147;
 
1130
 
 
1131
    t728 = y2 * (t87 + t127);
 
1132
    t729 = t531 + t728 + t530;
 
1133
    t730 = px2 * py3 * t2 * t142;
 
1134
 
 
1135
    t731 = px3 * py2 * t5 * t147;
 
1136
    t732 = px1 * t397;
 
1137
    t733 = t251 + t299 + t248;
 
1138
 
 
1139
    t734 = px2 * t2 * t142;
 
1140
    t735 = px3 * t5 * t147;
 
1141
    t736 = t735 + t240 + t734;
 
1142
 
 
1143
    t737 = t389 + t10 + t649;
 
1144
    t738 = t731 + t189 + t730;
 
1145
    t739 = px1 * t738;
 
1146
 
 
1147
    t740 = x2 * t165;
 
1148
    t741 = t740 + t204 + t203;
 
1149
    t742 = py1 * y2 * t741;
 
1150
    t743 = py1 * t736;
 
1151
 
 
1152
    t744 = px2 * py3 * t142;
 
1153
    t745 = px1 * t567;
 
1154
    t746 = t148 + t364 + t744;
 
1155
 
 
1156
    t747 = px3 * py2 * t5;
 
1157
    t748 = t389 + t95 + t747;
 
1158
    t749 = (t26 + t122) * y2;
 
1159
 
 
1160
    t750 = x2 * t280;
 
1161
    t751 = (t431 + t750) * y2;
 
1162
    t752 = - px3 * py2 * t5 * y3;
 
1163
 
 
1164
    t753 = x2 * (t322 + t143);
 
1165
    t754 = - px3 * t5;
 
1166
    t755 = t393 + t102 + t754;
 
1167
 
 
1168
    t756 = t128 + t292 + (t68 + t126) * y2;
 
1169
    t757 = x2 * t297;
 
1170
    t758 = x2 * (t204 + t327);
 
1171
 
 
1172
    t759 = t758 + t438 + (t436 + t757) * y2;
 
1173
    t760 = (t94 + t667) * y2;
 
1174
 
 
1175
    t761 = t203 + t249 + (t98 + t680) * y2;
 
1176
    t762 = px1 * (t140 + t253 * y2);
 
1177
 
 
1178
    t763 = - px3 * py2 * x2 * x3;
 
1179
    t764 = (t43 + t763) * y2;
 
1180
    t765 = - px2 * py3 * x2 * x3 * y3;
 
1181
 
 
1182
    t766 = px3 * x2 * x3;
 
1183
    t767 = px2 * x2 * x3 * y3;
 
1184
    t768 = t767 + t79 + (t78 + t766) * y2;
 
1185
 
 
1186
    t769 = px1 * (t120 + t700 + (t419 + t441 + t713) * y2);
 
1187
    t770 = t501 + t496 + t407;
 
1188
 
 
1189
    t771 = px3 * py2 * x3 * t142;
 
1190
    t772 = y2 * (t313 + t633 + t54);
 
1191
 
 
1192
    t773 = px2 * py3 * x2 * t147;
 
1193
    t774 = - px3 * py2 * t142;
 
1194
    t775 = t149 + t346 + t774;
 
1195
 
 
1196
    t776 = y2 * (t317 + t87 + t577);
 
1197
    t777 = t507 + t776 + t506;
 
1198
    t778 = px3 * t142;
 
1199
 
 
1200
    t779 = t177 + t354 + t778;
 
1201
    t780 = y2 * (t144 + t272);
 
1202
    t781 = y2 * (t203 + t292);
 
1203
 
 
1204
    t782 = t531 + t781 + t530;
 
1205
    t783 = px1 * (t336 + t258 + t333);
 
1206
    t784 = t690 + t94;
 
1207
 
 
1208
    t785 = x2 * t784;
 
1209
    t786 = (t431 + t785) * y2;
 
1210
    t787 = x2 * (t125 + t197);
 
1211
 
 
1212
    t788 = x2 * t629;
 
1213
    t789 = x2 * (t221 + t128);
 
1214
    t790 = t789 + t438 + (t436 + t788) * y2;
 
1215
 
 
1216
    t791 = - 2 * py2 * t2;
 
1217
    t792 = 2 * py3 * t2;
 
1218
    t793 = 2 * px2 * py3 * t2;
 
1219
 
 
1220
    t794 = (t793 + t10 + t42) * y2;
 
1221
    t795 = t5 * t324;
 
1222
    t796 = - 2 * px2 * t2;
 
1223
 
 
1224
    t797 = t5 * t329 + t80 + (t70 + t796 + t77) * y2;
 
1225
 
 
1226
    t798 = px1 * (t701 + t714 + (t400 + t399 + t699) * y2);
 
1227
 
 
1228
    t799 = px1 * (t5 * t259 + t401 * t142);
 
1229
    t800 = t429 * y2;
 
1230
    t801 = t503 + t800;
 
1231
 
 
1232
    t802 = t487 + t486;
 
1233
    t803 = t673 * t142;
 
1234
    t804 = - 2 * px2 * py3 * t147;
 
1235
 
 
1236
    t805 = x2 * (t804 + t148);
 
1237
    t806 = 2 * px2 * t147;
 
1238
 
 
1239
    t807 = x2 * (t178 + t806) + t728 + t574 * t142;
 
1240
    t808 = py1 * t755;
 
1241
    t809 = py1 * t779;
 
1242
 
 
1243
    t810 = y2 * (t58 + t144 + t273);
 
1244
    t811 = y2 * (t91 + t293 + t203);
 
1245
 
 
1246
    t812 = t507 + t811 + t506;
 
1247
    t813 = px1 * (t260 + t335 + t254);
 
1248
    t814 = 2 * py2 * t147;
 
1249
 
 
1250
    t815 = - 2 * py3 * t147;
 
1251
    t816 = (t389 + t42) * y2;
 
1252
    t817 = - py2 * py3 * t2;
 
1253
 
 
1254
    t818 = (t817 + py2 * py3 * x2 * x3) * y2;
 
1255
    t819 = - py2 * py3 * t5 * y3;
 
1256
 
 
1257
    t820 = py2 * py3 * x2 * x3 * y3;
 
1258
    t821 = px1 * (t820 + t819 + t818);
 
1259
    t822 = - py2 * py3 * t5;
 
1260
 
 
1261
    t823 = 2 * py2 * py3 * x2 * x3;
 
1262
    t824 = px1 * (t817 + t823 + t822);
 
1263
    t825 = (t431 + t9) * y2;
 
1264
 
 
1265
    t826 = py2 * py3 * t2;
 
1266
    t827 = (t826 - py2 * py3 * x2 * x3) * y2;
 
1267
    t828 = py2 * py3 * t5 * y3;
 
1268
 
 
1269
    t829 = - py2 * py3 * x2 * x3 * y3;
 
1270
    t830 = px1 * (t829 + t828 + t827);
 
1271
 
 
1272
    t831 = (py2 * py3 * x2 - py2 * py3 * x3) * y2;
 
1273
    t832 = - py2 * py3 * x2 * y3;
 
1274
 
 
1275
    t833 = py2 * py3 * x3 * y3;
 
1276
    t834 = px1 * (t833 + t832 + t831);
 
1277
 
 
1278
    t835 = (t690 + t94 + t122) * y2;
 
1279
    t836 = px1 * t693;
 
1280
    t837 = - py2 * t5 * y3;
 
1281
 
 
1282
    t838 = t560 + t837 + (t400 + t557) * y2;
 
1283
    t839 = x2 * t205;
 
1284
 
 
1285
    t840 = py1 * (t839 + x2 * t408 * y2);
 
1286
    t841 = (t20 + t51) * y2;
 
1287
    t842 = - py3 * x2;
 
1288
 
 
1289
    t843 = py2 * x2 * y3;
 
1290
    t844 = t135 + t843 + (t115 + t842) * y2;
 
1291
 
 
1292
    t845 = py1 * (t740 + t87 + t128 + (t98 + t68 + t63) * y2);
 
1293
    t846 = py2 * py3 * t5;
 
1294
 
 
1295
    t847 = - 2 * py2 * py3 * x2 * x3;
 
1296
    t848 = - py2 * x2 * x3;
 
1297
    t849 = - py3 * x2 * x3 * y3;
 
1298
 
 
1299
    t850 = t849 + t480 + (t419 + t848) * y2;
 
1300
    t851 = (py2 * py3 * x3 - py2 * py3 * x2) * y2;
 
1301
 
 
1302
    t852 = py2 * py3 * x2 * y3;
 
1303
    t853 = - py2 * py3 * x3 * y3;
 
1304
    t854 = x2 * t561;
 
1305
 
 
1306
    t855 = t854 + t136 + (t305 + t133 + t469) * y2;
 
1307
    t856 = py2 * py3 * t2 * t142;
 
1308
 
 
1309
    t857 = - 2 * py2 * py3 * x2 * x3 * y2 * y3;
 
1310
    t858 = py2 * py3 * t5 * t147;
 
1311
 
 
1312
    t859 = px1 * (t858 + t857 + t856);
 
1313
    t860 = - py2 * py3 * x3 * t142;
 
1314
 
 
1315
    t861 = y2 * (t833 + t852);
 
1316
    t862 = - py2 * py3 * x2 * t147;
 
1317
 
 
1318
    t863 = px1 * (t862 + t861 + t860);
 
1319
    t864 = - py2 * py3 * t2 * t142;
 
1320
 
 
1321
    t865 = 2 * py2 * py3 * x2 * x3 * y2 * y3;
 
1322
    t866 = - py2 * py3 * t5 * t147;
 
1323
 
 
1324
    t867 = py3 * t2 * t142;
 
1325
    t868 = py2 * t5 * t147;
 
1326
    t869 = t868 + t570 + t867;
 
1327
 
 
1328
    t870 = py2 * py3 * x3 * t142;
 
1329
    t871 = y2 * (t853 + t832);
 
1330
    t872 = py2 * py3 * x2 * t147;
 
1331
 
 
1332
    t873 = - py3 * x3 * t142;
 
1333
    t874 = - py2 * x2 * t147;
 
1334
    t875 = t874 + t521 + t873;
 
1335
 
 
1336
    t876 = py2 * x3 * t142;
 
1337
    t877 = py3 * x2 * t147;
 
1338
    t878 = t877 + t509 + t876;
 
1339
 
 
1340
    t879 = t287 * t142;
 
1341
    t880 = t596 + t87 + t128;
 
1342
    t881 = y2 * t880;
 
1343
    t882 = x2 * t207;
 
1344
 
 
1345
    t883 = py1 * (t882 + t881 + t879);
 
1346
    t884 = py1 * t662;
 
1347
 
 
1348
    t885 = px1 * (t826 + t847 + t846);
 
1349
    t886 = 2 * px3 * py2;
 
1350
 
 
1351
    t887 = (t94 + (t31 + t886) * x2) * y2;
 
1352
    t888 = px1 * (t853 + t852 + t851);
 
1353
 
 
1354
    t889 = py1 * t738;
 
1355
    t890 = px1 * (t866 + t865 + t864);
 
1356
 
 
1357
    t891 = px1 * (t872 + t871 + t870);
 
1358
    t892 = t656 * t142;
 
1359
    t893 = x2 * (t157 + t635);
 
1360
 
 
1361
    t894 = t221 + t577;
 
1362
    t895 = x2 * t253 * y2;
 
1363
    t896 = t701 + t895;
 
1364
    t897 = px1 * t896;
 
1365
 
 
1366
    t898 = (t20 + t279 + t122) * y2;
 
1367
 
 
1368
    t899 = py1 * (t596 + t204 + t203 + (t65 + t64 + t97) * y2);
 
1369
    t900 = t385 + t107;
 
1370
 
 
1371
    t901 = x2 * t900;
 
1372
    t902 = t901 + t136 + (t115 + t476 + t469) * y2;
 
1373
    t903 = px1 * t869;
 
1374
 
 
1375
    t904 = t874 + t612 + t873;
 
1376
    t905 = t408 * t142;
 
1377
    t906 = y2 * t741;
 
1378
    t907 = x2 * t168;
 
1379
 
 
1380
    t908 = py1 * (t907 + t906 + t905);
 
1381
    t909 = - py2 * py3 * t142;
 
1382
 
 
1383
    t910 = 2 * py2 * py3 * y2 * y3;
 
1384
    t911 = - py2 * py3 * t147;
 
1385
 
 
1386
    t912 = px1 * (t911 + t910 + t909);
 
1387
    t913 = t912 + py1 * t376;
 
1388
 
 
1389
    t914 = t481 + t117 + t428 + (t133 + (py3 + t709) * x2) * y2;
 
1390
    t915 = 2 * px3;
 
1391
 
 
1392
    t916 = t138 + t137 + t131 * y2;
 
1393
    t917 = px1 * t916;
 
1394
 
 
1395
    t918 = py1 * (t167 + t166 + t73 * t142);
 
1396
    t919 = py3 * t142;
 
1397
    t920 = t171 + t524 + t919;
 
1398
 
 
1399
    t921 = px1 * t920;
 
1400
    t922 = py2 * py3 * t142;
 
1401
    t923 = - 2 * py2 * py3 * y2 * y3;
 
1402
 
 
1403
    t924 = py2 * py3 * t147;
 
1404
    t925 = py1 * t513 + t924 + t923 + t922;
 
1405
    t926 = py1 * t420;
 
1406
 
 
1407
    t927 = py1 * t640;
 
1408
    t928 = t685 + t114;
 
1409
    t929 = x2 * (t172 + t814) + t621 + t928 * t142;
 
1410
 
 
1411
    t930 = px1 * (t924 + t923 + t922);
 
1412
    t931 = t930 + py1 * t347;
 
1413
 
 
1414
    t932 = py1 * t920 + t911 + t910 + t909;
 
1415
    t933 = t315 + t222;
 
1416
    t934 = py1 * t654;
 
1417
 
 
1418
    t935 = (t10 + t750) * y2;
 
1419
    t936 = t824 + py1 * t263;
 
1420
    t937 = py1 * t671;
 
1421
 
 
1422
    t938 = (t19 + t689) * y2;
 
1423
    t939 = (t10 + t785) * y2;
 
1424
    t940 = t296 + t314;
 
1425
 
 
1426
    t941 = py1 * (t78 + t436 + t62 * t5);
 
1427
    t942 = (t26 + t52 + t667) * y2;
 
1428
 
 
1429
    t943 = py1 * (t740 + t204 + t203 + t99 * y2);
 
1430
 
 
1431
    t944 = t611 + t118 + (t134 + t304 + t426) * y2;
 
1432
    t945 = (t431 + t541 + t42) * y2;
 
1433
 
 
1434
    t946 = t5 * t199;
 
1435
    t947 = t5 * t900 + t560 + (t419 + t791 + t557) * y2;
 
1436
    t948 = x2 * t287;
 
1437
 
 
1438
    t949 = t5 * t378;
 
1439
    t950 = py1 * (t289 + t949 + (t78 + t436 + t948) * y2);
 
1440
 
 
1441
    t951 = - py3 * t5;
 
1442
    t952 = t441 + t466 + t951;
 
1443
    t953 = py1 * t952 + t826 + t847 + t846;
 
1444
 
 
1445
    t954 = py3 * x2;
 
1446
    t955 = t117 + t620 + (t114 + t954) * y2;
 
1447
    t956 = py1 * t549;
 
1448
 
 
1449
    t957 = py3 * t5 * y3;
 
1450
    t958 = t585 + t957 + (t399 + t583) * y2;
 
1451
    t959 = (t389 + t763) * y2;
 
1452
 
 
1453
    t960 = (t309 + t19 + t667) * y2;
 
1454
    t961 = - 2 * px3;
 
1455
    t962 = px1 * t952;
 
1456
    t963 = x2 * t408;
 
1457
 
 
1458
    t964 = t5 * t165;
 
1459
    t965 = py1 * (t839 + t964 + (t70 + t393 + t963) * y2);
 
1460
 
 
1461
    t966 = t482 + t957 + (t399 + t478) * y2;
 
1462
    t967 = - 2 * px3 * py2;
 
1463
 
 
1464
    t968 = (t26 + (t16 + t967) * x2) * y2;
 
1465
 
 
1466
    t969 = t307 + t135 + t471 + (t134 + (t130 + t359) * x2) * y2;
 
1467
 
 
1468
    t970 = py1 * (t596 + t87 + t128 + t66 * y2);
 
1469
    t971 = t444 + t837 + (t400 + t440) * y2;
 
1470
 
 
1471
    t972 = t520 + t118 + (t685 + t114 + t426) * y2;
 
1472
    t973 = py1 * t405;
 
1473
 
 
1474
    t974 = t877 + t621 + t876;
 
1475
    t975 = - py2 * t2 * t142;
 
1476
    t976 = - py3 * t5 * t147;
 
1477
 
 
1478
    t977 = t976 + t588 + t975;
 
1479
    t978 = py1 * y2 * t880;
 
1480
    t979 = y2 * (t136 + t843);
 
1481
 
 
1482
    t980 = t522 + t979 + t518;
 
1483
    t981 = py1 * t276;
 
1484
    t982 = py1 * t572;
 
1485
    t983 = px1 * y2 * t334;
 
1486
 
 
1487
    t984 = px1 * t977;
 
1488
    t985 = (t94 + t51) * y2;
 
1489
    t986 = (t43 + t657) * y2;
 
1490
 
 
1491
    t987 = (t26 + t689) * y2;
 
1492
    t988 = t117 + t471 + (t134 + t954) * y2;
 
1493
 
 
1494
    t989 = py1 * (t740 + t287 * y2);
 
1495
    t990 = (t431 + t703) * y2;
 
1496
    t991 = - py3 * x2 * x3;
 
1497
 
 
1498
    t992 = - py2 * x2 * x3 * y3;
 
1499
    t993 = t992 + t480 + (t419 + t991) * y2;
 
1500
 
 
1501
    t994 = py1 * (t839 + t949 + (t78 + t436 + t963) * y2);
 
1502
    t995 = py3 * t5;
 
1503
 
 
1504
    t996 = t399 + t418 + t995;
 
1505
    t997 = t135 + t428 + (t133 + t842) * y2;
 
1506
    t998 = x2 * t928;
 
1507
 
 
1508
    t999 = x2 * (t118 + t609);
 
1509
    t1000 = t999 + t443 + (t441 + t998) * y2;
 
1510
 
 
1511
    t1001 = y2 * (t901 + t136 + t306);
 
1512
    t1002 = t510 + t1001 + t508;
 
1513
    t1003 = - py3 * t142;
 
1514
 
 
1515
    t1004 = t180 + t512 + t1003;
 
1516
    t1005 = y2 * (t117 + t428);
 
1517
    t1006 = t522 + t1005 + t518;
 
1518
 
 
1519
    t1007 = py1 * (t907 + t881 + t905);
 
1520
    t1008 = y2 * (t854 + t481 + t117);
 
1521
 
 
1522
    t1009 = t510 + t1008 + t508;
 
1523
    t1010 = 2 * px3 * t147;
 
1524
 
 
1525
    t1011 = py1 * (t5 * t207 + t394 * t142);
 
1526
    t1012 = t784 * t142;
 
1527
 
 
1528
    t1013 = 2 * px3 * py2 * t147;
 
1529
    t1014 = x2 * (t149 + t1013);
 
1530
 
 
1531
    t1015 = py1 * (t882 + t906 + t879);
 
1532
    t1016 = x2 * (t181 + t387) + t979 + t582 * t142;
 
1533
 
 
1534
    t1017 = (t43 + t674) * y2;
 
1535
    t1018 = x2 * t618;
 
1536
    t1019 = x2 * (t307 + t135);
 
1537
 
 
1538
    t1020 = t1019 + t443 + (t441 + t1018) * y2;
 
1539
    t1021 = - 2 * px3 * t2;
 
1540
 
 
1541
    t1022 = - 2 * px3 * py2 * t2;
 
1542
    t1023 = (t389 + t1022 + t9) * y2;
 
1543
    t1024 = t5 * t57;
 
1544
 
 
1545
    t1025 = t5 * t610 + t849 + (t400 + t687 + t848) * y2;
 
1546
 
 
1547
    t1026 = py1 * (t289 + t964 + (t70 + t393 + t948) * y2);
 
1548
    t1027 = px1 * t996;
 
1549
 
 
1550
    t1028 = px1 * t1004;
 
1551
    t1029 = x2 * t429 * y2;
 
1552
    t1030 = (t436 + t110) * y2;
 
1553
 
 
1554
    t1031 = (t441 + t557) * y2;
 
1555
    t1032 = (t393 + t77) * y2;
 
1556
    t1033 = (t399 + t848) * y2;
 
1557
 
 
1558
    t1034 = (t26 + t94 + t18) * y2;
 
1559
    t1035 = (t64 + t85) * y2;
 
1560
    t1036 = (t114 + t469) * y2;
 
1561
 
 
1562
    t1037 = (t98 + t628 + t126) * y2;
 
1563
    t1038 = (t134 + t304 + t842) * y2;
 
1564
 
 
1565
    t1039 = (t20 + t19 + t96) * y2;
 
1566
    t1040 = (t296 + t64 + t126) * y2;
 
1567
 
 
1568
    t1041 = (t685 + t114 + t842) * y2;
 
1569
    t1042 = (t98 + (t961 + px2) * x2) * y2;
 
1570
 
 
1571
    t1043 = t456 * t142;
 
1572
    t1044 = x2 * (t1010 + t166);
 
1573
 
 
1574
    t1045 = (t134 + (t710 + t105) * x2) * y2;
 
1575
    t1046 = t477 * t142;
 
1576
 
 
1577
    t1047 = x2 * (t815 + t171);
 
1578
    t1048 = t32 * t142;
 
1579
    t1049 = t171 + t526 + t524 + t919;
 
1580
 
 
1581
    t1050 = t536 + t166 + t366 + t365;
 
1582
    t1051 = (t389 + t10 + t430) * y2;
 
1583
 
 
1584
    t1052 = (t393 + t766) * y2;
 
1585
    t1053 = (t399 + t991) * y2;
 
1586
    t1054 = t17 * t5;
 
1587
 
 
1588
    t1055 = (t431 + t43 + t551) * y2;
 
1589
    t1056 = (t1021 + t436 + t77) * y2;
 
1590
    t1057 = t5 * t223;
 
1591
 
 
1592
    t1058 = (t792 + t441 + t848) * y2;
 
1593
    t1059 = t5 * t519;
 
1594
    t1060 = t338 * y2;
 
1595
 
 
1596
    t1061 = (t86 + t68 + t680) * y2;
 
1597
    t1062 = (t305 + t133 + t954) * y2;
 
1598
 
 
1599
    t1063 = (t115 + t426) * y2;
 
1600
    t1064 = (t400 + t1018) * y2;
 
1601
    t1065 = (t65 + t424) * y2;
 
1602
 
 
1603
    t1066 = (t70 + t788) * y2;
 
1604
    t1067 = (t70 + t757) * y2;
 
1605
    t1068 = (t400 + t998) * y2;
 
1606
 
 
1607
    t1069 = t21 * y2;
 
1608
    t1070 = (t68 + (t915 + t61) * x2) * y2;
 
1609
 
 
1610
    t1071 = (t133 + (t360 + py2) * x2) * y2;
 
1611
    t1072 = (t115 + t476 + t954) * y2;
 
1612
 
 
1613
    t1073 = (t65 + t314 + t680) * y2;
 
1614
 
 
1615
    trans->m[0][0]
 
1616
        = (x0 * (px0 * (x1 * (px1 * (y2 * (t388 + t387) + t142 * t386)
 
1617
                              + t383 + t372 + t371)
 
1618
                        + y1 * (x1 * (t369 + t382 + t156 + t346 + t345)
 
1619
                                + t337 + py1 * t301 + t285 + t283 + t281) + t381 + t151
 
1620
                        + t1 * (t141 + py1 * t92 + t58 + t54 + t53) + t146)
 
1621
                 + py0 * (y1 * (x1 * t380 + px1 * t332 + t219 + t218 + t217)
 
1622
                          + px1 * t234 + px1 * x1 * t379 + t1 * (px1 * t129 + t49 + t48 + t47))
 
1623
                 + y1 * (x1 * t377 + px1 * (t202 + t326 + t196) + t195) + px1 * t374
 
1624
                 + px1 * x1 * t373 + t1 * (px1 * (t125 + t124 + t123) + t269))
 
1625
           + y0 * (x0 * (px0 * (t261 + x1 * (t369 + t368 + t157 + t364 + t363) + py1 * t227
 
1626
                                + t202
 
1627
                                + y1
 
1628
                                * (x1
 
1629
                                   * (px1 * (t362 + t361 + (t360 + t359) * y2)
 
1630
                                      + t358 + t153 + t55 + t339)
 
1631
                                   + t320 + py1 * t294 + t144 + t273 + t272 + t271)
 
1632
                                + t201 + t196 + (t357 + t20 + t19 + t18) * t1)
 
1633
                         + py0 * (x1 * t356 + px1 * t252 + t194
 
1634
                                  + y1 * (px1 * t318 + px1 * x1 * t349 + t84 + t83 + t82)
 
1635
                                  + t193 + t192 + px1 * t99 * t1) + x1 * t348
 
1636
                         + px1 * (t247 + t246 + t244) + t278
 
1637
                         + y1 * (px1 * (t313 + t54 + t310) + t50 + px1 * x1 * t340)
 
1638
                         + px1 * t338 * t1)
 
1639
                   + px0 * (x1 * (t337 + py1 * t332 + t202 + t326 + t196)
 
1640
                            + t321 + px1 * t29 * t182 + t190
 
1641
                            + y1 * (x1 * (t320 + py1 * t318 + t313 + t54 + t310)
 
1642
                                    + px1 * (x2 * t308 + x2 * (t305 + t304) * y2) + t303 + t267
 
1643
                                    + t266) + t189 + t188 + (t302 + t10 + t27 + t25) * t1)
 
1644
                   + py0 * (x1 * (px1 * t301 + t194 + t193 + t192)
 
1645
                            + t295 + px1 * t29 * t179 + t186
 
1646
                            + y1 * (x1 * (px1 * t294 + t49 + t48 + t47) + px1 * t290) + t185 + t184
 
1647
                            + (t286 + t3 + t23 + t22) * t1)
 
1648
                   + x1 * (px1 * (t285 + t283 + t281) + t278) + t277 + t275 + px1 * t29 * t274
 
1649
                   + y1 * (x1 * (px1 * (t144 + t273 + t272 + t271) + t269) + px1 * t268)
 
1650
                   + (t264 + t262) * t1)
 
1651
           + px0 * (y1 * (x1 * (t261 + py1 * t252 + t247 + t246 + t244)
 
1652
                          + t243 + t238 + t237 + t236)
 
1653
                    + x1 * (t235 + t230 + t229) + px1 * t29 * t174
 
1654
                    + t1 * (t121 + py1 * t81 + t46 + t45 + t44))
 
1655
           + py0 * (y1 * (x1 * (px1 * t227 + t219 + t218 + t217)
 
1656
                          + px1 * t216 + t212 + t211 + t210)
 
1657
                    + px1 * t29 * t169 + px1 * x1 * t209 + t1 * (px1 * t113 + t40 + t39 + t38))
 
1658
           + y1 * (x1 * (px1 * (t202 + t201 + t196) + t195) + px1 * t191 + t187)
 
1659
           + px0 * t161 * t183 + px1 * t29 * t160 + px1 * x1 * t152
 
1660
           + t14 * (px0 * (x1 * (t141 + py1 * t129 + t125 + t124 + t123)
 
1661
                           + t121 + py1 * t113 + px1 * t29 * t109 + t13 + t12 + t11
 
1662
                           + (t104 + t43 + x1 * (t100 + t26 + t94 + t96) + t95 + t93) * y1)
 
1663
                    + py0 * (x1 * (px1 * t92 + t84 + t83 + t82)
 
1664
                             + px1 * t81 + px1 * t29 * t76 + t7 + t6 + t4
 
1665
                             + (px1 * t71 + t37 + px1 * x1 * t66 + t60 + t59) * y1)
 
1666
                    + x1 * (px1 * (t58 + t54 + t53) + t50) + px1 * (t46 + t45 + t44) + t41
 
1667
                    + px1 * t29 * t36 + (px1 * t28 + t24 + px1 * x1 * t21) * y1)
 
1668
           + t1 * (px1 * (t13 + t12 + t11) + t8));
 
1669
 
 
1670
    trans->m[0][1] =
 
1671
        (t161 * (px0 * (x1 * (t382 + t156 + t346 + t345)
 
1672
                        + py1 * t733 + t247
 
1673
                        + y1 * (t694 + x1 * (t358 + t153 + t55 + t339) + py1 * t681
 
1674
                                + t144 + t282 + t668) + t726 + t283 + t244
 
1675
                        + px1 * t646 * t1)
 
1676
                 + py0 * (x1 * (px1 * t601 + t343 + t342 + t341)
 
1677
                          + px1 * t729 + t219
 
1678
                          + y1 * (px1 * t692 + px1 * x1 * t76 + t49 + t48 + t47) + t218
 
1679
                          + t217 + px1 * t66 * t1) + x1 * (px1 * t746 + t375)
 
1680
                 + px1 * (t727 + t725 + t724) + t195
 
1681
                 + y1 * (px1 * (t325 + t143 + t691) + t269 + px1 * x1 * t36)
 
1682
                 + px1 * t21 * t1)
 
1683
         + x0 * (py0 * (t29 * t356 + t745 + t212
 
1684
                        + y1
 
1685
                        * (x1 * (px1 * t719 + t84 + t83 + t82)
 
1686
                           + px1 * t698 + t40 + t39 + t38) + px1 * x1 * y2 * t741
 
1687
                        + t211 + t210 + px1 * t650 * t1)
 
1688
                 + px0 * (t29 * (t602 + t148 + t364 + t744)
 
1689
                          + t743 + t722
 
1690
                          + y1 * (x1 * (t720 + py1 * t708 + t200 + t143 + t706)
 
1691
                                  + t702 + py1 * t684 + t676 + t659 + t675)
 
1692
                          + x1 * (t607 + px1 * y2 * (x2 * (t362 + t361) + t481 + t609) + t742)
 
1693
                          + t237 + t721 + px1 * t651 * t1) + t29 * t348 + t739 + t187
 
1694
                 + y1 * (x1 * (px1 * (t125 + t197 + t245 + t717) + t50)
 
1695
                         + px1 * (t696 + t13 + t695) + t8) + px1 * x1 * y2 * t606
 
1696
                 + px1 * t737 * t1)
 
1697
         + py0 * (x1 * (px1 * t736 + t186 + t185 + t184)
 
1698
                  + t29 * (px1 * t733 + t194 + t193 + t192)
 
1699
                  + y1 * (x1 * (px1 * t712 + t7 + t6 + t4) + t732) + px1 * t537 * t1)
 
1700
         + px0 * (x1 * (t568 + t731 + t189 + t730)
 
1701
                  + t29 * (py1 * t729 + t727 + t726 + t725 + t724)
 
1702
                  + y1 * (x1 * (t715 + py1 * t707 + t705 + t12 + t704) + t546 + t670 + t669)
 
1703
                  + px1 * t539 * t1) + x1 * (px1 * t723 + t275)
 
1704
         + t29 * (px1 * (t247 + t283 + t244) + t278)
 
1705
         + y0 * (x0 * (px0 * (x1 * (t720 + py1 * t719 + t125 + t197 + t245 + t717)
 
1706
                              + t715 + py1 * t712 + t29 * (t162 + t35 + t34 + t33) + t660
 
1707
                              + t659 + t658
 
1708
                              + (t688 + t711 + t431
 
1709
                                 + x1
 
1710
                                 * (px1 * (t305 + t304 + (t710 + t709) * x2)
 
1711
                                    + t100 + t26 + t94 + t96) + t27 + t661)
 
1712
                              * y1)
 
1713
                       + py0 * (x1 * (px1 * t708 + t49 + t48 + t47)
 
1714
                                + px1 * t707 + px1 * t29 * t349 + t40 + t39 + t38
 
1715
                                + (t286 + t3 + px1 * x1 * t99 + t23 + t22) * y1)
 
1716
                       + x1 * (px1 * (t200 + t143 + t706) + t269) + px1 * (t705 + t12 + t704)
 
1717
                       + t8 + px1 * t29 * t340 + (t264 + t262 + px1 * x1 * t338) * y1)
 
1718
                 + px0 * (x1 * (t702 + py1 * t698 + t696 + t13 + t695)
 
1719
                          + t29 * (t694 + py1 * t692 + t325 + t143 + t691) + t398 + t392 + t391
 
1720
                          + (x1 * (t688 + t104 + t43 + t95 + t93)
 
1721
                             + px1 * (x2 * (t558 + t687) + t5 * t686) + t412 + t665 + t664)
 
1722
                          * y1)
 
1723
                 + py0 * (x1 * (px1 * t684 + t7 + t6 + t4) + t682
 
1724
                          + t29 * (px1 * t681 + t84 + t83 + t82)
 
1725
                          + (px1 * t411 + x1 * t679) * y1)
 
1726
                 + x1 * (px1 * (t676 + t659 + t675) + t41) + t672
 
1727
                 + t29 * (px1 * (t144 + t282 + t668) + t50) + (px1 * t666 + x1 * t663) * y1)
 
1728
         + y1 * (x1 * (px1 * (t660 + t659 + t658) + t41) + t655) + px1 * t653 * t1
 
1729
         + px0 * t652 * t14)
 
1730
        ;
 
1731
 
 
1732
    trans->m[0][2] =
 
1733
        (x0 * (px0 * (y1 * (x1 * (t813 + py1 * t807 + t805 + t725 + t803)
 
1734
                            + t799 + t568 + t731 + t189 + t730)
 
1735
                      + x1 * (px1 * (x2 * y2 * (t815 + t814) + t142 * t308)
 
1736
                              + t235 + t230 + t229) + t29 * (t170 + t159 + t155)
 
1737
                      + t1 * (t769 + py1 * t759 + t753 + t752 + t751))
 
1738
               + py0 * (y1 * (x1 * (px1 * t812 + t194 + t193 + t192)
 
1739
                              + t295 + t186 + t185 + t184)
 
1740
                        + px1 * x1 * t234 + px1 * t29 * t379
 
1741
                        + t1 * (px1 * t768 + t7 + t6 + t4))
 
1742
               + y1 * (x1 * (px1 * (t773 + t810 + t771) + t278) + t277 + t275)
 
1743
               + px1 * x1 * t374 + px1 * t29 * t373
 
1744
               + t1 * (px1 * (t765 + t45 + t764) + t41))
 
1745
         + y0 * (x0 * (px0 * (x1 * (t813 + py1 * t812 + t773 + t810 + t771)
 
1746
                              + t495 + t29 * (t809 + t149 + t346 + t774)
 
1747
                              + y1
 
1748
                              * (x1
 
1749
                                 * (px1 * (x2 * t386 + t686 * y2)
 
1750
                                    + t598 + t503 + t800)
 
1751
                                 + t798 + py1 * t790 + t787 + t752 + t786)
 
1752
                              + (t808 + t389 + t95 + t747) * t1)
 
1753
                       + py0 * (x1 * (px1 * t807 + t219 + t218 + t217)
 
1754
                                + px1 * t490 + t29 * t380
 
1755
                                + y1 * (px1 * x1 * t597 + px1 * t797 + t7 + t6 + t4)
 
1756
                                + t679 * t1)
 
1757
                       + x1 * (px1 * (t805 + t725 + t803) + t195) + px1 * t802
 
1758
                       + t29 * t377
 
1759
                       + y1 * (px1 * x1 * t801 + px1 * (t795 + t46 + t794) + t41)
 
1760
                       + t663 * t1)
 
1761
                 + px0 * (x1 * (t799 + t243 + t238 + t237 + t236)
 
1762
                          + t29 * (t783 + py1 * t777 + t773 + t772 + t771)
 
1763
                          + y1 * (x1 * (t798 + py1 * t797 + t795 + t46 + t794)
 
1764
                                  + px1 * (t5 * (t481 + t609) + x2 * (t792 + t791) * y2)
 
1765
                                  + t546 + t670 + t669) + (t538 + t648 + t645) * t1)
 
1766
                 + py0 * (x1 * (t745 + t212 + t211 + t210)
 
1767
                          + t29 * (px1 * t782 + t219 + t218 + t217)
 
1768
                          + y1 * (x1 * (px1 * t790 + t40 + t39 + t38) + t682)
 
1769
                          + px1 * t411 * t1) + x1 * (t739 + t187)
 
1770
                 + t29 * (px1 * (t727 + t780 + t724) + t195)
 
1771
                 + y1 * (x1 * (px1 * (t787 + t752 + t786) + t8) + t672)
 
1772
                 + px1 * t666 * t1)
 
1773
         + t161 * (px0 * (y1
 
1774
                          * (t783 + x1 * (t368 + t157 + t364 + t363) + py1 * t782
 
1775
                             + t727 + t780 + t724)
 
1776
                          + x1 * (t383 + t372 + t371) + t500
 
1777
                          + t1 * (t762 + py1 * t756 + t125 + t272 + t749))
 
1778
                   + py0 * (y1
 
1779
                            * (x1 * (px1 * t779 + t352 + t351 + t350)
 
1780
                               + px1 * t777 + t194 + t193 + t192)
 
1781
                            + px1 * x1 * t169 + px1 * t209
 
1782
                            + t1 * (px1 * t761 + t84 + t83 + t82))
 
1783
                   + y1 * (x1 * (px1 * t775 + t344) + px1 * (t773 + t772 + t771)
 
1784
                           + t278) + px1 * x1 * t160
 
1785
                   + px1 * t152 + t1 * (px1 * (t144 + t245 + t760) + t50))
 
1786
         + px0 * t770
 
1787
         + t14 * (px0 * (x1 * (t769 + py1 * t768 + t765 + t45 + t764)
 
1788
                         + t29 * (t762 + py1 * t761 + t144 + t245 + t760) + t406
 
1789
                         + (t412 + x1 * (t711 + t431 + t27 + t661) + t665 + t664)
 
1790
                         * y1)
 
1791
                  + py0 * (x1 * (px1 * t759 + t40 + t39 + t38)
 
1792
                           + t732 + t29 * (px1 * t756 + t49 + t48 + t47)
 
1793
                           + (px1 * t537 + x1 * (px1 * t755 + t3 + t23 + t22)) * y1)
 
1794
                  + x1 * (px1 * (t753 + t752 + t751) + t8) + t655
 
1795
                  + t29 * (px1 * (t125 + t272 + t749) + t269)
 
1796
                  + (x1 * (px1 * t748 + t262) + px1 * t653) * y1));
 
1797
 
 
1798
    trans->m[1][0] = (x0 * (py0 * (x1 * (t516 + py1 * (y2 * (t631 + t806) + t142 * t933) + t372
 
1799
                                         + t371)
 
1800
                                   + y1 * (px1 * t929 + x1 * (t514 + t918 + t157 + t364 + t363) + t908
 
1801
                                           + t893 + t725 + t892) + t500 + t151
 
1802
                                   + t1 * (px1 * t855 + t845 + t325 + t125 + t835) + t146)
 
1803
                            + px0 * (y1 * (x1 * t932 + py1 * t904 + t872 + t871 + t870)
 
1804
                                     + py1 * x1 * t515 + py1 * t638
 
1805
                                     + t1 * (py1 * t844 + t833 + t832 + t831))
 
1806
                            + y1 * (x1 * t931 + t863 + py1 * (t247 + t810 + t244)) + py1 * t374
 
1807
                            + py1 * x1 * t373 + t1 * (t888 + py1 * (t54 + t282 + t841)))
 
1808
                      + y0 * (px0 * (x1 * (py1 * t929 + t862 + t861 + t860)
 
1809
                                     + t927 + py1 * t29 * t182 + t858
 
1810
                                     + y1 * (py1 * t896 + x1 * (py1 * t914 + t833 + t832 + t831)) + t857
 
1811
                                     + t856 + (t926 + t817 + t823 + t822) * t1)
 
1812
                              + x0 * (px0 * (x1 * t925 + py1 * t878 + t862
 
1813
                                             + y1
 
1814
                                             * (py1 * t902 + py1 * x1 * t916 + t853 + t852
 
1815
                                                + t851) + t861 + t860
 
1816
                                             + py1 * t693 * t1)
 
1817
                                      + py0 * (x1 * (t921 + t918 + t156 + t346 + t345)
 
1818
                                               + t883 + px1 * t875 + t247
 
1819
                                               + y1 * (x1 * (t917 + py1 * (t328 + t88 + (t915 + t718) * y2)
 
1820
                                                             + t153 + t55 + t339)
 
1821
                                                       + t899 + px1 * t914 + t322 + t143 + t245 + t887) + t772
 
1822
                                               + t244 + (t647 + t20 + t19 + t18) * t1) + x1 * t913 + t891
 
1823
                                      + py1 * (t202 + t780 + t196)
 
1824
                                      + y1 * (py1 * (t200 + t125 + t898) + t834 + py1 * x1 * t340)
 
1825
                                      + py1 * t338 * t1)
 
1826
                              + py0 * (x1 * (t908 + px1 * t904 + t247 + t810 + t244)
 
1827
                                       + t903 + py1 * t29 * t179 + t722
 
1828
                                       + y1 * (x1 * (px1 * t902 + t899 + t200 + t125 + t898)
 
1829
                                               + t897 + py1 * (x2 * t894 + x2 * (t86 + t628) * y2) + t267
 
1830
                                               + t266) + t237 + t721 + (t593 + t389 + t95 + t747) * t1)
 
1831
                              + x1 * (py1 * (t893 + t725 + t892) + t891) + t890 + t889 + py1 * t29 * t274
 
1832
                              + y1 * (x1 * (t888 + py1 * (t322 + t143 + t245 + t887)) + py1 * t268)
 
1833
                              + (t885 + t884) * t1)
 
1834
                      + py0 * (y1 * (x1 * (t883 + px1 * t878 + t202 + t780 + t196)
 
1835
                                     + t641 + t731 + t189 + t730)
 
1836
                               + x1 * (t639 + t230 + t229) + py1 * t29 * t169
 
1837
                               + t1 * (t840 + px1 * t850 + t13 + t752 + t825))
 
1838
                      + px0 * (y1 * (x1 * (py1 * t875 + t872 + t871 + t870)
 
1839
                                     + py1 * t869 + t866 + t865 + t864)
 
1840
                               + py1 * x1 * t499 + py1 * t29 * t174
 
1841
                               + t1 * (py1 * t838 + t829 + t828 + t827))
 
1842
                      + y1 * (x1 * (t863 + py1 * (t247 + t772 + t244)) + t859 + py1 * t723)
 
1843
                      + py0 * t161 * t183 + py1 * t29 * t160 + py1 * x1 * t152
 
1844
                      + t14 * (px0 * (x1 * (py1 * t855 + t853 + t852 + t851)
 
1845
                                      + py1 * t850 + py1 * t29 * t109 + t820 + t819 + t818
 
1846
                                      + (py1 * t592 + t826 + py1 * x1 * t646 + t847 + t846) * y1)
 
1847
                               + py0 * (x1 * (t845 + px1 * t844 + t54 + t282 + t841)
 
1848
                                        + t840 + px1 * t838 + py1 * t29 * t76 + t46 + t659 + t816
 
1849
                                        + (t421 + t431 + x1 * (t836 + t26 + t94 + t96) + t27 + t661) * y1)
 
1850
                               + x1 * (py1 * (t325 + t125 + t835) + t834) + t830
 
1851
                               + py1 * (t13 + t752 + t825) + py1 * t29 * t36
 
1852
                               + (t824 + py1 * t748 + py1 * x1 * t21) * y1)
 
1853
                      + t1 * (t821 + py1 * (t46 + t659 + t816)))
 
1854
        ;
 
1855
 
 
1856
    trans->m[1][1] = (t161 * (px0 * (x1 * (py1 * t603 + t911 + t910 + t909)
 
1857
                                     + py1 * t980 + t872
 
1858
                                     + y1 * (py1 * t944 + py1 * x1 * t109 + t833 + t832 + t831) + t871
 
1859
                                     + t870 + py1 * t646 * t1)
 
1860
                              + py0 * (x1 * (t514 + t157 + t364 + t363)
 
1861
                                       + px1 * t974 + t202
 
1862
                                       + y1 * (x1 * (t917 + t153 + t55 + t339)
 
1863
                                               + t943 + px1 * t955 + t143 + t124 + t938) + t978 + t725
 
1864
                                       + t196 + py1 * t66 * t1) + x1 * (t930 + py1 * t775) + t863
 
1865
                              + py1 * (t773 + t283 + t771)
 
1866
                              + y1 * (py1 * (t58 + t144 + t942) + t888 + py1 * x1 * t36)
 
1867
                              + py1 * t21 * t1)
 
1868
                      + x0 * (py0 * (t29 * (t604 + t149 + t346 + t774)
 
1869
                                     + t984 + t190
 
1870
                                     + y1 * (x1 * (px1 * t972 + t970 + t313 + t144 + t960)
 
1871
                                             + px1 * t958 + t950 + t787 + t12 + t939)
 
1872
                                     + x1 * (t607 + t983 + py1 * y2 * (x2 * (t328 + t88) + t293 + t327))
 
1873
                                     + t189 + t188 + py1 * t650 * t1)
 
1874
                              + px0 * (t29 * t925 + t982 + t866
 
1875
                                       + y1
 
1876
                                       * (x1 * (py1 * t969 + t853 + t852 + t851)
 
1877
                                          + py1 * t947 + t829 + t828 + t827)
 
1878
                                       + py1 * x1 * y2 * t334 + t865 + t864 + py1 * t651 * t1)
 
1879
                              + t29 * t913 + t859 + t981
 
1880
                              + y1 * (x1 * (t834 + py1 * (t633 + t54 + t272 + t968))
 
1881
                                      + py1 * (t946 + t46 + t945) + t821) + py1 * x1 * y2 * t606
 
1882
                              + py1 * t737 * t1)
 
1883
                      + py0 * (x1 * (t573 + t238 + t237 + t236)
 
1884
                               + t29 * (px1 * t980 + t773 + t978 + t283 + t771)
 
1885
                               + y1 * (x1 * (t965 + px1 * t971 + t765 + t659 + t959) + t550 + t670 + t669)
 
1886
                               + py1 * t537 * t1)
 
1887
                      + px0 * (x1 * (py1 * t977 + t858 + t857 + t856)
 
1888
                               + t29 * (py1 * t974 + t862 + t861 + t860)
 
1889
                               + y1 * (x1 * (py1 * t966 + t820 + t819 + t818) + t973) + py1 * t539 * t1)
 
1890
                      + x1 * (t890 + py1 * t191) + t29 * (t891 + py1 * (t202 + t725 + t196))
 
1891
                      + y0 * (x0 * (px0 * (x1 * (py1 * t972 + t833 + t832 + t831)
 
1892
                                           + py1 * t971 + py1 * t29 * t916 + t829 + t828 + t827
 
1893
                                           + (t926 + t817 + py1 * x1 * t693 + t823 + t822) * y1)
 
1894
                                    + py0 * (x1 * (t970 + px1 * t969 + t633 + t54 + t272 + t968)
 
1895
                                             + px1 * t966 + t965 + t29 * (t163 + t35 + t34 + t33) + t753 + t12
 
1896
                                             + t935
 
1897
                                             + (t962 + t941 + t43
 
1898
                                                + x1
 
1899
                                                * (t836 + py1 * (t86 + t628 + (t961 + t291) * x2)
 
1900
                                                   + t26 + t94 + t96) + t95 + t93)
 
1901
                                             * y1) + x1 * (py1 * (t313 + t144 + t960) + t888) + t821
 
1902
                                    + py1 * (t765 + t659 + t959) + py1 * t29 * t340
 
1903
                                    + (t885 + t884 + py1 * x1 * t338) * y1)
 
1904
                              + px0 * (x1 * (py1 * t958 + t820 + t819 + t818)
 
1905
                                       + t956 + t29 * (py1 * t955 + t853 + t852 + t851)
 
1906
                                       + (py1 * t415 + x1 * t953) * y1)
 
1907
                              + py0 * (x1 * (t950 + px1 * t947 + t946 + t46 + t945)
 
1908
                                       + t29 * (px1 * t944 + t943 + t58 + t144 + t942) + t406 + t392 + t391
 
1909
                                       + (x1 * (t421 + t941 + t431 + t27 + t661)
 
1910
                                          + t416 + py1 * (x2 * (t554 + t796) + t5 * t940) + t665 + t664)
 
1911
                                       * y1) + x1 * (py1 * (t787 + t12 + t939) + t830)
 
1912
                              + t29 * (t834 + py1 * (t143 + t124 + t938)) + t937
 
1913
                              + (x1 * t936 + py1 * t666) * y1)
 
1914
                      + y1 * (x1 * (py1 * (t753 + t12 + t935) + t830) + t934) + py1 * t653 * t1
 
1915
                      + py0 * t652 * t14)
 
1916
        ;
 
1917
 
 
1918
    trans->m[1][2] = (y0 * (x0 * (px0 * (x1 * (py1 * t1016 + t872 + t871 + t870)
 
1919
                                         + py1 * t494 + t29 * t932
 
1920
                                         + y1
 
1921
                                         * (py1 * t1025 + py1 * x1 * t599 + t820 + t819
 
1922
                                            + t818) + t953 * t1)
 
1923
                                  + py0 * (x1 * (t1015 + px1 * t1009 + t727 + t326 + t724)
 
1924
                                           + t29 * (t1028 + t148 + t364 + t744) + t491
 
1925
                                           + y1
 
1926
                                           * (x1
 
1927
                                              * (t600 + py1 * (x2 * t933 + t940 * y2) + t503
 
1928
                                                 + t800)
 
1929
                                              + px1 * t1020 + t1026 + t676 + t45 + t1017)
 
1930
                                           + (t1027 + t10 + t27 + t25) * t1)
 
1931
                                  + x1 * (py1 * (t1014 + t283 + t1012) + t863) + t29 * t931
 
1932
                                  + py1 * t802
 
1933
                                  + y1 * (py1 * x1 * t801 + py1 * (t1024 + t13 + t1023) + t830)
 
1934
                                  + t936 * t1)
 
1935
                            + py0 * (t29 * (t1007 + px1 * t1002 + t727 + t201 + t724)
 
1936
                                     + x1 * (t1011 + t641 + t731 + t189 + t730)
 
1937
                                     + y1 * (x1 * (t1026 + px1 * t1025 + t1024 + t13 + t1023)
 
1938
                                             + t550
 
1939
                                             + py1
 
1940
                                             * (t5 * (t293 + t327) + x2 * (t1021 + t697) * y2)
 
1941
                                             + t670 + t669) + (t540 + t648 + t645) * t1)
 
1942
                            + px0 * (x1 * (t982 + t866 + t865 + t864)
 
1943
                                     + t29 * (py1 * t1006 + t872 + t871 + t870)
 
1944
                                     + y1 * (x1 * (py1 * t1020 + t829 + t828 + t827) + t956)
 
1945
                                     + py1 * t415 * t1) + x1 * (t859 + t981)
 
1946
                            + t29 * (t863 + py1 * (t773 + t246 + t771))
 
1947
                            + y1 * (x1 * (py1 * (t676 + t45 + t1017) + t821) + t937)
 
1948
                            + py1 * t666 * t1)
 
1949
                      + x0 * (py0 * (y1 * (x1 * (px1 * t1016 + t1015 + t1014 + t283 + t1012)
 
1950
                                           + t1011 + t573 + t238 + t237 + t236)
 
1951
                                     + x1 * (t639
 
1952
                                             + py1 * (x2 * y2 * (t1010 + t300) + t142 * t894)
 
1953
                                             + t230 + t229) + t29 * (t175 + t159 + t155)
 
1954
                                     + t1 * (px1 * t1000 + t994 + t660 + t45 + t986))
 
1955
                              + px0 * (y1 * (x1 * (py1 * t1009 + t862 + t861 + t860)
 
1956
                                             + t927 + t858 + t857 + t856)
 
1957
                                       + py1 * t29 * t515 + py1 * x1 * t638
 
1958
                                       + t1 * (py1 * t993 + t820 + t819 + t818))
 
1959
                              + y1 * (x1 * (t891 + py1 * (t727 + t326 + t724)) + t890 + t889)
 
1960
                              + py1 * x1 * t374 + py1 * t29 * t373
 
1961
                              + t1 * (t830 + py1 * (t705 + t752 + t990)))
 
1962
                      + t161 * (py0 * (x1 * (t516 + t372 + t371)
 
1963
                                       + y1
 
1964
                                       * (x1 * (t921 + t156 + t346 + t345)
 
1965
                                          + t1007 + px1 * t1006 + t773 + t246 + t771) + t381
 
1966
                                       + t1 * (t989 + px1 * t997 + t54 + t245 + t985))
 
1967
                                + px0 * (y1
 
1968
                                         * (x1 * (py1 * t1004 + t924 + t923 + t922)
 
1969
                                            + py1 * t1002 + t862 + t861 + t860)
 
1970
                                         + py1 * t499 + py1 * x1 * t174
 
1971
                                         + t1 * (py1 * t988 + t853 + t852 + t851))
 
1972
                                + y1 * (x1 * (t912 + py1 * t746) + t891
 
1973
                                        + py1 * (t727 + t201 + t724))
 
1974
                                + py1 * x1 * t160 + py1 * t152
 
1975
                                + t1 * (t834 + py1 * (t143 + t272 + t987))) + py0 * t770
 
1976
                      + t14 * (px0 * (x1 * (py1 * t1000 + t829 + t828 + t827)
 
1977
                                      + t973 + t29 * (py1 * t997 + t833 + t832 + t831)
 
1978
                                      + (py1 * t539 + x1 * (py1 * t996 + t817 + t823 + t822))
 
1979
                                      * y1)
 
1980
                               + py0 * (x1 * (t994 + px1 * t993 + t705 + t752 + t990)
 
1981
                                        + t29 * (t989 + px1 * t988 + t143 + t272 + t987) + t398
 
1982
                                        + (t416 + x1 * (t962 + t43 + t95 + t93) + t665 + t664)
 
1983
                                        * y1) + x1 * (py1 * (t660 + t45 + t986) + t821)
 
1984
                               + t29 * (t888 + py1 * (t54 + t245 + t985)) + t934
 
1985
                               + (x1 * (t885 + py1 * t28) + py1 * t653) * y1));
 
1986
 
 
1987
    trans->m[2][0] = (x0 * (px0 * (y1 * (x1 * t617 + t586 + t877 + t1008 + t876)
 
1988
                                   + x1 * t515 + t637 + t1 * (t475 + t136 + t620 + t1036) + t636)
 
1989
                            + py0 * (y1 * (t579 + x1 * t627 + t251 + t811 + t248)
 
1990
                                     + x1 * t379 + t233 + t1 * (t455 + t87 + t298 + t1035) + t232)
 
1991
                            + x1 * (t516 + t383 + y2 * (t804 + t1013) + t142 * t595)
 
1992
                            + y1 * (px1 * (t1047 + t979 + t1046)
 
1993
                                    + x1 * (t921 + t368 + t157 + t156 + t1048)
 
1994
                                    + py1 * (t1044 + t728 + t1043) + t505 + t607 + t502) + t500 + t381
 
1995
                            + t1 * (px1 * (t611 + t135 + t1038) + py1 * (t330 + t128 + t1037) + t423
 
1996
                                    + t125 + t54 + t1034))
 
1997
                      + y0 * (x0 * (py0 * (x1 * t1050 + t533 + t226
 
1998
                                           + y1 * (t529 + t224 + x1 * t349 + t128 + t1040)
 
1999
                                           + t781 + t220 + t99 * t1)
 
2000
                                    + px0 * (t528 + x1 * t1049 + t874
 
2001
                                             + y1 * (t517 + x1 * t916 + t520 + t135 + t1041) + t1005
 
2002
                                             + t873 + t693 * t1)
 
2003
                                    + x1 * (t514 + t382 + t157 + t156 + t1048)
 
2004
                                    + px1 * (t877 + t1001 + t876) + py1 * (t251 + t776 + t248) + t608
 
2005
                                    + y1 * (x1 * (t917 + t358 + t56 + t323 + (t716 + t967) * y2)
 
2006
                                            + px1 * (t118 + t609 + t471 + t1045)
 
2007
                                            + py1 * (t204 + t327 + t249 + t1042) + t503 + t144 + t143
 
2008
                                            + t1039) + t504 + t605 + (t647 + t357) * t1)
 
2009
                              + px0 * (x1 * (t528 + t1047 + t979 + t1046)
 
2010
                                       + t643 + t29 * t182 + t571
 
2011
                                       + y1 * (x1 * (t475 + t118 + t609 + t471 + t1045) + t701 + t895)
 
2012
                                       + t570 + t569 + (t468 + t441 + t466 + t951) * t1)
 
2013
                              + py0 * (x1 * (t533 + t1044 + t728 + t1043)
 
2014
                                       + t642 + t29 * t179 + t566
 
2015
                                       + y1 * (x1 * (t455 + t204 + t327 + t249 + t1042) + t289 + t288)
 
2016
                                       + t214 + t565 + (t436 + t450 + t69 + t677) * t1)
 
2017
                              + x1 * (px1 * (t877 + t1008 + t876) + py1 * (t251 + t811 + t248) + t505
 
2018
                                      + t607 + t502) + t984 + t743
 
2019
                              + t29 * t274
 
2020
                              + y1 * (x1 * (px1 * (t520 + t135 + t1041) + py1 * (t224 + t128 + t1040)
 
2021
                                            + t503 + t144 + t143 + t1039)
 
2022
                                      + t897 + t303 + x2 * t634 + x2 * (t690 + t52) * y2)
 
2023
                              + (t1027 + t808) * t1)
 
2024
                      + py0 * (y1 * (x1 * (t579 + t251 + t776 + t248) + t581 + t735 + t240 + t734)
 
2025
                               + t29 * t169 + x1 * t209 + t1 * (t461 + t80 + t683 + t1032))
 
2026
                      + px0 * (y1 * (x1 * (t586 + t877 + t1001 + t876) + t591 + t976 + t588 + t975)
 
2027
                               + x1 * t499 + t29 * t174 + t1 * (t484 + t849 + t957 + t1033))
 
2028
                      + y1 * (x1 * (px1 * (t874 + t1005 + t873) + py1 * (t226 + t781 + t220) + t608
 
2029
                                    + t504 + t605)
 
2030
                              + t573 + t568) + t161 * t183 + x1 * (t639 + t235) + t29 * t160
 
2031
                      + t14 * (px0 * (x1 * (t517 + t611 + t135 + t1038)
 
2032
                                      + t615 + t29 * t109 + t560 + t443 + t1031
 
2033
                                      + (t399 + t613 + x1 * t646 + t418 + t995) * y1)
 
2034
                               + py0 * (x1 * (t529 + t330 + t128 + t1037)
 
2035
                                        + t625 + t29 * t76 + t112 + t438 + t1030
 
2036
                                        + (t623 + t393 + t102 + x1 * t66 + t754) * y1)
 
2037
                               + x1 * (px1 * (t136 + t620 + t1036) + py1 * (t87 + t298 + t1035) + t423
 
2038
                                       + t125 + t54 + t1034)
 
2039
                               + px1 * (t849 + t957 + t1033) + py1 * (t80 + t683 + t1032) + t434
 
2040
                               + t29 * t36 + t1029 + (t962 + t711 + x1 * (t836 + t100)) * y1)
 
2041
                      + t1 * (px1 * (t560 + t443 + t1031) + py1 * (t112 + t438 + t1030) + t434
 
2042
                              + t1029))
 
2043
        ;
 
2044
 
 
2045
    trans->m[2][1] = (t161 * (px0 * (x1 * (t616 + t180 + t512 + t1003)
 
2046
                                     + t586 + t510 + y1 * (t475 + t854 + x1 * t109 + t117 + t1062)
 
2047
                                     + t621 + t508 + t646 * t1)
 
2048
                              + py0 * (t579 + x1 * (t177 + t626 + t354 + t778) + t507
 
2049
                                       + y1 * (t455 + x1 * t76 + t91 + t203 + t1061) + t299
 
2050
                                       + t506 + t66 * t1) + x1 * (t921 + t368)
 
2051
                              + px1 * (t874 + t979 + t873) + py1 * (t226 + t728 + t220)
 
2052
                              + y1 * (x1 * (t917 + t358) + px1 * (t118 + t843 + t1063)
 
2053
                                      + py1 * (t204 + t127 + t1065) + t423 + t144
 
2054
                                      + t143 + t1060) + t504 + t21 * t1)
 
2055
                      + x0 * (py0 * (t29 * t1050 + t581 + t241
 
2056
                                     + y1
 
2057
                                     * (x1 * (t529 + t87 + t577 + t292 + t1070)
 
2058
                                        + t461 + t1057 + t80 + t1056) + x1 * y2 * t741 + t240
 
2059
                                     + t239 + t650 * t1)
 
2060
                              + px0 * (t591 + t29 * t1049 + t589
 
2061
                                       + y1 * (x1 * (t517 + t136 + t306 + t428 + t1071)
 
2062
                                               + t484 + t1059 + t849 + t1058) + x1 * y2 * t334 + t588
 
2063
                                       + t587 + t651 * t1) + t29 * (t1028 + t809) + t903 + t321
 
2064
                              + y1 * (x1 * (px1 * (t901 + t117 + t1072) + py1 * (t317 + t203 + t1073)
 
2065
                                            + t503 + t125 + t54 + t1069)
 
2066
                                      + px1 * (t1019 + t837 + t1064) + py1 * (t789 + t111 + t1066) + t267
 
2067
                                      + t433 + t1055)
 
2068
                              + x1 * (y2 * (x2 * (t56 + t323) + t322 + t273) + t983 + t742) + t737 * t1)
 
2069
                      + py0 * (x1 * (t642 + t215 + t214 + t213) + t29 * (t533 + t226 + t728 + t220)
 
2070
                               + y1
 
2071
                               * (x1 * (t625 + t758 + t111 + t1067)
 
2072
                                  + t396 + t395) + t537 * t1)
 
2073
                      + px0 * (x1 * (t643 + t868 + t570 + t867) + t29 * (t528 + t874 + t979 + t873)
 
2074
                               + y1
 
2075
                               * (x1 * (t615 + t999 + t837 + t1068)
 
2076
                                  + t404 + t402) + t539 * t1)
 
2077
                      + x1 * (t641 + t243)
 
2078
                      + t29 * (px1 * (t510 + t621 + t508) + py1 * (t507 + t299 + t506) + t504)
 
2079
                      + y0 * (x0 * (py0 * (x1 * (t455 + t317 + t203 + t1073)
 
2080
                                           + t461 + t29 * t349 + t767 + t683 + t1052
 
2081
                                           + (t436 + t450 + x1 * t99 + t69 + t677) * y1)
 
2082
                                    + px0 * (x1 * (t475 + t901 + t117 + t1072)
 
2083
                                             + t484 + t29 * t916 + t992 + t957 + t1053
 
2084
                                             + (t468 + t441 + t466 + x1 * t693 + t951) * y1)
 
2085
                                    + x1 * (px1 * (t136 + t306 + t428 + t1071)
 
2086
                                            + py1 * (t87 + t577 + t292 + t1070) + t503 + t125 + t54 + t1069)
 
2087
                                    + px1 * (t999 + t837 + t1068) + py1 * (t758 + t111 + t1067)
 
2088
                                    + t29 * (t163 + t162) + t434 + t553 + t1051
 
2089
                                    + (t421 + t104 + t431 + t43
 
2090
                                       + x1 * (t836 + t100 + t690 + t52 + (t270 + t886) * x2)
 
2091
                                       + t1054)
 
2092
                                    * y1)
 
2093
                              + py0 * (x1 * (t625 + t789 + t111 + t1066)
 
2094
                                       + t29 * (t529 + t204 + t127 + t1065) + t544 + t543
 
2095
                                       + (x1 * t624 + t410 + t409) * y1)
 
2096
                              + px0 * (x1 * (t615 + t1019 + t837 + t1064)
 
2097
                                       + t29 * (t517 + t118 + t843 + t1063) + t548 + t547
 
2098
                                       + (x1 * t614 + t414 + t413) * y1)
 
2099
                              + t29 * (px1 * (t854 + t117 + t1062) + py1 * (t91 + t203 + t1061) + t423
 
2100
                                       + t144 + t143 + t1060)
 
2101
                              + x1 * (px1 * (t1059 + t849 + t1058) + py1 * (t1057 + t80 + t1056) + t267
 
2102
                                      + t433 + t1055) + t406 + t398
 
2103
                              + (t416 + x1 * (t962 + t711 + t431 + t43 + t1054) + t412
 
2104
                                 + x2 * (t793 + t1022) + t5 * t594)
 
2105
                              * y1)
 
2106
                      + y1 * (x1 * (px1 * (t992 + t957 + t1053) + py1 * (t767 + t683 + t1052) + t434
 
2107
                                    + t553 + t1051)
 
2108
                              + t550 + t546) + t653 * t1 + t652 * t14)
 
2109
        ;
 
2110
    trans->m[2][2] = t644;
 
2111
}
 
2112
 
 
2113
static void
 
2114
print_trans (const char *header, struct pixman_f_transform *trans)
 
2115
{
 
2116
    int i, j;
 
2117
    double max;
 
2118
 
 
2119
    max = 0;
 
2120
 
 
2121
    printf ("%s\n", header);
 
2122
 
 
2123
    for (i = 0; i < 3; ++i)
 
2124
    {
 
2125
        for (j = 0; j < 3; ++j)
 
2126
        {
 
2127
            double a = fabs (trans->m[i][j]);
 
2128
 
 
2129
            if (a > max)
 
2130
                max = a;
 
2131
        }
 
2132
    }
 
2133
 
 
2134
    if (max == 0.0)
 
2135
        max = 1.0;
 
2136
 
 
2137
    for (i = 0; i < 3; ++i)
 
2138
    {
 
2139
        printf ("{ ");
 
2140
        for (j = 0; j < 3; ++j)
 
2141
        {
 
2142
            printf ("D2F (%.5f)%s", 16384 * (trans->m[i][j] / max), j == 2 ? "" : ", ");
 
2143
        }
 
2144
 
 
2145
        printf ("},\n");
 
2146
    }
 
2147
}
 
2148
 
 
2149
int
 
2150
main ()
 
2151
{
 
2152
    struct pixman_f_transform t;
 
2153
 
 
2154
#if 0
 
2155
    quad_to_quad (75, 200,
 
2156
                  325, 200,
 
2157
                  450, 335,
 
2158
                  -50, 335,
 
2159
 
 
2160
                  0, 0,
 
2161
                  400, 0,
 
2162
                  400, 400,
 
2163
                  0, 400,
 
2164
 
 
2165
                  &t);
 
2166
#endif
 
2167
    quad_to_quad (
 
2168
        1, 0,
 
2169
        1, 2,
 
2170
        2, 2,
 
2171
        2, 0,
 
2172
 
 
2173
        1, 0,
 
2174
        1, 112,
 
2175
        2, 2,
 
2176
        2, 0,
 
2177
 
 
2178
        &t);
 
2179
 
 
2180
    print_trans ("0->0", &t);
 
2181
 
 
2182
    return 0;
 
2183
}