#include #include #include #include #include char *name = "math.abs on complex"; double f_if(float re, float im) { float r16355 = re; float r16356 = r16355 * r16355; float r16357 = im; float r16358 = r16357 * r16357; float r16359 = r16356 + r16358; float r16360 = sqrt(r16359); return r16360; } double f_id(double re, double im) { double r16361 = re; double r16362 = r16361 * r16361; double r16363 = im; double r16364 = r16363 * r16363; double r16365 = r16362 + r16364; double r16366 = sqrt(r16365); return r16366; } double f_of(float re, float im) { float r16367 = re; float r16368 = -6.854525568260506e+16f; bool r16369 = r16367 <= r16368; float r16370 = -r16367; float r16371 = 96025305088.0f; bool r16372 = r16367 <= r16371; float r16373 = r16367 * r16367; float r16374 = im; float r16375 = r16374 * r16374; float r16376 = r16373 + r16375; float r16377 = sqrt(r16376); float r16378 = 0.5f; float r16379 = r16378 * r16374; float r16380 = r16367 / r16374; float r16381 = r16379 / r16380; float r16382 = r16367 + r16381; float r16383 = r16372 ? r16377 : r16382; float r16384 = r16369 ? r16370 : r16383; return r16384; } double f_od(double re, double im) { double r16385 = re; double r16386 = -6.854525568260506e+16; bool r16387 = r16385 <= r16386; double r16388 = -r16385; double r16389 = 96025305088.0; bool r16390 = r16385 <= r16389; double r16391 = r16385 * r16385; double r16392 = im; double r16393 = r16392 * r16392; double r16394 = r16391 + r16393; double r16395 = sqrt(r16394); double r16396 = 0.5; double r16397 = r16396 * r16392; double r16398 = r16385 / r16392; double r16399 = r16397 / r16398; double r16400 = r16385 + r16399; double r16401 = r16390 ? r16395 : r16400; double r16402 = r16387 ? r16388 : r16401; return r16402; } 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 r16403, r16404, r16405, r16406, r16407, r16408; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16403); mpfr_init(r16404); mpfr_init(r16405); mpfr_init(r16406); mpfr_init(r16407); mpfr_init(r16408); } double f_im(double re, double im) { mpfr_set_d(r16403, re, MPFR_RNDN); mpfr_mul(r16404, r16403, r16403, MPFR_RNDN); mpfr_set_d(r16405, im, MPFR_RNDN); mpfr_mul(r16406, r16405, r16405, MPFR_RNDN); mpfr_add(r16407, r16404, r16406, MPFR_RNDN); mpfr_sqrt(r16408, r16407, MPFR_RNDN); return mpfr_get_d(r16408, MPFR_RNDN); } static mpfr_t r16409, r16410, r16411, r16412, r16413, r16414, r16415, r16416, r16417, r16418, r16419, r16420, r16421, r16422, r16423, r16424, r16425, r16426; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16409); mpfr_init_set_str(r16410, "-6.8545256f+16", 10, MPFR_RNDN); mpfr_init(r16411); mpfr_init(r16412); mpfr_init_set_str(r16413, "9.6025305f+10", 10, MPFR_RNDN); mpfr_init(r16414); mpfr_init(r16415); mpfr_init(r16416); mpfr_init(r16417); mpfr_init(r16418); mpfr_init(r16419); mpfr_init_set_str(r16420, "1/2", 10, MPFR_RNDN); mpfr_init(r16421); mpfr_init(r16422); mpfr_init(r16423); mpfr_init(r16424); mpfr_init(r16425); mpfr_init(r16426); } double f_fm(double re, double im) { mpfr_set_d(r16409, re, MPFR_RNDN); ; mpfr_set_si(r16411, mpfr_cmp(r16409, r16410) <= 0, MPFR_RNDN); mpfr_neg(r16412, r16409, MPFR_RNDN); ; mpfr_set_si(r16414, mpfr_cmp(r16409, r16413) <= 0, MPFR_RNDN); mpfr_sqr(r16415, r16409, MPFR_RNDN); mpfr_set_d(r16416, im, MPFR_RNDN); mpfr_mul(r16417, r16416, r16416, MPFR_RNDN); mpfr_add(r16418, r16415, r16417, MPFR_RNDN); mpfr_sqrt(r16419, r16418, MPFR_RNDN); ; mpfr_mul(r16421, r16420, r16416, MPFR_RNDN); mpfr_div(r16422, r16409, r16416, MPFR_RNDN); mpfr_div(r16423, r16421, r16422, MPFR_RNDN); mpfr_add(r16424, r16409, r16423, MPFR_RNDN); if (mpfr_get_si(r16414, MPFR_RNDN)) { mpfr_set(r16425, r16419, MPFR_RNDN); } else { mpfr_set(r16425, r16424, MPFR_RNDN); }; if (mpfr_get_si(r16411, MPFR_RNDN)) { mpfr_set(r16426, r16412, MPFR_RNDN); } else { mpfr_set(r16426, r16425, MPFR_RNDN); }; return mpfr_get_d(r16426, MPFR_RNDN); } static mpfr_t r16427, r16428, r16429, r16430, r16431, r16432, r16433, r16434, r16435, r16436, r16437, r16438, r16439, r16440, r16441, r16442, r16443, r16444; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16427); mpfr_init_set_str(r16428, "-6.8545256f+16", 10, MPFR_RNDN); mpfr_init(r16429); mpfr_init(r16430); mpfr_init_set_str(r16431, "9.6025305f+10", 10, MPFR_RNDN); mpfr_init(r16432); mpfr_init(r16433); mpfr_init(r16434); mpfr_init(r16435); mpfr_init(r16436); mpfr_init(r16437); mpfr_init_set_str(r16438, "1/2", 10, MPFR_RNDN); mpfr_init(r16439); mpfr_init(r16440); mpfr_init(r16441); mpfr_init(r16442); mpfr_init(r16443); mpfr_init(r16444); } double f_dm(double re, double im) { mpfr_set_d(r16427, re, MPFR_RNDN); ; mpfr_set_si(r16429, mpfr_cmp(r16427, r16428) <= 0, MPFR_RNDN); mpfr_neg(r16430, r16427, MPFR_RNDN); ; mpfr_set_si(r16432, mpfr_cmp(r16427, r16431) <= 0, MPFR_RNDN); mpfr_sqr(r16433, r16427, MPFR_RNDN); mpfr_set_d(r16434, im, MPFR_RNDN); mpfr_mul(r16435, r16434, r16434, MPFR_RNDN); mpfr_add(r16436, r16433, r16435, MPFR_RNDN); mpfr_sqrt(r16437, r16436, MPFR_RNDN); ; mpfr_mul(r16439, r16438, r16434, MPFR_RNDN); mpfr_div(r16440, r16427, r16434, MPFR_RNDN); mpfr_div(r16441, r16439, r16440, MPFR_RNDN); mpfr_add(r16442, r16427, r16441, MPFR_RNDN); if (mpfr_get_si(r16432, MPFR_RNDN)) { mpfr_set(r16443, r16437, MPFR_RNDN); } else { mpfr_set(r16443, r16442, MPFR_RNDN); }; if (mpfr_get_si(r16429, MPFR_RNDN)) { mpfr_set(r16444, r16430, MPFR_RNDN); } else { mpfr_set(r16444, r16443, MPFR_RNDN); }; return mpfr_get_d(r16444, MPFR_RNDN); }