#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 r27706 = 0.5;
        float r27707 = re;
        float r27708 = cos(r27707);
        float r27709 = r27706 * r27708;
        float r27710 = 0;
        float r27711 = im;
        float r27712 = r27710 - r27711;
        float r27713 = exp(r27712);
        float r27714 = exp(r27711);
        float r27715 = r27713 - r27714;
        float r27716 = r27709 * r27715;
        return r27716;
}

double f_id(double re, double im) {
        double r27717 = 0.5;
        double r27718 = re;
        double r27719 = cos(r27718);
        double r27720 = r27717 * r27719;
        double r27721 = 0;
        double r27722 = im;
        double r27723 = r27721 - r27722;
        double r27724 = exp(r27723);
        double r27725 = exp(r27722);
        double r27726 = r27724 - r27725;
        double r27727 = r27720 * r27726;
        return r27727;
}


double f_of(float re, float im) {
        float r27728 = im;
        float r27729 = 1/3;
        float r27730 = r27729 * r27728;
        float r27731 = 2;
        float r27732 = fma(r27730, r27728, r27731);
        float r27733 = 5;
        float r27734 = pow(r27728, r27733);
        float r27735 = 1/60;
        float r27736 = r27734 * r27735;
        float r27737 = fma(r27728, r27732, r27736);
        float r27738 = re;
        float r27739 = cos(r27738);
        float r27740 = 0.5;
        float r27741 = -r27740;
        float r27742 = r27739 * r27741;
        float r27743 = r27737 * r27742;
        return r27743;
}

double f_od(double re, double im) {
        double r27744 = im;
        double r27745 = 1/3;
        double r27746 = r27745 * r27744;
        double r27747 = 2;
        double r27748 = fma(r27746, r27744, r27747);
        double r27749 = 5;
        double r27750 = pow(r27744, r27749);
        double r27751 = 1/60;
        double r27752 = r27750 * r27751;
        double r27753 = fma(r27744, r27748, r27752);
        double r27754 = re;
        double r27755 = cos(r27754);
        double r27756 = 0.5;
        double r27757 = -r27756;
        double r27758 = r27755 * r27757;
        double r27759 = r27753 * r27758;
        return r27759;
}

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 r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767, r27768, r27769, r27770;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r27760, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27761);
        mpfr_init(r27762);
        mpfr_init(r27763);
        mpfr_init_set_str(r27764, "0", 10, MPFR_RNDN);
        mpfr_init(r27765);
        mpfr_init(r27766);
        mpfr_init(r27767);
        mpfr_init(r27768);
        mpfr_init(r27769);
        mpfr_init(r27770);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r27761, re, MPFR_RNDN);
        mpfr_cos(r27762, r27761, MPFR_RNDN);
        mpfr_mul(r27763, r27760, r27762, MPFR_RNDN);
        ;
        mpfr_set_d(r27765, im, MPFR_RNDN);
        mpfr_sub(r27766, r27764, r27765, MPFR_RNDN);
        mpfr_exp(r27767, r27766, MPFR_RNDN);
        mpfr_exp(r27768, r27765, MPFR_RNDN);
        mpfr_sub(r27769, r27767, r27768, MPFR_RNDN);
        mpfr_mul(r27770, r27763, r27769, MPFR_RNDN);
        return mpfr_get_d(r27770, MPFR_RNDN);
}

static mpfr_t r27771, r27772, r27773, r27774, r27775, r27776, r27777, r27778, r27779, r27780, r27781, r27782, r27783, r27784, r27785, r27786;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27771);
        mpfr_init_set_str(r27772, "1/3", 10, MPFR_RNDN);
        mpfr_init(r27773);
        mpfr_init_set_str(r27774, "2", 10, MPFR_RNDN);
        mpfr_init(r27775);
        mpfr_init_set_str(r27776, "5", 10, MPFR_RNDN);
        mpfr_init(r27777);
        mpfr_init_set_str(r27778, "1/60", 10, MPFR_RNDN);
        mpfr_init(r27779);
        mpfr_init(r27780);
        mpfr_init(r27781);
        mpfr_init(r27782);
        mpfr_init_set_str(r27783, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27784);
        mpfr_init(r27785);
        mpfr_init(r27786);
}

double f_fm(double re, double im) {
        mpfr_set_d(r27771, im, MPFR_RNDN);
        ;
        mpfr_mul(r27773, r27772, r27771, MPFR_RNDN);
        ;
        mpfr_fma(r27775, r27773, r27771, r27774, MPFR_RNDN);
        ;
        mpfr_pow(r27777, r27771, r27776, MPFR_RNDN);
        ;
        mpfr_mul(r27779, r27777, r27778, MPFR_RNDN);
        mpfr_fma(r27780, r27771, r27775, r27779, MPFR_RNDN);
        mpfr_set_d(r27781, re, MPFR_RNDN);
        mpfr_cos(r27782, r27781, MPFR_RNDN);
        ;
        mpfr_neg(r27784, r27783, MPFR_RNDN);
        mpfr_mul(r27785, r27782, r27784, MPFR_RNDN);
        mpfr_mul(r27786, r27780, r27785, MPFR_RNDN);
        return mpfr_get_d(r27786, MPFR_RNDN);
}

static mpfr_t r27787, r27788, r27789, r27790, r27791, r27792, r27793, r27794, r27795, r27796, r27797, r27798, r27799, r27800, r27801, r27802;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27787);
        mpfr_init_set_str(r27788, "1/3", 10, MPFR_RNDN);
        mpfr_init(r27789);
        mpfr_init_set_str(r27790, "2", 10, MPFR_RNDN);
        mpfr_init(r27791);
        mpfr_init_set_str(r27792, "5", 10, MPFR_RNDN);
        mpfr_init(r27793);
        mpfr_init_set_str(r27794, "1/60", 10, MPFR_RNDN);
        mpfr_init(r27795);
        mpfr_init(r27796);
        mpfr_init(r27797);
        mpfr_init(r27798);
        mpfr_init_set_str(r27799, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27800);
        mpfr_init(r27801);
        mpfr_init(r27802);
}

double f_dm(double re, double im) {
        mpfr_set_d(r27787, im, MPFR_RNDN);
        ;
        mpfr_mul(r27789, r27788, r27787, MPFR_RNDN);
        ;
        mpfr_fma(r27791, r27789, r27787, r27790, MPFR_RNDN);
        ;
        mpfr_pow(r27793, r27787, r27792, MPFR_RNDN);
        ;
        mpfr_mul(r27795, r27793, r27794, MPFR_RNDN);
        mpfr_fma(r27796, r27787, r27791, r27795, MPFR_RNDN);
        mpfr_set_d(r27797, re, MPFR_RNDN);
        mpfr_cos(r27798, r27797, MPFR_RNDN);
        ;
        mpfr_neg(r27800, r27799, MPFR_RNDN);
        mpfr_mul(r27801, r27798, r27800, MPFR_RNDN);
        mpfr_mul(r27802, r27796, r27801, MPFR_RNDN);
        return mpfr_get_d(r27802, MPFR_RNDN);
}