#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "powComplex, real part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r19531 = x_re;
        float r19532 = r19531 * r19531;
        float r19533 = x_im;
        float r19534 = r19533 * r19533;
        float r19535 = r19532 + r19534;
        float r19536 = sqrt(r19535);
        float r19537 = log(r19536);
        float r19538 = y_re;
        float r19539 = r19537 * r19538;
        float r19540 = atan2(r19533, r19531);
        float r19541 = y_im;
        float r19542 = r19540 * r19541;
        float r19543 = r19539 - r19542;
        float r19544 = exp(r19543);
        float r19545 = r19537 * r19541;
        float r19546 = r19540 * r19538;
        float r19547 = r19545 + r19546;
        float r19548 = cos(r19547);
        float r19549 = r19544 * r19548;
        return r19549;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r19550 = x_re;
        double r19551 = r19550 * r19550;
        double r19552 = x_im;
        double r19553 = r19552 * r19552;
        double r19554 = r19551 + r19553;
        double r19555 = sqrt(r19554);
        double r19556 = log(r19555);
        double r19557 = y_re;
        double r19558 = r19556 * r19557;
        double r19559 = atan2(r19552, r19550);
        double r19560 = y_im;
        double r19561 = r19559 * r19560;
        double r19562 = r19558 - r19561;
        double r19563 = exp(r19562);
        double r19564 = r19556 * r19560;
        double r19565 = r19559 * r19557;
        double r19566 = r19564 + r19565;
        double r19567 = cos(r19566);
        double r19568 = r19563 * r19567;
        return r19568;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r19569 = x_im;
        float r19570 = -5.4922881955228e-310;
        bool r19571 = r19569 <= r19570;
        float r19572 = y_re;
        float r19573 = x_re;
        float r19574 = atan2(r19569, r19573);
        float r19575 = r19572 * r19574;
        float r19576 = y_im;
        float r19577 = -r19569;
        float r19578 = log(r19577);
        float r19579 = r19576 * r19578;
        float r19580 = r19575 + r19579;
        float r19581 = cos(r19580);
        float r19582 = exp(r19576);
        float r19583 = pow(r19582, r19574);
        float r19584 = r19569 * r19569;
        float r19585 = r19573 * r19573;
        float r19586 = r19584 + r19585;
        float r19587 = sqrt(r19586);
        float r19588 = pow(r19587, r19572);
        float r19589 = r19583 / r19588;
        float r19590 = r19581 / r19589;
        float r19591 = log(r19569);
        float r19592 = r19576 * r19591;
        float r19593 = r19574 * r19572;
        float r19594 = r19592 + r19593;
        float r19595 = cos(r19594);
        float r19596 = r19595 / r19589;
        float r19597 = r19571 ? r19590 : r19596;
        return r19597;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r19598 = x_im;
        double r19599 = -5.4922881955228e-310;
        bool r19600 = r19598 <= r19599;
        double r19601 = y_re;
        double r19602 = x_re;
        double r19603 = atan2(r19598, r19602);
        double r19604 = r19601 * r19603;
        double r19605 = y_im;
        double r19606 = -r19598;
        double r19607 = log(r19606);
        double r19608 = r19605 * r19607;
        double r19609 = r19604 + r19608;
        double r19610 = cos(r19609);
        double r19611 = exp(r19605);
        double r19612 = pow(r19611, r19603);
        double r19613 = r19598 * r19598;
        double r19614 = r19602 * r19602;
        double r19615 = r19613 + r19614;
        double r19616 = sqrt(r19615);
        double r19617 = pow(r19616, r19601);
        double r19618 = r19612 / r19617;
        double r19619 = r19610 / r19618;
        double r19620 = log(r19598);
        double r19621 = r19605 * r19620;
        double r19622 = r19603 * r19601;
        double r19623 = r19621 + r19622;
        double r19624 = cos(r19623);
        double r19625 = r19624 / r19618;
        double r19626 = r19600 ? r19619 : r19625;
        return r19626;
}

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 r19627, r19628, r19629, r19630, r19631, r19632, r19633, r19634, r19635, r19636, r19637, r19638, r19639, r19640, r19641, r19642, r19643, r19644, r19645;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r19627);
        mpfr_init(r19628);
        mpfr_init(r19629);
        mpfr_init(r19630);
        mpfr_init(r19631);
        mpfr_init(r19632);
        mpfr_init(r19633);
        mpfr_init(r19634);
        mpfr_init(r19635);
        mpfr_init(r19636);
        mpfr_init(r19637);
        mpfr_init(r19638);
        mpfr_init(r19639);
        mpfr_init(r19640);
        mpfr_init(r19641);
        mpfr_init(r19642);
        mpfr_init(r19643);
        mpfr_init(r19644);
        mpfr_init(r19645);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r19627, x_re, MPFR_RNDN);
        mpfr_mul(r19628, r19627, r19627, MPFR_RNDN);
        mpfr_set_d(r19629, x_im, MPFR_RNDN);
        mpfr_mul(r19630, r19629, r19629, MPFR_RNDN);
        mpfr_add(r19631, r19628, r19630, MPFR_RNDN);
        mpfr_sqrt(r19632, r19631, MPFR_RNDN);
        mpfr_log(r19633, r19632, MPFR_RNDN);
        mpfr_set_d(r19634, y_re, MPFR_RNDN);
        mpfr_mul(r19635, r19633, r19634, MPFR_RNDN);
        mpfr_atan2(r19636, r19629, r19627, MPFR_RNDN);
        mpfr_set_d(r19637, y_im, MPFR_RNDN);
        mpfr_mul(r19638, r19636, r19637, MPFR_RNDN);
        mpfr_sub(r19639, r19635, r19638, MPFR_RNDN);
        mpfr_exp(r19640, r19639, MPFR_RNDN);
        mpfr_mul(r19641, r19633, r19637, MPFR_RNDN);
        mpfr_mul(r19642, r19636, r19634, MPFR_RNDN);
        mpfr_add(r19643, r19641, r19642, MPFR_RNDN);
        mpfr_cos(r19644, r19643, MPFR_RNDN);
        mpfr_mul(r19645, r19640, r19644, MPFR_RNDN);
        return mpfr_get_d(r19645, MPFR_RNDN);
}

