#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 r20529 = x_re; float r20530 = r20529 * r20529; float r20531 = x_im; float r20532 = r20531 * r20531; float r20533 = r20530 + r20532; float r20534 = sqrt(r20533); float r20535 = log(r20534); float r20536 = y_re; float r20537 = r20535 * r20536; float r20538 = atan2(r20531, r20529); float r20539 = y_im; float r20540 = r20538 * r20539; float r20541 = r20537 - r20540; float r20542 = exp(r20541); float r20543 = r20535 * r20539; float r20544 = r20538 * r20536; float r20545 = r20543 + r20544; float r20546 = sin(r20545); float r20547 = r20542 * r20546; return r20547; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r20548 = x_re; double r20549 = r20548 * r20548; double r20550 = x_im; double r20551 = r20550 * r20550; double r20552 = r20549 + r20551; double r20553 = sqrt(r20552); double r20554 = log(r20553); double r20555 = y_re; double r20556 = r20554 * r20555; double r20557 = atan2(r20550, r20548); double r20558 = y_im; double r20559 = r20557 * r20558; double r20560 = r20556 - r20559; double r20561 = exp(r20560); double r20562 = r20554 * r20558; double r20563 = r20557 * r20555; double r20564 = r20562 + r20563; double r20565 = sin(r20564); double r20566 = r20561 * r20565; return r20566; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r20567 = x_im; float r20568 = x_re; float r20569 = hypot(r20567, r20568); float r20570 = log(r20569); float r20571 = y_im; float r20572 = atan2(r20567, r20568); float r20573 = y_re; float r20574 = r20572 * r20573; float r20575 = fma(r20570, r20571, r20574); float r20576 = sin(r20575); float r20577 = expm1(r20576); float r20578 = log1p(r20577); float r20579 = r20571 * r20572; float r20580 = r20570 * r20573; float r20581 = r20579 - r20580; float r20582 = exp(r20581); float r20583 = r20578 / r20582; return r20583; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r20584 = x_im; double r20585 = x_re; double r20586 = hypot(r20584, r20585); double r20587 = log(r20586); double r20588 = y_im; double r20589 = atan2(r20584, r20585); double r20590 = y_re; double r20591 = r20589 * r20590; double r20592 = fma(r20587, r20588, r20591); double r20593 = sin(r20592); double r20594 = expm1(r20593); double r20595 = log1p(r20594); double r20596 = r20588 * r20589; double r20597 = r20587 * r20590; double r20598 = r20596 - r20597; double r20599 = exp(r20598); double r20600 = r20595 / r20599; return r20600; } 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 r20601, r20602, r20603, r20604, r20605, r20606, r20607, r20608, r20609, r20610, r20611, r20612, r20613, r20614, r20615, r20616, r20617, r20618, r20619; void setup_mpfr_f_im() { mpfr_set_default_prec(2192); mpfr_init(r20601); mpfr_init(r20602); mpfr_init(r20603); mpfr_init(r20604); mpfr_init(r20605); mpfr_init(r20606); mpfr_init(r20607); mpfr_init(r20608); mpfr_init(r20609); mpfr_init(r20610); mpfr_init(r20611); mpfr_init(r20612); mpfr_init(r20613); mpfr_init(r20614); mpfr_init(r20615); mpfr_init(r20616); mpfr_init(r20617); mpfr_init(r20618); mpfr_init(r20619); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r20601, x_re, MPFR_RNDN); mpfr_mul(r20602, r20601, r20601, MPFR_RNDN); mpfr_set_d(r20603, x_im, MPFR_RNDN); mpfr_mul(r20604, r20603, r20603, MPFR_RNDN); mpfr_add(r20605, r20602, r20604, MPFR_RNDN); mpfr_sqrt(r20606, r20605, MPFR_RNDN); mpfr_log(r20607, r20606, MPFR_RNDN); mpfr_set_d(r20608, y_re, MPFR_RNDN); mpfr_mul(r20609, r20607, r20608, MPFR_RNDN); mpfr_atan2(r20610, r20603, r20601, MPFR_RNDN); mpfr_set_d(r20611, y_im, MPFR_RNDN); mpfr_mul(r20612, r20610, r20611, MPFR_RNDN); mpfr_sub(r20613, r20609, r20612, MPFR_RNDN); mpfr_exp(r20614, r20613, MPFR_RNDN); mpfr_mul(r20615, r20607, r20611, MPFR_RNDN); mpfr_mul(r20616, r20610, r20608, MPFR_RNDN); mpfr_add(r20617, r20615, r20616, MPFR_RNDN); mpfr_sin(r20618, r20617, MPFR_RNDN); mpfr_mul(r20619, r20614, r20618, MPFR_RNDN); return mpfr_get_d(r20619, MPFR_RNDN); } static mpfr_t r20620, r20621, r20622, r20623, r20624, r20625, r20626, r20627, r20628, r20629, r20630, r20631, r20632, r20633, r20634, r20635, r20636; void setup_mpfr_f_fm() { mpfr_set_default_prec(2192); 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_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r20620, x_im, MPFR_RNDN); mpfr_set_d(r20621, x_re, MPFR_RNDN); mpfr_hypot(r20622, r20620, r20621, MPFR_RNDN); mpfr_log(r20623, r20622, MPFR_RNDN); mpfr_set_d(r20624, y_im, MPFR_RNDN); mpfr_atan2(r20625, r20620, r20621, MPFR_RNDN); mpfr_set_d(r20626, y_re, MPFR_RNDN); mpfr_mul(r20627, r20625, r20626, MPFR_RNDN); mpfr_fma(r20628, r20623, r20624, r20627, MPFR_RNDN); mpfr_sin(r20629, r20628, MPFR_RNDN); mpfr_expm1(r20630, r20629, MPFR_RNDN); mpfr_log1p(r20631, r20630, MPFR_RNDN); mpfr_mul(r20632, r20624, r20625, MPFR_RNDN); mpfr_mul(r20633, r20623, r20626, MPFR_RNDN); mpfr_sub(r20634, r20632, r20633, MPFR_RNDN); mpfr_exp(r20635, r20634, MPFR_RNDN); mpfr_div(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_dm() { 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_dm(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); }