diff options
| author | Søren Sandmann Pedersen <ssp@redhat.com> | 2011-10-12 04:48:33 -0400 |
|---|---|---|
| committer | Søren Sandmann Pedersen <ssp@redhat.com> | 2012-03-24 16:29:27 -0400 |
| commit | 41863fbabb6dd08871aed47beba5d08ae2ae3cf6 (patch) | |
| tree | 91ad9144add216cef41be60cf76ebf489af26694 /demos | |
| parent | cf0d0d63645bcb6425a1e2d7b5d9f1e26e205247 (diff) | |
demos: Add quad2quad program
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/
Diffstat (limited to 'demos')
| -rw-r--r-- | demos/quad2quad.c | 2183 |
1 files changed, 2183 insertions, 0 deletions
diff --git a/demos/quad2quad.c b/demos/quad2quad.c new file mode 100644 index 0000000..66b838f --- /dev/null +++ b/demos/quad2quad.c @@ -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 +quad_to_quad (double x0, double y0, + double x1, double y1, + double x2, double y2, + double x3, double y3, + + double px0, double py0, + double px1, double py1, + double px2, double py2, + double px3, double py3, + + struct pixman_f_transform *trans) +{ + 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; + + printf ("%s\n", header); + + 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 + quad_to_quad (75, 200, + 325, 200, + 450, 335, + -50, 335, + + 0, 0, + 400, 0, + 400, 400, + 0, 400, + + &t); +#endif + quad_to_quad ( + 1, 0, + 1, 2, + 2, 2, + 2, 0, + + 1, 0, + 1, 112, + 2, 2, + 2, 0, + + &t); + + print_trans ("0->0", &t); + + return 0; +} |