static mpfr_t r19646, r19647, r19648, r19649, r19650, r19651, r19652, r19653, r19654, r19655, r19656, r19657, r19658, r19659, r19660, r19661, r19662, r19663, r19664, r19665, r19666, r19667, r19668, r19669, r19670, r19671, r19672, r19673, r19674;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r19646);
        mpfr_init_set_str(r19647, "-5.4922881955228e-310", 10, MPFR_RNDN);
        mpfr_init(r19648);
        mpfr_init(r19649);
        mpfr_init(r19650);
        mpfr_init(r19651);
        mpfr_init(r19652);
        mpfr_init(r19653);
        mpfr_init(r19654);
        mpfr_init(r19655);
        mpfr_init(r19656);
        mpfr_init(r19657);
        mpfr_init(r19658);
        mpfr_init(r19659);
        mpfr_init(r19660);
        mpfr_init(r19661);
        mpfr_init(r19662);
        mpfr_init(r19663);
        mpfr_init(r19664);
        mpfr_init(r19665);
        mpfr_init(r19666);
        mpfr_init(r19667);
        mpfr_init(r19668);
        mpfr_init(r19669);
        mpfr_init(r19670);
        mpfr_init(r19671);
        mpfr_init(r19672);
        mpfr_init(r19673);
        mpfr_init(r19674);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r19646, x_im, MPFR_RNDN);
        ;
        mpfr_set_si(r19648, mpfr_cmp(r19646, r19647) <= 0, MPFR_RNDN);
        mpfr_set_d(r19649, y_re, MPFR_RNDN);
        mpfr_set_d(r19650, x_re, MPFR_RNDN);
        mpfr_atan2(r19651, r19646, r19650, MPFR_RNDN);
        mpfr_mul(r19652, r19649, r19651, MPFR_RNDN);
        mpfr_set_d(r19653, y_im, MPFR_RNDN);
        mpfr_neg(r19654, r19646, MPFR_RNDN);
        mpfr_log(r19655, r19654, MPFR_RNDN);
        mpfr_mul(r19656, r19653, r19655, MPFR_RNDN);
        mpfr_add(r19657, r19652, r19656, MPFR_RNDN);
        mpfr_cos(r19658, r19657, MPFR_RNDN);
        mpfr_exp(r19659, r19653, MPFR_RNDN);
        mpfr_pow(r19660, r19659, r19651, MPFR_RNDN);
        mpfr_mul(r19661, r19646, r19646, MPFR_RNDN);
        mpfr_mul(r19662, r19650, r19650, MPFR_RNDN);
        mpfr_add(r19663, r19661, r19662, MPFR_RNDN);
        mpfr_sqrt(r19664, r19663, MPFR_RNDN);
        mpfr_pow(r19665, r19664, r19649, MPFR_RNDN);
        mpfr_div(r19666, r19660, r19665, MPFR_RNDN);
        mpfr_div(r19667, r19658, r19666, MPFR_RNDN);
        mpfr_log(r19668, r19646, MPFR_RNDN);
        mpfr_mul(r19669, r19653, r19668, MPFR_RNDN);
        mpfr_mul(r19670, r19651, r19649, MPFR_RNDN);
        mpfr_add(r19671, r19669, r19670, MPFR_RNDN);
        mpfr_cos(r19672, r19671, MPFR_RNDN);
        mpfr_div(r19673, r19672, r19666, MPFR_RNDN);
        if (mpfr_get_si(r19648, MPFR_RNDN)) { mpfr_set(r19674, r19667, MPFR_RNDN); } else { mpfr_set(r19674, r19673, MPFR_RNDN); };
        return mpfr_get_d(r19674, MPFR_RNDN);
}

