qemu/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 }