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
quad_to_quad(double x0,double y0,double x1,double y1,double x2,double y2,double x3,double y3,double px0,double py0,double px1,double py1,double px2,double py2,double px3,double py3,struct pixman_f_transform * trans)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
print_trans(const char * header,struct pixman_f_transform * trans)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
main()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 }
2184