#include #include #include #include #include char *name = "Linear.Matrix:det33 from linear-1.19.1.3"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i, float j) { float r17129 = x; float r17130 = y; float r17131 = z; float r17132 = r17130 * r17131; float r17133 = t; float r17134 = a; float r17135 = r17133 * r17134; float r17136 = r17132 - r17135; float r17137 = r17129 * r17136; float r17138 = b; float r17139 = c; float r17140 = r17139 * r17131; float r17141 = i; float r17142 = r17141 * r17134; float r17143 = r17140 - r17142; float r17144 = r17138 * r17143; float r17145 = r17137 - r17144; float r17146 = j; float r17147 = r17139 * r17133; float r17148 = r17141 * r17130; float r17149 = r17147 - r17148; float r17150 = r17146 * r17149; float r17151 = r17145 + r17150; return r17151; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j) { double r17152 = x; double r17153 = y; double r17154 = z; double r17155 = r17153 * r17154; double r17156 = t; double r17157 = a; double r17158 = r17156 * r17157; double r17159 = r17155 - r17158; double r17160 = r17152 * r17159; double r17161 = b; double r17162 = c; double r17163 = r17162 * r17154; double r17164 = i; double r17165 = r17164 * r17157; double r17166 = r17163 - r17165; double r17167 = r17161 * r17166; double r17168 = r17160 - r17167; double r17169 = j; double r17170 = r17162 * r17156; double r17171 = r17164 * r17153; double r17172 = r17170 - r17171; double r17173 = r17169 * r17172; double r17174 = r17168 + r17173; return r17174; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j) { float r17175 = x; float r17176 = y; float r17177 = z; float r17178 = r17176 * r17177; float r17179 = t; float r17180 = a; float r17181 = r17179 * r17180; float r17182 = r17178 - r17181; float r17183 = r17175 * r17182; float r17184 = b; float r17185 = c; float r17186 = r17185 * r17177; float r17187 = i; float r17188 = r17187 * r17180; float r17189 = r17186 - r17188; float r17190 = r17184 * r17189; float r17191 = cbrt(r17190); float r17192 = r17191 * (r17191 * r17191); float r17193 = r17183 - r17192; float r17194 = j; float r17195 = r17185 * r17179; float r17196 = r17187 * r17176; float r17197 = r17195 - r17196; float r17198 = r17194 * r17197; float r17199 = r17193 + r17198; return r17199; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j) { double r17200 = x; double r17201 = y; double r17202 = z; double r17203 = r17201 * r17202; double r17204 = t; double r17205 = a; double r17206 = r17204 * r17205; double r17207 = r17203 - r17206; double r17208 = r17200 * r17207; double r17209 = b; double r17210 = c; double r17211 = r17210 * r17202; double r17212 = i; double r17213 = r17212 * r17205; double r17214 = r17211 - r17213; double r17215 = r17209 * r17214; double r17216 = cbrt(r17215); double r17217 = r17216 * (r17216 * r17216); double r17218 = r17208 - r17217; double r17219 = j; double r17220 = r17210 * r17204; double r17221 = r17212 * r17201; double r17222 = r17220 - r17221; double r17223 = r17219 * r17222; double r17224 = r17218 + r17223; return r17224; } void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) { mpfr_fmod(r, n, d, rmd); if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd); } static mpfr_t r17225, r17226, r17227, r17228, r17229, r17230, r17231, r17232, r17233, r17234, r17235, r17236, r17237, r17238, r17239, r17240, r17241, r17242, r17243, r17244, r17245, r17246, r17247; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17225); mpfr_init(r17226); mpfr_init(r17227); mpfr_init(r17228); mpfr_init(r17229); mpfr_init(r17230); mpfr_init(r17231); mpfr_init(r17232); mpfr_init(r17233); mpfr_init(r17234); mpfr_init(r17235); mpfr_init(r17236); mpfr_init(r17237); mpfr_init(r17238); mpfr_init(r17239); mpfr_init(r17240); mpfr_init(r17241); mpfr_init(r17242); mpfr_init(r17243); mpfr_init(r17244); mpfr_init(r17245); mpfr_init(r17246); mpfr_init(r17247); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r17225, x, MPFR_RNDN); mpfr_set_d(r17226, y, MPFR_RNDN); mpfr_set_d(r17227, z, MPFR_RNDN); mpfr_mul(r17228, r17226, r17227, MPFR_RNDN); mpfr_set_d(r17229, t, MPFR_RNDN); mpfr_set_d(r17230, a, MPFR_RNDN); mpfr_mul(r17231, r17229, r17230, MPFR_RNDN); mpfr_sub(r17232, r17228, r17231, MPFR_RNDN); mpfr_mul(r17233, r17225, r17232, MPFR_RNDN); mpfr_set_d(r17234, b, MPFR_RNDN); mpfr_set_d(r17235, c, MPFR_RNDN); mpfr_mul(r17236, r17235, r17227, MPFR_RNDN); mpfr_set_d(r17237, i, MPFR_RNDN); mpfr_mul(r17238, r17237, r17230, MPFR_RNDN); mpfr_sub(r17239, r17236, r17238, MPFR_RNDN); mpfr_mul(r17240, r17234, r17239, MPFR_RNDN); mpfr_sub(r17241, r17233, r17240, MPFR_RNDN); mpfr_set_d(r17242, j, MPFR_RNDN); mpfr_mul(r17243, r17235, r17229, MPFR_RNDN); mpfr_mul(r17244, r17237, r17226, MPFR_RNDN); mpfr_sub(r17245, r17243, r17244, MPFR_RNDN); mpfr_mul(r17246, r17242, r17245, MPFR_RNDN); mpfr_add(r17247, r17241, r17246, MPFR_RNDN); return mpfr_get_d(r17247, MPFR_RNDN); } static mpfr_t r17248, r17249, r17250, r17251, r17252, r17253, r17254, r17255, r17256, r17257, r17258, r17259, r17260, r17261, r17262, r17263, r17264, r17265, r17266, r17267, r17268, r17269, r17270, r17271, r17272; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17248); mpfr_init(r17249); mpfr_init(r17250); mpfr_init(r17251); mpfr_init(r17252); mpfr_init(r17253); mpfr_init(r17254); mpfr_init(r17255); mpfr_init(r17256); mpfr_init(r17257); mpfr_init(r17258); mpfr_init(r17259); mpfr_init(r17260); mpfr_init(r17261); mpfr_init(r17262); mpfr_init(r17263); mpfr_init(r17264); mpfr_init(r17265); mpfr_init(r17266); mpfr_init(r17267); mpfr_init(r17268); mpfr_init(r17269); mpfr_init(r17270); mpfr_init(r17271); mpfr_init(r17272); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r17248, x, MPFR_RNDN); mpfr_set_d(r17249, y, MPFR_RNDN); mpfr_set_d(r17250, z, MPFR_RNDN); mpfr_mul(r17251, r17249, r17250, MPFR_RNDN); mpfr_set_d(r17252, t, MPFR_RNDN); mpfr_set_d(r17253, a, MPFR_RNDN); mpfr_mul(r17254, r17252, r17253, MPFR_RNDN); mpfr_sub(r17255, r17251, r17254, MPFR_RNDN); mpfr_mul(r17256, r17248, r17255, MPFR_RNDN); mpfr_set_d(r17257, b, MPFR_RNDN); mpfr_set_d(r17258, c, MPFR_RNDN); mpfr_mul(r17259, r17258, r17250, MPFR_RNDN); mpfr_set_d(r17260, i, MPFR_RNDN); mpfr_mul(r17261, r17260, r17253, MPFR_RNDN); mpfr_sub(r17262, r17259, r17261, MPFR_RNDN); mpfr_mul(r17263, r17257, r17262, MPFR_RNDN); mpfr_cbrt(r17264, r17263, MPFR_RNDN); mpfr_mul(r17265, r17264, r17264, MPFR_RNDN); mpfr_mul(r17265, r17265, r17264, MPFR_RNDN); mpfr_sub(r17266, r17256, r17265, MPFR_RNDN); mpfr_set_d(r17267, j, MPFR_RNDN); mpfr_mul(r17268, r17258, r17252, MPFR_RNDN); mpfr_mul(r17269, r17260, r17249, MPFR_RNDN); mpfr_sub(r17270, r17268, r17269, MPFR_RNDN); mpfr_mul(r17271, r17267, r17270, MPFR_RNDN); mpfr_add(r17272, r17266, r17271, MPFR_RNDN); return mpfr_get_d(r17272, MPFR_RNDN); } static mpfr_t r17273, r17274, r17275, r17276, r17277, r17278, r17279, r17280, r17281, r17282, r17283, r17284, r17285, r17286, r17287, r17288, r17289, r17290, r17291, r17292, r17293, r17294, r17295, r17296, r17297; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17273); mpfr_init(r17274); mpfr_init(r17275); mpfr_init(r17276); mpfr_init(r17277); mpfr_init(r17278); mpfr_init(r17279); mpfr_init(r17280); mpfr_init(r17281); mpfr_init(r17282); mpfr_init(r17283); mpfr_init(r17284); mpfr_init(r17285); mpfr_init(r17286); mpfr_init(r17287); mpfr_init(r17288); mpfr_init(r17289); mpfr_init(r17290); mpfr_init(r17291); mpfr_init(r17292); mpfr_init(r17293); mpfr_init(r17294); mpfr_init(r17295); mpfr_init(r17296); mpfr_init(r17297); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r17273, x, MPFR_RNDN); mpfr_set_d(r17274, y, MPFR_RNDN); mpfr_set_d(r17275, z, MPFR_RNDN); mpfr_mul(r17276, r17274, r17275, MPFR_RNDN); mpfr_set_d(r17277, t, MPFR_RNDN); mpfr_set_d(r17278, a, MPFR_RNDN); mpfr_mul(r17279, r17277, r17278, MPFR_RNDN); mpfr_sub(r17280, r17276, r17279, MPFR_RNDN); mpfr_mul(r17281, r17273, r17280, MPFR_RNDN); mpfr_set_d(r17282, b, MPFR_RNDN); mpfr_set_d(r17283, c, MPFR_RNDN); mpfr_mul(r17284, r17283, r17275, MPFR_RNDN); mpfr_set_d(r17285, i, MPFR_RNDN); mpfr_mul(r17286, r17285, r17278, MPFR_RNDN); mpfr_sub(r17287, r17284, r17286, MPFR_RNDN); mpfr_mul(r17288, r17282, r17287, MPFR_RNDN); mpfr_cbrt(r17289, r17288, MPFR_RNDN); mpfr_mul(r17290, r17289, r17289, MPFR_RNDN); mpfr_mul(r17290, r17290, r17289, MPFR_RNDN); mpfr_sub(r17291, r17281, r17290, MPFR_RNDN); mpfr_set_d(r17292, j, MPFR_RNDN); mpfr_mul(r17293, r17283, r17277, MPFR_RNDN); mpfr_mul(r17294, r17285, r17274, MPFR_RNDN); mpfr_sub(r17295, r17293, r17294, MPFR_RNDN); mpfr_mul(r17296, r17292, r17295, MPFR_RNDN); mpfr_add(r17297, r17291, r17296, MPFR_RNDN); return mpfr_get_d(r17297, MPFR_RNDN); }