#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 r19149 = 0.5f;
        float r19150 = re;
        float r19151 = sin(r19150);
        float r19152 = r19149 * r19151;
        float r19153 = im;
        float r19154 = -r19153;
        float r19155 = exp(r19154);
        float r19156 = exp(r19153);
        float r19157 = r19155 - r19156;
        float r19158 = r19152 * r19157;
        return r19158;
}

double f_id(double re, double im) {
        double r19159 = 0.5;
        double r19160 = re;
        double r19161 = sin(r19160);
        double r19162 = r19159 * r19161;
        double r19163 = im;
        double r19164 = -r19163;
        double r19165 = exp(r19164);
        double r19166 = exp(r19163);
        double r19167 = r19165 - r19166;
        double r19168 = r19162 * r19167;
        return r19168;
}


double f_of(float re, float im) {
        float r19169 = 0.0005208333333333333f;
        float r19170 = im;
        float r19171 = 5.0f;
        float r19172 = pow(r19170, r19171);
        float r19173 = r19170 * (r19170 * r19170);
        float r19174 = 0.041666666666666664f;
        float r19175 = fma(r19173, r19174, r19170);
        float r19176 = fma(r19169, r19172, r19175);
        float r19177 = -r19176;
        float r19178 = re;
        float r19179 = sin(r19178);
        float r19180 = 0.5f;
        float r19181 = r19179 * r19180;
        float r19182 = -r19170;
        float r19183 = exp(r19182);
        float r19184 = sqrt(r19183);
        float r19185 = exp(r19170);
        float r19186 = sqrt(r19185);
        float r19187 = r19186 * r19181;
        float r19188 = fma(r19181, r19184, r19187);
        float r19189 = r19177 * r19188;
        return r19189;
}

double f_od(double re, double im) {
        double r19190 = 0.0005208333333333333;
        double r19191 = im;
        double r19192 = 5.0;
        double r19193 = pow(r19191, r19192);
        double r19194 = r19191 * (r19191 * r19191);
        double r19195 = 0.041666666666666664;
        double r19196 = fma(r19194, r19195, r19191);
        double r19197 = fma(r19190, r19193, r19196);
        double r19198 = -r19197;
        double r19199 = re;
        double r19200 = sin(r19199);
        double r19201 = 0.5;
        double r19202 = r19200 * r19201;
        double r19203 = -r19191;
        double r19204 = exp(r19203);
        double r19205 = sqrt(r19204);
        double r19206 = exp(r19191);
        double r19207 = sqrt(r19206);
        double r19208 = r19207 * r19202;
        double r19209 = fma(r19202, r19205, r19208);
        double r19210 = r19198 * r19209;
        return r19210;
}

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 r19211, r19212, r19213, r19214, r19215, r19216, r19217, r19218, r19219, r19220;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19211, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19212);
        mpfr_init(r19213);
        mpfr_init(r19214);
        mpfr_init(r19215);
        mpfr_init(r19216);
        mpfr_init(r19217);
        mpfr_init(r19218);
        mpfr_init(r19219);
        mpfr_init(r19220);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r19212, re, MPFR_RNDN);
        mpfr_sin(r19213, r19212, MPFR_RNDN);
        mpfr_mul(r19214, r19211, r19213, MPFR_RNDN);
        mpfr_set_d(r19215, im, MPFR_RNDN);
        mpfr_neg(r19216, r19215, MPFR_RNDN);
        mpfr_exp(r19217, r19216, MPFR_RNDN);
        mpfr_exp(r19218, r19215, MPFR_RNDN);
        mpfr_sub(r19219, r19217, r19218, MPFR_RNDN);
        mpfr_mul(r19220, r19214, r19219, MPFR_RNDN);
        return mpfr_get_d(r19220, MPFR_RNDN);
}

static mpfr_t r19221, r19222, r19223, r19224, r19225, r19226, r19227, r19228, r19229, r19230, r19231, r19232, r19233, r19234, r19235, r19236, r19237, r19238, r19239, r19240, r19241;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19221, "1/1920", 10, MPFR_RNDN);
        mpfr_init(r19222);
        mpfr_init_set_str(r19223, "5", 10, MPFR_RNDN);
        mpfr_init(r19224);
        mpfr_init(r19225);
        mpfr_init_set_str(r19226, "1/24", 10, MPFR_RNDN);
        mpfr_init(r19227);
        mpfr_init(r19228);
        mpfr_init(r19229);
        mpfr_init(r19230);
        mpfr_init(r19231);
        mpfr_init_set_str(r19232, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19233);
        mpfr_init(r19234);
        mpfr_init(r19235);
        mpfr_init(r19236);
        mpfr_init(r19237);
        mpfr_init(r19238);
        mpfr_init(r19239);
        mpfr_init(r19240);
        mpfr_init(r19241);
}

