#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 r19185 = 0.5f;
        float r19186 = re;
        float r19187 = sin(r19186);
        float r19188 = r19185 * r19187;
        float r19189 = im;
        float r19190 = -r19189;
        float r19191 = exp(r19190);
        float r19192 = exp(r19189);
        float r19193 = r19191 - r19192;
        float r19194 = r19188 * r19193;
        return r19194;
}

double f_id(double re, double im) {
        double r19195 = 0.5;
        double r19196 = re;
        double r19197 = sin(r19196);
        double r19198 = r19195 * r19197;
        double r19199 = im;
        double r19200 = -r19199;
        double r19201 = exp(r19200);
        double r19202 = exp(r19199);
        double r19203 = r19201 - r19202;
        double r19204 = r19198 * r19203;
        return r19204;
}


double f_of(float re, float im) {
        float r19205 = im;
        float r19206 = r19205 * (r19205 * r19205);
        float r19207 = 0.3333333333333333f;
        float r19208 = 5.0f;
        float r19209 = pow(r19205, r19208);
        float r19210 = 0.016666666666666666f;
        float r19211 = 2.0f;
        float r19212 = r19205 * r19211;
        float r19213 = fma(r19209, r19210, r19212);
        float r19214 = fma(r19206, r19207, r19213);
        float r19215 = re;
        float r19216 = sin(r19215);
        float r19217 = 0.5f;
        float r19218 = -r19217;
        float r19219 = r19216 * r19218;
        float r19220 = r19214 * r19219;
        return r19220;
}

double f_od(double re, double im) {
        double r19221 = im;
        double r19222 = r19221 * (r19221 * r19221);
        double r19223 = 0.3333333333333333;
        double r19224 = 5.0;
        double r19225 = pow(r19221, r19224);
        double r19226 = 0.016666666666666666;
        double r19227 = 2.0;
        double r19228 = r19221 * r19227;
        double r19229 = fma(r19225, r19226, r19228);
        double r19230 = fma(r19222, r19223, r19229);
        double r19231 = re;
        double r19232 = sin(r19231);
        double r19233 = 0.5;
        double r19234 = -r19233;
        double r19235 = r19232 * r19234;
        double r19236 = r19230 * r19235;
        return r19236;
}

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 r19237, r19238, r19239, r19240, r19241, r19242, r19243, r19244, r19245, r19246;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19237, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19238);
        mpfr_init(r19239);
        mpfr_init(r19240);
        mpfr_init(r19241);
        mpfr_init(r19242);
        mpfr_init(r19243);
        mpfr_init(r19244);
        mpfr_init(r19245);
        mpfr_init(r19246);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r19238, re, MPFR_RNDN);
        mpfr_sin(r19239, r19238, MPFR_RNDN);
        mpfr_mul(r19240, r19237, r19239, MPFR_RNDN);
        mpfr_set_d(r19241, im, MPFR_RNDN);
        mpfr_neg(r19242, r19241, MPFR_RNDN);
        mpfr_exp(r19243, r19242, MPFR_RNDN);
        mpfr_exp(r19244, r19241, MPFR_RNDN);
        mpfr_sub(r19245, r19243, r19244, MPFR_RNDN);
        mpfr_mul(r19246, r19240, r19245, MPFR_RNDN);
        return mpfr_get_d(r19246, MPFR_RNDN);
}

static mpfr_t r19247, r19248, r19249, r19250, r19251, r19252, r19253, r19254, r19255, r19256, r19257, r19258, r19259, r19260, r19261, r19262;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19247);
        mpfr_init(r19248);
        mpfr_init_set_str(r19249, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r19250, "5", 10, MPFR_RNDN);
        mpfr_init(r19251);
        mpfr_init_set_str(r19252, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r19253, "2", 10, MPFR_RNDN);
        mpfr_init(r19254);
        mpfr_init(r19255);
        mpfr_init(r19256);
        mpfr_init(r19257);
        mpfr_init(r19258);
        mpfr_init_set_str(r19259, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19260);
        mpfr_init(r19261);
        mpfr_init(r19262);
}

double f_fm(double re, double im) {
        mpfr_set_d(r19247, im, MPFR_RNDN);
        mpfr_mul(r19248, r19247, r19247, MPFR_RNDN); mpfr_mul(r19248, r19248, r19247, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r19251, r19247, r19250, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r19254, r19247, r19253, MPFR_RNDN);
        mpfr_fma(r19255, r19251, r19252, r19254, MPFR_RNDN);
        mpfr_fma(r19256, r19248, r19249, r19255, MPFR_RNDN);
        mpfr_set_d(r19257, re, MPFR_RNDN);
        mpfr_sin(r19258, r19257, MPFR_RNDN);
        ;
        mpfr_neg(r19260, r19259, MPFR_RNDN);
        mpfr_mul(r19261, r19258, r19260, MPFR_RNDN);
        mpfr_mul(r19262, r19256, r19261, MPFR_RNDN);
        return mpfr_get_d(r19262, MPFR_RNDN);
}

static mpfr_t r19263, r19264, r19265, r19266, r19267, r19268, r19269, r19270, r19271, r19272, r19273, r19274, r19275, r19276, r19277, r19278;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19263);
        mpfr_init(r19264);
        mpfr_init_set_str(r19265, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r19266, "5", 10, MPFR_RNDN);
        mpfr_init(r19267);
        mpfr_init_set_str(r19268, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r19269, "2", 10, MPFR_RNDN);
        mpfr_init(r19270);
        mpfr_init(r19271);
        mpfr_init(r19272);
        mpfr_init(r19273);
        mpfr_init(r19274);
        mpfr_init_set_str(r19275, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19276);
        mpfr_init(r19277);
        mpfr_init(r19278);
}

double f_dm(double re, double im) {
        mpfr_set_d(r19263, im, MPFR_RNDN);
        mpfr_mul(r19264, r19263, r19263, MPFR_RNDN); mpfr_mul(r19264, r19264, r19263, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r19267, r19263, r19266, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r19270, r19263, r19269, MPFR_RNDN);
        mpfr_fma(r19271, r19267, r19268, r19270, MPFR_RNDN);
        mpfr_fma(r19272, r19264, r19265, r19271, MPFR_RNDN);
        mpfr_set_d(r19273, re, MPFR_RNDN);
        mpfr_sin(r19274, r19273, MPFR_RNDN);
        ;
        mpfr_neg(r19276, r19275, MPFR_RNDN);
        mpfr_mul(r19277, r19274, r19276, MPFR_RNDN);
        mpfr_mul(r19278, r19272, r19277, MPFR_RNDN);
        return mpfr_get_d(r19278, MPFR_RNDN);
}