static mpfr_t r19675, r19676, r19677, r19678, r19679, r19680, r19681, r19682, r19683, r19684, r19685, r19686, r19687, r19688, r19689, r19690, r19691, r19692, r19693, r19694, r19695, r19696, r19697, r19698, r19699, r19700, r19701, r19702, r19703;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r19675);
        mpfr_init_set_str(r19676, "-5.4922881955228e-310", 10, MPFR_RNDN);
        mpfr_init(r19677);
        mpfr_init(r19678);
        mpfr_init(r19679);
        mpfr_init(r19680);
        mpfr_init(r19681);
        mpfr_init(r19682);
        mpfr_init(r19683);
        mpfr_init(r19684);
        mpfr_init(r19685);
        mpfr_init(r19686);
        mpfr_init(r19687);
        mpfr_init(r19688);
        mpfr_init(r19689);
        mpfr_init(r19690);
        mpfr_init(r19691);
        mpfr_init(r19692);
        mpfr_init(r19693);
        mpfr_init(r19694);
        mpfr_init(r19695);
        mpfr_init(r19696);
        mpfr_init(r19697);
        mpfr_init(r19698);
        mpfr_init(r19699);
        mpfr_init(r19700);
        mpfr_init(r19701);
        mpfr_init(r19702);
        mpfr_init(r19703);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r19675, x_im, MPFR_RNDN);
        ;
        mpfr_set_si(r19677, mpfr_cmp(r19675, r19676) <= 0, MPFR_RNDN);
        mpfr_set_d(r19678, y_re, MPFR_RNDN);
        mpfr_set_d(r19679, x_re, MPFR_RNDN);
        mpfr_atan2(r19680, r19675, r19679, MPFR_RNDN);
        mpfr_mul(r19681, r19678, r19680, MPFR_RNDN);
        mpfr_set_d(r19682, y_im, MPFR_RNDN);
        mpfr_neg(r19683, r19675, MPFR_RNDN);
        mpfr_log(r19684, r19683, MPFR_RNDN);
        mpfr_mul(r19685, r19682, r19684, MPFR_RNDN);
        mpfr_add(r19686, r19681, r19685, MPFR_RNDN);
        mpfr_cos(r19687, r19686, MPFR_RNDN);
        mpfr_exp(r19688, r19682, MPFR_RNDN);
        mpfr_pow(r19689, r19688, r19680, MPFR_RNDN);
        mpfr_mul(r19690, r19675, r19675, MPFR_RNDN);
        mpfr_mul(r19691, r19679, r19679, MPFR_RNDN);
        mpfr_add(r19692, r19690, r19691, MPFR_RNDN);
        mpfr_sqrt(r19693, r19692, MPFR_RNDN);
        mpfr_pow(r19694, r19693, r19678, MPFR_RNDN);
        mpfr_div(r19695, r19689, r19694, MPFR_RNDN);
        mpfr_div(r19696, r19687, r19695, MPFR_RNDN);
        mpfr_log(r19697, r19675, MPFR_RNDN);
        mpfr_mul(r19698, r19682, r19697, MPFR_RNDN);
        mpfr_mul(r19699, r19680, r19678, MPFR_RNDN);
        mpfr_add(r19700, r19698, r19699, MPFR_RNDN);
        mpfr_cos(r19701, r19700, MPFR_RNDN);
        mpfr_div(r19702, r19701, r19695, MPFR_RNDN);
        if (mpfr_get_si(r19677, MPFR_RNDN)) { mpfr_set(r19703, r19696, MPFR_RNDN); } else { mpfr_set(r19703, r19702, MPFR_RNDN); };
        return mpfr_get_d(r19703, MPFR_RNDN);
}