#include #include #include #include #include char *name = "math.sqrt on complex, real part"; double f_if(float re, float im) { float r25146 = 0.5; float r25147 = 2.0; float r25148 = re; float r25149 = r25148 * r25148; float r25150 = im; float r25151 = r25150 * r25150; float r25152 = r25149 + r25151; float r25153 = sqrt(r25152); float r25154 = r25153 + r25148; float r25155 = r25147 * r25154; float r25156 = sqrt(r25155); float r25157 = r25146 * r25156; return r25157; } double f_id(double re, double im) { double r25158 = 0.5; double r25159 = 2.0; double r25160 = re; double r25161 = r25160 * r25160; double r25162 = im; double r25163 = r25162 * r25162; double r25164 = r25161 + r25163; double r25165 = sqrt(r25164); double r25166 = r25165 + r25160; double r25167 = r25159 * r25166; double r25168 = sqrt(r25167); double r25169 = r25158 * r25168; return r25169; } double f_of(float re, float im) { float r25170 = 0.5; float r25171 = re; float r25172 = im; float r25173 = hypot(r25171, r25172); float r25174 = 2.0; float r25175 = r25171 * r25174; float r25176 = fma(r25173, r25174, r25175); float r25177 = cbrt(r25176); float r25178 = r25177 * r25177; float r25179 = r25178 * r25177; float r25180 = sqrt(r25179); float r25181 = r25170 * r25180; float r25182 = 0.0; bool r25183 = r25181 <= r25182; float r25184 = 1.0; float r25185 = -1; float r25186 = r25185 / r25172; float r25187 = r25184 / r25186; float r25188 = r25185 / r25171; float r25189 = r25188 / r25186; float r25190 = r25187 * r25189; float r25191 = sqrt(r25190); float r25192 = r25191 * r25170; float r25193 = sqrt(r25176); float r25194 = r25170 * r25193; float r25195 = r25183 ? r25192 : r25194; return r25195; } double f_od(double re, double im) { double r25196 = 0.5; double r25197 = re; double r25198 = im; double r25199 = hypot(r25197, r25198); double r25200 = 2.0; double r25201 = r25197 * r25200; double r25202 = fma(r25199, r25200, r25201); double r25203 = cbrt(r25202); double r25204 = r25203 * r25203; double r25205 = r25204 * r25203; double r25206 = sqrt(r25205); double r25207 = r25196 * r25206; double r25208 = 0.0; bool r25209 = r25207 <= r25208; double r25210 = 1.0; double r25211 = -1; double r25212 = r25211 / r25198; double r25213 = r25210 / r25212; double r25214 = r25211 / r25197; double r25215 = r25214 / r25212; double r25216 = r25213 * r25215; double r25217 = sqrt(r25216); double r25218 = r25217 * r25196; double r25219 = sqrt(r25202); double r25220 = r25196 * r25219; double r25221 = r25209 ? r25218 : r25220; return r25221; } 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 r25222, r25223, r25224, r25225, r25226, r25227, r25228, r25229, r25230, r25231, r25232, r25233; void setup_mpfr_f_im() { mpfr_set_default_prec(3472); mpfr_init_set_str(r25222, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r25223, "2.0", 10, MPFR_RNDN); mpfr_init(r25224); mpfr_init(r25225); mpfr_init(r25226); mpfr_init(r25227); mpfr_init(r25228); mpfr_init(r25229); mpfr_init(r25230); mpfr_init(r25231); mpfr_init(r25232); mpfr_init(r25233); } double f_im(double re, double im) { ; ; mpfr_set_d(r25224, re, MPFR_RNDN); mpfr_mul(r25225, r25224, r25224, MPFR_RNDN); mpfr_set_d(r25226, im, MPFR_RNDN); mpfr_mul(r25227, r25226, r25226, MPFR_RNDN); mpfr_add(r25228, r25225, r25227, MPFR_RNDN); mpfr_sqrt(r25229, r25228, MPFR_RNDN); mpfr_add(r25230, r25229, r25224, MPFR_RNDN); mpfr_mul(r25231, r25223, r25230, MPFR_RNDN); mpfr_sqrt(r25232, r25231, MPFR_RNDN); mpfr_mul(r25233, r25222, r25232, MPFR_RNDN); return mpfr_get_d(r25233, MPFR_RNDN); } static mpfr_t r25234, r25235, r25236, r25237, r25238, r25239, r25240, r25241, r25242, r25243, r25244, r25245, r25246, r25247, r25248, r25249, r25250, r25251, r25252, r25253, r25254, r25255, r25256, r25257, r25258, r25259; void setup_mpfr_f_fm() { mpfr_set_default_prec(3472); mpfr_init_set_str(r25234, "0.5", 10, MPFR_RNDN); mpfr_init(r25235); mpfr_init(r25236); mpfr_init(r25237); mpfr_init_set_str(r25238, "2.0", 10, MPFR_RNDN); mpfr_init(r25239); mpfr_init(r25240); mpfr_init(r25241); mpfr_init(r25242); mpfr_init(r25243); mpfr_init(r25244); mpfr_init(r25245); mpfr_init_set_str(r25246, "0.0", 10, MPFR_RNDN); mpfr_init(r25247); mpfr_init_set_str(r25248, "1.0", 10, MPFR_RNDN); mpfr_init_set_str(r25249, "-1", 10, MPFR_RNDN); mpfr_init(r25250); mpfr_init(r25251); mpfr_init(r25252); mpfr_init(r25253); mpfr_init(r25254); mpfr_init(r25255); mpfr_init(r25256); mpfr_init(r25257); mpfr_init(r25258); mpfr_init(r25259); } double f_fm(double re, double im) { ; mpfr_set_d(r25235, re, MPFR_RNDN); mpfr_set_d(r25236, im, MPFR_RNDN); mpfr_hypot(r25237, r25235, r25236, MPFR_RNDN); ; mpfr_mul(r25239, r25235, r25238, MPFR_RNDN); mpfr_fma(r25240, r25237, r25238, r25239, MPFR_RNDN); mpfr_cbrt(r25241, r25240, MPFR_RNDN); mpfr_mul(r25242, r25241, r25241, MPFR_RNDN); mpfr_mul(r25243, r25242, r25241, MPFR_RNDN); mpfr_sqrt(r25244, r25243, MPFR_RNDN); mpfr_mul(r25245, r25234, r25244, MPFR_RNDN); ; mpfr_set_si(r25247, mpfr_cmp(r25245, r25246) <= 0, MPFR_RNDN); ; ; mpfr_div(r25250, r25249, r25236, MPFR_RNDN); mpfr_div(r25251, r25248, r25250, MPFR_RNDN); mpfr_div(r25252, r25249, r25235, MPFR_RNDN); mpfr_div(r25253, r25252, r25250, MPFR_RNDN); mpfr_mul(r25254, r25251, r25253, MPFR_RNDN); mpfr_sqrt(r25255, r25254, MPFR_RNDN); mpfr_mul(r25256, r25255, r25234, MPFR_RNDN); mpfr_sqrt(r25257, r25240, MPFR_RNDN); mpfr_mul(r25258, r25234, r25257, MPFR_RNDN); if (mpfr_get_si(r25247, MPFR_RNDN)) { mpfr_set(r25259, r25256, MPFR_RNDN); } else { mpfr_set(r25259, r25258, MPFR_RNDN); }; return mpfr_get_d(r25259, MPFR_RNDN); } static mpfr_t r25260, r25261, r25262, r25263, r25264, r25265, r25266, r25267, r25268, r25269, r25270, r25271, r25272, r25273, r25274, r25275, r25276, r25277, r25278, r25279, r25280, r25281, r25282, r25283, r25284, r25285; void setup_mpfr_f_dm() { mpfr_set_default_prec(3472); mpfr_init_set_str(r25260, "0.5", 10, MPFR_RNDN); mpfr_init(r25261); mpfr_init(r25262); mpfr_init(r25263); mpfr_init_set_str(r25264, "2.0", 10, MPFR_RNDN); mpfr_init(r25265); mpfr_init(r25266); mpfr_init(r25267); mpfr_init(r25268); mpfr_init(r25269); mpfr_init(r25270); mpfr_init(r25271); mpfr_init_set_str(r25272, "0.0", 10, MPFR_RNDN); mpfr_init(r25273); mpfr_init_set_str(r25274, "1.0", 10, MPFR_RNDN); mpfr_init_set_str(r25275, "-1", 10, MPFR_RNDN); mpfr_init(r25276); mpfr_init(r25277); mpfr_init(r25278); mpfr_init(r25279); mpfr_init(r25280); mpfr_init(r25281); mpfr_init(r25282); mpfr_init(r25283); mpfr_init(r25284); mpfr_init(r25285); } double f_dm(double re, double im) { ; mpfr_set_d(r25261, re, MPFR_RNDN); mpfr_set_d(r25262, im, MPFR_RNDN); mpfr_hypot(r25263, r25261, r25262, MPFR_RNDN); ; mpfr_mul(r25265, r25261, r25264, MPFR_RNDN); mpfr_fma(r25266, r25263, r25264, r25265, MPFR_RNDN); mpfr_cbrt(r25267, r25266, MPFR_RNDN); mpfr_mul(r25268, r25267, r25267, MPFR_RNDN); mpfr_mul(r25269, r25268, r25267, MPFR_RNDN); mpfr_sqrt(r25270, r25269, MPFR_RNDN); mpfr_mul(r25271, r25260, r25270, MPFR_RNDN); ; mpfr_set_si(r25273, mpfr_cmp(r25271, r25272) <= 0, MPFR_RNDN); ; ; mpfr_div(r25276, r25275, r25262, MPFR_RNDN); mpfr_div(r25277, r25274, r25276, MPFR_RNDN); mpfr_div(r25278, r25275, r25261, MPFR_RNDN); mpfr_div(r25279, r25278, r25276, MPFR_RNDN); mpfr_mul(r25280, r25277, r25279, MPFR_RNDN); mpfr_sqrt(r25281, r25280, MPFR_RNDN); mpfr_mul(r25282, r25281, r25260, MPFR_RNDN); mpfr_sqrt(r25283, r25266, MPFR_RNDN); mpfr_mul(r25284, r25260, r25283, MPFR_RNDN); if (mpfr_get_si(r25273, MPFR_RNDN)) { mpfr_set(r25285, r25282, MPFR_RNDN); } else { mpfr_set(r25285, r25284, MPFR_RNDN); }; return mpfr_get_d(r25285, MPFR_RNDN); }