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

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

double f_if(float re, float im) {
        float r18337 = 0.5f;
        float r18338 = re;
        float r18339 = sin(r18338);
        float r18340 = r18337 * r18339;
        float r18341 = im;
        float r18342 = -r18341;
        float r18343 = exp(r18342);
        float r18344 = exp(r18341);
        float r18345 = r18343 - r18344;
        float r18346 = r18340 * r18345;
        return r18346;
}

double f_id(double re, double im) {
        double r18347 = 0.5;
        double r18348 = re;
        double r18349 = sin(r18348);
        double r18350 = r18347 * r18349;
        double r18351 = im;
        double r18352 = -r18351;
        double r18353 = exp(r18352);
        double r18354 = exp(r18351);
        double r18355 = r18353 - r18354;
        double r18356 = r18350 * r18355;
        return r18356;
}


double f_of(float re, float im) {
        float r18357 = im;
        float r18358 = r18357 * (r18357 * r18357);
        float r18359 = 0.3333333432674408f;
        float r18360 = 5.0f;
        float r18361 = pow(r18357, r18360);
        float r18362 = 0.01666666753590107f;
        float r18363 = 2.0f;
        float r18364 = r18357 * r18363;
        float r18365 = fma(r18361, r18362, r18364);
        float r18366 = fma(r18358, r18359, r18365);
        float r18367 = re;
        float r18368 = sin(r18367);
        float r18369 = 0.5f;
        float r18370 = -r18369;
        float r18371 = r18368 * r18370;
        float r18372 = r18366 * r18371;
        return r18372;
}

double f_od(double re, double im) {
        double r18373 = im;
        double r18374 = r18373 * (r18373 * r18373);
        double r18375 = 0.3333333432674408;
        double r18376 = 5.0;
        double r18377 = pow(r18373, r18376);
        double r18378 = 0.01666666753590107;
        double r18379 = 2.0;
        double r18380 = r18373 * r18379;
        double r18381 = fma(r18377, r18378, r18380);
        double r18382 = fma(r18374, r18375, r18381);
        double r18383 = re;
        double r18384 = sin(r18383);
        double r18385 = 0.5;
        double r18386 = -r18385;
        double r18387 = r18384 * r18386;
        double r18388 = r18382 * r18387;
        return r18388;
}

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 r18389, r18390, r18391, r18392, r18393, r18394, r18395, r18396, r18397, r18398;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18389, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18390);
        mpfr_init(r18391);
        mpfr_init(r18392);
        mpfr_init(r18393);
        mpfr_init(r18394);
        mpfr_init(r18395);
        mpfr_init(r18396);
        mpfr_init(r18397);
        mpfr_init(r18398);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r18390, re, MPFR_RNDN);
        mpfr_sin(r18391, r18390, MPFR_RNDN);
        mpfr_mul(r18392, r18389, r18391, MPFR_RNDN);
        mpfr_set_d(r18393, im, MPFR_RNDN);
        mpfr_neg(r18394, r18393, MPFR_RNDN);
        mpfr_exp(r18395, r18394, MPFR_RNDN);
        mpfr_exp(r18396, r18393, MPFR_RNDN);
        mpfr_sub(r18397, r18395, r18396, MPFR_RNDN);
        mpfr_mul(r18398, r18392, r18397, MPFR_RNDN);
        return mpfr_get_d(r18398, MPFR_RNDN);
}

static mpfr_t r18399, r18400, r18401, r18402, r18403, r18404, r18405, r18406, r18407, r18408, r18409, r18410, r18411, r18412, r18413, r18414;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18399);
        mpfr_init(r18400);
        mpfr_init_set_str(r18401, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r18402, "5", 10, MPFR_RNDN);
        mpfr_init(r18403);
        mpfr_init_set_str(r18404, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r18405, "2", 10, MPFR_RNDN);
        mpfr_init(r18406);
        mpfr_init(r18407);
        mpfr_init(r18408);
        mpfr_init(r18409);
        mpfr_init(r18410);
        mpfr_init_set_str(r18411, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18412);
        mpfr_init(r18413);
        mpfr_init(r18414);
}

double f_fm(double re, double im) {
        mpfr_set_d(r18399, im, MPFR_RNDN);
        mpfr_mul(r18400, r18399, r18399, MPFR_RNDN); mpfr_mul(r18400, r18400, r18399, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r18403, r18399, r18402, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18406, r18399, r18405, MPFR_RNDN);
        mpfr_fma(r18407, r18403, r18404, r18406, MPFR_RNDN);
        mpfr_fma(r18408, r18400, r18401, r18407, MPFR_RNDN);
        mpfr_set_d(r18409, re, MPFR_RNDN);
        mpfr_sin(r18410, r18409, MPFR_RNDN);
        ;
        mpfr_neg(r18412, r18411, MPFR_RNDN);
        mpfr_mul(r18413, r18410, r18412, MPFR_RNDN);
        mpfr_mul(r18414, r18408, r18413, MPFR_RNDN);
        return mpfr_get_d(r18414, MPFR_RNDN);
}

static mpfr_t r18415, r18416, r18417, r18418, r18419, r18420, r18421, r18422, r18423, r18424, r18425, r18426, r18427, r18428, r18429, r18430;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18415);
        mpfr_init(r18416);
        mpfr_init_set_str(r18417, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r18418, "5", 10, MPFR_RNDN);
        mpfr_init(r18419);
        mpfr_init_set_str(r18420, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r18421, "2", 10, MPFR_RNDN);
        mpfr_init(r18422);
        mpfr_init(r18423);
        mpfr_init(r18424);
        mpfr_init(r18425);
        mpfr_init(r18426);
        mpfr_init_set_str(r18427, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18428);
        mpfr_init(r18429);
        mpfr_init(r18430);
}

double f_dm(double re, double im) {
        mpfr_set_d(r18415, im, MPFR_RNDN);
        mpfr_mul(r18416, r18415, r18415, MPFR_RNDN); mpfr_mul(r18416, r18416, r18415, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r18419, r18415, r18418, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18422, r18415, r18421, MPFR_RNDN);
        mpfr_fma(r18423, r18419, r18420, r18422, MPFR_RNDN);
        mpfr_fma(r18424, r18416, r18417, r18423, MPFR_RNDN);
        mpfr_set_d(r18425, re, MPFR_RNDN);
        mpfr_sin(r18426, r18425, MPFR_RNDN);
        ;
        mpfr_neg(r18428, r18427, MPFR_RNDN);
        mpfr_mul(r18429, r18426, r18428, MPFR_RNDN);
        mpfr_mul(r18430, r18424, r18429, MPFR_RNDN);
        return mpfr_get_d(r18430, MPFR_RNDN);
}