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

char *name = "powComplex, imaginary part";

double f_if(float x_re, float x_im, float y_re, float y_im) {
        float r7590 = x_re;
        float r7591 = r7590 * r7590;
        float r7592 = x_im;
        float r7593 = r7592 * r7592;
        float r7594 = r7591 + r7593;
        float r7595 = sqrt(r7594);
        float r7596 = log(r7595);
        float r7597 = y_re;
        float r7598 = r7596 * r7597;
        float r7599 = atan2(r7592, r7590);
        float r7600 = y_im;
        float r7601 = r7599 * r7600;
        float r7602 = r7598 - r7601;
        float r7603 = exp(r7602);
        float r7604 = r7596 * r7600;
        float r7605 = r7599 * r7597;
        float r7606 = r7604 + r7605;
        float r7607 = sin(r7606);
        float r7608 = r7603 * r7607;
        return r7608;
}

double f_id(double x_re, double x_im, double y_re, double y_im) {
        double r7609 = x_re;
        double r7610 = r7609 * r7609;
        double r7611 = x_im;
        double r7612 = r7611 * r7611;
        double r7613 = r7610 + r7612;
        double r7614 = sqrt(r7613);
        double r7615 = log(r7614);
        double r7616 = y_re;
        double r7617 = r7615 * r7616;
        double r7618 = atan2(r7611, r7609);
        double r7619 = y_im;
        double r7620 = r7618 * r7619;
        double r7621 = r7617 - r7620;
        double r7622 = exp(r7621);
        double r7623 = r7615 * r7619;
        double r7624 = r7618 * r7616;
        double r7625 = r7623 + r7624;
        double r7626 = sin(r7625);
        double r7627 = r7622 * r7626;
        return r7627;
}


double f_of(float x_re, float x_im, float y_re, float y_im) {
        float r7628 = y_im;
        float r7629 = x_re;
        float r7630 = x_im;
        float r7631 = hypot(r7629, r7630);
        float r7632 = log(r7631);
        float r7633 = atan2(r7630, r7629);
        float r7634 = y_re;
        float r7635 = r7633 * r7634;
        float r7636 = fma(r7628, r7632, r7635);
        float r7637 = sin(r7636);
        float r7638 = expm1(r7637);
        float r7639 = log1p(r7638);
        float r7640 = r7633 * r7628;
        float r7641 = r7634 * r7632;
        float r7642 = r7640 - r7641;
        float r7643 = exp(r7642);
        float r7644 = r7639 / r7643;
        return r7644;
}

double f_od(double x_re, double x_im, double y_re, double y_im) {
        double r7645 = y_im;
        double r7646 = x_re;
        double r7647 = x_im;
        double r7648 = hypot(r7646, r7647);
        double r7649 = log(r7648);
        double r7650 = atan2(r7647, r7646);
        double r7651 = y_re;
        double r7652 = r7650 * r7651;
        double r7653 = fma(r7645, r7649, r7652);
        double r7654 = sin(r7653);
        double r7655 = expm1(r7654);
        double r7656 = log1p(r7655);
        double r7657 = r7650 * r7645;
        double r7658 = r7651 * r7649;
        double r7659 = r7657 - r7658;
        double r7660 = exp(r7659);
        double r7661 = r7656 / r7660;
        return r7661;
}

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 r7662, r7663, r7664, r7665, r7666, r7667, r7668, r7669, r7670, r7671, r7672, r7673, r7674, r7675, r7676, r7677, r7678, r7679, r7680;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2128);
        mpfr_init(r7662);
        mpfr_init(r7663);
        mpfr_init(r7664);
        mpfr_init(r7665);
        mpfr_init(r7666);
        mpfr_init(r7667);
        mpfr_init(r7668);
        mpfr_init(r7669);
        mpfr_init(r7670);
        mpfr_init(r7671);
        mpfr_init(r7672);
        mpfr_init(r7673);
        mpfr_init(r7674);
        mpfr_init(r7675);
        mpfr_init(r7676);
        mpfr_init(r7677);
        mpfr_init(r7678);
        mpfr_init(r7679);
        mpfr_init(r7680);
}

double f_im(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r7662, x_re, MPFR_RNDN);
        mpfr_mul(r7663, r7662, r7662, MPFR_RNDN);
        mpfr_set_d(r7664, x_im, MPFR_RNDN);
        mpfr_mul(r7665, r7664, r7664, MPFR_RNDN);
        mpfr_add(r7666, r7663, r7665, MPFR_RNDN);
        mpfr_sqrt(r7667, r7666, MPFR_RNDN);
        mpfr_log(r7668, r7667, MPFR_RNDN);
        mpfr_set_d(r7669, y_re, MPFR_RNDN);
        mpfr_mul(r7670, r7668, r7669, MPFR_RNDN);
        mpfr_atan2(r7671, r7664, r7662, MPFR_RNDN);
        mpfr_set_d(r7672, y_im, MPFR_RNDN);
        mpfr_mul(r7673, r7671, r7672, MPFR_RNDN);
        mpfr_sub(r7674, r7670, r7673, MPFR_RNDN);
        mpfr_exp(r7675, r7674, MPFR_RNDN);
        mpfr_mul(r7676, r7668, r7672, MPFR_RNDN);
        mpfr_mul(r7677, r7671, r7669, MPFR_RNDN);
        mpfr_add(r7678, r7676, r7677, MPFR_RNDN);
        mpfr_sin(r7679, r7678, MPFR_RNDN);
        mpfr_mul(r7680, r7675, r7679, MPFR_RNDN);
        return mpfr_get_d(r7680, MPFR_RNDN);
}

