Søren Sandmann sandmann at cs.au.dk
Sat Mar 24 13:24:59 PDT 2012

```From: Søren Sandmann Pedersen <ssp at redhat.com>

This program can compute the projective transformation that transforms
one quadrilateral into another. The code is basically maxima[1] output
translated into C.

[1] http://maxima.sourceforge.net/
---
1 files changed, 2183 insertions(+), 0 deletions(-)

new file mode 100644
index 0000000..a8c05b1
--- /dev/null
@@ -0,0 +1,2183 @@
+#include <math.h>
+#include <stdio.h>
+#include <pixman.h>
+
+/* This code is basically the output of Maxima translated into C.
+ *
+ * See http://maxima.sourceforge.net/
+ */
+static void
+	      double x1, double y1,
+	      double x2, double y2,
+	      double x3, double y3,
+
+	      double px0, double py0,
+	      double px1, double py1,
+	      double px2, double py2,
+	      double px3, double py3,
+
+	      struct pixman_f_transform *trans)
+{
+    double
+	t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15, t16, t17, t18,
+	t19, t20, t21, t22, t23, t24, t25, t26, t27, t28, t29, t30, t31, t32, t33, t34,
+	t35, t36, t37, t38, t39, t40, t41, t42, t43, t44, t45, t46, t47, t48, t49, t50,
+	t51, t52, t53, t54, t55, t56, t57, t58, t59, t60, t61, t62, t63, t64, t65, t66,
+	t67, t68, t69, t70, t71, t72, t73, t74, t75, t76, t77, t78, t79, t80, t81, t82,
+	t83, t84, t85, t86, t87, t88, t89, t90, t91, t92, t93, t94, t95, t96, t97, t98,
+	t99, t100, t101, t102, t103, t104, t105, t106, t107, t108, t109, t110, t111,
+	t112, t113, t114, t115, t116, t117, t118, t119, t120, t121, t122, t123,
+	t124, t125, t126, t127, t128, t129, t130, t131, t132, t133, t134, t135,
+	t136, t137, t138, t139, t140, t141, t142, t143, t144, t145, t146, t147,
+	t148, t149, t150, t151, t152, t153, t154, t155, t156, t157, t158, t159,
+	t160, t161, t162, t163, t164, t165, t166, t167, t168, t169, t170, t171,
+	t172, t173, t174, t175, t176, t177, t178, t179, t180, t181, t182, t183,
+	t184, t185, t186, t187, t188, t189, t190, t191, t192, t193, t194, t195,
+	t196, t197, t198, t199, t200, t201, t202, t203, t204, t205, t206, t207,
+	t208, t209, t210, t211, t212, t213, t214, t215, t216, t217, t218, t219,
+	t220, t221, t222, t223, t224, t225, t226, t227, t228, t229, t230, t231,
+	t232, t233, t234, t235, t236, t237, t238, t239, t240, t241, t242, t243,
+	t244, t245, t246, t247, t248, t249, t250, t251, t252, t253, t254, t255,
+	t256, t257, t258, t259, t260, t261, t262, t263, t264, t265, t266, t267,
+	t268, t269, t270, t271, t272, t273, t274, t275, t276, t277, t278, t279,
+	t280, t281, t282, t283, t284, t285, t286, t287, t288, t289, t290, t291,
+	t292, t293, t294, t295, t296, t297, t298, t299, t300, t301, t302, t303,
+	t304, t305, t306, t307, t308, t309, t310, t311, t312, t313, t314, t315,
+	t316, t317, t318, t319, t320, t321, t322, t323, t324, t325, t326, t327,
+	t328, t329, t330, t331, t332, t333, t334, t335, t336, t337, t338, t339,
+	t340, t341, t342, t343, t344, t345, t346, t347, t348, t349, t350, t351,
+	t352, t353, t354, t355, t356, t357, t358, t359, t360, t361, t362, t363,
+	t364, t365, t366, t367, t368, t369, t370, t371, t372, t373, t374, t375,
+	t376, t377, t378, t379, t380, t381, t382, t383, t384, t385, t386, t387,
+	t388, t389, t390, t391, t392, t393, t394, t395, t396, t397, t398, t399,
+	t400, t401, t402, t403, t404, t405, t406, t407, t408, t409, t410, t411,
+	t412, t413, t414, t415, t416, t417, t418, t419, t420, t421, t422, t423,
+	t424, t425, t426, t427, t428, t429, t430, t431, t432, t433, t434, t435,
+	t436, t437, t438, t439, t440, t441, t442, t443, t444, t445, t446, t447,
+	t448, t449, t450, t451, t452, t453, t454, t455, t456, t457, t458, t459,
+	t460, t461, t462, t463, t464, t465, t466, t467, t468, t469, t470, t471,
+	t472, t473, t474, t475, t476, t477, t478, t479, t480, t481, t482, t483,
+	t484, t485, t486, t487, t488, t489, t490, t491, t492, t493, t494, t495,
+	t496, t497, t498, t499, t500, t501, t502, t503, t504, t505, t506, t507,
+	t508, t509, t510, t511, t512, t513, t514, t515, t516, t517, t518, t519,
+	t520, t521, t522, t523, t524, t525, t526, t527, t528, t529, t530, t531,
+	t532, t533, t534, t535, t536, t537, t538, t539, t540, t541, t542, t543,
+	t544, t545, t546, t547, t548, t549, t550, t551, t552, t553, t554, t555,
+	t556, t557, t558, t559, t560, t561, t562, t563, t564, t565, t566, t567,
+	t568, t569, t570, t571, t572, t573, t574, t575, t576, t577, t578, t579,
+	t580, t581, t582, t583, t584, t585, t586, t587, t588, t589, t590, t591,
+	t592, t593, t594, t595, t596, t597, t598, t599, t600, t601, t602, t603,
+	t604, t605, t606, t607, t608, t609, t610, t611, t612, t613, t614, t615,
+	t616, t617, t618, t619, t620, t621, t622, t623, t624, t625, t626, t627,
+	t628, t629, t630, t631, t632, t633, t634, t635, t636, t637, t638, t639,
+	t640, t641, t642, t643, t644, t645, t646, t647, t648, t649, t650, t651,
+	t652, t653, t654, t655, t656, t657, t658, t659, t660, t661, t662, t663,
+	t664, t665, t666, t667, t668, t669, t670, t671, t672, t673, t674, t675,
+	t676, t677, t678, t679, t680, t681, t682, t683, t684, t685, t686, t687,
+	t688, t689, t690, t691, t692, t693, t694, t695, t696, t697, t698, t699,
+	t700, t701, t702, t703, t704, t705, t706, t707, t708, t709, t710, t711,
+	t712, t713, t714, t715, t716, t717, t718, t719, t720, t721, t722, t723,
+	t724, t725, t726, t727, t728, t729, t730, t731, t732, t733, t734, t735,
+	t736, t737, t738, t739, t740, t741, t742, t743, t744, t745, t746, t747,
+	t748, t749, t750, t751, t752, t753, t754, t755, t756, t757, t758, t759,
+	t760, t761, t762, t763, t764, t765, t766, t767, t768, t769, t770, t771,
+	t772, t773, t774, t775, t776, t777, t778, t779, t780, t781, t782, t783,
+	t784, t785, t786, t787, t788, t789, t790, t791, t792, t793, t794, t795,
+	t796, t797, t798, t799, t800, t801, t802, t803, t804, t805, t806, t807,
+	t808, t809, t810, t811, t812, t813, t814, t815, t816, t817, t818, t819,
+	t820, t821, t822, t823, t824, t825, t826, t827, t828, t829, t830, t831,
+	t832, t833, t834, t835, t836, t837, t838, t839, t840, t841, t842, t843,
+	t844, t845, t846, t847, t848, t849, t850, t851, t852, t853, t854, t855,
+	t856, t857, t858, t859, t860, t861, t862, t863, t864, t865, t866, t867,
+	t868, t869, t870, t871, t872, t873, t874, t875, t876, t877, t878, t879,
+	t880, t881, t882, t883, t884, t885, t886, t887, t888, t889, t890, t891,
+	t892, t893, t894, t895, t896, t897, t898, t899, t900, t901, t902, t903,
+	t904, t905, t906, t907, t908, t909, t910, t911, t912, t913, t914, t915,
+	t916, t917, t918, t919, t920, t921, t922, t923, t924, t925, t926, t927,
+	t928, t929, t930, t931, t932, t933, t934, t935, t936, t937, t938, t939,
+	t940, t941, t942, t943, t944, t945, t946, t947, t948, t949, t950, t951,
+	t952, t953, t954, t955, t956, t957, t958, t959, t960, t961, t962, t963,
+	t964, t965, t966, t967, t968, t969, t970, t971, t972, t973, t974, t975,
+	t976, t977, t978, t979, t980, t981, t982, t983, t984, t985, t986, t987,
+	t988, t989, t990, t991, t992, t993, t994, t995, t996, t997, t998, t999,
+	t1000, t1001, t1002, t1003, t1004, t1005, t1006, t1007, t1008, t1009,
+	t1010, t1011, t1012, t1013, t1014, t1015, t1016, t1017, t1018, t1019,
+	t1020, t1021, t1022, t1023, t1024, t1025, t1026, t1027, t1028, t1029,
+	t1030, t1031, t1032, t1033, t1034, t1035, t1036, t1037, t1038, t1039,
+	t1040, t1041, t1042, t1043, t1044, t1045, t1046, t1047, t1048, t1049,
+	t1050, t1051, t1052, t1053, t1054, t1055, t1056, t1057, t1058, t1059,
+	t1060, t1061, t1062, t1063, t1064, t1065, t1066, t1067, t1068, t1069,
+	t1070, t1071, t1072, t1073;
+
+    t1 = y1 * y1;
+    t2 = x3 * x3;
+    t3 = px2 * px3 * t2;
+    t4 = (t3 - px2 * px3 * x2 * x3) * y2;
+    t5 = x2 * x2;
+    t6 = px2 * px3 * t5 * y3;
+
+    t7 = - px2 * px3 * x2 * x3 * y3;
+    t8 = py1 * (t7 + t6 + t4);
+    t9 = px3 * py2 * x2 * x3;
+
+    t10 = - px3 * py2 * t2;
+    t11 = (t10 + t9) * y2;
+    t12 = - px2 * py3 * t5 * y3;
+
+    t13 = px2 * py3 * x2 * x3 * y3;
+    t14 = y0 * y0;
+    t15 = - px3 * py2;
+    t16 = px2 * py3;
+
+    t17 = t16 + t15;
+    t18 = t17 * x2;
+    t19 = px3 * py2 * x3;
+    t20 = - px2 * py3 * x3;
+
+    t21 = t20 + t19 + t18;
+    t22 = px2 * px3 * t5;
+    t23 = - 2 * px2 * px3 * x2 * x3;
+
+    t24 = py1 * (t3 + t23 + t22);
+    t25 = - px2 * py3 * t5;
+    t26 = px2 * py3 * x3;
+
+    t27 = x2 * (t26 + t19);
+    t28 = t10 + t27 + t25;
+    t29 = x1 * x1;
+    t30 = px3 * py2;
+
+    t31 = - px2 * py3;
+    t32 = t31 + t30;
+    t33 = t32 * y2;
+    t34 = - px3 * py2 * y3;
+
+    t35 = px2 * py3 * y3;
+    t36 = t35 + t34 + t33;
+    t37 = - px2 * px3 * t2;
+
+    t38 = (t37 + px2 * px3 * x2 * x3) * y2;
+    t39 = - px2 * px3 * t5 * y3;
+
+    t40 = px2 * px3 * x2 * x3 * y3;
+    t41 = py1 * (t40 + t39 + t38);
+    t42 = - px2 * py3 * x2 * x3;
+
+    t43 = px3 * py2 * t2;
+    t44 = (t43 + t42) * y2;
+    t45 = px2 * py3 * t5 * y3;
+
+    t46 = - px3 * py2 * x2 * x3 * y3;
+    t47 = (px2 * px3 * x3 - px2 * px3 * x2) * y2;
+
+    t48 = px2 * px3 * x2 * y3;
+    t49 = - px2 * px3 * x3 * y3;
+    t50 = py1 * (t49 + t48 + t47);
+
+    t51 = px2 * py3 * x2;
+    t52 = - 2 * px3 * py2 * x3;
+    t53 = (t26 + t52 + t51) * y2;
+
+    t54 = px3 * py2 * x3 * y3;
+    t55 = px3 * py2 * y3;
+    t56 = - 2 * px2 * py3 * y3;
+    t57 = t56 + t55;
+
+    t58 = x2 * t57;
+    t59 = - px2 * px3 * t5;
+    t60 = 2 * px2 * px3 * x2 * x3;
+    t61 = - px2;
+
+    t62 = px3 + t61;
+    t63 = t62 * x2;
+    t64 = px2 * x3;
+    t65 = - px3 * x3;
+    t66 = t65 + t64 + t63;
+
+    t67 = px2 * t5;
+    t68 = - px2 * x3;
+    t69 = x2 * (t65 + t68);
+    t70 = px3 * t2;
+
+    t71 = t70 + t69 + t67;
+    t72 = - px3;
+    t73 = t72 + px2;
+    t74 = - px2 * y3;
+    t75 = px3 * y3;
+
+    t76 = t75 + t74 + t73 * y2;
+    t77 = px2 * x2 * x3;
+    t78 = - px3 * t2;
+    t79 = - px2 * t5 * y3;
+
+    t80 = px3 * x2 * x3 * y3;
+    t81 = t80 + t79 + (t78 + t77) * y2;
+
+    t82 = (px2 * px3 * x2 - px2 * px3 * x3) * y2;
+    t83 = - px2 * px3 * x2 * y3;
+
+    t84 = px2 * px3 * x3 * y3;
+    t85 = - px2 * x2;
+    t86 = 2 * px3 * x3;
+    t87 = - px3 * x3 * y3;
+
+    t88 = 2 * px2 * y3;
+    t89 = - px3 * y3;
+    t90 = t89 + t88;
+    t91 = x2 * t90;
+
+    t92 = t91 + t87 + (t86 + t68 + t85) * y2;
+    t93 = px2 * py3 * t5;
+    t94 = - px3 * py2 * x3;
+
+    t95 = x2 * (t20 + t94);
+    t96 = t32 * x2;
+    t97 = t73 * x2;
+    t98 = px3 * x3;
+
+    t99 = t98 + t68 + t97;
+    t100 = py1 * t99;
+    t101 = - px2 * t5;
+    t102 = x2 * (t98 + t64);
+
+    t103 = t78 + t102 + t101;
+    t104 = py1 * t103;
+    t105 = - py2;
+    t106 = py3 + t105;
+
+    t107 = py2 * y3;
+    t108 = - py3 * y3;
+    t109 = t108 + t107 + t106 * y2;
+    t110 = - px3 * x2 * x3;
+
+    t111 = px2 * t5 * y3;
+    t112 = - px2 * x2 * x3 * y3;
+    t113 = t112 + t111 + (t70 + t110) * y2;
+
+    t114 = - py2 * x3;
+    t115 = py3 * x3;
+    t116 = t115 + t114;
+    t117 = py2 * x3 * y3;
+
+    t118 = - py3 * x3 * y3;
+    t119 = t118 + t117;
+    t120 = x2 * t119;
+
+    t121 = px1 * (t120 + x2 * t116 * y2);
+    t122 = - px3 * py2 * x2;
+    t123 = (t19 + t122) * y2;
+
+    t124 = px2 * py3 * x2 * y3;
+    t125 = - px2 * py3 * x3 * y3;
+    t126 = px3 * x2;
+
+    t127 = - px2 * x2 * y3;
+    t128 = px2 * x3 * y3;
+    t129 = t128 + t127 + (t65 + t126) * y2;
+
+    t130 = - py3;
+    t131 = t130 + py2;
+    t132 = t131 * x2;
+    t133 = py2 * x3;
+    t134 = - py3 * x3;
+
+    t135 = - py2 * x3 * y3;
+    t136 = py3 * x3 * y3;
+    t137 = - py2 * y3;
+    t138 = py3 * y3;
+
+    t139 = t138 + t137;
+    t140 = x2 * t139;
+
+    t141 = px1 * (t140 + t136 + t135 + (t134 + t133 + t132) * y2);
+    t142 = y2 * y2;
+
+    t143 = - px3 * py2 * x3 * y3;
+    t144 = px2 * py3 * x3 * y3;
+    t145 = t144 + t143;
+
+    t146 = t142 * t145;
+    t147 = y3 * y3;
+    t148 = px3 * py2 * t147;
+    t149 = - px2 * py3 * t147;
+
+    t150 = t149 + t148;
+    t151 = x2 * y2 * t150;
+    t152 = t151 + t146;
+    t153 = - px2 * py3 * y3;
+
+    t154 = t153 + t55;
+    t155 = t142 * t154;
+    t156 = - px3 * py2 * t147;
+
+    t157 = px2 * py3 * t147;
+    t158 = t157 + t156;
+    t159 = y2 * t158;
+    t160 = t159 + t155;
+
+    t161 = x0 * x0;
+    t162 = py1 * t76;
+    t163 = px1 * t109;
+    t164 = px2 * y3;
+    t165 = t89 + t164;
+
+    t166 = - px2 * t147;
+    t167 = px3 * t147;
+    t168 = t167 + t166;
+
+    t169 = y2 * t168 + t142 * t165;
+    t170 = py1 * t169;
+    t171 = py2 * t147;
+
+    t172 = - py3 * t147;
+    t173 = t172 + t171;
+    t174 = y2 * t173 + t142 * t139;
+
+    t175 = px1 * t174;
+    t176 = t17 * t142;
+    t177 = px2 * t147;
+    t178 = - px3 * t147;
+
+    t179 = t178 + t177 + t62 * t142;
+    t180 = - py2 * t147;
+    t181 = py3 * t147;
+
+    t182 = t181 + t180 + t131 * t142;
+
+    t183 = y1 * (px1 * t182 + py1 * t179 + t149 + t148 + t176)
+	+ t175 + t170 + t159 + t1 * (t163 + t162 + t35 + t34 + t33) + t155;
+
+    t184 = - px2 * px3 * t2 * t142;
+    t185 = 2 * px2 * px3 * x2 * x3 * y2 * y3;
+
+    t186 = - px2 * px3 * t5 * t147;
+    t187 = py1 * (t186 + t185 + t184);
+
+    t188 = px3 * py2 * t2 * t142;
+    t189 = x2 * y2 * (t125 + t143);
+    t190 = px2 * py3 * t5 * t147;
+
+    t191 = t190 + t189 + t188;
+    t192 = px2 * px3 * x3 * t142;
+    t193 = y2 * (t49 + t83);
+
+    t194 = px2 * px3 * x2 * t147;
+    t195 = py1 * (t194 + t193 + t192);
+
+    t196 = - px3 * py2 * x3 * t142;
+    t197 = 2 * px3 * py2 * x3 * y3;
+    t198 = 2 * px2 * py3 * y3;
+
+    t199 = t198 + t34;
+    t200 = x2 * t199;
+    t201 = y2 * (t200 + t125 + t197);
+
+    t202 = - px2 * py3 * x2 * t147;
+    t203 = - px2 * x3 * y3;
+    t204 = px3 * x3 * y3;
+
+    t205 = t204 + t203;
+    t206 = t142 * t205;
+    t207 = t178 + t177;
+    t208 = x2 * y2 * t207;
+
+    t209 = t208 + t206;
+    t210 = px2 * px3 * t2 * t142;
+    t211 = - 2 * px2 * px3 * x2 * x3 * y2 * y3;
+
+    t212 = px2 * px3 * t5 * t147;
+    t213 = - px3 * t2 * t142;
+    t214 = x2 * y2 * (t204 + t128);
+
+    t215 = - px2 * t5 * t147;
+    t216 = t215 + t214 + t213;
+    t217 = - px2 * px3 * x3 * t142;
+
+    t218 = y2 * (t84 + t48);
+    t219 = - px2 * px3 * x2 * t147;
+    t220 = px3 * x3 * t142;
+
+    t221 = - 2 * px3 * x3 * y3;
+    t222 = - 2 * px2 * y3;
+    t223 = t75 + t222;
+    t224 = x2 * t223;
+
+    t225 = y2 * (t224 + t221 + t128);
+    t226 = px2 * x2 * t147;
+    t227 = t226 + t225 + t220;
+
+    t228 = t125 + t54;
+    t229 = t142 * t228;
+    t230 = x2 * y2 * t158;
+    t231 = t87 + t128;
+
+    t232 = t142 * t231;
+    t233 = x2 * y2 * t168;
+    t234 = t233 + t232;
+    t235 = py1 * t234;
+
+    t236 = - px3 * py2 * t2 * t142;
+    t237 = x2 * y2 * (t144 + t54);
+
+    t238 = - px2 * py3 * t5 * t147;
+    t239 = px3 * t2 * t142;
+    t240 = x2 * y2 * (t87 + t203);
+
+    t241 = px2 * t5 * t147;
+    t242 = t241 + t240 + t239;
+    t243 = py1 * t242;
+
+    t244 = px2 * py3 * x3 * t142;
+    t245 = - px2 * py3 * x2 * y3;
+    t246 = y2 * (t143 + t245);
+
+    t247 = px3 * py2 * x2 * t147;
+    t248 = - px2 * x3 * t142;
+    t249 = px2 * x2 * y3;
+
+    t250 = y2 * (t204 + t249);
+    t251 = - px3 * x2 * t147;
+    t252 = t251 + t250 + t248;
+
+    t253 = t134 + t133;
+    t254 = t253 * t142;
+    t255 = t108 + t107;
+    t256 = x2 * t255;
+
+    t257 = t256 + t136 + t135;
+    t258 = y2 * t257;
+    t259 = t181 + t180;
+    t260 = x2 * t259;
+
+    t261 = px1 * (t260 + t258 + t254);
+    t262 = py1 * (t37 + t60 + t59);
+
+    t263 = t43 + t95 + t93;
+    t264 = px1 * t263;
+    t265 = t26 + t94;
+    t266 = x2 * t265 * y2;
+
+    t267 = x2 * t228;
+    t268 = t267 + t266;
+    t269 = py1 * (t84 + t83 + t82);
+
+    t270 = - 2 * px2 * py3;
+    t271 = (t26 + (t270 + t30) * x2) * y2;
+    t272 = px3 * py2 * x2 * y3;
+
+    t273 = - 2 * px3 * py2 * x3 * y3;
+    t274 = t149 + t148 + t176;
+
+    t275 = py1 * (t212 + t211 + t210);
+    t276 = t238 + t237 + t236;
+    t277 = px1 * t276;
+
+    t278 = py1 * (t219 + t218 + t217);
+    t279 = 2 * px3 * py2 * x3;
+    t280 = t20 + t279;
+
+    t281 = t280 * t142;
+    t282 = - px3 * py2 * x2 * y3;
+    t283 = y2 * (t125 + t282);
+
+    t284 = 2 * px2 * py3 * t147;
+    t285 = x2 * (t284 + t156);
+    t286 = px1 * t103;
+
+    t287 = t98 + t68;
+    t288 = x2 * t287 * y2;
+    t289 = x2 * t231;
+    t290 = t289 + t288;
+
+    t291 = 2 * px2;
+    t292 = - px3 * x2 * y3;
+    t293 = 2 * px3 * x3 * y3;
+
+    t294 = t293 + t203 + t292 + (t68 + (t72 + t291) * x2) * y2;
+    t295 = px1 * t242;
+
+    t296 = - 2 * px3 * x3;
+    t297 = t296 + t64;
+    t298 = px3 * x2 * y3;
+    t299 = y2 * (t128 + t298);
+
+    t300 = - 2 * px2 * t147;
+    t301 = x2 * (t167 + t300) + t299 + t297 * t142;
+    t302 = py1 * t71;
+
+    t303 = py1 * t290;
+    t304 = 2 * py2 * x3;
+    t305 = - 2 * py3 * x3;
+    t306 = - 2 * py2 * x3 * y3;
+
+    t307 = 2 * py3 * x3 * y3;
+    t308 = t307 + t306;
+    t309 = - 2 * px2 * py3 * x3;
+
+    t310 = (t309 + t19 + t51) * y2;
+    t311 = - 2 * px3 * py2 * y3;
+    t312 = t35 + t311;
+
+    t313 = x2 * t312;
+    t314 = 2 * px2 * x3;
+    t315 = 2 * px3 * y3;
+    t316 = t315 + t74;
+
+    t317 = x2 * t316;
+    t318 = t317 + t87 + (t65 + t314 + t85) * y2;
+    t319 = t106 * x2;
+
+    t320 = px1 * (t256 + t118 + t117 + (t115 + t114 + t319) * y2);
+    t321 = py1 * t216;
+
+    t322 = 2 * px2 * py3 * x3 * y3;
+    t323 = 2 * px3 * py2 * y3;
+    t324 = t153 + t323;
+
+    t325 = x2 * t324;
+    t326 = y2 * (t325 + t322 + t143);
+    t327 = - 2 * px2 * x3 * y3;
+
+    t328 = - 2 * px3 * y3;
+    t329 = t328 + t164;
+    t330 = x2 * t329;
+
+    t331 = y2 * (t330 + t204 + t327);
+    t332 = t226 + t331 + t220;
+    t333 = t116 * t142;
+
+    t334 = t140 + t118 + t117;
+    t335 = y2 * t334;
+    t336 = x2 * t173;
+
+    t337 = px1 * (t336 + t335 + t333);
+    t338 = t26 + t94 + t96;
+    t339 = t17 * y2;
+
+    t340 = t153 + t55 + t339;
+    t341 = px2 * px3 * t142;
+    t342 = - 2 * px2 * px3 * y2 * y3;
+
+    t343 = px2 * px3 * t147;
+    t344 = py1 * (t343 + t342 + t341);
+    t345 = - px2 * py3 * t142;
+
+    t346 = y2 * (t35 + t55);
+    t347 = t156 + t346 + t345;
+    t348 = px1 * t347 + t344;
+
+    t349 = t89 + t164 + t62 * y2;
+    t350 = - px2 * px3 * t142;
+    t351 = 2 * px2 * px3 * y2 * y3;
+
+    t352 = - px2 * px3 * t147;
+    t353 = px2 * t142;
+    t354 = y2 * (t89 + t74);
+
+    t355 = t167 + t354 + t353;
+    t356 = px1 * t355 + t352 + t351 + t350;
+    t357 = py1 * t66;
+
+    t358 = py1 * t349;
+    t359 = 2 * py2;
+    t360 = - 2 * py3;
+    t361 = - 2 * py2 * y3;
+
+    t362 = 2 * py3 * y3;
+    t363 = px3 * py2 * t142;
+    t364 = y2 * (t153 + t34);
+
+    t365 = - px3 * t142;
+    t366 = y2 * (t75 + t164);
+    t367 = t166 + t366 + t365;
+
+    t368 = py1 * t367;
+    t369 = px1 * (t172 + t171 + t106 * t142);
+    t370 = t35 + t34;
+
+    t371 = t142 * t370;
+    t372 = y2 * t150;
+    t373 = t372 + t371;
+    t374 = t230 + t229;
+
+    t375 = py1 * (t352 + t351 + t350);
+    t376 = t157 + t364 + t363;
+    t377 = px1 * t376 + t375;
+
+    t378 = t75 + t74;
+    t379 = y2 * t207 + t142 * t378;
+    t380 = px1 * t367 + t343 + t342 + t341;
+
+    t381 = py1 * t209;
+    t382 = py1 * t355;
+    t383 = py1 * t379;
+    t384 = 2 * py2 * y3;
+
+    t385 = - 2 * py3 * y3;
+    t386 = t385 + t384;
+    t387 = - 2 * py2 * t147;
+    t388 = 2 * py3 * t147;
+
+    t389 = px2 * py3 * t2;
+    t390 = t389 + t10;
+    t391 = x2 * t390 * y2;
+    t392 = t5 * t228;
+
+    t393 = - px2 * t2;
+    t394 = t70 + t393;
+    t395 = x2 * t394 * y2;
+    t396 = t5 * t231;
+
+    t397 = t396 + t395;
+    t398 = py1 * t397;
+    t399 = py2 * t2;
+    t400 = - py3 * t2;
+
+    t401 = t400 + t399;
+    t402 = x2 * t401 * y2;
+    t403 = t136 + t135;
+    t404 = t5 * t403;
+
+    t405 = t404 + t402;
+    t406 = px1 * t405;
+    t407 = t1 * (t406 + t398 + t392 + t391);
+
+    t408 = t65 + t64;
+    t409 = t5 * t408;
+    t410 = x2 * t394;
+    t411 = t410 + t409;
+
+    t412 = py1 * t411;
+    t413 = t5 * t116;
+    t414 = x2 * t401;
+    t415 = t414 + t413;
+
+    t416 = px1 * t415;
+    t417 = py2 * t5;
+    t418 = x2 * (t134 + t114);
+    t419 = py3 * t2;
+
+    t420 = t419 + t418 + t417;
+    t421 = px1 * t420;
+    t422 = t265 * y2;
+    t423 = x2 * t154;
+
+    t424 = px2 * x2;
+    t425 = (t68 + t424) * y2;
+    t426 = - py2 * x2;
+    t427 = (t133 + t426) * y2;
+
+    t428 = py3 * x2 * y3;
+    t429 = t20 + t19;
+    t430 = x2 * t429;
+    t431 = - px2 * py3 * t2;
+
+    t432 = (t431 + t43 + t430) * y2;
+    t433 = t5 * t370;
+    t434 = x2 * t145;
+
+    t435 = - px2 * x2 * x3;
+    t436 = px2 * t2;
+    t437 = (t436 + t435) * y2;
+    t438 = px3 * t5 * y3;
+
+    t439 = - px3 * x2 * x3 * y3;
+    t440 = py2 * x2 * x3;
+    t441 = - py2 * t2;
+
+    t442 = (t441 + t440) * y2;
+    t443 = - py3 * t5 * y3;
+    t444 = py3 * x2 * x3 * y3;
+
+    t445 = t5 * t287;
+    t446 = t78 + t436;
+    t447 = x2 * t446;
+    t448 = - t2;
+
+    t449 = t448 + 2 * x2 * x3 - t5;
+    t450 = px1 * t449;
+    t451 = (t98 + t85) * y2;
+    t452 = - x2 * y3;
+
+    t453 = x3 * y3;
+    t454 = t453 + t452 + (x2 - x3) * y2;
+    t455 = px1 * t454;
+    t456 = t65 + t314;
+
+    t457 = x2 * t456;
+    t458 = (t78 + t457) * y2;
+    t459 = x2 * (t293 + t203);
+
+    t460 = - x2 * x3 * y3 + t5 * y3 + (t2 - x2 * x3) * y2;
+    t461 = px1 * t460;
+    t462 = t5 * t253;
+
+    t463 = t419 + t441;
+    t464 = x2 * t463;
+    t465 = - py2 * t5;
+    t466 = x2 * (t115 + t133);
+
+    t467 = t2 - 2 * x2 * x3 + t5;
+    t468 = py1 * t467;
+    t469 = py2 * x2;
+    t470 = (t134 + t469) * y2;
+
+    t471 = - py2 * x2 * y3;
+    t472 = x2 * y3;
+    t473 = - x3 * y3;
+    t474 = t473 + t472 + (x3 - x2) * y2;
+
+    t475 = py1 * t474;
+    t476 = - 2 * py2 * x3;
+    t477 = t115 + t476;
+    t478 = x2 * t477;
+
+    t479 = (t419 + t478) * y2;
+    t480 = py2 * t5 * y3;
+    t481 = - 2 * py3 * x3 * y3;
+
+    t482 = x2 * (t481 + t117);
+    t483 = x2 * x3 * y3 - t5 * y3 + (t448 + x2 * x3) * y2;
+
+    t484 = py1 * t483;
+    t485 = t431 + t43;
+    t486 = t485 * t142;
+    t487 = t5 * t158;
+
+    t488 = t446 * t142;
+    t489 = t5 * t168;
+    t490 = t489 + t488;
+    t491 = py1 * t490;
+
+    t492 = t463 * t142;
+    t493 = t5 * t173;
+    t494 = t493 + t492;
+    t495 = px1 * t494;
+
+    t496 = x1 * y1 * (t495 + t491 + t487 + t486);
+    t497 = t142 * t119;
+    t498 = x2 * y2 * t259;
+
+    t499 = t498 + t497;
+    t500 = px1 * t499;
+    t501 = t29 * (t500 + t381 + t151 + t146);
+
+    t502 = t429 * t142;
+    t503 = x2 * t370;
+    t504 = y2 * (t503 + t125 + t54);
+    t505 = x2 * t158;
+
+    t506 = - px3 * x3 * t142;
+    t507 = - px2 * x2 * t147;
+    t508 = py3 * x3 * t142;
+
+    t509 = y2 * (t118 + t471);
+    t510 = py2 * x2 * t147;
+    t511 = - py2 * t142;
+
+    t512 = y2 * (t138 + t107);
+    t513 = t172 + t512 + t511;
+    t514 = px1 * t513;
+
+    t515 = y2 * t259 + t142 * t255;
+    t516 = px1 * t515;
+    t517 = py1 * t454;
+
+    t518 = - py2 * x3 * t142;
+    t519 = t108 + t384;
+    t520 = x2 * t519;
+
+    t521 = y2 * (t520 + t307 + t135);
+    t522 = - py3 * x2 * t147;
+    t523 = py2 * t142;
+
+    t524 = y2 * (t108 + t137);
+    t525 = - t147 + 2 * y2 * y3 - t142;
+    t526 = py1 * t525;
+
+    t527 = x2 * t147 + y2 * (t473 + t452) + x3 * t142;
+    t528 = py1 * t527;
+    t529 = px1 * t474;
+
+    t530 = px2 * x3 * t142;
+    t531 = px3 * x2 * t147;
+
+    t532 = - x2 * t147 + y2 * (t453 + t472) - x3 * t142;
+    t533 = px1 * t532;
+
+    t534 = - px2 * t142;
+    t535 = t147 - 2 * y2 * y3 + t142;
+    t536 = px1 * t535;
+
+    t537 = t447 + t445;
+    t538 = py1 * t537;
+    t539 = t464 + t462;
+    t540 = px1 * t539;
+
+    t541 = 2 * px3 * py2 * t2;
+    t542 = - 2 * px2 * py3 * t2;
+    t543 = x2 * t446 * y2;
+
+    t544 = t5 * t205;
+    t545 = t544 + t543;
+    t546 = py1 * t545;
+    t547 = x2 * t463 * y2;
+
+    t548 = t5 * t119;
+    t549 = t548 + t547;
+    t550 = px1 * t549;
+    t551 = x2 * t265;
+
+    t552 = (t389 + t10 + t551) * y2;
+    t553 = t5 * t154;
+    t554 = 2 * px3 * t2;
+
+    t555 = (t554 + t393 + t110) * y2;
+    t556 = t5 * t90;
+    t557 = py3 * x2 * x3;
+
+    t558 = - 2 * py3 * t2;
+    t559 = (t558 + t399 + t557) * y2;
+    t560 = py2 * x2 * x3 * y3;
+
+    t561 = t138 + t361;
+    t562 = t5 * t561;
+    t563 = t390 * t142;
+    t564 = t5 * t150;
+
+    t565 = - px2 * t2 * t142;
+    t566 = - px3 * t5 * t147;
+    t567 = t566 + t214 + t565;
+
+    t568 = py1 * t567;
+    t569 = py2 * t2 * t142;
+    t570 = x2 * y2 * (t118 + t135);
+
+    t571 = py3 * t5 * t147;
+    t572 = t571 + t570 + t569;
+    t573 = px1 * t572;
+    t574 = t86 + t68;
+
+    t575 = x2 * t574;
+    t576 = (t78 + t575) * y2;
+    t577 = 2 * px2 * x3 * y3;
+
+    t578 = x2 * (t87 + t577);
+    t579 = px1 * t527;
+
+    t580 = - t5 * t147 + 2 * x2 * x3 * y2 * y3 - t2 * t142;
+    t581 = px1 * t580;
+    t582 = t305 + t133;
+
+    t583 = x2 * t582;
+    t584 = (t419 + t583) * y2;
+    t585 = x2 * (t136 + t306);
+
+    t586 = py1 * t532;
+    t587 = - py3 * t2 * t142;
+    t588 = x2 * y2 * (t136 + t117);
+
+    t589 = - py2 * t5 * t147;
+    t590 = t5 * t147 - 2 * x2 * x3 * y2 * y3 + t2 * t142;
+
+    t591 = py1 * t590;
+    t592 = t400 + t466 + t465;
+    t593 = px1 * t592;
+    t594 = t309 + t279;
+
+    t595 = t198 + t311;
+    t596 = x2 * t378;
+    t597 = t596 + t408 * y2;
+    t598 = py1 * t597;
+
+    t599 = t256 + t116 * y2;
+    t600 = px1 * t599;
+    t601 = t178 + t366 + t534;
+
+    t602 = py1 * t601;
+    t603 = t181 + t524 + t523;
+    t604 = px1 * t603;
+    t605 = t265 * t142;
+
+    t606 = t423 + t144 + t143;
+    t607 = y2 * t606;
+    t608 = x2 * t150;
+    t609 = 2 * py2 * x3 * y3;
+
+    t610 = t362 + t137;
+    t611 = x2 * t610;
+    t612 = y2 * (t611 + t118 + t609);
+
+    t613 = py1 * t449;
+    t614 = t419 + t613 + t418 + t417;
+    t615 = py1 * t460;
+
+    t616 = py1 * t535;
+    t617 = t616 + t172 + t512 + t511;
+    t618 = t134 + t304;
+
+    t619 = t618 * t142;
+    t620 = - py3 * x2 * y3;
+    t621 = y2 * (t135 + t620);
+
+    t622 = x2 * (t388 + t180);
+    t623 = px1 * t467;
+    t624 = t623 + t78 + t102 + t101;
+
+    t625 = px1 * t483;
+    t626 = px1 * t525;
+    t627 = t167 + t626 + t354 + t353;
+
+    t628 = - 2 * px2 * x3;
+    t629 = t98 + t628;
+    t630 = t629 * t142;
+    t631 = - 2 * px3 * t147;
+
+    t632 = x2 * (t631 + t177);
+    t633 = - 2 * px2 * py3 * x3 * y3;
+    t634 = t633 + t197;
+
+    t635 = - 2 * px3 * py2 * t147;
+    t636 = t142 * t403;
+    t637 = x2 * y2 * t173;
+
+    t638 = t637 + t636;
+    t639 = px1 * t638;
+    t640 = t589 + t588 + t587;
+    t641 = px1 * t640;
+
+    t642 = px1 * t590;
+    t643 = py1 * t580;
+
+    t644 = (x0 * (px0 * (y1 * (x1 * (t528 + t522 + t612 + t518)
+				 + t643 + t571 + t570 + t569)
+			   + t29 * t515 + x1 * t638 + t1 * (t615 + t444 + t443 + t442))
+		    + py0 * (y1 * (x1 * (t533 + t531 + t331 + t530)
+				   + t642 + t566 + t214 + t565)
+			     + x1 * t234 + t29 * t379 + t1 * (t625 + t439 + t438 + t437))
+		    + y1 * (x1 * (px1 * (t622 + t621 + t619) + py1 * (t632 + t299 + t630)
+				  + t608 + t607 + t605)
+			    + t641 + t243 + t564 + t563)
+		    + x1 * (t639 + t235 + x2 * y2 * (t284 + t635) + t142 * t634)
+		    + t29 * (t175 + t170)
+		    + t1 * (px1 * (t482 + t480 + t479) + py1 * (t459 + t79 + t458) + t434
+			    + t433 + t432))
+	      + y0 * (x0 * (py0 * (x1 * (t579 + t632 + t299 + t630)
+				   + t489 + t29 * t627
+				   + y1 * (x1 * t597 + t625 + t556 + t112 + t555) + t488
+				   + t624 * t1)
+			    + px0 * (x1 * (t586 + t622 + t621 + t619)
+				     + t29 * t617 + t493
+				     + y1 * (x1 * t599 + t615 + t562 + t560 + t559) + t492
+				     + t614 * t1)
+			    + x1 * (px1 * (t522 + t612 + t518) + py1 * (t531 + t331 + t530)
+				    + t608 + t607 + t605)
+			    + t29 * (t604 + t602) + t487
+			    + y1 * (x1 * (t600 + t598 + x2 * t595 + t594 * y2)
+				    + px1 * (t585 + t480 + t584) + py1 * (t578 + t79 + t576) + t267
+				    + t553 + t552) + t486 + (t593 + t302) * t1)
+		      + px0 * (x1 * (t591 + t589 + t588 + t587)
+			       + t29 * (t586 + t510 + t509 + t508)
+			       + y1 * (x1 * (t484 + t585 + t480 + t584) + t548 + t547) + t415 * t1)
+		      + py0 * (x1 * (t581 + t241 + t240 + t239)
+			       + t29 * (t579 + t507 + t250 + t506)
+			       + y1 * (x1 * (t461 + t578 + t79 + t576) + t544 + t543) + t411 * t1)
+		      + x1 * (t573 + t568 + t564 + t563)
+		      + t29 * (px1 * (t522 + t521 + t518) + py1 * (t531 + t225 + t530) + t505
+			       + t504 + t502)
+		      + y1 * (x1 * (px1 * (t562 + t560 + t559) + py1 * (t556 + t112 + t555)
+				    + t267 + t553 + t552)
+			      + t550 + t546 + t5 * (t322 + t273) + x2 * (t542 + t541) * y2)
+		      + (t540 + t538) * t1)
+	      + t161 * (py0 * (y1 * (x1 * (t536 + t178 + t366 + t534)
+				     + t533 + t531 + t225 + t530)
+			       + x1 * t169 + t208 + t1 * (t529 + t204 + t292 + t425) + t206)
+			+ px0 * (y1 * (t528 + x1 * (t181 + t526 + t524 + t523) + t522 + t521
+				       + t518)
+				 + x1 * t174 + t498 + t1 * (t517 + t118 + t428 + t427) + t497)
+			+ x1 * (t516 + t383)
+			+ y1 * (x1 * (t514 + t382) + px1 * (t510 + t509 + t508)
+				+ py1 * (t507 + t250 + t506) + t505 + t504
+				+ t502) + t151
+			+ t1 * (px1 * (t136 + t471 + t470) + py1 * (t87 + t249 + t451) + t423
+				+ t422) + t146) + t501 + t496
+	      + t14 * (px0 * (x1 * (t484 + t482 + t480 + t479)
+			      + t29 * (t475 + t136 + t471 + t470) + t404 + t402
+			      + (x1 * (t468 + t400 + t466 + t465) + t464 + t462) * y1)
+		       + py0 * (x1 * (t461 + t459 + t79 + t458)
+				+ t29 * (t455 + t87 + t249 + t451) + t396 + t395
+				+ (x1 * (t70 + t450 + t69 + t67) + t447 + t445) * y1)
+		       + x1 * (px1 * (t444 + t443 + t442) + py1 * (t439 + t438 + t437) + t434
+			       + t433 + t432)
+		       + t29 * (px1 * (t118 + t428 + t427) + py1 * (t204 + t292 + t425) + t423
+				+ t422) + t392 + t391
+		       + (x1 * (t421 + t104) + t416 + t412) * y1) + t407);
+    t645 = t5 * t265;
+
+    t646 = t115 + t114 + t132;
+    t647 = px1 * t646;
+    t648 = x2 * t485;
+    t649 = t32 * t5;
+
+    t650 = t70 + t393 + t73 * t5;
+    t651 = t400 + t399 + t106 * t5;
+
+    t652 = t540 + x1 * (px1 * t651 + py1 * t650 + t389 + t10 + t649) + t538 + t648
+	+ t29 * (t647 + t357 + t20 + t19 + t18) + t645;
+    t653 = t648 + t645;
+
+    t654 = t392 + t391;
+    t655 = px1 * t654;
+    t656 = t309 + t19;
+    t657 = x2 * t656;
+
+    t658 = (t389 + t657) * y2;
+    t659 = px3 * py2 * t5 * y3;
+    t660 = x2 * (t144 + t273);
+
+    t661 = - px3 * py2 * t5;
+    t662 = t431 + t27 + t661;
+    t663 = px1 * t662 + t24;
+
+    t664 = t5 * t429;
+    t665 = x2 * t390;
+    t666 = t665 + t664;
+    t667 = px3 * py2 * x2;
+
+    t668 = (t20 + t667) * y2;
+    t669 = x2 * t485 * y2;
+    t670 = t5 * t145;
+    t671 = t670 + t669;
+
+    t672 = px1 * t671;
+    t673 = t26 + t52;
+    t674 = x2 * t673;
+    t675 = (t389 + t674) * y2;
+
+    t676 = x2 * (t633 + t54);
+    t677 = px3 * t5;
+    t678 = t436 + t69 + t677;
+
+    t679 = px1 * t678 + t37 + t60 + t59;
+    t680 = - px3 * x2;
+
+    t681 = t203 + t298 + (t64 + t680) * y2;
+    t682 = px1 * t545;
+    t683 = - px3 * t5 * y3;
+
+    t684 = t578 + t683 + (t393 + t575) * y2;
+    t685 = 2 * py3 * x3;
+    t686 = t685 + t476;
+
+    t687 = 2 * py2 * t2;
+    t688 = px1 * (t419 + t441 + t131 * t5);
+    t689 = - px2 * py3 * x2;
+
+    t690 = 2 * px2 * py3 * x3;
+    t691 = (t690 + t94 + t689) * y2;
+
+    t692 = t330 + t204 + (t98 + t628 + t424) * y2;
+    t693 = t134 + t133 + t319;
+
+    t694 = px1 * (t140 + t118 + t117 + t693 * y2);
+    t695 = (t542 + t43 + t9) * y2;
+
+    t696 = t5 * t312;
+    t697 = 2 * px2 * t2;
+    t698 = t5 * t316 + t112 + (t78 + t697 + t110) * y2;
+
+    t699 = x2 * t253;
+    t700 = t5 * t255;
+    t701 = x2 * t403;
+
+    t702 = px1 * (t701 + t700 + (t419 + t441 + t699) * y2);
+    t703 = px2 * py3 * x2 * x3;
+
+    t704 = (t10 + t703) * y2;
+    t705 = px3 * py2 * x2 * x3 * y3;
+    t706 = (t20 + t279 + t689) * y2;
+
+    t707 = t439 + t111 + (t70 + t435) * y2;
+    t708 = t224 + t204 + (t296 + t64 + t424) * y2;
+
+    t709 = - 2 * py2;
+    t710 = 2 * py3;
+    t711 = py1 * t678;
+
+    t712 = t459 + t683 + (t393 + t457) * y2;
+    t713 = x2 * t116;
+    t714 = t5 * t139;
+
+    t715 = px1 * (t120 + t714 + (t400 + t399 + t713) * y2);
+    t716 = 2 * px2 * py3;
+
+    t717 = (t94 + (t716 + t15) * x2) * y2;
+    t718 = - 2 * px2;
+
+    t719 = t221 + t128 + t249 + (t98 + (px3 + t718) * x2) * y2;
+
+    t720 = px1 * (t256 + t136 + t135 + t646 * y2);
+    t721 = - px2 * py3 * t2 * t142;
+
+    t722 = - px3 * py2 * t5 * t147;
+    t723 = t722 + t237 + t721;
+    t724 = - px2 * py3 * x3 * t142;
+
+    t725 = y2 * (t54 + t124);
+    t726 = px1 * y2 * t257;
+    t727 = - px3 * py2 * x2 * t147;
+
+    t728 = y2 * (t87 + t127);
+    t729 = t531 + t728 + t530;
+    t730 = px2 * py3 * t2 * t142;
+
+    t731 = px3 * py2 * t5 * t147;
+    t732 = px1 * t397;
+    t733 = t251 + t299 + t248;
+
+    t734 = px2 * t2 * t142;
+    t735 = px3 * t5 * t147;
+    t736 = t735 + t240 + t734;
+
+    t737 = t389 + t10 + t649;
+    t738 = t731 + t189 + t730;
+    t739 = px1 * t738;
+
+    t740 = x2 * t165;
+    t741 = t740 + t204 + t203;
+    t742 = py1 * y2 * t741;
+    t743 = py1 * t736;
+
+    t744 = px2 * py3 * t142;
+    t745 = px1 * t567;
+    t746 = t148 + t364 + t744;
+
+    t747 = px3 * py2 * t5;
+    t748 = t389 + t95 + t747;
+    t749 = (t26 + t122) * y2;
+
+    t750 = x2 * t280;
+    t751 = (t431 + t750) * y2;
+    t752 = - px3 * py2 * t5 * y3;
+
+    t753 = x2 * (t322 + t143);
+    t754 = - px3 * t5;
+    t755 = t393 + t102 + t754;
+
+    t756 = t128 + t292 + (t68 + t126) * y2;
+    t757 = x2 * t297;
+    t758 = x2 * (t204 + t327);
+
+    t759 = t758 + t438 + (t436 + t757) * y2;
+    t760 = (t94 + t667) * y2;
+
+    t761 = t203 + t249 + (t98 + t680) * y2;
+    t762 = px1 * (t140 + t253 * y2);
+
+    t763 = - px3 * py2 * x2 * x3;
+    t764 = (t43 + t763) * y2;
+    t765 = - px2 * py3 * x2 * x3 * y3;
+
+    t766 = px3 * x2 * x3;
+    t767 = px2 * x2 * x3 * y3;
+    t768 = t767 + t79 + (t78 + t766) * y2;
+
+    t769 = px1 * (t120 + t700 + (t419 + t441 + t713) * y2);
+    t770 = t501 + t496 + t407;
+
+    t771 = px3 * py2 * x3 * t142;
+    t772 = y2 * (t313 + t633 + t54);
+
+    t773 = px2 * py3 * x2 * t147;
+    t774 = - px3 * py2 * t142;
+    t775 = t149 + t346 + t774;
+
+    t776 = y2 * (t317 + t87 + t577);
+    t777 = t507 + t776 + t506;
+    t778 = px3 * t142;
+
+    t779 = t177 + t354 + t778;
+    t780 = y2 * (t144 + t272);
+    t781 = y2 * (t203 + t292);
+
+    t782 = t531 + t781 + t530;
+    t783 = px1 * (t336 + t258 + t333);
+    t784 = t690 + t94;
+
+    t785 = x2 * t784;
+    t786 = (t431 + t785) * y2;
+    t787 = x2 * (t125 + t197);
+
+    t788 = x2 * t629;
+    t789 = x2 * (t221 + t128);
+    t790 = t789 + t438 + (t436 + t788) * y2;
+
+    t791 = - 2 * py2 * t2;
+    t792 = 2 * py3 * t2;
+    t793 = 2 * px2 * py3 * t2;
+
+    t794 = (t793 + t10 + t42) * y2;
+    t795 = t5 * t324;
+    t796 = - 2 * px2 * t2;
+
+    t797 = t5 * t329 + t80 + (t70 + t796 + t77) * y2;
+
+    t798 = px1 * (t701 + t714 + (t400 + t399 + t699) * y2);
+
+    t799 = px1 * (t5 * t259 + t401 * t142);
+    t800 = t429 * y2;
+    t801 = t503 + t800;
+
+    t802 = t487 + t486;
+    t803 = t673 * t142;
+    t804 = - 2 * px2 * py3 * t147;
+
+    t805 = x2 * (t804 + t148);
+    t806 = 2 * px2 * t147;
+
+    t807 = x2 * (t178 + t806) + t728 + t574 * t142;
+    t808 = py1 * t755;
+    t809 = py1 * t779;
+
+    t810 = y2 * (t58 + t144 + t273);
+    t811 = y2 * (t91 + t293 + t203);
+
+    t812 = t507 + t811 + t506;
+    t813 = px1 * (t260 + t335 + t254);
+    t814 = 2 * py2 * t147;
+
+    t815 = - 2 * py3 * t147;
+    t816 = (t389 + t42) * y2;
+    t817 = - py2 * py3 * t2;
+
+    t818 = (t817 + py2 * py3 * x2 * x3) * y2;
+    t819 = - py2 * py3 * t5 * y3;
+
+    t820 = py2 * py3 * x2 * x3 * y3;
+    t821 = px1 * (t820 + t819 + t818);
+    t822 = - py2 * py3 * t5;
+
+    t823 = 2 * py2 * py3 * x2 * x3;
+    t824 = px1 * (t817 + t823 + t822);
+    t825 = (t431 + t9) * y2;
+
+    t826 = py2 * py3 * t2;
+    t827 = (t826 - py2 * py3 * x2 * x3) * y2;
+    t828 = py2 * py3 * t5 * y3;
+
+    t829 = - py2 * py3 * x2 * x3 * y3;
+    t830 = px1 * (t829 + t828 + t827);
+
+    t831 = (py2 * py3 * x2 - py2 * py3 * x3) * y2;
+    t832 = - py2 * py3 * x2 * y3;
+
+    t833 = py2 * py3 * x3 * y3;
+    t834 = px1 * (t833 + t832 + t831);
+
+    t835 = (t690 + t94 + t122) * y2;
+    t836 = px1 * t693;
+    t837 = - py2 * t5 * y3;
+
+    t838 = t560 + t837 + (t400 + t557) * y2;
+    t839 = x2 * t205;
+
+    t840 = py1 * (t839 + x2 * t408 * y2);
+    t841 = (t20 + t51) * y2;
+    t842 = - py3 * x2;
+
+    t843 = py2 * x2 * y3;
+    t844 = t135 + t843 + (t115 + t842) * y2;
+
+    t845 = py1 * (t740 + t87 + t128 + (t98 + t68 + t63) * y2);
+    t846 = py2 * py3 * t5;
+
+    t847 = - 2 * py2 * py3 * x2 * x3;
+    t848 = - py2 * x2 * x3;
+    t849 = - py3 * x2 * x3 * y3;
+
+    t850 = t849 + t480 + (t419 + t848) * y2;
+    t851 = (py2 * py3 * x3 - py2 * py3 * x2) * y2;
+
+    t852 = py2 * py3 * x2 * y3;
+    t853 = - py2 * py3 * x3 * y3;
+    t854 = x2 * t561;
+
+    t855 = t854 + t136 + (t305 + t133 + t469) * y2;
+    t856 = py2 * py3 * t2 * t142;
+
+    t857 = - 2 * py2 * py3 * x2 * x3 * y2 * y3;
+    t858 = py2 * py3 * t5 * t147;
+
+    t859 = px1 * (t858 + t857 + t856);
+    t860 = - py2 * py3 * x3 * t142;
+
+    t861 = y2 * (t833 + t852);
+    t862 = - py2 * py3 * x2 * t147;
+
+    t863 = px1 * (t862 + t861 + t860);
+    t864 = - py2 * py3 * t2 * t142;
+
+    t865 = 2 * py2 * py3 * x2 * x3 * y2 * y3;
+    t866 = - py2 * py3 * t5 * t147;
+
+    t867 = py3 * t2 * t142;
+    t868 = py2 * t5 * t147;
+    t869 = t868 + t570 + t867;
+
+    t870 = py2 * py3 * x3 * t142;
+    t871 = y2 * (t853 + t832);
+    t872 = py2 * py3 * x2 * t147;
+
+    t873 = - py3 * x3 * t142;
+    t874 = - py2 * x2 * t147;
+    t875 = t874 + t521 + t873;
+
+    t876 = py2 * x3 * t142;
+    t877 = py3 * x2 * t147;
+    t878 = t877 + t509 + t876;
+
+    t879 = t287 * t142;
+    t880 = t596 + t87 + t128;
+    t881 = y2 * t880;
+    t882 = x2 * t207;
+
+    t883 = py1 * (t882 + t881 + t879);
+    t884 = py1 * t662;
+
+    t885 = px1 * (t826 + t847 + t846);
+    t886 = 2 * px3 * py2;
+
+    t887 = (t94 + (t31 + t886) * x2) * y2;
+    t888 = px1 * (t853 + t852 + t851);
+
+    t889 = py1 * t738;
+    t890 = px1 * (t866 + t865 + t864);
+
+    t891 = px1 * (t872 + t871 + t870);
+    t892 = t656 * t142;
+    t893 = x2 * (t157 + t635);
+
+    t894 = t221 + t577;
+    t895 = x2 * t253 * y2;
+    t896 = t701 + t895;
+    t897 = px1 * t896;
+
+    t898 = (t20 + t279 + t122) * y2;
+
+    t899 = py1 * (t596 + t204 + t203 + (t65 + t64 + t97) * y2);
+    t900 = t385 + t107;
+
+    t901 = x2 * t900;
+    t902 = t901 + t136 + (t115 + t476 + t469) * y2;
+    t903 = px1 * t869;
+
+    t904 = t874 + t612 + t873;
+    t905 = t408 * t142;
+    t906 = y2 * t741;
+    t907 = x2 * t168;
+
+    t908 = py1 * (t907 + t906 + t905);
+    t909 = - py2 * py3 * t142;
+
+    t910 = 2 * py2 * py3 * y2 * y3;
+    t911 = - py2 * py3 * t147;
+
+    t912 = px1 * (t911 + t910 + t909);
+    t913 = t912 + py1 * t376;
+
+    t914 = t481 + t117 + t428 + (t133 + (py3 + t709) * x2) * y2;
+    t915 = 2 * px3;
+
+    t916 = t138 + t137 + t131 * y2;
+    t917 = px1 * t916;
+
+    t918 = py1 * (t167 + t166 + t73 * t142);
+    t919 = py3 * t142;
+    t920 = t171 + t524 + t919;
+
+    t921 = px1 * t920;
+    t922 = py2 * py3 * t142;
+    t923 = - 2 * py2 * py3 * y2 * y3;
+
+    t924 = py2 * py3 * t147;
+    t925 = py1 * t513 + t924 + t923 + t922;
+    t926 = py1 * t420;
+
+    t927 = py1 * t640;
+    t928 = t685 + t114;
+    t929 = x2 * (t172 + t814) + t621 + t928 * t142;
+
+    t930 = px1 * (t924 + t923 + t922);
+    t931 = t930 + py1 * t347;
+
+    t932 = py1 * t920 + t911 + t910 + t909;
+    t933 = t315 + t222;
+    t934 = py1 * t654;
+
+    t935 = (t10 + t750) * y2;
+    t936 = t824 + py1 * t263;
+    t937 = py1 * t671;
+
+    t938 = (t19 + t689) * y2;
+    t939 = (t10 + t785) * y2;
+    t940 = t296 + t314;
+
+    t941 = py1 * (t78 + t436 + t62 * t5);
+    t942 = (t26 + t52 + t667) * y2;
+
+    t943 = py1 * (t740 + t204 + t203 + t99 * y2);
+
+    t944 = t611 + t118 + (t134 + t304 + t426) * y2;
+    t945 = (t431 + t541 + t42) * y2;
+
+    t946 = t5 * t199;
+    t947 = t5 * t900 + t560 + (t419 + t791 + t557) * y2;
+    t948 = x2 * t287;
+
+    t949 = t5 * t378;
+    t950 = py1 * (t289 + t949 + (t78 + t436 + t948) * y2);
+
+    t951 = - py3 * t5;
+    t952 = t441 + t466 + t951;
+    t953 = py1 * t952 + t826 + t847 + t846;
+
+    t954 = py3 * x2;
+    t955 = t117 + t620 + (t114 + t954) * y2;
+    t956 = py1 * t549;
+
+    t957 = py3 * t5 * y3;
+    t958 = t585 + t957 + (t399 + t583) * y2;
+    t959 = (t389 + t763) * y2;
+
+    t960 = (t309 + t19 + t667) * y2;
+    t961 = - 2 * px3;
+    t962 = px1 * t952;
+    t963 = x2 * t408;
+
+    t964 = t5 * t165;
+    t965 = py1 * (t839 + t964 + (t70 + t393 + t963) * y2);
+
+    t966 = t482 + t957 + (t399 + t478) * y2;
+    t967 = - 2 * px3 * py2;
+
+    t968 = (t26 + (t16 + t967) * x2) * y2;
+
+    t969 = t307 + t135 + t471 + (t134 + (t130 + t359) * x2) * y2;
+
+    t970 = py1 * (t596 + t87 + t128 + t66 * y2);
+    t971 = t444 + t837 + (t400 + t440) * y2;
+
+    t972 = t520 + t118 + (t685 + t114 + t426) * y2;
+    t973 = py1 * t405;
+
+    t974 = t877 + t621 + t876;
+    t975 = - py2 * t2 * t142;
+    t976 = - py3 * t5 * t147;
+
+    t977 = t976 + t588 + t975;
+    t978 = py1 * y2 * t880;
+    t979 = y2 * (t136 + t843);
+
+    t980 = t522 + t979 + t518;
+    t981 = py1 * t276;
+    t982 = py1 * t572;
+    t983 = px1 * y2 * t334;
+
+    t984 = px1 * t977;
+    t985 = (t94 + t51) * y2;
+    t986 = (t43 + t657) * y2;
+
+    t987 = (t26 + t689) * y2;
+    t988 = t117 + t471 + (t134 + t954) * y2;
+
+    t989 = py1 * (t740 + t287 * y2);
+    t990 = (t431 + t703) * y2;
+    t991 = - py3 * x2 * x3;
+
+    t992 = - py2 * x2 * x3 * y3;
+    t993 = t992 + t480 + (t419 + t991) * y2;
+
+    t994 = py1 * (t839 + t949 + (t78 + t436 + t963) * y2);
+    t995 = py3 * t5;
+
+    t996 = t399 + t418 + t995;
+    t997 = t135 + t428 + (t133 + t842) * y2;
+    t998 = x2 * t928;
+
+    t999 = x2 * (t118 + t609);
+    t1000 = t999 + t443 + (t441 + t998) * y2;
+
+    t1001 = y2 * (t901 + t136 + t306);
+    t1002 = t510 + t1001 + t508;
+    t1003 = - py3 * t142;
+
+    t1004 = t180 + t512 + t1003;
+    t1005 = y2 * (t117 + t428);
+    t1006 = t522 + t1005 + t518;
+
+    t1007 = py1 * (t907 + t881 + t905);
+    t1008 = y2 * (t854 + t481 + t117);
+
+    t1009 = t510 + t1008 + t508;
+    t1010 = 2 * px3 * t147;
+
+    t1011 = py1 * (t5 * t207 + t394 * t142);
+    t1012 = t784 * t142;
+
+    t1013 = 2 * px3 * py2 * t147;
+    t1014 = x2 * (t149 + t1013);
+
+    t1015 = py1 * (t882 + t906 + t879);
+    t1016 = x2 * (t181 + t387) + t979 + t582 * t142;
+
+    t1017 = (t43 + t674) * y2;
+    t1018 = x2 * t618;
+    t1019 = x2 * (t307 + t135);
+
+    t1020 = t1019 + t443 + (t441 + t1018) * y2;
+    t1021 = - 2 * px3 * t2;
+
+    t1022 = - 2 * px3 * py2 * t2;
+    t1023 = (t389 + t1022 + t9) * y2;
+    t1024 = t5 * t57;
+
+    t1025 = t5 * t610 + t849 + (t400 + t687 + t848) * y2;
+
+    t1026 = py1 * (t289 + t964 + (t70 + t393 + t948) * y2);
+    t1027 = px1 * t996;
+
+    t1028 = px1 * t1004;
+    t1029 = x2 * t429 * y2;
+    t1030 = (t436 + t110) * y2;
+
+    t1031 = (t441 + t557) * y2;
+    t1032 = (t393 + t77) * y2;
+    t1033 = (t399 + t848) * y2;
+
+    t1034 = (t26 + t94 + t18) * y2;
+    t1035 = (t64 + t85) * y2;
+    t1036 = (t114 + t469) * y2;
+
+    t1037 = (t98 + t628 + t126) * y2;
+    t1038 = (t134 + t304 + t842) * y2;
+
+    t1039 = (t20 + t19 + t96) * y2;
+    t1040 = (t296 + t64 + t126) * y2;
+
+    t1041 = (t685 + t114 + t842) * y2;
+    t1042 = (t98 + (t961 + px2) * x2) * y2;
+
+    t1043 = t456 * t142;
+    t1044 = x2 * (t1010 + t166);
+
+    t1045 = (t134 + (t710 + t105) * x2) * y2;
+    t1046 = t477 * t142;
+
+    t1047 = x2 * (t815 + t171);
+    t1048 = t32 * t142;
+    t1049 = t171 + t526 + t524 + t919;
+
+    t1050 = t536 + t166 + t366 + t365;
+    t1051 = (t389 + t10 + t430) * y2;
+
+    t1052 = (t393 + t766) * y2;
+    t1053 = (t399 + t991) * y2;
+    t1054 = t17 * t5;
+
+    t1055 = (t431 + t43 + t551) * y2;
+    t1056 = (t1021 + t436 + t77) * y2;
+    t1057 = t5 * t223;
+
+    t1058 = (t792 + t441 + t848) * y2;
+    t1059 = t5 * t519;
+    t1060 = t338 * y2;
+
+    t1061 = (t86 + t68 + t680) * y2;
+    t1062 = (t305 + t133 + t954) * y2;
+
+    t1063 = (t115 + t426) * y2;
+    t1064 = (t400 + t1018) * y2;
+    t1065 = (t65 + t424) * y2;
+
+    t1066 = (t70 + t788) * y2;
+    t1067 = (t70 + t757) * y2;
+    t1068 = (t400 + t998) * y2;
+
+    t1069 = t21 * y2;
+    t1070 = (t68 + (t915 + t61) * x2) * y2;
+
+    t1071 = (t133 + (t360 + py2) * x2) * y2;
+    t1072 = (t115 + t476 + t954) * y2;
+
+    t1073 = (t65 + t314 + t680) * y2;
+
+    trans->m[0][0]
+	= (x0 * (px0 * (x1 * (px1 * (y2 * (t388 + t387) + t142 * t386)
+			      + t383 + t372 + t371)
+			+ y1 * (x1 * (t369 + t382 + t156 + t346 + t345)
+				+ t337 + py1 * t301 + t285 + t283 + t281) + t381 + t151
+			+ t1 * (t141 + py1 * t92 + t58 + t54 + t53) + t146)
+		 + py0 * (y1 * (x1 * t380 + px1 * t332 + t219 + t218 + t217)
+			  + px1 * t234 + px1 * x1 * t379 + t1 * (px1 * t129 + t49 + t48 + t47))
+		 + y1 * (x1 * t377 + px1 * (t202 + t326 + t196) + t195) + px1 * t374
+		 + px1 * x1 * t373 + t1 * (px1 * (t125 + t124 + t123) + t269))
+	   + y0 * (x0 * (px0 * (t261 + x1 * (t369 + t368 + t157 + t364 + t363) + py1 * t227
+				+ t202
+				+ y1
+				* (x1
+				   * (px1 * (t362 + t361 + (t360 + t359) * y2)
+				      + t358 + t153 + t55 + t339)
+				   + t320 + py1 * t294 + t144 + t273 + t272 + t271)
+				+ t201 + t196 + (t357 + t20 + t19 + t18) * t1)
+			 + py0 * (x1 * t356 + px1 * t252 + t194
+				  + y1 * (px1 * t318 + px1 * x1 * t349 + t84 + t83 + t82)
+				  + t193 + t192 + px1 * t99 * t1) + x1 * t348
+			 + px1 * (t247 + t246 + t244) + t278
+			 + y1 * (px1 * (t313 + t54 + t310) + t50 + px1 * x1 * t340)
+			 + px1 * t338 * t1)
+		   + px0 * (x1 * (t337 + py1 * t332 + t202 + t326 + t196)
+			    + t321 + px1 * t29 * t182 + t190
+			    + y1 * (x1 * (t320 + py1 * t318 + t313 + t54 + t310)
+				    + px1 * (x2 * t308 + x2 * (t305 + t304) * y2) + t303 + t267
+				    + t266) + t189 + t188 + (t302 + t10 + t27 + t25) * t1)
+		   + py0 * (x1 * (px1 * t301 + t194 + t193 + t192)
+			    + t295 + px1 * t29 * t179 + t186
+			    + y1 * (x1 * (px1 * t294 + t49 + t48 + t47) + px1 * t290) + t185 + t184
+			    + (t286 + t3 + t23 + t22) * t1)
+		   + x1 * (px1 * (t285 + t283 + t281) + t278) + t277 + t275 + px1 * t29 * t274
+		   + y1 * (x1 * (px1 * (t144 + t273 + t272 + t271) + t269) + px1 * t268)
+		   + (t264 + t262) * t1)
+	   + px0 * (y1 * (x1 * (t261 + py1 * t252 + t247 + t246 + t244)
+			  + t243 + t238 + t237 + t236)
+		    + x1 * (t235 + t230 + t229) + px1 * t29 * t174
+		    + t1 * (t121 + py1 * t81 + t46 + t45 + t44))
+	   + py0 * (y1 * (x1 * (px1 * t227 + t219 + t218 + t217)
+			  + px1 * t216 + t212 + t211 + t210)
+		    + px1 * t29 * t169 + px1 * x1 * t209 + t1 * (px1 * t113 + t40 + t39 + t38))
+	   + y1 * (x1 * (px1 * (t202 + t201 + t196) + t195) + px1 * t191 + t187)
+	   + px0 * t161 * t183 + px1 * t29 * t160 + px1 * x1 * t152
+	   + t14 * (px0 * (x1 * (t141 + py1 * t129 + t125 + t124 + t123)
+			   + t121 + py1 * t113 + px1 * t29 * t109 + t13 + t12 + t11
+			   + (t104 + t43 + x1 * (t100 + t26 + t94 + t96) + t95 + t93) * y1)
+		    + py0 * (x1 * (px1 * t92 + t84 + t83 + t82)
+			     + px1 * t81 + px1 * t29 * t76 + t7 + t6 + t4
+			     + (px1 * t71 + t37 + px1 * x1 * t66 + t60 + t59) * y1)
+		    + x1 * (px1 * (t58 + t54 + t53) + t50) + px1 * (t46 + t45 + t44) + t41
+		    + px1 * t29 * t36 + (px1 * t28 + t24 + px1 * x1 * t21) * y1)
+	   + t1 * (px1 * (t13 + t12 + t11) + t8));
+
+    trans->m[0][1] =
+	(t161 * (px0 * (x1 * (t382 + t156 + t346 + t345)
+			+ py1 * t733 + t247
+			+ y1 * (t694 + x1 * (t358 + t153 + t55 + t339) + py1 * t681
+				+ t144 + t282 + t668) + t726 + t283 + t244
+			+ px1 * t646 * t1)
+		 + py0 * (x1 * (px1 * t601 + t343 + t342 + t341)
+			  + px1 * t729 + t219
+			  + y1 * (px1 * t692 + px1 * x1 * t76 + t49 + t48 + t47) + t218
+			  + t217 + px1 * t66 * t1) + x1 * (px1 * t746 + t375)
+		 + px1 * (t727 + t725 + t724) + t195
+		 + y1 * (px1 * (t325 + t143 + t691) + t269 + px1 * x1 * t36)
+		 + px1 * t21 * t1)
+	 + x0 * (py0 * (t29 * t356 + t745 + t212
+			+ y1
+			* (x1 * (px1 * t719 + t84 + t83 + t82)
+			   + px1 * t698 + t40 + t39 + t38) + px1 * x1 * y2 * t741
+			  + t211 + t210 + px1 * t650 * t1)
+		   + px0 * (t29 * (t602 + t148 + t364 + t744)
+			    + t743 + t722
+			    + y1 * (x1 * (t720 + py1 * t708 + t200 + t143 + t706)
+				    + t702 + py1 * t684 + t676 + t659 + t675)
+			    + x1 * (t607 + px1 * y2 * (x2 * (t362 + t361) + t481 + t609) + t742)
+			    + t237 + t721 + px1 * t651 * t1) + t29 * t348 + t739 + t187
+		   + y1 * (x1 * (px1 * (t125 + t197 + t245 + t717) + t50)
+			   + px1 * (t696 + t13 + t695) + t8) + px1 * x1 * y2 * t606
+		   + px1 * t737 * t1)
+	   + py0 * (x1 * (px1 * t736 + t186 + t185 + t184)
+		    + t29 * (px1 * t733 + t194 + t193 + t192)
+		    + y1 * (x1 * (px1 * t712 + t7 + t6 + t4) + t732) + px1 * t537 * t1)
+	   + px0 * (x1 * (t568 + t731 + t189 + t730)
+		    + t29 * (py1 * t729 + t727 + t726 + t725 + t724)
+		    + y1 * (x1 * (t715 + py1 * t707 + t705 + t12 + t704) + t546 + t670 + t669)
+		    + px1 * t539 * t1) + x1 * (px1 * t723 + t275)
+	   + t29 * (px1 * (t247 + t283 + t244) + t278)
+	   + y0 * (x0 * (px0 * (x1 * (t720 + py1 * t719 + t125 + t197 + t245 + t717)
+				+ t715 + py1 * t712 + t29 * (t162 + t35 + t34 + t33) + t660
+				+ t659 + t658
+				+ (t688 + t711 + t431
+				   + x1
+				   * (px1 * (t305 + t304 + (t710 + t709) * x2)
+				      + t100 + t26 + t94 + t96) + t27 + t661)
+				* y1)
+			 + py0 * (x1 * (px1 * t708 + t49 + t48 + t47)
+				  + px1 * t707 + px1 * t29 * t349 + t40 + t39 + t38
+				  + (t286 + t3 + px1 * x1 * t99 + t23 + t22) * y1)
+			 + x1 * (px1 * (t200 + t143 + t706) + t269) + px1 * (t705 + t12 + t704)
+			 + t8 + px1 * t29 * t340 + (t264 + t262 + px1 * x1 * t338) * y1)
+		   + px0 * (x1 * (t702 + py1 * t698 + t696 + t13 + t695)
+			    + t29 * (t694 + py1 * t692 + t325 + t143 + t691) + t398 + t392 + t391
+			    + (x1 * (t688 + t104 + t43 + t95 + t93)
+			       + px1 * (x2 * (t558 + t687) + t5 * t686) + t412 + t665 + t664)
+			    * y1)
+		   + py0 * (x1 * (px1 * t684 + t7 + t6 + t4) + t682
+			    + t29 * (px1 * t681 + t84 + t83 + t82)
+			    + (px1 * t411 + x1 * t679) * y1)
+		   + x1 * (px1 * (t676 + t659 + t675) + t41) + t672
+		   + t29 * (px1 * (t144 + t282 + t668) + t50) + (px1 * t666 + x1 * t663) * y1)
+	   + y1 * (x1 * (px1 * (t660 + t659 + t658) + t41) + t655) + px1 * t653 * t1
+	   + px0 * t652 * t14)
+	;
+
+    trans->m[0][2] =
+	(x0 * (px0 * (y1 * (x1 * (t813 + py1 * t807 + t805 + t725 + t803)
+				     + t799 + t568 + t731 + t189 + t730)
+			       + x1 * (px1 * (x2 * y2 * (t815 + t814) + t142 * t308)
+				       + t235 + t230 + t229) + t29 * (t170 + t159 + t155)
+			       + t1 * (t769 + py1 * t759 + t753 + t752 + t751))
+			+ py0 * (y1 * (x1 * (px1 * t812 + t194 + t193 + t192)
+				       + t295 + t186 + t185 + t184)
+				 + px1 * x1 * t234 + px1 * t29 * t379
+				 + t1 * (px1 * t768 + t7 + t6 + t4))
+			+ y1 * (x1 * (px1 * (t773 + t810 + t771) + t278) + t277 + t275)
+			+ px1 * x1 * t374 + px1 * t29 * t373
+			+ t1 * (px1 * (t765 + t45 + t764) + t41))
+		  + y0 * (x0 * (px0 * (x1 * (t813 + py1 * t812 + t773 + t810 + t771)
+				       + t495 + t29 * (t809 + t149 + t346 + t774)
+				       + y1
+				       * (x1
+					  * (px1 * (x2 * t386 + t686 * y2)
+					     + t598 + t503 + t800)
+					  + t798 + py1 * t790 + t787 + t752 + t786)
+				       + (t808 + t389 + t95 + t747) * t1)
+				+ py0 * (x1 * (px1 * t807 + t219 + t218 + t217)
+					 + px1 * t490 + t29 * t380
+					 + y1 * (px1 * x1 * t597 + px1 * t797 + t7 + t6 + t4)
+					 + t679 * t1)
+				+ x1 * (px1 * (t805 + t725 + t803) + t195) + px1 * t802
+				+ t29 * t377
+				+ y1 * (px1 * x1 * t801 + px1 * (t795 + t46 + t794) + t41)
+				+ t663 * t1)
+			  + px0 * (x1 * (t799 + t243 + t238 + t237 + t236)
+				   + t29 * (t783 + py1 * t777 + t773 + t772 + t771)
+				   + y1 * (x1 * (t798 + py1 * t797 + t795 + t46 + t794)
+					   + px1 * (t5 * (t481 + t609) + x2 * (t792 + t791) * y2)
+					   + t546 + t670 + t669) + (t538 + t648 + t645) * t1)
+			  + py0 * (x1 * (t745 + t212 + t211 + t210)
+				   + t29 * (px1 * t782 + t219 + t218 + t217)
+				   + y1 * (x1 * (px1 * t790 + t40 + t39 + t38) + t682)
+				   + px1 * t411 * t1) + x1 * (t739 + t187)
+			  + t29 * (px1 * (t727 + t780 + t724) + t195)
+			  + y1 * (x1 * (px1 * (t787 + t752 + t786) + t8) + t672)
+			  + px1 * t666 * t1)
+		  + t161 * (px0 * (y1
+				   * (t783 + x1 * (t368 + t157 + t364 + t363) + py1 * t782
+				      + t727 + t780 + t724)
+				   + x1 * (t383 + t372 + t371) + t500
+				   + t1 * (t762 + py1 * t756 + t125 + t272 + t749))
+			    + py0 * (y1
+				     * (x1 * (px1 * t779 + t352 + t351 + t350)
+					+ px1 * t777 + t194 + t193 + t192)
+				     + px1 * x1 * t169 + px1 * t209
+				     + t1 * (px1 * t761 + t84 + t83 + t82))
+			    + y1 * (x1 * (px1 * t775 + t344) + px1 * (t773 + t772 + t771)
+				    + t278) + px1 * x1 * t160
+			    + px1 * t152 + t1 * (px1 * (t144 + t245 + t760) + t50))
+		  + px0 * t770
+		  + t14 * (px0 * (x1 * (t769 + py1 * t768 + t765 + t45 + t764)
+				  + t29 * (t762 + py1 * t761 + t144 + t245 + t760) + t406
+				  + (t412 + x1 * (t711 + t431 + t27 + t661) + t665 + t664)
+				  * y1)
+			   + py0 * (x1 * (px1 * t759 + t40 + t39 + t38)
+				    + t732 + t29 * (px1 * t756 + t49 + t48 + t47)
+				    + (px1 * t537 + x1 * (px1 * t755 + t3 + t23 + t22)) * y1)
+			   + x1 * (px1 * (t753 + t752 + t751) + t8) + t655
+			   + t29 * (px1 * (t125 + t272 + t749) + t269)
+			   + (x1 * (px1 * t748 + t262) + px1 * t653) * y1));
+
+    trans->m[1][0] = (x0 * (py0 * (x1 * (t516 + py1 * (y2 * (t631 + t806) + t142 * t933) + t372
+			      + t371)
+			+ y1 * (px1 * t929 + x1 * (t514 + t918 + t157 + t364 + t363) + t908
+				+ t893 + t725 + t892) + t500 + t151
+			+ t1 * (px1 * t855 + t845 + t325 + t125 + t835) + t146)
+		 + px0 * (y1 * (x1 * t932 + py1 * t904 + t872 + t871 + t870)
+			  + py1 * x1 * t515 + py1 * t638
+			  + t1 * (py1 * t844 + t833 + t832 + t831))
+		 + y1 * (x1 * t931 + t863 + py1 * (t247 + t810 + t244)) + py1 * t374
+		 + py1 * x1 * t373 + t1 * (t888 + py1 * (t54 + t282 + t841)))
+	   + y0 * (px0 * (x1 * (py1 * t929 + t862 + t861 + t860)
+			  + t927 + py1 * t29 * t182 + t858
+			  + y1 * (py1 * t896 + x1 * (py1 * t914 + t833 + t832 + t831)) + t857
+			  + t856 + (t926 + t817 + t823 + t822) * t1)
+		   + x0 * (px0 * (x1 * t925 + py1 * t878 + t862
+				  + y1
+				  * (py1 * t902 + py1 * x1 * t916 + t853 + t852
+				     + t851) + t861 + t860
+				  + py1 * t693 * t1)
+			   + py0 * (x1 * (t921 + t918 + t156 + t346 + t345)
+				    + t883 + px1 * t875 + t247
+				    + y1 * (x1 * (t917 + py1 * (t328 + t88 + (t915 + t718) * y2)
+						  + t153 + t55 + t339)
+					    + t899 + px1 * t914 + t322 + t143 + t245 + t887) + t772
+				    + t244 + (t647 + t20 + t19 + t18) * t1) + x1 * t913 + t891
+			   + py1 * (t202 + t780 + t196)
+			   + y1 * (py1 * (t200 + t125 + t898) + t834 + py1 * x1 * t340)
+			   + py1 * t338 * t1)
+		   + py0 * (x1 * (t908 + px1 * t904 + t247 + t810 + t244)
+			    + t903 + py1 * t29 * t179 + t722
+			    + y1 * (x1 * (px1 * t902 + t899 + t200 + t125 + t898)
+				    + t897 + py1 * (x2 * t894 + x2 * (t86 + t628) * y2) + t267
+				    + t266) + t237 + t721 + (t593 + t389 + t95 + t747) * t1)
+		   + x1 * (py1 * (t893 + t725 + t892) + t891) + t890 + t889 + py1 * t29 * t274
+		   + y1 * (x1 * (t888 + py1 * (t322 + t143 + t245 + t887)) + py1 * t268)
+		   + (t885 + t884) * t1)
+	   + py0 * (y1 * (x1 * (t883 + px1 * t878 + t202 + t780 + t196)
+			  + t641 + t731 + t189 + t730)
+		    + x1 * (t639 + t230 + t229) + py1 * t29 * t169
+		    + t1 * (t840 + px1 * t850 + t13 + t752 + t825))
+	   + px0 * (y1 * (x1 * (py1 * t875 + t872 + t871 + t870)
+			  + py1 * t869 + t866 + t865 + t864)
+		    + py1 * x1 * t499 + py1 * t29 * t174
+		    + t1 * (py1 * t838 + t829 + t828 + t827))
+	   + y1 * (x1 * (t863 + py1 * (t247 + t772 + t244)) + t859 + py1 * t723)
+	   + py0 * t161 * t183 + py1 * t29 * t160 + py1 * x1 * t152
+	   + t14 * (px0 * (x1 * (py1 * t855 + t853 + t852 + t851)
+			   + py1 * t850 + py1 * t29 * t109 + t820 + t819 + t818
+			   + (py1 * t592 + t826 + py1 * x1 * t646 + t847 + t846) * y1)
+		    + py0 * (x1 * (t845 + px1 * t844 + t54 + t282 + t841)
+			     + t840 + px1 * t838 + py1 * t29 * t76 + t46 + t659 + t816
+			     + (t421 + t431 + x1 * (t836 + t26 + t94 + t96) + t27 + t661) * y1)
+		    + x1 * (py1 * (t325 + t125 + t835) + t834) + t830
+		    + py1 * (t13 + t752 + t825) + py1 * t29 * t36
+		    + (t824 + py1 * t748 + py1 * x1 * t21) * y1)
+	   + t1 * (t821 + py1 * (t46 + t659 + t816)))
+	;
+
+    trans->m[1][1] = (t161 * (px0 * (x1 * (py1 * t603 + t911 + t910 + t909)
+			  + py1 * t980 + t872
+			  + y1 * (py1 * t944 + py1 * x1 * t109 + t833 + t832 + t831) + t871
+			  + t870 + py1 * t646 * t1)
+		   + py0 * (x1 * (t514 + t157 + t364 + t363)
+			    + px1 * t974 + t202
+			    + y1 * (x1 * (t917 + t153 + t55 + t339)
+				    + t943 + px1 * t955 + t143 + t124 + t938) + t978 + t725
+			    + t196 + py1 * t66 * t1) + x1 * (t930 + py1 * t775) + t863
+		   + py1 * (t773 + t283 + t771)
+		   + y1 * (py1 * (t58 + t144 + t942) + t888 + py1 * x1 * t36)
+		   + py1 * t21 * t1)
+	   + x0 * (py0 * (t29 * (t604 + t149 + t346 + t774)
+			  + t984 + t190
+			  + y1 * (x1 * (px1 * t972 + t970 + t313 + t144 + t960)
+				  + px1 * t958 + t950 + t787 + t12 + t939)
+			  + x1 * (t607 + t983 + py1 * y2 * (x2 * (t328 + t88) + t293 + t327))
+			  + t189 + t188 + py1 * t650 * t1)
+		   + px0 * (t29 * t925 + t982 + t866
+			    + y1
+			    * (x1 * (py1 * t969 + t853 + t852 + t851)
+			       + py1 * t947 + t829 + t828 + t827)
+			    + py1 * x1 * y2 * t334 + t865 + t864 + py1 * t651 * t1)
+		   + t29 * t913 + t859 + t981
+		   + y1 * (x1 * (t834 + py1 * (t633 + t54 + t272 + t968))
+			   + py1 * (t946 + t46 + t945) + t821) + py1 * x1 * y2 * t606
+		   + py1 * t737 * t1)
+	   + py0 * (x1 * (t573 + t238 + t237 + t236)
+		    + t29 * (px1 * t980 + t773 + t978 + t283 + t771)
+		    + y1 * (x1 * (t965 + px1 * t971 + t765 + t659 + t959) + t550 + t670 + t669)
+		    + py1 * t537 * t1)
+	   + px0 * (x1 * (py1 * t977 + t858 + t857 + t856)
+		    + t29 * (py1 * t974 + t862 + t861 + t860)
+		    + y1 * (x1 * (py1 * t966 + t820 + t819 + t818) + t973) + py1 * t539 * t1)
+	   + x1 * (t890 + py1 * t191) + t29 * (t891 + py1 * (t202 + t725 + t196))
+	   + y0 * (x0 * (px0 * (x1 * (py1 * t972 + t833 + t832 + t831)
+				+ py1 * t971 + py1 * t29 * t916 + t829 + t828 + t827
+				+ (t926 + t817 + py1 * x1 * t693 + t823 + t822) * y1)
+			 + py0 * (x1 * (t970 + px1 * t969 + t633 + t54 + t272 + t968)
+				  + px1 * t966 + t965 + t29 * (t163 + t35 + t34 + t33) + t753 + t12
+				  + t935
+				  + (t962 + t941 + t43
+				     + x1
+				     * (t836 + py1 * (t86 + t628 + (t961 + t291) * x2)
+					+ t26 + t94 + t96) + t95 + t93)
+				  * y1) + x1 * (py1 * (t313 + t144 + t960) + t888) + t821
+			 + py1 * (t765 + t659 + t959) + py1 * t29 * t340
+			 + (t885 + t884 + py1 * x1 * t338) * y1)
+		   + px0 * (x1 * (py1 * t958 + t820 + t819 + t818)
+			    + t956 + t29 * (py1 * t955 + t853 + t852 + t851)
+			    + (py1 * t415 + x1 * t953) * y1)
+		   + py0 * (x1 * (t950 + px1 * t947 + t946 + t46 + t945)
+			    + t29 * (px1 * t944 + t943 + t58 + t144 + t942) + t406 + t392 + t391
+			    + (x1 * (t421 + t941 + t431 + t27 + t661)
+			       + t416 + py1 * (x2 * (t554 + t796) + t5 * t940) + t665 + t664)
+			    * y1) + x1 * (py1 * (t787 + t12 + t939) + t830)
+		   + t29 * (t834 + py1 * (t143 + t124 + t938)) + t937
+		   + (x1 * t936 + py1 * t666) * y1)
+	   + y1 * (x1 * (py1 * (t753 + t12 + t935) + t830) + t934) + py1 * t653 * t1
+	   + py0 * t652 * t14)
+	;
+
+    trans->m[1][2] = (y0 * (x0 * (px0 * (x1 * (py1 * t1016 + t872 + t871 + t870)
+				     + py1 * t494 + t29 * t932
+				     + y1
+				     * (py1 * t1025 + py1 * x1 * t599 + t820 + t819
+					+ t818) + t953 * t1)
+			      + py0 * (x1 * (t1015 + px1 * t1009 + t727 + t326 + t724)
+				       + t29 * (t1028 + t148 + t364 + t744) + t491
+				       + y1
+				       * (x1
+					  * (t600 + py1 * (x2 * t933 + t940 * y2) + t503
+					     + t800)
+					  + px1 * t1020 + t1026 + t676 + t45 + t1017)
+				       + (t1027 + t10 + t27 + t25) * t1)
+			      + x1 * (py1 * (t1014 + t283 + t1012) + t863) + t29 * t931
+			      + py1 * t802
+			      + y1 * (py1 * x1 * t801 + py1 * (t1024 + t13 + t1023) + t830)
+			      + t936 * t1)
+			+ py0 * (t29 * (t1007 + px1 * t1002 + t727 + t201 + t724)
+				 + x1 * (t1011 + t641 + t731 + t189 + t730)
+				 + y1 * (x1 * (t1026 + px1 * t1025 + t1024 + t13 + t1023)
+					 + t550
+					 + py1
+					 * (t5 * (t293 + t327) + x2 * (t1021 + t697) * y2)
+					 + t670 + t669) + (t540 + t648 + t645) * t1)
+			+ px0 * (x1 * (t982 + t866 + t865 + t864)
+				 + t29 * (py1 * t1006 + t872 + t871 + t870)
+				 + y1 * (x1 * (py1 * t1020 + t829 + t828 + t827) + t956)
+				 + py1 * t415 * t1) + x1 * (t859 + t981)
+			+ t29 * (t863 + py1 * (t773 + t246 + t771))
+			+ y1 * (x1 * (py1 * (t676 + t45 + t1017) + t821) + t937)
+			+ py1 * t666 * t1)
+		  + x0 * (py0 * (y1 * (x1 * (px1 * t1016 + t1015 + t1014 + t283 + t1012)
+				       + t1011 + t573 + t238 + t237 + t236)
+				 + x1 * (t639
+					 + py1 * (x2 * y2 * (t1010 + t300) + t142 * t894)
+					 + t230 + t229) + t29 * (t175 + t159 + t155)
+				 + t1 * (px1 * t1000 + t994 + t660 + t45 + t986))
+			  + px0 * (y1 * (x1 * (py1 * t1009 + t862 + t861 + t860)
+					 + t927 + t858 + t857 + t856)
+				   + py1 * t29 * t515 + py1 * x1 * t638
+				   + t1 * (py1 * t993 + t820 + t819 + t818))
+			  + y1 * (x1 * (t891 + py1 * (t727 + t326 + t724)) + t890 + t889)
+			  + py1 * x1 * t374 + py1 * t29 * t373
+			  + t1 * (t830 + py1 * (t705 + t752 + t990)))
+		  + t161 * (py0 * (x1 * (t516 + t372 + t371)
+				   + y1
+				   * (x1 * (t921 + t156 + t346 + t345)
+				      + t1007 + px1 * t1006 + t773 + t246 + t771) + t381
+				   + t1 * (t989 + px1 * t997 + t54 + t245 + t985))
+			    + px0 * (y1
+				     * (x1 * (py1 * t1004 + t924 + t923 + t922)
+					+ py1 * t1002 + t862 + t861 + t860)
+				     + py1 * t499 + py1 * x1 * t174
+				     + t1 * (py1 * t988 + t853 + t852 + t851))
+			    + y1 * (x1 * (t912 + py1 * t746) + t891
+				    + py1 * (t727 + t201 + t724))
+			    + py1 * x1 * t160 + py1 * t152
+			    + t1 * (t834 + py1 * (t143 + t272 + t987))) + py0 * t770
+		  + t14 * (px0 * (x1 * (py1 * t1000 + t829 + t828 + t827)
+				  + t973 + t29 * (py1 * t997 + t833 + t832 + t831)
+				  + (py1 * t539 + x1 * (py1 * t996 + t817 + t823 + t822))
+				  * y1)
+			   + py0 * (x1 * (t994 + px1 * t993 + t705 + t752 + t990)
+				    + t29 * (t989 + px1 * t988 + t143 + t272 + t987) + t398
+				    + (t416 + x1 * (t962 + t43 + t95 + t93) + t665 + t664)
+				    * y1) + x1 * (py1 * (t660 + t45 + t986) + t821)
+			   + t29 * (t888 + py1 * (t54 + t245 + t985)) + t934
+			   + (x1 * (t885 + py1 * t28) + py1 * t653) * y1));
+
+    trans->m[2][0] = (x0 * (px0 * (y1 * (x1 * t617 + t586 + t877 + t1008 + t876)
+			+ x1 * t515 + t637 + t1 * (t475 + t136 + t620 + t1036) + t636)
+		 + py0 * (y1 * (t579 + x1 * t627 + t251 + t811 + t248)
+			  + x1 * t379 + t233 + t1 * (t455 + t87 + t298 + t1035) + t232)
+		 + x1 * (t516 + t383 + y2 * (t804 + t1013) + t142 * t595)
+		 + y1 * (px1 * (t1047 + t979 + t1046)
+			 + x1 * (t921 + t368 + t157 + t156 + t1048)
+			 + py1 * (t1044 + t728 + t1043) + t505 + t607 + t502) + t500 + t381
+		 + t1 * (px1 * (t611 + t135 + t1038) + py1 * (t330 + t128 + t1037) + t423
+			 + t125 + t54 + t1034))
+	   + y0 * (x0 * (py0 * (x1 * t1050 + t533 + t226
+				+ y1 * (t529 + t224 + x1 * t349 + t128 + t1040)
+				+ t781 + t220 + t99 * t1)
+			 + px0 * (t528 + x1 * t1049 + t874
+				  + y1 * (t517 + x1 * t916 + t520 + t135 + t1041) + t1005
+				  + t873 + t693 * t1)
+			 + x1 * (t514 + t382 + t157 + t156 + t1048)
+			 + px1 * (t877 + t1001 + t876) + py1 * (t251 + t776 + t248) + t608
+			 + y1 * (x1 * (t917 + t358 + t56 + t323 + (t716 + t967) * y2)
+				 + px1 * (t118 + t609 + t471 + t1045)
+				 + py1 * (t204 + t327 + t249 + t1042) + t503 + t144 + t143
+				 + t1039) + t504 + t605 + (t647 + t357) * t1)
+		   + px0 * (x1 * (t528 + t1047 + t979 + t1046)
+			    + t643 + t29 * t182 + t571
+			    + y1 * (x1 * (t475 + t118 + t609 + t471 + t1045) + t701 + t895)
+			    + t570 + t569 + (t468 + t441 + t466 + t951) * t1)
+		   + py0 * (x1 * (t533 + t1044 + t728 + t1043)
+			    + t642 + t29 * t179 + t566
+			    + y1 * (x1 * (t455 + t204 + t327 + t249 + t1042) + t289 + t288)
+			    + t214 + t565 + (t436 + t450 + t69 + t677) * t1)
+		   + x1 * (px1 * (t877 + t1008 + t876) + py1 * (t251 + t811 + t248) + t505
+			   + t607 + t502) + t984 + t743
+		   + t29 * t274
+		   + y1 * (x1 * (px1 * (t520 + t135 + t1041) + py1 * (t224 + t128 + t1040)
+				 + t503 + t144 + t143 + t1039)
+			   + t897 + t303 + x2 * t634 + x2 * (t690 + t52) * y2)
+		   + (t1027 + t808) * t1)
+	   + py0 * (y1 * (x1 * (t579 + t251 + t776 + t248) + t581 + t735 + t240 + t734)
+		    + t29 * t169 + x1 * t209 + t1 * (t461 + t80 + t683 + t1032))
+	   + px0 * (y1 * (x1 * (t586 + t877 + t1001 + t876) + t591 + t976 + t588 + t975)
+		    + x1 * t499 + t29 * t174 + t1 * (t484 + t849 + t957 + t1033))
+	   + y1 * (x1 * (px1 * (t874 + t1005 + t873) + py1 * (t226 + t781 + t220) + t608
+			 + t504 + t605)
+		   + t573 + t568) + t161 * t183 + x1 * (t639 + t235) + t29 * t160
+	   + t14 * (px0 * (x1 * (t517 + t611 + t135 + t1038)
+			   + t615 + t29 * t109 + t560 + t443 + t1031
+			   + (t399 + t613 + x1 * t646 + t418 + t995) * y1)
+		    + py0 * (x1 * (t529 + t330 + t128 + t1037)
+			     + t625 + t29 * t76 + t112 + t438 + t1030
+			     + (t623 + t393 + t102 + x1 * t66 + t754) * y1)
+		    + x1 * (px1 * (t136 + t620 + t1036) + py1 * (t87 + t298 + t1035) + t423
+			    + t125 + t54 + t1034)
+		    + px1 * (t849 + t957 + t1033) + py1 * (t80 + t683 + t1032) + t434
+		    + t29 * t36 + t1029 + (t962 + t711 + x1 * (t836 + t100)) * y1)
+	   + t1 * (px1 * (t560 + t443 + t1031) + py1 * (t112 + t438 + t1030) + t434
+		   + t1029))
+	;
+
+    trans->m[2][1] = (t161 * (px0 * (x1 * (t616 + t180 + t512 + t1003)
+			  + t586 + t510 + y1 * (t475 + t854 + x1 * t109 + t117 + t1062)
+			  + t621 + t508 + t646 * t1)
+		   + py0 * (t579 + x1 * (t177 + t626 + t354 + t778) + t507
+			    + y1 * (t455 + x1 * t76 + t91 + t203 + t1061) + t299
+			    + t506 + t66 * t1) + x1 * (t921 + t368)
+		   + px1 * (t874 + t979 + t873) + py1 * (t226 + t728 + t220)
+		   + y1 * (x1 * (t917 + t358) + px1 * (t118 + t843 + t1063)
+			   + py1 * (t204 + t127 + t1065) + t423 + t144
+			   + t143 + t1060) + t504 + t21 * t1)
+	   + x0 * (py0 * (t29 * t1050 + t581 + t241
+			  + y1
+			  * (x1 * (t529 + t87 + t577 + t292 + t1070)
+			     + t461 + t1057 + t80 + t1056) + x1 * y2 * t741 + t240
+			  + t239 + t650 * t1)
+		   + px0 * (t591 + t29 * t1049 + t589
+			    + y1 * (x1 * (t517 + t136 + t306 + t428 + t1071)
+				    + t484 + t1059 + t849 + t1058) + x1 * y2 * t334 + t588
+			    + t587 + t651 * t1) + t29 * (t1028 + t809) + t903 + t321
+		   + y1 * (x1 * (px1 * (t901 + t117 + t1072) + py1 * (t317 + t203 + t1073)
+				 + t503 + t125 + t54 + t1069)
+			   + px1 * (t1019 + t837 + t1064) + py1 * (t789 + t111 + t1066) + t267
+			   + t433 + t1055)
+		   + x1 * (y2 * (x2 * (t56 + t323) + t322 + t273) + t983 + t742) + t737 * t1)
+	   + py0 * (x1 * (t642 + t215 + t214 + t213) + t29 * (t533 + t226 + t728 + t220)
+		    + y1
+		    * (x1 * (t625 + t758 + t111 + t1067)
+		       + t396 + t395) + t537 * t1)
+	   + px0 * (x1 * (t643 + t868 + t570 + t867) + t29 * (t528 + t874 + t979 + t873)
+		    + y1
+		    * (x1 * (t615 + t999 + t837 + t1068)
+		       + t404 + t402) + t539 * t1)
+	   + x1 * (t641 + t243)
+	   + t29 * (px1 * (t510 + t621 + t508) + py1 * (t507 + t299 + t506) + t504)
+	   + y0 * (x0 * (py0 * (x1 * (t455 + t317 + t203 + t1073)
+				+ t461 + t29 * t349 + t767 + t683 + t1052
+				+ (t436 + t450 + x1 * t99 + t69 + t677) * y1)
+			 + px0 * (x1 * (t475 + t901 + t117 + t1072)
+				  + t484 + t29 * t916 + t992 + t957 + t1053
+				  + (t468 + t441 + t466 + x1 * t693 + t951) * y1)
+			 + x1 * (px1 * (t136 + t306 + t428 + t1071)
+				 + py1 * (t87 + t577 + t292 + t1070) + t503 + t125 + t54 + t1069)
+			 + px1 * (t999 + t837 + t1068) + py1 * (t758 + t111 + t1067)
+			 + t29 * (t163 + t162) + t434 + t553 + t1051
+			 + (t421 + t104 + t431 + t43
+			    + x1 * (t836 + t100 + t690 + t52 + (t270 + t886) * x2)
+			    + t1054)
+			 * y1)
+		   + py0 * (x1 * (t625 + t789 + t111 + t1066)
+			    + t29 * (t529 + t204 + t127 + t1065) + t544 + t543
+			    + (x1 * t624 + t410 + t409) * y1)
+		   + px0 * (x1 * (t615 + t1019 + t837 + t1064)
+			    + t29 * (t517 + t118 + t843 + t1063) + t548 + t547
+			    + (x1 * t614 + t414 + t413) * y1)
+		   + t29 * (px1 * (t854 + t117 + t1062) + py1 * (t91 + t203 + t1061) + t423
+			    + t144 + t143 + t1060)
+		   + x1 * (px1 * (t1059 + t849 + t1058) + py1 * (t1057 + t80 + t1056) + t267
+			   + t433 + t1055) + t406 + t398
+		   + (t416 + x1 * (t962 + t711 + t431 + t43 + t1054) + t412
+		      + x2 * (t793 + t1022) + t5 * t594)
+		   * y1)
+	   + y1 * (x1 * (px1 * (t992 + t957 + t1053) + py1 * (t767 + t683 + t1052) + t434
+			 + t553 + t1051)
+		   + t550 + t546) + t653 * t1 + t652 * t14)
+	;
+    trans->m[2][2] = t644;
+}
+
+static void
+print_trans (const char *header, struct pixman_f_transform *trans)
+{
+    int i, j;
+    double max;
+
+    max = 0;
+
+
+    for (i = 0; i < 3; ++i)
+    {
+	for (j = 0; j < 3; ++j)
+	{
+	    double a = fabs (trans->m[i][j]);
+
+	    if (a > max)
+		max = a;
+	}
+    }
+
+    if (max == 0.0)
+	max = 1.0;
+
+    for (i = 0; i < 3; ++i)
+    {
+	printf ("{ ");
+	for (j = 0; j < 3; ++j)
+	{
+	    printf ("D2F (%.5f)%s", 16384 * (trans->m[i][j] / max), j == 2 ? "" : ", ");
+	}
+
+	printf ("},\n");
+    }
+}
+
+int
+main ()
+{
+    struct pixman_f_transform t;
+
+#if 0
+		  325, 200,
+		  450, 335,
+		  -50, 335,
+
+		  0, 0,
+		  400, 0,
+		  400, 400,
+		  0, 400,
+
+		  &t);
+#endif
+	1, 0,
+	1, 2,
+	2, 2,
+	2, 0,
+
+	1, 0,
+	1, 112,
+	2, 2,
+	2, 0,
+
+	&t);
+
+    print_trans ("0->0", &t);
+
+    return 0;
+}
--
1.6.0.6

```