#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 r18479 = 0.5f;
        float r18480 = re;
        float r18481 = cos(r18480);
        float r18482 = r18479 * r18481;
        float r18483 = 0.0f;
        float r18484 = im;
        float r18485 = r18483 - r18484;
        float r18486 = exp(r18485);
        float r18487 = exp(r18484);
        float r18488 = r18486 - r18487;
        float r18489 = r18482 * r18488;
        return r18489;
}

double f_id(double re, double im) {
        double r18490 = 0.5;
        double r18491 = re;
        double r18492 = cos(r18491);
        double r18493 = r18490 * r18492;
        double r18494 = 0.0;
        double r18495 = im;
        double r18496 = r18494 - r18495;
        double r18497 = exp(r18496);
        double r18498 = exp(r18495);
        double r18499 = r18497 - r18498;
        double r18500 = r18493 * r18499;
        return r18500;
}


double f_of(float re, float im) {
        float r18501 = re;
        float r18502 = cos(r18501);
        float r18503 = 0.5f;
        float r18504 = r18502 * r18503;
        float r18505 = im;
        float r18506 = 5.0f;
        float r18507 = pow(r18505, r18506);
        float r18508 = 0.0005208333604969084f;
        float r18509 = 0.0416666679084301f;
        float r18510 = r18505 * (r18505 * r18505);
        float r18511 = fma(r18509, r18510, r18505);
        float r18512 = fma(r18507, r18508, r18511);
        float r18513 = -r18512;
        float r18514 = exp(r18505);
        float r18515 = sqrt(r18514);
        float r18516 = -r18505;
        float r18517 = exp(r18516);
        float r18518 = sqrt(r18517);
        float r18519 = r18515 + r18518;
        float r18520 = r18513 * r18519;
        float r18521 = r18504 * r18520;
        return r18521;
}

double f_od(double re, double im) {
        double r18522 = re;
        double r18523 = cos(r18522);
        double r18524 = 0.5;
        double r18525 = r18523 * r18524;
        double r18526 = im;
        double r18527 = 5.0;
        double r18528 = pow(r18526, r18527);
        double r18529 = 0.0005208333604969084;
        double r18530 = 0.0416666679084301;
        double r18531 = r18526 * (r18526 * r18526);
        double r18532 = fma(r18530, r18531, r18526);
        double r18533 = fma(r18528, r18529, r18532);
        double r18534 = -r18533;
        double r18535 = exp(r18526);
        double r18536 = sqrt(r18535);
        double r18537 = -r18526;
        double r18538 = exp(r18537);
        double r18539 = sqrt(r18538);
        double r18540 = r18536 + r18539;
        double r18541 = r18534 * r18540;
        double r18542 = r18525 * r18541;
        return r18542;
}

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 r18543, r18544, r18545, r18546, r18547, r18548, r18549, r18550, r18551, r18552, r18553;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18543, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18544);
        mpfr_init(r18545);
        mpfr_init(r18546);
        mpfr_init_set_str(r18547, "0", 10, MPFR_RNDN);
        mpfr_init(r18548);
        mpfr_init(r18549);
        mpfr_init(r18550);
        mpfr_init(r18551);
        mpfr_init(r18552);
        mpfr_init(r18553);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r18544, re, MPFR_RNDN);
        mpfr_cos(r18545, r18544, MPFR_RNDN);
        mpfr_mul(r18546, r18543, r18545, MPFR_RNDN);
        ;
        mpfr_set_d(r18548, im, MPFR_RNDN);
        mpfr_sub(r18549, r18547, r18548, MPFR_RNDN);
        mpfr_exp(r18550, r18549, MPFR_RNDN);
        mpfr_exp(r18551, r18548, MPFR_RNDN);
        mpfr_sub(r18552, r18550, r18551, MPFR_RNDN);
        mpfr_mul(r18553, r18546, r18552, MPFR_RNDN);
        return mpfr_get_d(r18553, MPFR_RNDN);
}

static mpfr_t r18554, r18555, r18556, r18557, r18558, r18559, r18560, r18561, r18562, r18563, r18564, r18565, r18566, r18567, r18568, r18569, r18570, r18571, r18572, r18573, r18574;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18554);
        mpfr_init(r18555);
        mpfr_init_set_str(r18556, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18557);
        mpfr_init(r18558);
        mpfr_init_set_str(r18559, "5", 10, MPFR_RNDN);
        mpfr_init(r18560);
        mpfr_init_set_str(r18561, "1/1920", 10, MPFR_RNDN);
        mpfr_init_set_str(r18562, "1/24", 10, MPFR_RNDN);
        mpfr_init(r18563);
        mpfr_init(r18564);
        mpfr_init(r18565);
        mpfr_init(r18566);
        mpfr_init(r18567);
        mpfr_init(r18568);
        mpfr_init(r18569);
        mpfr_init(r18570);
        mpfr_init(r18571);
        mpfr_init(r18572);
        mpfr_init(r18573);
        mpfr_init(r18574);
}

