#include #include #include #include #include char *name = "math.sqrt on complex, real part"; double f_if(float re, float im) { float r25410 = 0.5; float r25411 = 2.0; float r25412 = re; float r25413 = r25412 * r25412; float r25414 = im; float r25415 = r25414 * r25414; float r25416 = r25413 + r25415; float r25417 = sqrt(r25416); float r25418 = r25417 + r25412; float r25419 = r25411 * r25418; float r25420 = sqrt(r25419); float r25421 = r25410 * r25420; return r25421; } double f_id(double re, double im) { double r25422 = 0.5; double r25423 = 2.0; double r25424 = re; double r25425 = r25424 * r25424; double r25426 = im; double r25427 = r25426 * r25426; double r25428 = r25425 + r25427; double r25429 = sqrt(r25428); double r25430 = r25429 + r25424; double r25431 = r25423 * r25430; double r25432 = sqrt(r25431); double r25433 = r25422 * r25432; return r25433; } double f_of(float re, float im) { float r25434 = 0.5; float r25435 = 2.0; float r25436 = re; float r25437 = im; float r25438 = hypot(r25436, r25437); float r25439 = r25436 + r25438; float r25440 = log1p(r25439); float r25441 = expm1(r25440); float r25442 = r25435 * r25441; float r25443 = sqrt(r25442); float r25444 = r25434 * r25443; float r25445 = 1.434644772517729e-150; bool r25446 = r25444 <= r25445; float r25447 = r25437 * r25437; float r25448 = r25435 * r25447; float r25449 = sqrt(r25448); float r25450 = r25438 - r25436; float r25451 = sqrt(r25450); float r25452 = r25449 / r25451; float r25453 = r25434 * r25452; float r25454 = r25436 * r25435; float r25455 = fma(r25438, r25435, r25454); float r25456 = sqrt(r25455); float r25457 = r25434 * r25456; float r25458 = r25446 ? r25453 : r25457; return r25458; } double f_od(double re, double im) { double r25459 = 0.5; double r25460 = 2.0; double r25461 = re; double r25462 = im; double r25463 = hypot(r25461, r25462); double r25464 = r25461 + r25463; double r25465 = log1p(r25464); double r25466 = expm1(r25465); double r25467 = r25460 * r25466; double r25468 = sqrt(r25467); double r25469 = r25459 * r25468; double r25470 = 1.434644772517729e-150; bool r25471 = r25469 <= r25470; double r25472 = r25462 * r25462; double r25473 = r25460 * r25472; double r25474 = sqrt(r25473); double r25475 = r25463 - r25461; double r25476 = sqrt(r25475); double r25477 = r25474 / r25476; double r25478 = r25459 * r25477; double r25479 = r25461 * r25460; double r25480 = fma(r25463, r25460, r25479); double r25481 = sqrt(r25480); double r25482 = r25459 * r25481; double r25483 = r25471 ? r25478 : r25482; return r25483; } 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 r25484, r25485, r25486, r25487, r25488, r25489, r25490, r25491, r25492, r25493, r25494, r25495; void setup_mpfr_f_im() { mpfr_set_default_prec(3664); mpfr_init_set_str(r25484, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r25485, "2.0", 10, MPFR_RNDN); mpfr_init(r25486); mpfr_init(r25487); mpfr_init(r25488); mpfr_init(r25489); mpfr_init(r25490); mpfr_init(r25491); mpfr_init(r25492); mpfr_init(r25493); mpfr_init(r25494); mpfr_init(r25495); } double f_im(double re, double im) { ; ; mpfr_set_d(r25486, re, MPFR_RNDN); mpfr_mul(r25487, r25486, r25486, MPFR_RNDN); mpfr_set_d(r25488, im, MPFR_RNDN); mpfr_mul(r25489, r25488, r25488, MPFR_RNDN); mpfr_add(r25490, r25487, r25489, MPFR_RNDN); mpfr_sqrt(r25491, r25490, MPFR_RNDN); mpfr_add(r25492, r25491, r25486, MPFR_RNDN); mpfr_mul(r25493, r25485, r25492, MPFR_RNDN); mpfr_sqrt(r25494, r25493, MPFR_RNDN); mpfr_mul(r25495, r25484, r25494, MPFR_RNDN); return mpfr_get_d(r25495, MPFR_RNDN); } static mpfr_t r25496, r25497, r25498, r25499, r25500, r25501, r25502, r25503, r25504, r25505, r25506, r25507, r25508, r25509, r25510, r25511, r25512, r25513, r25514, r25515, r25516, r25517, r25518, r25519, r25520; void setup_mpfr_f_fm() { mpfr_set_default_prec(3664); mpfr_init_set_str(r25496, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r25497, "2.0", 10, MPFR_RNDN); mpfr_init(r25498); mpfr_init(r25499); mpfr_init(r25500); mpfr_init(r25501); mpfr_init(r25502); mpfr_init(r25503); mpfr_init(r25504); mpfr_init(r25505); mpfr_init(r25506); mpfr_init_set_str(r25507, "1.434644772517729e-150", 10, MPFR_RNDN); mpfr_init(r25508); mpfr_init(r25509); mpfr_init(r25510); mpfr_init(r25511); mpfr_init(r25512); mpfr_init(r25513); mpfr_init(r25514); mpfr_init(r25515); mpfr_init(r25516); mpfr_init(r25517); mpfr_init(r25518); mpfr_init(r25519); mpfr_init(r25520); } double f_fm(double re, double im) { ; ; mpfr_set_d(r25498, re, MPFR_RNDN); mpfr_set_d(r25499, im, MPFR_RNDN); mpfr_hypot(r25500, r25498, r25499, MPFR_RNDN); mpfr_add(r25501, r25498, r25500, MPFR_RNDN); mpfr_log1p(r25502, r25501, MPFR_RNDN); mpfr_expm1(r25503, r25502, MPFR_RNDN); mpfr_mul(r25504, r25497, r25503, MPFR_RNDN); mpfr_sqrt(r25505, r25504, MPFR_RNDN); mpfr_mul(r25506, r25496, r25505, MPFR_RNDN); ; mpfr_set_si(r25508, mpfr_cmp(r25506, r25507) <= 0, MPFR_RNDN); mpfr_mul(r25509, r25499, r25499, MPFR_RNDN); mpfr_mul(r25510, r25497, r25509, MPFR_RNDN); mpfr_sqrt(r25511, r25510, MPFR_RNDN); mpfr_sub(r25512, r25500, r25498, MPFR_RNDN); mpfr_sqrt(r25513, r25512, MPFR_RNDN); mpfr_div(r25514, r25511, r25513, MPFR_RNDN); mpfr_mul(r25515, r25496, r25514, MPFR_RNDN); mpfr_mul(r25516, r25498, r25497, MPFR_RNDN); mpfr_fma(r25517, r25500, r25497, r25516, MPFR_RNDN); mpfr_sqrt(r25518, r25517, MPFR_RNDN); mpfr_mul(r25519, r25496, r25518, MPFR_RNDN); if (mpfr_get_si(r25508, MPFR_RNDN)) { mpfr_set(r25520, r25515, MPFR_RNDN); } else { mpfr_set(r25520, r25519, MPFR_RNDN); }; return mpfr_get_d(r25520, MPFR_RNDN); } static mpfr_t r25521, r25522, r25523, r25524, r25525, r25526, r25527, r25528, r25529, r25530, r25531, r25532, r25533, r25534, r25535, r25536, r25537, r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545; void setup_mpfr_f_dm() { mpfr_set_default_prec(3664); mpfr_init_set_str(r25521, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r25522, "2.0", 10, MPFR_RNDN); mpfr_init(r25523); mpfr_init(r25524); mpfr_init(r25525); mpfr_init(r25526); mpfr_init(r25527); mpfr_init(r25528); mpfr_init(r25529); mpfr_init(r25530); mpfr_init(r25531); mpfr_init_set_str(r25532, "1.434644772517729e-150", 10, MPFR_RNDN); mpfr_init(r25533); mpfr_init(r25534); mpfr_init(r25535); mpfr_init(r25536); mpfr_init(r25537); mpfr_init(r25538); mpfr_init(r25539); mpfr_init(r25540); mpfr_init(r25541); mpfr_init(r25542); mpfr_init(r25543); mpfr_init(r25544); mpfr_init(r25545); } double f_dm(double re, double im) { ; ; mpfr_set_d(r25523, re, MPFR_RNDN); mpfr_set_d(r25524, im, MPFR_RNDN); mpfr_hypot(r25525, r25523, r25524, MPFR_RNDN); mpfr_add(r25526, r25523, r25525, MPFR_RNDN); mpfr_log1p(r25527, r25526, MPFR_RNDN); mpfr_expm1(r25528, r25527, MPFR_RNDN); mpfr_mul(r25529, r25522, r25528, MPFR_RNDN); mpfr_sqrt(r25530, r25529, MPFR_RNDN); mpfr_mul(r25531, r25521, r25530, MPFR_RNDN); ; mpfr_set_si(r25533, mpfr_cmp(r25531, r25532) <= 0, MPFR_RNDN); mpfr_mul(r25534, r25524, r25524, MPFR_RNDN); mpfr_mul(r25535, r25522, r25534, MPFR_RNDN); mpfr_sqrt(r25536, r25535, MPFR_RNDN); mpfr_sub(r25537, r25525, r25523, MPFR_RNDN); mpfr_sqrt(r25538, r25537, MPFR_RNDN); mpfr_div(r25539, r25536, r25538, MPFR_RNDN); mpfr_mul(r25540, r25521, r25539, MPFR_RNDN); mpfr_mul(r25541, r25523, r25522, MPFR_RNDN); mpfr_fma(r25542, r25525, r25522, r25541, MPFR_RNDN); mpfr_sqrt(r25543, r25542, MPFR_RNDN); mpfr_mul(r25544, r25521, r25543, MPFR_RNDN); if (mpfr_get_si(r25533, MPFR_RNDN)) { mpfr_set(r25545, r25540, MPFR_RNDN); } else { mpfr_set(r25545, r25544, MPFR_RNDN); }; return mpfr_get_d(r25545, MPFR_RNDN); }