123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347 |
- t2 = dax*(1.0/2.0);
- t3 = dax_b*(1.0/2.0);
- t4 = t2-t3;
- t5 = day*(1.0/2.0);
- t6 = day_b*(1.0/2.0);
- t7 = t5-t6;
- t8 = daz*(1.0/2.0);
- t9 = daz_b*(1.0/2.0);
- t10 = t8-t9;
- t11 = q2*t4*(1.0/2.0);
- t12 = q1*t7*(1.0/2.0);
- t13 = q0*t10*(1.0/2.0);
- t14 = q2*(1.0/2.0);
- t15 = q3*t4*(1.0/2.0);
- t16 = q1*t10*(1.0/2.0);
- t17 = q1*(1.0/2.0);
- t18 = q0*t4*(1.0/2.0);
- t19 = q3*t7*(1.0/2.0);
- t20 = q0*(1.0/2.0);
- t21 = q2*t7*(1.0/2.0);
- t22 = q3*t10*(1.0/2.0);
- t23 = q0*t7*(1.0/2.0);
- t24 = q3*(1.0/2.0);
- t25 = q1*t4*(1.0/2.0);
- t26 = q2*t10*(1.0/2.0);
- t27 = t11+t12+t13-t24;
- t28 = t14+t15+t16-t23;
- t29 = q2*t28*2.0;
- t30 = t17+t18+t19-t26;
- t31 = q1*t30*2.0;
- t32 = t20+t21+t22-t25;
- t33 = q0*t32*2.0;
- t40 = q3*t27*2.0;
- t34 = t29+t31+t33-t40;
- t35 = q0*q0;
- t36 = q1*q1;
- t37 = q2*q2;
- t38 = q3*q3;
- t39 = t35+t36+t37+t38;
- t41 = t11+t12-t13+t24;
- t42 = t14-t15+t16+t23;
- t43 = q1*t42*2.0;
- t44 = -t17+t18+t19+t26;
- t45 = q2*t44*2.0;
- t46 = t20-t21+t22+t25;
- t47 = q3*t46*2.0;
- t57 = q0*t41*2.0;
- t48 = t43+t45+t47-t57;
- t49 = -t14+t15+t16+t23;
- t50 = q0*t49*2.0;
- t51 = t11-t12+t13+t24;
- t52 = t20+t21-t22+t25;
- t53 = q2*t52*2.0;
- t54 = t17-t18+t19+t26;
- t55 = q3*t54*2.0;
- t58 = q1*t51*2.0;
- t56 = t50+t53+t55-t58;
- t59 = OP[0][0]*t34;
- t60 = OP[1][0]*t48;
- t66 = OP[6][0]*t39;
- t67 = OP[2][0]*t56;
- t61 = t59+t60-t66-t67;
- t62 = OP[0][1]*t34;
- t63 = OP[1][1]*t48;
- t64 = OP[0][2]*t34;
- t65 = OP[1][2]*t48;
- t71 = OP[6][1]*t39;
- t72 = OP[2][1]*t56;
- t68 = t62+t63-t71-t72;
- t79 = OP[6][2]*t39;
- t80 = OP[2][2]*t56;
- t69 = t64+t65-t79-t80;
- t70 = t35+t36-t37-t38;
- t73 = q0*q2*2.0;
- t74 = q1*q3*2.0;
- t75 = t73+t74;
- t76 = q0*q3*2.0;
- t78 = q1*q2*2.0;
- t77 = t76-t78;
- t81 = t35-t36+t37-t38;
- t82 = q0*q1*2.0;
- t86 = q2*q3*2.0;
- t83 = t82-t86;
- t84 = t76+t78;
- t85 = t35-t36-t37+t38;
- t87 = t82+t86;
- t88 = t73-t74;
- t89 = OP[0][6]*t34;
- t90 = OP[1][6]*t48;
- t265 = OP[6][6]*t39;
- t91 = t89+t90-t265-OP[2][6]*t56;
- t92 = OP[0][7]*t34;
- t93 = OP[1][7]*t48;
- t266 = OP[6][7]*t39;
- t94 = t92+t93-t266-OP[2][7]*t56;
- t95 = OP[0][8]*t34;
- t96 = OP[1][8]*t48;
- t267 = OP[6][8]*t39;
- t97 = t95+t96-t267-OP[2][8]*t56;
- t98 = q0*t27*2.0;
- t99 = q1*t28*2.0;
- t100 = q3*t32*2.0;
- t110 = q2*t30*2.0;
- t101 = t98+t99+t100-t110;
- t102 = q0*t46*2.0;
- t103 = q2*t42*2.0;
- t104 = q3*t41*2.0;
- t111 = q1*t44*2.0;
- t105 = t102+t103+t104-t111;
- t106 = q1*t52*2.0;
- t107 = q2*t51*2.0;
- t108 = q3*t49*2.0;
- t112 = q0*t54*2.0;
- t109 = t106+t107+t108-t112;
- t113 = OP[7][0]*t39;
- t114 = OP[0][0]*t101;
- t123 = OP[1][0]*t105;
- t124 = OP[2][0]*t109;
- t115 = t113+t114-t123-t124;
- t116 = OP[7][1]*t39;
- t117 = OP[0][1]*t101;
- t129 = OP[1][1]*t105;
- t130 = OP[2][1]*t109;
- t118 = t116+t117-t129-t130;
- t119 = OP[7][2]*t39;
- t120 = OP[0][2]*t101;
- t135 = OP[1][2]*t105;
- t136 = OP[2][2]*t109;
- t121 = t119+t120-t135-t136;
- t122 = t39*t39;
- t125 = q1*t27*2.0;
- t126 = q2*t32*2.0;
- t127 = q3*t30*2.0;
- t170 = q0*t28*2.0;
- t128 = t125+t126+t127-t170;
- t131 = q0*t44*2.0;
- t132 = q1*t46*2.0;
- t133 = q2*t41*2.0;
- t171 = q3*t42*2.0;
- t134 = t131+t132+t133-t171;
- t137 = q0*t52*2.0;
- t138 = q1*t54*2.0;
- t139 = q3*t51*2.0;
- t172 = q2*t49*2.0;
- t140 = t137+t138+t139-t172;
- t141 = dvy*t70;
- t142 = dvx*t77;
- t143 = t141+t142;
- t144 = dvx*t75;
- t145 = dvy*t75;
- t146 = dvz*t77;
- t147 = t145+t146;
- t148 = dvx*t81;
- t188 = dvy*t84;
- t149 = t148-t188;
- t150 = dvz*t81;
- t151 = dvy*t83;
- t152 = t150+t151;
- t153 = dvx*t83;
- t154 = dvz*t84;
- t155 = t153+t154;
- t156 = dvx*t85;
- t157 = dvz*t88;
- t158 = t156+t157;
- t159 = dvy*t85;
- t189 = dvz*t87;
- t160 = t159-t189;
- t161 = dvx*t87;
- t162 = dvy*t88;
- t163 = t161+t162;
- t164 = OP[7][6]*t39;
- t165 = OP[0][6]*t101;
- t166 = OP[7][7]*t39;
- t167 = OP[0][7]*t101;
- t168 = OP[7][8]*t39;
- t169 = OP[0][8]*t101;
- t173 = OP[8][0]*t39;
- t174 = OP[1][0]*t134;
- t182 = OP[0][0]*t128;
- t183 = OP[2][0]*t140;
- t175 = t173+t174-t182-t183;
- t176 = OP[8][1]*t39;
- t177 = OP[1][1]*t134;
- t184 = OP[0][1]*t128;
- t185 = OP[2][1]*t140;
- t178 = t176+t177-t184-t185;
- t179 = OP[8][2]*t39;
- t180 = OP[1][2]*t134;
- t186 = OP[0][2]*t128;
- t187 = OP[2][2]*t140;
- t181 = t179+t180-t186-t187;
- t190 = OP[8][6]*t39;
- t191 = OP[1][6]*t134;
- t192 = OP[8][7]*t39;
- t193 = OP[1][7]*t134;
- t194 = OP[8][8]*t39;
- t195 = OP[1][8]*t134;
- t197 = dvz*t70;
- t196 = t144-t197;
- t198 = OP[0][0]*t147;
- t204 = OP[2][0]*t143;
- t205 = OP[1][0]*t196;
- t199 = OP[3][0]+t198-t204-t205;
- t200 = OP[0][1]*t147;
- t206 = OP[2][1]*t143;
- t207 = OP[1][1]*t196;
- t201 = OP[3][1]+t200-t206-t207;
- t202 = OP[0][2]*t147;
- t208 = OP[2][2]*t143;
- t209 = OP[1][2]*t196;
- t203 = OP[3][2]+t202-t208-t209;
- t210 = -t144+t197;
- t211 = OP[1][0]*t210;
- t212 = OP[3][0]+t198-t204+t211;
- t213 = OP[1][1]*t210;
- t214 = OP[3][1]+t200-t206+t213;
- t215 = OP[1][2]*t210;
- t216 = OP[3][2]+t202-t208+t215;
- t217 = OP[0][6]*t147;
- t218 = OP[0][7]*t147;
- t219 = OP[0][8]*t147;
- t220 = OP[1][0]*t155;
- t221 = OP[2][0]*t149;
- t229 = OP[0][0]*t152;
- t222 = OP[4][0]+t220+t221-t229;
- t223 = OP[1][1]*t155;
- t224 = OP[2][1]*t149;
- t230 = OP[0][1]*t152;
- t225 = OP[4][1]+t223+t224-t230;
- t226 = OP[1][2]*t155;
- t227 = OP[2][2]*t149;
- t231 = OP[0][2]*t152;
- t228 = OP[4][2]+t226+t227-t231;
- t232 = dvxNoise*t70*t84;
- t233 = OP[1][6]*t155;
- t234 = OP[2][6]*t149;
- t235 = OP[4][6]+t233+t234-OP[0][6]*t152;
- t236 = OP[1][7]*t155;
- t237 = OP[2][7]*t149;
- t238 = OP[4][7]+t236+t237-OP[0][7]*t152;
- t239 = OP[1][8]*t155;
- t240 = OP[2][8]*t149;
- t241 = OP[4][8]+t239+t240-OP[0][8]*t152;
- t242 = OP[2][0]*t163;
- t243 = OP[0][0]*t160;
- t251 = OP[1][0]*t158;
- t244 = OP[5][0]+t242+t243-t251;
- t245 = OP[2][1]*t163;
- t246 = OP[0][1]*t160;
- t252 = OP[1][1]*t158;
- t247 = OP[5][1]+t245+t246-t252;
- t248 = OP[2][2]*t163;
- t249 = OP[0][2]*t160;
- t253 = OP[1][2]*t158;
- t250 = OP[5][2]+t248+t249-t253;
- t254 = dvzNoise*t75*t85;
- t255 = dvyNoise*t81*t87;
- t256 = OP[2][6]*t163;
- t257 = OP[0][6]*t160;
- t258 = OP[5][6]+t256+t257-OP[1][6]*t158;
- t259 = OP[2][7]*t163;
- t260 = OP[0][7]*t160;
- t261 = OP[5][7]+t259+t260-OP[1][7]*t158;
- t262 = OP[2][8]*t163;
- t263 = OP[0][8]*t160;
- t264 = OP[5][8]+t262+t263-OP[1][8]*t158;
- A0[0][0] = daxNoise*t122+t34*t61+t48*t68-t56*t69-t39*t91;
- A0[0][0] = -t39*t94-t61*t101+t68*t105+t69*t109;
- A0[0][1] = -t39*t97-t68*t134+t69*t140+t128*(t59+t60-t66-t67);
- A0[0][2] = OP[0][3]*t34+OP[1][3]*t48-OP[2][3]*t56-OP[6][3]*t39-t69*t143+t147*(t59+t60-t66-t67)-t196*(t62+t63-t71-t72);
- A0[0][3] = OP[0][4]*t34+OP[1][4]*t48-OP[2][4]*t56-OP[6][4]*t39-t61*t152+t69*t149+t155*(t62+t63-t71-t72);
- A0[0][4] = OP[0][5]*t34+OP[1][5]*t48-OP[2][5]*t56-OP[6][5]*t39+t61*t160-t68*t158+t163*(t64+t65-t79-t80);
- A0[0][5] = t91;
- A0[0][6] = t94;
- A0[0][7] = t97;
- A0[0][0] = -t34*t115-t48*t118+t56*t121+t39*(t164+t165-OP[1][6]*t105-OP[2][6]*t109);
- A0[0][0] = dayNoise*t122+t101*t115-t105*t118-t109*t121+t39*(t166+t167-OP[1][7]*t105-OP[2][7]*t109);
- A0[0][1] = -t115*t128+t118*t134-t121*t140+t39*(t168+t169-OP[1][8]*t105-OP[2][8]*t109);
- A0[0][2] = -OP[0][3]*t101-OP[7][3]*t39+OP[1][3]*t105+OP[2][3]*t109-t115*t147+t121*t143+t118*t196;
- A0[0][3] = -OP[0][4]*t101-OP[7][4]*t39+OP[1][4]*t105+OP[2][4]*t109+t115*t152-t121*t149-t118*t155;
- A0[0][4] = -OP[0][5]*t101-OP[7][5]*t39+OP[1][5]*t105+OP[2][5]*t109-t115*t160+t118*t158-t121*t163;
- A0[0][5] = -t164-t165+OP[1][6]*t105+OP[2][6]*t109;
- A0[0][6] = -t166-t167+OP[1][7]*t105+OP[2][7]*t109;
- A0[0][7] = -t168-t169+OP[1][8]*t105+OP[2][8]*t109;
- A0[1][0] = -t34*t175-t48*t178+t56*t181+t39*(t190+t191-OP[0][6]*t128-OP[2][6]*t140);
- A0[1][0] = t101*t175-t105*t178-t109*t181+t39*(t192+t193-OP[0][7]*t128-OP[2][7]*t140);
- A0[1][1] = dazNoise*t122-t128*t175+t134*t178-t140*t181+t39*(t194+t195-OP[0][8]*t128-OP[2][8]*t140);
- A0[1][2] = -OP[8][3]*t39+OP[0][3]*t128-OP[1][3]*t134+OP[2][3]*t140-t147*t175+t143*t181+t178*t196;
- A0[1][3] = -OP[8][4]*t39+OP[0][4]*t128-OP[1][4]*t134+OP[2][4]*t140+t152*t175-t149*t181-t155*t178;
- A0[1][4] = -OP[8][5]*t39+OP[0][5]*t128-OP[1][5]*t134+OP[2][5]*t140-t160*t175+t158*t178-t163*t181;
- A0[1][5] = -t190-t191+OP[0][6]*t128+OP[2][6]*t140;
- A0[1][6] = -t192-t193+OP[0][7]*t128+OP[2][7]*t140;
- A0[1][7] = -t194-t195+OP[0][8]*t128+OP[2][8]*t140;
- A0[2][0] = -t39*(OP[3][6]+t217-OP[2][6]*t143-OP[1][6]*t196)+t34*t199+t48*t201-t56*t203;
- A0[2][0] = -t39*(OP[3][7]+t218-OP[2][7]*t143-OP[1][7]*t196)-t101*t199+t105*t201+t109*t203;
- A0[2][1] = -t39*(OP[3][8]+t219-OP[2][8]*t143-OP[1][8]*t196)+t128*t199-t134*t201+t140*t203;
- A0[2][2] = OP[3][3]+OP[0][3]*t147-OP[2][3]*t143+OP[1][3]*t210-t143*t203+t147*t212+t210*t214+dvxNoise*(t70*t70)+dvyNoise*(t77*t77)+dvzNoise*(t75*t75);
- A0[2][3] = OP[3][4]+t232+OP[0][4]*t147-OP[2][4]*t143+OP[1][4]*t210-t152*t212+t149*t216+t155*t214-dvyNoise*t77*t81-dvzNoise*t75*t83;
- A0[2][4] = OP[3][5]+t254+OP[0][5]*t147-OP[2][5]*t143+OP[1][5]*t210-t158*t214+t160*t212+t163*t216-dvxNoise*t70*t88-dvyNoise*t77*t87;
- A0[2][5] = OP[3][6]+t217-OP[2][6]*t143+OP[1][6]*t210;
- A0[2][6] = OP[3][7]+t218-OP[2][7]*t143+OP[1][7]*t210;
- A0[2][7] = OP[3][8]+t219-OP[2][8]*t143+OP[1][8]*t210;
- A0[3][0] = t34*t222+t48*t225-t39*t235-t56*t228;
- A0[3][0] = -t39*t238-t101*t222+t105*t225+t109*t228;
- A0[3][1] = -t39*t241+t128*t222-t134*t225+t140*t228;
- A0[3][2] = OP[4][3]+t232-OP[0][3]*t152+OP[1][3]*t155+OP[2][3]*t149+t147*t222-t143*t228+t210*t225-dvyNoise*t77*t81-dvzNoise*t75*t83;
- A0[3][3] = OP[4][4]-OP[0][4]*t152+OP[1][4]*t155+OP[2][4]*t149-t152*t222+t149*t228+t155*t225+dvxNoise*(t84*t84)+dvyNoise*(t81*t81)+dvzNoise*(t83*t83);
- A0[3][4] = OP[4][5]+t255-OP[0][5]*t152+OP[1][5]*t155+OP[2][5]*t149+t160*t222-t158*t225+t163*t228-dvxNoise*t84*t88-dvzNoise*t83*t85;
- A0[3][5] = t235;
- A0[3][6] = t238;
- A0[3][7] = t241;
- A0[4][0] = t34*t244+t48*t247-t39*t258-t56*t250;
- A0[4][0] = -t39*t261-t101*t244+t105*t247+t109*t250;
- A0[4][1] = -t39*t264+t128*t244-t134*t247+t140*t250;
- A0[4][2] = OP[5][3]+t254+OP[0][3]*t160-OP[1][3]*t158+OP[2][3]*t163+t147*t244-t143*t250+t210*t247-dvxNoise*t70*t88-dvyNoise*t77*t87;
- A0[4][3] = OP[5][4]+t255+OP[0][4]*t160-OP[1][4]*t158+OP[2][4]*t163-t152*t244+t149*t250+t155*t247-dvxNoise*t84*t88-dvzNoise*t83*t85;
- A0[4][4] = OP[5][5]+OP[0][5]*t160-OP[1][5]*t158+OP[2][5]*t163+t160*t244-t158*t247+t163*t250+dvxNoise*(t88*t88)+dvyNoise*(t87*t87)+dvzNoise*(t85*t85);
- A0[4][5] = t258;
- A0[4][6] = t261;
- A0[4][7] = t264;
- A0[5][0] = -t265+OP[6][0]*t34+OP[6][1]*t48-OP[6][2]*t56;
- A0[5][0] = -t266-OP[6][0]*t101+OP[6][1]*t105+OP[6][2]*t109;
- A0[5][1] = -t267+OP[6][0]*t128-OP[6][1]*t134+OP[6][2]*t140;
- A0[5][2] = OP[6][3]-OP[6][2]*t143+OP[6][0]*t147+OP[6][1]*t210;
- A0[5][3] = OP[6][4]+OP[6][2]*t149-OP[6][0]*t152+OP[6][1]*t155;
- A0[5][4] = OP[6][5]-OP[6][1]*t158+OP[6][0]*t160+OP[6][2]*t163;
- A0[5][5] = OP[6][6];
- A0[5][6] = OP[6][7];
- A0[5][7] = OP[6][8];
- A0[6][0] = -t164+OP[7][0]*t34+OP[7][1]*t48-OP[7][2]*t56;
- A0[6][0] = -t166-OP[7][0]*t101+OP[7][1]*t105+OP[7][2]*t109;
- A0[6][1] = -t168+OP[7][0]*t128-OP[7][1]*t134+OP[7][2]*t140;
- A0[6][2] = OP[7][3]-OP[7][2]*t143+OP[7][0]*t147+OP[7][1]*t210;
- A0[6][3] = OP[7][4]+OP[7][2]*t149-OP[7][0]*t152+OP[7][1]*t155;
- A0[6][4] = OP[7][5]-OP[7][1]*t158+OP[7][0]*t160+OP[7][2]*t163;
- A0[6][5] = OP[7][6];
- A0[6][6] = OP[7][7];
- A0[6][7] = OP[7][8];
- A0[7][0] = -t190+OP[8][0]*t34+OP[8][1]*t48-OP[8][2]*t56;
- A0[7][0] = -t192-OP[8][0]*t101+OP[8][1]*t105+OP[8][2]*t109;
- A0[7][1] = -t194+OP[8][0]*t128-OP[8][1]*t134+OP[8][2]*t140;
- A0[7][2] = OP[8][3]-OP[8][2]*t143+OP[8][0]*t147+OP[8][1]*t210;
- A0[7][3] = OP[8][4]+OP[8][2]*t149-OP[8][0]*t152+OP[8][1]*t155;
- A0[7][4] = OP[8][5]-OP[8][1]*t158+OP[8][0]*t160+OP[8][2]*t163;
- A0[7][5] = OP[8][6];
- A0[7][6] = OP[8][7];
- A0[7][7] = OP[8][8];
|