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

char *name = "math.sin on complex, imaginary part";

double f_if(float re, float im) {
        float r18687 = 0.5f;
        float r18688 = re;
        float r18689 = cos(r18688);
        float r18690 = r18687 * r18689;
        float r18691 = 0.0f;
        float r18692 = im;
        float r18693 = r18691 - r18692;
        float r18694 = exp(r18693);
        float r18695 = exp(r18692);
        float r18696 = r18694 - r18695;
        float r18697 = r18690 * r18696;
        return r18697;
}

double f_id(double re, double im) {
        double r18698 = 0.5;
        double r18699 = re;
        double r18700 = cos(r18699);
        double r18701 = r18698 * r18700;
        double r18702 = 0.0;
        double r18703 = im;
        double r18704 = r18702 - r18703;
        double r18705 = exp(r18704);
        double r18706 = exp(r18703);
        double r18707 = r18705 - r18706;
        double r18708 = r18701 * r18707;
        return r18708;
}


double f_of(float re, float im) {
        float r18709 = re;
        float r18710 = cos(r18709);
        float r18711 = 0.5f;
        float r18712 = r18710 * r18711;
        float r18713 = im;
        float r18714 = 5.0f;
        float r18715 = pow(r18713, r18714);
        float r18716 = 0.0005208333604969084f;
        float r18717 = 0.0416666679084301f;
        float r18718 = r18713 * (r18713 * r18713);
        float r18719 = fma(r18717, r18718, r18713);
        float r18720 = fma(r18715, r18716, r18719);
        float r18721 = -r18720;
        float r18722 = exp(r18713);
        float r18723 = sqrt(r18722);
        float r18724 = -r18713;
        float r18725 = exp(r18724);
        float r18726 = sqrt(r18725);
        float r18727 = r18723 + r18726;
        float r18728 = r18721 * r18727;
        float r18729 = r18712 * r18728;
        return r18729;
}

double f_od(double re, double im) {
        double r18730 = re;
        double r18731 = cos(r18730);
        double r18732 = 0.5;
        double r18733 = r18731 * r18732;
        double r18734 = im;
        double r18735 = 5.0;
        double r18736 = pow(r18734, r18735);
        double r18737 = 0.0005208333604969084;
        double r18738 = 0.0416666679084301;
        double r18739 = r18734 * (r18734 * r18734);
        double r18740 = fma(r18738, r18739, r18734);
        double r18741 = fma(r18736, r18737, r18740);
        double r18742 = -r18741;
        double r18743 = exp(r18734);
        double r18744 = sqrt(r18743);
        double r18745 = -r18734;
        double r18746 = exp(r18745);
        double r18747 = sqrt(r18746);
        double r18748 = r18744 + r18747;
        double r18749 = r18742 * r18748;
        double r18750 = r18733 * r18749;
        return r18750;
}

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 r18751, r18752, r18753, r18754, r18755, r18756, r18757, r18758, r18759, r18760, r18761;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18751, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18752);
        mpfr_init(r18753);
        mpfr_init(r18754);
        mpfr_init_set_str(r18755, "0", 10, MPFR_RNDN);
        mpfr_init(r18756);
        mpfr_init(r18757);
        mpfr_init(r18758);
        mpfr_init(r18759);
        mpfr_init(r18760);
        mpfr_init(r18761);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r18752, re, MPFR_RNDN);
        mpfr_cos(r18753, r18752, MPFR_RNDN);
        mpfr_mul(r18754, r18751, r18753, MPFR_RNDN);
        ;
        mpfr_set_d(r18756, im, MPFR_RNDN);
        mpfr_sub(r18757, r18755, r18756, MPFR_RNDN);
        mpfr_exp(r18758, r18757, MPFR_RNDN);
        mpfr_exp(r18759, r18756, MPFR_RNDN);
        mpfr_sub(r18760, r18758, r18759, MPFR_RNDN);
        mpfr_mul(r18761, r18754, r18760, MPFR_RNDN);
        return mpfr_get_d(r18761, MPFR_RNDN);
}