double f_fm(double re, double im) {
        ;
        mpfr_set_d(r19222, im, MPFR_RNDN);
        ;
        mpfr_pow(r19224, r19222, r19223, MPFR_RNDN);
        mpfr_mul(r19225, r19222, r19222, MPFR_RNDN); mpfr_mul(r19225, r19225, r19222, MPFR_RNDN);
        ;
        mpfr_fma(r19227, r19225, r19226, r19222, MPFR_RNDN);
        mpfr_fma(r19228, r19221, r19224, r19227, MPFR_RNDN);
        mpfr_neg(r19229, r19228, MPFR_RNDN);
        mpfr_set_d(r19230, re, MPFR_RNDN);
        mpfr_sin(r19231, r19230, MPFR_RNDN);
        ;
        mpfr_mul(r19233, r19231, r19232, MPFR_RNDN);
        mpfr_neg(r19234, r19222, MPFR_RNDN);
        mpfr_exp(r19235, r19234, MPFR_RNDN);
        mpfr_sqrt(r19236, r19235, MPFR_RNDN);
        mpfr_exp(r19237, r19222, MPFR_RNDN);
        mpfr_sqrt(r19238, r19237, MPFR_RNDN);
        mpfr_mul(r19239, r19238, r19233, MPFR_RNDN);
        mpfr_fma(r19240, r19233, r19236, r19239, MPFR_RNDN);
        mpfr_mul(r19241, r19229, r19240, MPFR_RNDN);
        return mpfr_get_d(r19241, MPFR_RNDN);
}

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

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19242, "1/1920", 10, MPFR_RNDN);
        mpfr_init(r19243);
        mpfr_init_set_str(r19244, "5", 10, MPFR_RNDN);
        mpfr_init(r19245);
        mpfr_init(r19246);
        mpfr_init_set_str(r19247, "1/24", 10, MPFR_RNDN);
        mpfr_init(r19248);
        mpfr_init(r19249);
        mpfr_init(r19250);
        mpfr_init(r19251);
        mpfr_init(r19252);
        mpfr_init_set_str(r19253, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19254);
        mpfr_init(r19255);
        mpfr_init(r19256);
        mpfr_init(r19257);
        mpfr_init(r19258);
        mpfr_init(r19259);
        mpfr_init(r19260);
        mpfr_init(r19261);
        mpfr_init(r19262);
}

double f_dm(double re, double im) {
        ;
        mpfr_set_d(r19243, im, MPFR_RNDN);
        ;
        mpfr_pow(r19245, r19243, r19244, MPFR_RNDN);
        mpfr_mul(r19246, r19243, r19243, MPFR_RNDN); mpfr_mul(r19246, r19246, r19243, MPFR_RNDN);
        ;
        mpfr_fma(r19248, r19246, r19247, r19243, MPFR_RNDN);
        mpfr_fma(r19249, r19242, r19245, r19248, MPFR_RNDN);
        mpfr_neg(r19250, r19249, MPFR_RNDN);
        mpfr_set_d(r19251, re, MPFR_RNDN);
        mpfr_sin(r19252, r19251, MPFR_RNDN);
        ;
        mpfr_mul(r19254, r19252, r19253, MPFR_RNDN);
        mpfr_neg(r19255, r19243, MPFR_RNDN);
        mpfr_exp(r19256, r19255, MPFR_RNDN);
        mpfr_sqrt(r19257, r19256, MPFR_RNDN);
        mpfr_exp(r19258, r19243, MPFR_RNDN);
        mpfr_sqrt(r19259, r19258, MPFR_RNDN);
        mpfr_mul(r19260, r19259, r19254, MPFR_RNDN);
        mpfr_fma(r19261, r19254, r19257, r19260, MPFR_RNDN);
        mpfr_mul(r19262, r19250, r19261, MPFR_RNDN);
        return mpfr_get_d(r19262, MPFR_RNDN);
}