double f_fm(double re, double im) {
        mpfr_set_d(r18554, re, MPFR_RNDN);
        mpfr_cos(r18555, r18554, MPFR_RNDN);
        ;
        mpfr_mul(r18557, r18555, r18556, MPFR_RNDN);
        mpfr_set_d(r18558, im, MPFR_RNDN);
        ;
        mpfr_pow(r18560, r18558, r18559, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18563, r18558, r18558, MPFR_RNDN); mpfr_mul(r18563, r18563, r18558, MPFR_RNDN);
        mpfr_fma(r18564, r18562, r18563, r18558, MPFR_RNDN);
        mpfr_fma(r18565, r18560, r18561, r18564, MPFR_RNDN);
        mpfr_neg(r18566, r18565, MPFR_RNDN);
        mpfr_exp(r18567, r18558, MPFR_RNDN);
        mpfr_sqrt(r18568, r18567, MPFR_RNDN);
        mpfr_neg(r18569, r18558, MPFR_RNDN);
        mpfr_exp(r18570, r18569, MPFR_RNDN);
        mpfr_sqrt(r18571, r18570, MPFR_RNDN);
        mpfr_add(r18572, r18568, r18571, MPFR_RNDN);
        mpfr_mul(r18573, r18566, r18572, MPFR_RNDN);
        mpfr_mul(r18574, r18557, r18573, MPFR_RNDN);
        return mpfr_get_d(r18574, MPFR_RNDN);
}

static mpfr_t r18575, r18576, r18577, r18578, r18579, r18580, r18581, r18582, r18583, r18584, r18585, r18586, r18587, r18588, r18589, r18590, r18591, r18592, r18593, r18594, r18595;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18575);
        mpfr_init(r18576);
        mpfr_init_set_str(r18577, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18578);
        mpfr_init(r18579);
        mpfr_init_set_str(r18580, "5", 10, MPFR_RNDN);
        mpfr_init(r18581);
        mpfr_init_set_str(r18582, "1/1920", 10, MPFR_RNDN);
        mpfr_init_set_str(r18583, "1/24", 10, MPFR_RNDN);
        mpfr_init(r18584);
        mpfr_init(r18585);
        mpfr_init(r18586);
        mpfr_init(r18587);
        mpfr_init(r18588);
        mpfr_init(r18589);
        mpfr_init(r18590);
        mpfr_init(r18591);
        mpfr_init(r18592);
        mpfr_init(r18593);
        mpfr_init(r18594);
        mpfr_init(r18595);
}

double f_dm(double re, double im) {
        mpfr_set_d(r18575, re, MPFR_RNDN);
        mpfr_cos(r18576, r18575, MPFR_RNDN);
        ;
        mpfr_mul(r18578, r18576, r18577, MPFR_RNDN);
        mpfr_set_d(r18579, im, MPFR_RNDN);
        ;
        mpfr_pow(r18581, r18579, r18580, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18584, r18579, r18579, MPFR_RNDN); mpfr_mul(r18584, r18584, r18579, MPFR_RNDN);
        mpfr_fma(r18585, r18583, r18584, r18579, MPFR_RNDN);
        mpfr_fma(r18586, r18581, r18582, r18585, MPFR_RNDN);
        mpfr_neg(r18587, r18586, MPFR_RNDN);
        mpfr_exp(r18588, r18579, MPFR_RNDN);
        mpfr_sqrt(r18589, r18588, MPFR_RNDN);
        mpfr_neg(r18590, r18579, MPFR_RNDN);
        mpfr_exp(r18591, r18590, MPFR_RNDN);
        mpfr_sqrt(r18592, r18591, MPFR_RNDN);
        mpfr_add(r18593, r18589, r18592, MPFR_RNDN);
        mpfr_mul(r18594, r18587, r18593, MPFR_RNDN);
        mpfr_mul(r18595, r18578, r18594, MPFR_RNDN);
        return mpfr_get_d(r18595, MPFR_RNDN);
}