static mpfr_t r18762, r18763, r18764, r18765, r18766, r18767, r18768, r18769, r18770, r18771, r18772, r18773, r18774, r18775, r18776, r18777, r18778, r18779, r18780, r18781, r18782;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18762);
        mpfr_init(r18763);
        mpfr_init_set_str(r18764, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18765);
        mpfr_init(r18766);
        mpfr_init_set_str(r18767, "5", 10, MPFR_RNDN);
        mpfr_init(r18768);
        mpfr_init_set_str(r18769, "1/1920", 10, MPFR_RNDN);
        mpfr_init_set_str(r18770, "1/24", 10, MPFR_RNDN);
        mpfr_init(r18771);
        mpfr_init(r18772);
        mpfr_init(r18773);
        mpfr_init(r18774);
        mpfr_init(r18775);
        mpfr_init(r18776);
        mpfr_init(r18777);
        mpfr_init(r18778);
        mpfr_init(r18779);
        mpfr_init(r18780);
        mpfr_init(r18781);
        mpfr_init(r18782);
}

double f_fm(double re, double im) {
        mpfr_set_d(r18762, re, MPFR_RNDN);
        mpfr_cos(r18763, r18762, MPFR_RNDN);
        ;
        mpfr_mul(r18765, r18763, r18764, MPFR_RNDN);
        mpfr_set_d(r18766, im, MPFR_RNDN);
        ;
        mpfr_pow(r18768, r18766, r18767, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18771, r18766, r18766, MPFR_RNDN); mpfr_mul(r18771, r18771, r18766, MPFR_RNDN);
        mpfr_fma(r18772, r18770, r18771, r18766, MPFR_RNDN);
        mpfr_fma(r18773, r18768, r18769, r18772, MPFR_RNDN);
        mpfr_neg(r18774, r18773, MPFR_RNDN);
        mpfr_exp(r18775, r18766, MPFR_RNDN);
        mpfr_sqrt(r18776, r18775, MPFR_RNDN);
        mpfr_neg(r18777, r18766, MPFR_RNDN);
        mpfr_exp(r18778, r18777, MPFR_RNDN);
        mpfr_sqrt(r18779, r18778, MPFR_RNDN);
        mpfr_add(r18780, r18776, r18779, MPFR_RNDN);
        mpfr_mul(r18781, r18774, r18780, MPFR_RNDN);
        mpfr_mul(r18782, r18765, r18781, MPFR_RNDN);
        return mpfr_get_d(r18782, MPFR_RNDN);
}

static mpfr_t r18783, r18784, r18785, r18786, r18787, r18788, r18789, r18790, r18791, r18792, r18793, r18794, r18795, r18796, r18797, r18798, r18799, r18800, r18801, r18802, r18803;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18783);
        mpfr_init(r18784);
        mpfr_init_set_str(r18785, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18786);
        mpfr_init(r18787);
        mpfr_init_set_str(r18788, "5", 10, MPFR_RNDN);
        mpfr_init(r18789);
        mpfr_init_set_str(r18790, "1/1920", 10, MPFR_RNDN);
        mpfr_init_set_str(r18791, "1/24", 10, MPFR_RNDN);
        mpfr_init(r18792);
        mpfr_init(r18793);
        mpfr_init(r18794);
        mpfr_init(r18795);
        mpfr_init(r18796);
        mpfr_init(r18797);
        mpfr_init(r18798);
        mpfr_init(r18799);
        mpfr_init(r18800);
        mpfr_init(r18801);
        mpfr_init(r18802);
        mpfr_init(r18803);
}

double f_dm(double re, double im) {
        mpfr_set_d(r18783, re, MPFR_RNDN);
        mpfr_cos(r18784, r18783, MPFR_RNDN);
        ;
        mpfr_mul(r18786, r18784, r18785, MPFR_RNDN);
        mpfr_set_d(r18787, im, MPFR_RNDN);
        ;
        mpfr_pow(r18789, r18787, r18788, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18792, r18787, r18787, MPFR_RNDN); mpfr_mul(r18792, r18792, r18787, MPFR_RNDN);
        mpfr_fma(r18793, r18791, r18792, r18787, MPFR_RNDN);
        mpfr_fma(r18794, r18789, r18790, r18793, MPFR_RNDN);
        mpfr_neg(r18795, r18794, MPFR_RNDN);
        mpfr_exp(r18796, r18787, MPFR_RNDN);
        mpfr_sqrt(r18797, r18796, MPFR_RNDN);
        mpfr_neg(r18798, r18787, MPFR_RNDN);
        mpfr_exp(r18799, r18798, MPFR_RNDN);
        mpfr_sqrt(r18800, r18799, MPFR_RNDN);
        mpfr_add(r18801, r18797, r18800, MPFR_RNDN);
        mpfr_mul(r18802, r18795, r18801, MPFR_RNDN);
        mpfr_mul(r18803, r18786, r18802, MPFR_RNDN);
        return mpfr_get_d(r18803, MPFR_RNDN);
}