#include #include #include #include #include char *name = "math.log/1 on complex, real part"; double f_if(float re, float im) { float r22490 = re; float r22491 = r22490 * r22490; float r22492 = im; float r22493 = r22492 * r22492; float r22494 = r22491 + r22493; float r22495 = sqrt(r22494); float r22496 = log(r22495); return r22496; } double f_id(double re, double im) { double r22497 = re; double r22498 = r22497 * r22497; double r22499 = im; double r22500 = r22499 * r22499; double r22501 = r22498 + r22500; double r22502 = sqrt(r22501); double r22503 = log(r22502); return r22503; } double f_of(float re, float im) { float r22504 = re; float r22505 = -r22504; float r22506 = -1.1311462515526873e+108; bool r22507 = r22505 <= r22506; float r22508 = log(r22504); float r22509 = -1.3345219303414604e-129; bool r22510 = r22505 <= r22509; float r22511 = r22504 * r22504; float r22512 = im; float r22513 = r22512 * r22512; float r22514 = r22511 + r22513; float r22515 = sqrt(r22514); float r22516 = log(r22515); float r22517 = 7.28790255135744e-188; bool r22518 = r22505 <= r22517; float r22519 = log(r22512); float r22520 = 6.881791039338204e+19; bool r22521 = r22505 <= r22520; float r22522 = log(r22505); float r22523 = r22521 ? r22516 : r22522; float r22524 = r22518 ? r22519 : r22523; float r22525 = r22510 ? r22516 : r22524; float r22526 = r22507 ? r22508 : r22525; return r22526; } double f_od(double re, double im) { double r22527 = re; double r22528 = -r22527; double r22529 = -1.1311462515526873e+108; bool r22530 = r22528 <= r22529; double r22531 = log(r22527); double r22532 = -1.3345219303414604e-129; bool r22533 = r22528 <= r22532; double r22534 = r22527 * r22527; double r22535 = im; double r22536 = r22535 * r22535; double r22537 = r22534 + r22536; double r22538 = sqrt(r22537); double r22539 = log(r22538); double r22540 = 7.28790255135744e-188; bool r22541 = r22528 <= r22540; double r22542 = log(r22535); double r22543 = 6.881791039338204e+19; bool r22544 = r22528 <= r22543; double r22545 = log(r22528); double r22546 = r22544 ? r22539 : r22545; double r22547 = r22541 ? r22542 : r22546; double r22548 = r22533 ? r22539 : r22547; double r22549 = r22530 ? r22531 : r22548; return r22549; } 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 r22550, r22551, r22552, r22553, r22554, r22555, r22556; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22550); mpfr_init(r22551); mpfr_init(r22552); mpfr_init(r22553); mpfr_init(r22554); mpfr_init(r22555); mpfr_init(r22556); } double f_im(double re, double im) { mpfr_set_d(r22550, re, MPFR_RNDN); mpfr_mul(r22551, r22550, r22550, MPFR_RNDN); mpfr_set_d(r22552, im, MPFR_RNDN); mpfr_mul(r22553, r22552, r22552, MPFR_RNDN); mpfr_add(r22554, r22551, r22553, MPFR_RNDN); mpfr_sqrt(r22555, r22554, MPFR_RNDN); mpfr_log(r22556, r22555, MPFR_RNDN); return mpfr_get_d(r22556, MPFR_RNDN); } static mpfr_t r22557, r22558, r22559, r22560, r22561, r22562, r22563, r22564, r22565, r22566, r22567, r22568, r22569, r22570, r22571, r22572, r22573, r22574, r22575, r22576, r22577, r22578, r22579; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22557); mpfr_init(r22558); mpfr_init_set_str(r22559, "-1.1311462515526873e+108", 10, MPFR_RNDN); mpfr_init(r22560); mpfr_init(r22561); mpfr_init_set_str(r22562, "-1.3345219303414604e-129", 10, MPFR_RNDN); mpfr_init(r22563); mpfr_init(r22564); mpfr_init(r22565); mpfr_init(r22566); mpfr_init(r22567); mpfr_init(r22568); mpfr_init(r22569); mpfr_init_set_str(r22570, "7.28790255135744e-188", 10, MPFR_RNDN); mpfr_init(r22571); mpfr_init(r22572); mpfr_init_set_str(r22573, "6.881791039338204e+19", 10, MPFR_RNDN); mpfr_init(r22574); mpfr_init(r22575); mpfr_init(r22576); mpfr_init(r22577); mpfr_init(r22578); mpfr_init(r22579); } double f_fm(double re, double im) { mpfr_set_d(r22557, re, MPFR_RNDN); mpfr_neg(r22558, r22557, MPFR_RNDN); ; mpfr_set_si(r22560, mpfr_cmp(r22558, r22559) <= 0, MPFR_RNDN); mpfr_log(r22561, r22557, MPFR_RNDN); ; mpfr_set_si(r22563, mpfr_cmp(r22558, r22562) <= 0, MPFR_RNDN); mpfr_mul(r22564, r22557, r22557, MPFR_RNDN); mpfr_set_d(r22565, im, MPFR_RNDN); mpfr_mul(r22566, r22565, r22565, MPFR_RNDN); mpfr_add(r22567, r22564, r22566, MPFR_RNDN); mpfr_sqrt(r22568, r22567, MPFR_RNDN); mpfr_log(r22569, r22568, MPFR_RNDN); ; mpfr_set_si(r22571, mpfr_cmp(r22558, r22570) <= 0, MPFR_RNDN); mpfr_log(r22572, r22565, MPFR_RNDN); ; mpfr_set_si(r22574, mpfr_cmp(r22558, r22573) <= 0, MPFR_RNDN); mpfr_log(r22575, r22558, MPFR_RNDN); if (mpfr_get_si(r22574, MPFR_RNDN)) { mpfr_set(r22576, r22569, MPFR_RNDN); } else { mpfr_set(r22576, r22575, MPFR_RNDN); }; if (mpfr_get_si(r22571, MPFR_RNDN)) { mpfr_set(r22577, r22572, MPFR_RNDN); } else { mpfr_set(r22577, r22576, MPFR_RNDN); }; if (mpfr_get_si(r22563, MPFR_RNDN)) { mpfr_set(r22578, r22569, MPFR_RNDN); } else { mpfr_set(r22578, r22577, MPFR_RNDN); }; if (mpfr_get_si(r22560, MPFR_RNDN)) { mpfr_set(r22579, r22561, MPFR_RNDN); } else { mpfr_set(r22579, r22578, MPFR_RNDN); }; return mpfr_get_d(r22579, MPFR_RNDN); } static mpfr_t r22580, r22581, r22582, r22583, r22584, r22585, r22586, r22587, r22588, r22589, r22590, r22591, r22592, r22593, r22594, r22595, r22596, r22597, r22598, r22599, r22600, r22601, r22602; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22580); mpfr_init(r22581); mpfr_init_set_str(r22582, "-1.1311462515526873e+108", 10, MPFR_RNDN); mpfr_init(r22583); mpfr_init(r22584); mpfr_init_set_str(r22585, "-1.3345219303414604e-129", 10, MPFR_RNDN); mpfr_init(r22586); mpfr_init(r22587); mpfr_init(r22588); mpfr_init(r22589); mpfr_init(r22590); mpfr_init(r22591); mpfr_init(r22592); mpfr_init_set_str(r22593, "7.28790255135744e-188", 10, MPFR_RNDN); mpfr_init(r22594); mpfr_init(r22595); mpfr_init_set_str(r22596, "6.881791039338204e+19", 10, MPFR_RNDN); mpfr_init(r22597); mpfr_init(r22598); mpfr_init(r22599); mpfr_init(r22600); mpfr_init(r22601); mpfr_init(r22602); } double f_dm(double re, double im) { mpfr_set_d(r22580, re, MPFR_RNDN); mpfr_neg(r22581, r22580, MPFR_RNDN); ; mpfr_set_si(r22583, mpfr_cmp(r22581, r22582) <= 0, MPFR_RNDN); mpfr_log(r22584, r22580, MPFR_RNDN); ; mpfr_set_si(r22586, mpfr_cmp(r22581, r22585) <= 0, MPFR_RNDN); mpfr_mul(r22587, r22580, r22580, MPFR_RNDN); mpfr_set_d(r22588, im, MPFR_RNDN); mpfr_mul(r22589, r22588, r22588, MPFR_RNDN); mpfr_add(r22590, r22587, r22589, MPFR_RNDN); mpfr_sqrt(r22591, r22590, MPFR_RNDN); mpfr_log(r22592, r22591, MPFR_RNDN); ; mpfr_set_si(r22594, mpfr_cmp(r22581, r22593) <= 0, MPFR_RNDN); mpfr_log(r22595, r22588, MPFR_RNDN); ; mpfr_set_si(r22597, mpfr_cmp(r22581, r22596) <= 0, MPFR_RNDN); mpfr_log(r22598, r22581, MPFR_RNDN); if (mpfr_get_si(r22597, MPFR_RNDN)) { mpfr_set(r22599, r22592, MPFR_RNDN); } else { mpfr_set(r22599, r22598, MPFR_RNDN); }; if (mpfr_get_si(r22594, MPFR_RNDN)) { mpfr_set(r22600, r22595, MPFR_RNDN); } else { mpfr_set(r22600, r22599, MPFR_RNDN); }; if (mpfr_get_si(r22586, MPFR_RNDN)) { mpfr_set(r22601, r22592, MPFR_RNDN); } else { mpfr_set(r22601, r22600, MPFR_RNDN); }; if (mpfr_get_si(r22583, MPFR_RNDN)) { mpfr_set(r22602, r22584, MPFR_RNDN); } else { mpfr_set(r22602, r22601, MPFR_RNDN); }; return mpfr_get_d(r22602, MPFR_RNDN); }