#include #include #include #include #include char *name = "math.sqrt on complex, imaginary part, im greater than 0 branch"; double f_if(float re, float im) { float r19019 = 0.5f; float r19020 = 2.0f; float r19021 = re; float r19022 = r19021 * r19021; float r19023 = im; float r19024 = r19023 * r19023; float r19025 = r19022 - r19024; float r19026 = sqrt(r19025); float r19027 = r19026 + r19021; float r19028 = r19020 * r19027; float r19029 = sqrt(r19028); float r19030 = r19019 * r19029; return r19030; } double f_id(double re, double im) { double r19031 = 0.5; double r19032 = 2.0; double r19033 = re; double r19034 = r19033 * r19033; double r19035 = im; double r19036 = r19035 * r19035; double r19037 = r19034 - r19036; double r19038 = sqrt(r19037); double r19039 = r19038 + r19033; double r19040 = r19032 * r19039; double r19041 = sqrt(r19040); double r19042 = r19031 * r19041; return r19042; } double f_of(float re, float im) { float r19043 = re; float r19044 = 3.280266654840657e-252f; bool r19045 = r19043 <= r19044; float r19046 = 0.5f; float r19047 = 2.0f; float r19048 = -1.0f; float r19049 = r19048 * r19043; float r19050 = r19049 + r19043; float r19051 = r19047 * r19050; float r19052 = sqrt(r19051); float r19053 = r19046 * r19052; float r19054 = im; float r19055 = r19043 + r19054; float r19056 = sqrt(r19055); float r19057 = r19043 - r19054; float r19058 = sqrt(r19057); float r19059 = r19056 * r19058; float r19060 = r19059 + r19043; float r19061 = r19047 * r19060; float r19062 = sqrt(r19061); float r19063 = r19046 * r19062; float r19064 = r19045 ? r19053 : r19063; return r19064; } double f_od(double re, double im) { double r19065 = re; double r19066 = 3.280266654840657e-252; bool r19067 = r19065 <= r19066; double r19068 = 0.5; double r19069 = 2.0; double r19070 = -1.0; double r19071 = r19070 * r19065; double r19072 = r19071 + r19065; double r19073 = r19069 * r19072; double r19074 = sqrt(r19073); double r19075 = r19068 * r19074; double r19076 = im; double r19077 = r19065 + r19076; double r19078 = sqrt(r19077); double r19079 = r19065 - r19076; double r19080 = sqrt(r19079); double r19081 = r19078 * r19080; double r19082 = r19081 + r19065; double r19083 = r19069 * r19082; double r19084 = sqrt(r19083); double r19085 = r19068 * r19084; double r19086 = r19067 ? r19075 : r19085; return r19086; } 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 r19087, r19088, r19089, r19090, r19091, r19092, r19093, r19094, r19095, r19096, r19097, r19098; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19087, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r19088, "2.0", 10, MPFR_RNDN); mpfr_init(r19089); mpfr_init(r19090); mpfr_init(r19091); mpfr_init(r19092); mpfr_init(r19093); mpfr_init(r19094); mpfr_init(r19095); mpfr_init(r19096); mpfr_init(r19097); mpfr_init(r19098); } double f_im(double re, double im) { ; ; mpfr_set_d(r19089, re, MPFR_RNDN); mpfr_mul(r19090, r19089, r19089, MPFR_RNDN); mpfr_set_d(r19091, im, MPFR_RNDN); mpfr_mul(r19092, r19091, r19091, MPFR_RNDN); mpfr_sub(r19093, r19090, r19092, MPFR_RNDN); mpfr_sqrt(r19094, r19093, MPFR_RNDN); mpfr_add(r19095, r19094, r19089, MPFR_RNDN); mpfr_mul(r19096, r19088, r19095, MPFR_RNDN); mpfr_sqrt(r19097, r19096, MPFR_RNDN); mpfr_mul(r19098, r19087, r19097, MPFR_RNDN); return mpfr_get_d(r19098, MPFR_RNDN); } static mpfr_t r19099, r19100, r19101, r19102, r19103, r19104, r19105, r19106, r19107, r19108, r19109, r19110, r19111, r19112, r19113, r19114, r19115, r19116, r19117, r19118, r19119, r19120; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19099); mpfr_init_set_str(r19100, "3.280266654840657e-252", 10, MPFR_RNDN); mpfr_init(r19101); mpfr_init_set_str(r19102, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r19103, "2.0", 10, MPFR_RNDN); mpfr_init_set_str(r19104, "-1", 10, MPFR_RNDN); mpfr_init(r19105); mpfr_init(r19106); mpfr_init(r19107); mpfr_init(r19108); mpfr_init(r19109); mpfr_init(r19110); mpfr_init(r19111); mpfr_init(r19112); mpfr_init(r19113); mpfr_init(r19114); mpfr_init(r19115); mpfr_init(r19116); mpfr_init(r19117); mpfr_init(r19118); mpfr_init(r19119); mpfr_init(r19120); } double f_fm(double re, double im) { mpfr_set_d(r19099, re, MPFR_RNDN); ; mpfr_set_si(r19101, mpfr_cmp(r19099, r19100) <= 0, MPFR_RNDN); ; ; ; mpfr_mul(r19105, r19104, r19099, MPFR_RNDN); mpfr_add(r19106, r19105, r19099, MPFR_RNDN); mpfr_mul(r19107, r19103, r19106, MPFR_RNDN); mpfr_sqrt(r19108, r19107, MPFR_RNDN); mpfr_mul(r19109, r19102, r19108, MPFR_RNDN); mpfr_set_d(r19110, im, MPFR_RNDN); mpfr_add(r19111, r19099, r19110, MPFR_RNDN); mpfr_sqrt(r19112, r19111, MPFR_RNDN); mpfr_sub(r19113, r19099, r19110, MPFR_RNDN); mpfr_sqrt(r19114, r19113, MPFR_RNDN); mpfr_mul(r19115, r19112, r19114, MPFR_RNDN); mpfr_add(r19116, r19115, r19099, MPFR_RNDN); mpfr_mul(r19117, r19103, r19116, MPFR_RNDN); mpfr_sqrt(r19118, r19117, MPFR_RNDN); mpfr_mul(r19119, r19102, r19118, MPFR_RNDN); if (mpfr_get_si(r19101, MPFR_RNDN)) { mpfr_set(r19120, r19109, MPFR_RNDN); } else { mpfr_set(r19120, r19119, MPFR_RNDN); }; return mpfr_get_d(r19120, MPFR_RNDN); } static mpfr_t r19121, r19122, r19123, r19124, r19125, r19126, r19127, r19128, r19129, r19130, r19131, r19132, r19133, r19134, r19135, r19136, r19137, r19138, r19139, r19140, r19141, r19142; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19121); mpfr_init_set_str(r19122, "3.280266654840657e-252", 10, MPFR_RNDN); mpfr_init(r19123); mpfr_init_set_str(r19124, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r19125, "2.0", 10, MPFR_RNDN); mpfr_init_set_str(r19126, "-1", 10, MPFR_RNDN); mpfr_init(r19127); mpfr_init(r19128); mpfr_init(r19129); mpfr_init(r19130); mpfr_init(r19131); mpfr_init(r19132); mpfr_init(r19133); mpfr_init(r19134); mpfr_init(r19135); mpfr_init(r19136); mpfr_init(r19137); mpfr_init(r19138); mpfr_init(r19139); mpfr_init(r19140); mpfr_init(r19141); mpfr_init(r19142); } double f_dm(double re, double im) { mpfr_set_d(r19121, re, MPFR_RNDN); ; mpfr_set_si(r19123, mpfr_cmp(r19121, r19122) <= 0, MPFR_RNDN); ; ; ; mpfr_mul(r19127, r19126, r19121, MPFR_RNDN); mpfr_add(r19128, r19127, r19121, MPFR_RNDN); mpfr_mul(r19129, r19125, r19128, MPFR_RNDN); mpfr_sqrt(r19130, r19129, MPFR_RNDN); mpfr_mul(r19131, r19124, r19130, MPFR_RNDN); mpfr_set_d(r19132, im, MPFR_RNDN); mpfr_add(r19133, r19121, r19132, MPFR_RNDN); mpfr_sqrt(r19134, r19133, MPFR_RNDN); mpfr_sub(r19135, r19121, r19132, MPFR_RNDN); mpfr_sqrt(r19136, r19135, MPFR_RNDN); mpfr_mul(r19137, r19134, r19136, MPFR_RNDN); mpfr_add(r19138, r19137, r19121, MPFR_RNDN); mpfr_mul(r19139, r19125, r19138, MPFR_RNDN); mpfr_sqrt(r19140, r19139, MPFR_RNDN); mpfr_mul(r19141, r19124, r19140, MPFR_RNDN); if (mpfr_get_si(r19123, MPFR_RNDN)) { mpfr_set(r19142, r19131, MPFR_RNDN); } else { mpfr_set(r19142, r19141, MPFR_RNDN); }; return mpfr_get_d(r19142, MPFR_RNDN); }