#include #include #include #include #include char *name = "math.abs on complex"; double f_if(float re, float im) { float r16295 = re; float r16296 = r16295 * r16295; float r16297 = im; float r16298 = r16297 * r16297; float r16299 = r16296 + r16298; float r16300 = sqrt(r16299); return r16300; } double f_id(double re, double im) { double r16301 = re; double r16302 = r16301 * r16301; double r16303 = im; double r16304 = r16303 * r16303; double r16305 = r16302 + r16304; double r16306 = sqrt(r16305); return r16306; } double f_of(float re, float im) { float r16307 = re; float r16308 = -1.80096568655872e+17f; bool r16309 = r16307 <= r16308; float r16310 = -r16307; float r16311 = 8580637917184.0f; bool r16312 = r16307 <= r16311; float r16313 = r16307 * r16307; float r16314 = im; float r16315 = r16314 * r16314; float r16316 = r16313 + r16315; float r16317 = sqrt(r16316); float r16318 = 0.5f; float r16319 = r16318 * r16314; float r16320 = r16307 / r16314; float r16321 = r16319 / r16320; float r16322 = r16307 + r16321; float r16323 = r16312 ? r16317 : r16322; float r16324 = r16309 ? r16310 : r16323; return r16324; } double f_od(double re, double im) { double r16325 = re; double r16326 = -1.80096568655872e+17; bool r16327 = r16325 <= r16326; double r16328 = -r16325; double r16329 = 8580637917184.0; bool r16330 = r16325 <= r16329; double r16331 = r16325 * r16325; double r16332 = im; double r16333 = r16332 * r16332; double r16334 = r16331 + r16333; double r16335 = sqrt(r16334); double r16336 = 0.5; double r16337 = r16336 * r16332; double r16338 = r16325 / r16332; double r16339 = r16337 / r16338; double r16340 = r16325 + r16339; double r16341 = r16330 ? r16335 : r16340; double r16342 = r16327 ? r16328 : r16341; return r16342; } 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 r16343, r16344, r16345, r16346, r16347, r16348; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16343); mpfr_init(r16344); mpfr_init(r16345); mpfr_init(r16346); mpfr_init(r16347); mpfr_init(r16348); } double f_im(double re, double im) { mpfr_set_d(r16343, re, MPFR_RNDN); mpfr_mul(r16344, r16343, r16343, MPFR_RNDN); mpfr_set_d(r16345, im, MPFR_RNDN); mpfr_mul(r16346, r16345, r16345, MPFR_RNDN); mpfr_add(r16347, r16344, r16346, MPFR_RNDN); mpfr_sqrt(r16348, r16347, MPFR_RNDN); return mpfr_get_d(r16348, MPFR_RNDN); } static mpfr_t r16349, r16350, r16351, r16352, r16353, r16354, r16355, r16356, r16357, r16358, r16359, r16360, r16361, r16362, r16363, r16364, r16365, r16366; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16349); mpfr_init_set_str(r16350, "-1.8009657f+17", 10, MPFR_RNDN); mpfr_init(r16351); mpfr_init(r16352); mpfr_init_set_str(r16353, "8.580638f+12", 10, MPFR_RNDN); mpfr_init(r16354); mpfr_init(r16355); mpfr_init(r16356); mpfr_init(r16357); mpfr_init(r16358); mpfr_init(r16359); mpfr_init_set_str(r16360, "1/2", 10, MPFR_RNDN); mpfr_init(r16361); mpfr_init(r16362); mpfr_init(r16363); mpfr_init(r16364); mpfr_init(r16365); mpfr_init(r16366); } double f_fm(double re, double im) { mpfr_set_d(r16349, re, MPFR_RNDN); ; mpfr_set_si(r16351, mpfr_cmp(r16349, r16350) <= 0, MPFR_RNDN); mpfr_neg(r16352, r16349, MPFR_RNDN); ; mpfr_set_si(r16354, mpfr_cmp(r16349, r16353) <= 0, MPFR_RNDN); mpfr_sqr(r16355, r16349, MPFR_RNDN); mpfr_set_d(r16356, im, MPFR_RNDN); mpfr_mul(r16357, r16356, r16356, MPFR_RNDN); mpfr_add(r16358, r16355, r16357, MPFR_RNDN); mpfr_sqrt(r16359, r16358, MPFR_RNDN); ; mpfr_mul(r16361, r16360, r16356, MPFR_RNDN); mpfr_div(r16362, r16349, r16356, MPFR_RNDN); mpfr_div(r16363, r16361, r16362, MPFR_RNDN); mpfr_add(r16364, r16349, r16363, MPFR_RNDN); if (mpfr_get_si(r16354, MPFR_RNDN)) { mpfr_set(r16365, r16359, MPFR_RNDN); } else { mpfr_set(r16365, r16364, MPFR_RNDN); }; if (mpfr_get_si(r16351, MPFR_RNDN)) { mpfr_set(r16366, r16352, MPFR_RNDN); } else { mpfr_set(r16366, r16365, MPFR_RNDN); }; return mpfr_get_d(r16366, MPFR_RNDN); } static mpfr_t r16367, r16368, r16369, r16370, r16371, r16372, r16373, r16374, r16375, r16376, r16377, r16378, r16379, r16380, r16381, r16382, r16383, r16384; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16367); mpfr_init_set_str(r16368, "-1.8009657f+17", 10, MPFR_RNDN); mpfr_init(r16369); mpfr_init(r16370); mpfr_init_set_str(r16371, "8.580638f+12", 10, MPFR_RNDN); mpfr_init(r16372); mpfr_init(r16373); mpfr_init(r16374); mpfr_init(r16375); mpfr_init(r16376); mpfr_init(r16377); mpfr_init_set_str(r16378, "1/2", 10, MPFR_RNDN); mpfr_init(r16379); mpfr_init(r16380); mpfr_init(r16381); mpfr_init(r16382); mpfr_init(r16383); mpfr_init(r16384); } double f_dm(double re, double im) { mpfr_set_d(r16367, re, MPFR_RNDN); ; mpfr_set_si(r16369, mpfr_cmp(r16367, r16368) <= 0, MPFR_RNDN); mpfr_neg(r16370, r16367, MPFR_RNDN); ; mpfr_set_si(r16372, mpfr_cmp(r16367, r16371) <= 0, MPFR_RNDN); mpfr_sqr(r16373, r16367, MPFR_RNDN); mpfr_set_d(r16374, im, MPFR_RNDN); mpfr_mul(r16375, r16374, r16374, MPFR_RNDN); mpfr_add(r16376, r16373, r16375, MPFR_RNDN); mpfr_sqrt(r16377, r16376, MPFR_RNDN); ; mpfr_mul(r16379, r16378, r16374, MPFR_RNDN); mpfr_div(r16380, r16367, r16374, MPFR_RNDN); mpfr_div(r16381, r16379, r16380, MPFR_RNDN); mpfr_add(r16382, r16367, r16381, MPFR_RNDN); if (mpfr_get_si(r16372, MPFR_RNDN)) { mpfr_set(r16383, r16377, MPFR_RNDN); } else { mpfr_set(r16383, r16382, MPFR_RNDN); }; if (mpfr_get_si(r16369, MPFR_RNDN)) { mpfr_set(r16384, r16370, MPFR_RNDN); } else { mpfr_set(r16384, r16383, MPFR_RNDN); }; return mpfr_get_d(r16384, MPFR_RNDN); }