#include #include #include #include #include char *name = "powComplex, imaginary part"; double f_if(float x_re, float x_im, float y_re, float y_im) { float r20546 = x_re; float r20547 = r20546 * r20546; float r20548 = x_im; float r20549 = r20548 * r20548; float r20550 = r20547 + r20549; float r20551 = sqrt(r20550); float r20552 = log(r20551); float r20553 = y_re; float r20554 = r20552 * r20553; float r20555 = atan2(r20548, r20546); float r20556 = y_im; float r20557 = r20555 * r20556; float r20558 = r20554 - r20557; float r20559 = exp(r20558); float r20560 = r20552 * r20556; float r20561 = r20555 * r20553; float r20562 = r20560 + r20561; float r20563 = sin(r20562); float r20564 = r20559 * r20563; return r20564; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r20565 = x_re; double r20566 = r20565 * r20565; double r20567 = x_im; double r20568 = r20567 * r20567; double r20569 = r20566 + r20568; double r20570 = sqrt(r20569); double r20571 = log(r20570); double r20572 = y_re; double r20573 = r20571 * r20572; double r20574 = atan2(r20567, r20565); double r20575 = y_im; double r20576 = r20574 * r20575; double r20577 = r20573 - r20576; double r20578 = exp(r20577); double r20579 = r20571 * r20575; double r20580 = r20574 * r20572; double r20581 = r20579 + r20580; double r20582 = sin(r20581); double r20583 = r20578 * r20582; return r20583; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r20584 = x_im; float r20585 = x_re; float r20586 = hypot(r20584, r20585); float r20587 = log(r20586); float r20588 = y_im; float r20589 = atan2(r20584, r20585); float r20590 = y_re; float r20591 = r20589 * r20590; float r20592 = fma(r20587, r20588, r20591); float r20593 = sin(r20592); float r20594 = expm1(r20593); float r20595 = log1p(r20594); float r20596 = r20588 * r20589; float r20597 = r20587 * r20590; float r20598 = r20596 - r20597; float r20599 = exp(r20598); float r20600 = r20595 / r20599; return r20600; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r20601 = x_im; double r20602 = x_re; double r20603 = hypot(r20601, r20602); double r20604 = log(r20603); double r20605 = y_im; double r20606 = atan2(r20601, r20602); double r20607 = y_re; double r20608 = r20606 * r20607; double r20609 = fma(r20604, r20605, r20608); double r20610 = sin(r20609); double r20611 = expm1(r20610); double r20612 = log1p(r20611); double r20613 = r20605 * r20606; double r20614 = r20604 * r20607; double r20615 = r20613 - r20614; double r20616 = exp(r20615); double r20617 = r20612 / r20616; return r20617; } 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 r20618, r20619, r20620, r20621, r20622, r20623, r20624, r20625, r20626, r20627, r20628, r20629, r20630, r20631, r20632, r20633, r20634, r20635, r20636; void setup_mpfr_f_im() { mpfr_set_default_prec(2192); mpfr_init(r20618); mpfr_init(r20619); mpfr_init(r20620); mpfr_init(r20621); mpfr_init(r20622); mpfr_init(r20623); mpfr_init(r20624); mpfr_init(r20625); mpfr_init(r20626); mpfr_init(r20627); mpfr_init(r20628); mpfr_init(r20629); mpfr_init(r20630); mpfr_init(r20631); mpfr_init(r20632); mpfr_init(r20633); mpfr_init(r20634); mpfr_init(r20635); mpfr_init(r20636); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r20618, x_re, MPFR_RNDN); mpfr_mul(r20619, r20618, r20618, MPFR_RNDN); mpfr_set_d(r20620, x_im, MPFR_RNDN); mpfr_mul(r20621, r20620, r20620, MPFR_RNDN); mpfr_add(r20622, r20619, r20621, MPFR_RNDN); mpfr_sqrt(r20623, r20622, MPFR_RNDN); mpfr_log(r20624, r20623, MPFR_RNDN); mpfr_set_d(r20625, y_re, MPFR_RNDN); mpfr_mul(r20626, r20624, r20625, MPFR_RNDN); mpfr_atan2(r20627, r20620, r20618, MPFR_RNDN); mpfr_set_d(r20628, y_im, MPFR_RNDN); mpfr_mul(r20629, r20627, r20628, MPFR_RNDN); mpfr_sub(r20630, r20626, r20629, MPFR_RNDN); mpfr_exp(r20631, r20630, MPFR_RNDN); mpfr_mul(r20632, r20624, r20628, MPFR_RNDN); mpfr_mul(r20633, r20627, r20625, MPFR_RNDN); mpfr_add(r20634, r20632, r20633, MPFR_RNDN); mpfr_sin(r20635, r20634, MPFR_RNDN); mpfr_mul(r20636, r20631, r20635, MPFR_RNDN); return mpfr_get_d(r20636, MPFR_RNDN); } static mpfr_t r20637, r20638, r20639, r20640, r20641, r20642, r20643, r20644, r20645, r20646, r20647, r20648, r20649, r20650, r20651, r20652, r20653; void setup_mpfr_f_fm() { mpfr_set_default_prec(2192); mpfr_init(r20637); mpfr_init(r20638); mpfr_init(r20639); mpfr_init(r20640); mpfr_init(r20641); mpfr_init(r20642); mpfr_init(r20643); mpfr_init(r20644); mpfr_init(r20645); mpfr_init(r20646); mpfr_init(r20647); mpfr_init(r20648); mpfr_init(r20649); mpfr_init(r20650); mpfr_init(r20651); mpfr_init(r20652); mpfr_init(r20653); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r20637, x_im, MPFR_RNDN); mpfr_set_d(r20638, x_re, MPFR_RNDN); mpfr_hypot(r20639, r20637, r20638, MPFR_RNDN); mpfr_log(r20640, r20639, MPFR_RNDN); mpfr_set_d(r20641, y_im, MPFR_RNDN); mpfr_atan2(r20642, r20637, r20638, MPFR_RNDN); mpfr_set_d(r20643, y_re, MPFR_RNDN); mpfr_mul(r20644, r20642, r20643, MPFR_RNDN); mpfr_fma(r20645, r20640, r20641, r20644, MPFR_RNDN); mpfr_sin(r20646, r20645, MPFR_RNDN); mpfr_expm1(r20647, r20646, MPFR_RNDN); mpfr_log1p(r20648, r20647, MPFR_RNDN); mpfr_mul(r20649, r20641, r20642, MPFR_RNDN); mpfr_mul(r20650, r20640, r20643, MPFR_RNDN); mpfr_sub(r20651, r20649, r20650, MPFR_RNDN); mpfr_exp(r20652, r20651, MPFR_RNDN); mpfr_div(r20653, r20648, r20652, MPFR_RNDN); return mpfr_get_d(r20653, MPFR_RNDN); } static mpfr_t r20654, r20655, r20656, r20657, r20658, r20659, r20660, r20661, r20662, r20663, r20664, r20665, r20666, r20667, r20668, r20669, r20670; void setup_mpfr_f_dm() { mpfr_set_default_prec(2192); mpfr_init(r20654); mpfr_init(r20655); mpfr_init(r20656); mpfr_init(r20657); mpfr_init(r20658); mpfr_init(r20659); mpfr_init(r20660); mpfr_init(r20661); mpfr_init(r20662); mpfr_init(r20663); mpfr_init(r20664); mpfr_init(r20665); mpfr_init(r20666); mpfr_init(r20667); mpfr_init(r20668); mpfr_init(r20669); mpfr_init(r20670); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r20654, x_im, MPFR_RNDN); mpfr_set_d(r20655, x_re, MPFR_RNDN); mpfr_hypot(r20656, r20654, r20655, MPFR_RNDN); mpfr_log(r20657, r20656, MPFR_RNDN); mpfr_set_d(r20658, y_im, MPFR_RNDN); mpfr_atan2(r20659, r20654, r20655, MPFR_RNDN); mpfr_set_d(r20660, y_re, MPFR_RNDN); mpfr_mul(r20661, r20659, r20660, MPFR_RNDN); mpfr_fma(r20662, r20657, r20658, r20661, MPFR_RNDN); mpfr_sin(r20663, r20662, MPFR_RNDN); mpfr_expm1(r20664, r20663, MPFR_RNDN); mpfr_log1p(r20665, r20664, MPFR_RNDN); mpfr_mul(r20666, r20658, r20659, MPFR_RNDN); mpfr_mul(r20667, r20657, r20660, MPFR_RNDN); mpfr_sub(r20668, r20666, r20667, MPFR_RNDN); mpfr_exp(r20669, r20668, MPFR_RNDN); mpfr_div(r20670, r20665, r20669, MPFR_RNDN); return mpfr_get_d(r20670, MPFR_RNDN); }