static mpfr_t r7681, r7682, r7683, r7684, r7685, r7686, r7687, r7688, r7689, r7690, r7691, r7692, r7693, r7694, r7695, r7696, r7697;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2128);
        mpfr_init(r7681);
        mpfr_init(r7682);
        mpfr_init(r7683);
        mpfr_init(r7684);
        mpfr_init(r7685);
        mpfr_init(r7686);
        mpfr_init(r7687);
        mpfr_init(r7688);
        mpfr_init(r7689);
        mpfr_init(r7690);
        mpfr_init(r7691);
        mpfr_init(r7692);
        mpfr_init(r7693);
        mpfr_init(r7694);
        mpfr_init(r7695);
        mpfr_init(r7696);
        mpfr_init(r7697);
}

double f_fm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r7681, y_im, MPFR_RNDN);
        mpfr_set_d(r7682, x_re, MPFR_RNDN);
        mpfr_set_d(r7683, x_im, MPFR_RNDN);
        mpfr_hypot(r7684, r7682, r7683, MPFR_RNDN);
        mpfr_log(r7685, r7684, MPFR_RNDN);
        mpfr_atan2(r7686, r7683, r7682, MPFR_RNDN);
        mpfr_set_d(r7687, y_re, MPFR_RNDN);
        mpfr_mul(r7688, r7686, r7687, MPFR_RNDN);
        mpfr_fma(r7689, r7681, r7685, r7688, MPFR_RNDN);
        mpfr_sin(r7690, r7689, MPFR_RNDN);
        mpfr_expm1(r7691, r7690, MPFR_RNDN);
        mpfr_log1p(r7692, r7691, MPFR_RNDN);
        mpfr_mul(r7693, r7686, r7681, MPFR_RNDN);
        mpfr_mul(r7694, r7687, r7685, MPFR_RNDN);
        mpfr_sub(r7695, r7693, r7694, MPFR_RNDN);
        mpfr_exp(r7696, r7695, MPFR_RNDN);
        mpfr_div(r7697, r7692, r7696, MPFR_RNDN);
        return mpfr_get_d(r7697, MPFR_RNDN);
}

static mpfr_t r7698, r7699, r7700, r7701, r7702, r7703, r7704, r7705, r7706, r7707, r7708, r7709, r7710, r7711, r7712, r7713, r7714;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2128);
        mpfr_init(r7698);
        mpfr_init(r7699);
        mpfr_init(r7700);
        mpfr_init(r7701);
        mpfr_init(r7702);
        mpfr_init(r7703);
        mpfr_init(r7704);
        mpfr_init(r7705);
        mpfr_init(r7706);
        mpfr_init(r7707);
        mpfr_init(r7708);
        mpfr_init(r7709);
        mpfr_init(r7710);
        mpfr_init(r7711);
        mpfr_init(r7712);
        mpfr_init(r7713);
        mpfr_init(r7714);
}

double f_dm(double x_re, double x_im, double y_re, double y_im) {
        mpfr_set_d(r7698, y_im, MPFR_RNDN);
        mpfr_set_d(r7699, x_re, MPFR_RNDN);
        mpfr_set_d(r7700, x_im, MPFR_RNDN);
        mpfr_hypot(r7701, r7699, r7700, MPFR_RNDN);
        mpfr_log(r7702, r7701, MPFR_RNDN);
        mpfr_atan2(r7703, r7700, r7699, MPFR_RNDN);
        mpfr_set_d(r7704, y_re, MPFR_RNDN);
        mpfr_mul(r7705, r7703, r7704, MPFR_RNDN);
        mpfr_fma(r7706, r7698, r7702, r7705, MPFR_RNDN);
        mpfr_sin(r7707, r7706, MPFR_RNDN);
        mpfr_expm1(r7708, r7707, MPFR_RNDN);
        mpfr_log1p(r7709, r7708, MPFR_RNDN);
        mpfr_mul(r7710, r7703, r7698, MPFR_RNDN);
        mpfr_mul(r7711, r7704, r7702, MPFR_RNDN);
        mpfr_sub(r7712, r7710, r7711, MPFR_RNDN);
        mpfr_exp(r7713, r7712, MPFR_RNDN);
        mpfr_div(r7714, r7709, r7713, MPFR_RNDN);
        return mpfr_get_d(r7714, MPFR_RNDN);
}