#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 r20151 = 0.5f;
        float r20152 = re;
        float r20153 = cos(r20152);
        float r20154 = r20151 * r20153;
        float r20155 = 0.0f;
        float r20156 = im;
        float r20157 = r20155 - r20156;
        float r20158 = exp(r20157);
        float r20159 = exp(r20156);
        float r20160 = r20158 - r20159;
        float r20161 = r20154 * r20160;
        return r20161;
}

double f_id(double re, double im) {
        double r20162 = 0.5;
        double r20163 = re;
        double r20164 = cos(r20163);
        double r20165 = r20162 * r20164;
        double r20166 = 0.0;
        double r20167 = im;
        double r20168 = r20166 - r20167;
        double r20169 = exp(r20168);
        double r20170 = exp(r20167);
        double r20171 = r20169 - r20170;
        double r20172 = r20165 * r20171;
        return r20172;
}


double f_of(float re, float im) {
        float r20173 = re;
        float r20174 = cos(r20173);
        float r20175 = 0.5f;
        float r20176 = r20174 * r20175;
        float r20177 = 0.0416666679084301f;
        float r20178 = -r20177;
        float r20179 = im;
        float r20180 = r20179 * (r20179 * r20179);
        float r20181 = r20178 * r20180;
        float r20182 = -r20179;
        float r20183 = r20181 + r20182;
        float r20184 = 0.0005208333604969084f;
        float r20185 = -r20184;
        float r20186 = 5.0f;
        float r20187 = pow(r20179, r20186);
        float r20188 = r20185 * r20187;
        float r20189 = r20183 + r20188;
        float r20190 = exp(r20179);
        float r20191 = sqrt(r20190);
        float r20192 = exp(r20182);
        float r20193 = sqrt(r20192);
        float r20194 = r20191 + r20193;
        float r20195 = r20189 * r20194;
        float r20196 = r20176 * r20195;
        return r20196;
}

double f_od(double re, double im) {
        double r20197 = re;
        double r20198 = cos(r20197);
        double r20199 = 0.5;
        double r20200 = r20198 * r20199;
        double r20201 = 0.0416666679084301;
        double r20202 = -r20201;
        double r20203 = im;
        double r20204 = r20203 * (r20203 * r20203);
        double r20205 = r20202 * r20204;
        double r20206 = -r20203;
        double r20207 = r20205 + r20206;
        double r20208 = 0.0005208333604969084;
        double r20209 = -r20208;
        double r20210 = 5.0;
        double r20211 = pow(r20203, r20210);
        double r20212 = r20209 * r20211;
        double r20213 = r20207 + r20212;
        double r20214 = exp(r20203);
        double r20215 = sqrt(r20214);
        double r20216 = exp(r20206);
        double r20217 = sqrt(r20216);
        double r20218 = r20215 + r20217;
        double r20219 = r20213 * r20218;
        double r20220 = r20200 * r20219;
        return r20220;
}

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 r20221, r20222, r20223, r20224, r20225, r20226, r20227, r20228, r20229, r20230, r20231;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r20221, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20222);
        mpfr_init(r20223);
        mpfr_init(r20224);
        mpfr_init_set_str(r20225, "0", 10, MPFR_RNDN);
        mpfr_init(r20226);
        mpfr_init(r20227);
        mpfr_init(r20228);
        mpfr_init(r20229);
        mpfr_init(r20230);
        mpfr_init(r20231);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r20222, re, MPFR_RNDN);
        mpfr_cos(r20223, r20222, MPFR_RNDN);
        mpfr_mul(r20224, r20221, r20223, MPFR_RNDN);
        ;
        mpfr_set_d(r20226, im, MPFR_RNDN);
        mpfr_sub(r20227, r20225, r20226, MPFR_RNDN);
        mpfr_exp(r20228, r20227, MPFR_RNDN);
        mpfr_exp(r20229, r20226, MPFR_RNDN);
        mpfr_sub(r20230, r20228, r20229, MPFR_RNDN);
        mpfr_mul(r20231, r20224, r20230, MPFR_RNDN);
        return mpfr_get_d(r20231, MPFR_RNDN);
}

static mpfr_t r20232, r20233, r20234, r20235, r20236, r20237, r20238, r20239, r20240, r20241, r20242, r20243, r20244, r20245, r20246, r20247, r20248, r20249, r20250, r20251, r20252, r20253, r20254, r20255;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20232);
        mpfr_init(r20233);
        mpfr_init_set_str(r20234, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20235);
        mpfr_init_set_str(r20236, "1/24", 10, MPFR_RNDN);
        mpfr_init(r20237);
        mpfr_init(r20238);
        mpfr_init(r20239);
        mpfr_init(r20240);
        mpfr_init(r20241);
        mpfr_init(r20242);
        mpfr_init_set_str(r20243, "1/1920", 10, MPFR_RNDN);
        mpfr_init(r20244);
        mpfr_init_set_str(r20245, "5", 10, MPFR_RNDN);
        mpfr_init(r20246);
        mpfr_init(r20247);
        mpfr_init(r20248);
        mpfr_init(r20249);
        mpfr_init(r20250);
        mpfr_init(r20251);
        mpfr_init(r20252);
        mpfr_init(r20253);
        mpfr_init(r20254);
        mpfr_init(r20255);
}

double f_fm(double re, double im) {
        mpfr_set_d(r20232, re, MPFR_RNDN);
        mpfr_cos(r20233, r20232, MPFR_RNDN);
        ;
        mpfr_mul(r20235, r20233, r20234, MPFR_RNDN);
        ;
        mpfr_neg(r20237, r20236, MPFR_RNDN);
        mpfr_set_d(r20238, im, MPFR_RNDN);
        mpfr_mul(r20239, r20238, r20238, MPFR_RNDN); mpfr_mul(r20239, r20239, r20238, MPFR_RNDN);
        mpfr_mul(r20240, r20237, r20239, MPFR_RNDN);
        mpfr_neg(r20241, r20238, MPFR_RNDN);
        mpfr_add(r20242, r20240, r20241, MPFR_RNDN);
        ;
        mpfr_neg(r20244, r20243, MPFR_RNDN);
        ;
        mpfr_pow(r20246, r20238, r20245, MPFR_RNDN);
        mpfr_mul(r20247, r20244, r20246, MPFR_RNDN);
        mpfr_add(r20248, r20242, r20247, MPFR_RNDN);
        mpfr_exp(r20249, r20238, MPFR_RNDN);
        mpfr_sqrt(r20250, r20249, MPFR_RNDN);
        mpfr_exp(r20251, r20241, MPFR_RNDN);
        mpfr_sqrt(r20252, r20251, MPFR_RNDN);
        mpfr_add(r20253, r20250, r20252, MPFR_RNDN);
        mpfr_mul(r20254, r20248, r20253, MPFR_RNDN);
        mpfr_mul(r20255, r20235, r20254, MPFR_RNDN);
        return mpfr_get_d(r20255, MPFR_RNDN);
}

static mpfr_t r20256, r20257, r20258, r20259, r20260, r20261, r20262, r20263, r20264, r20265, r20266, r20267, r20268, r20269, r20270, r20271, r20272, r20273, r20274, r20275, r20276, r20277, r20278, r20279;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20256);
        mpfr_init(r20257);
        mpfr_init_set_str(r20258, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20259);
        mpfr_init_set_str(r20260, "1/24", 10, MPFR_RNDN);
        mpfr_init(r20261);
        mpfr_init(r20262);
        mpfr_init(r20263);
        mpfr_init(r20264);
        mpfr_init(r20265);
        mpfr_init(r20266);
        mpfr_init_set_str(r20267, "1/1920", 10, MPFR_RNDN);
        mpfr_init(r20268);
        mpfr_init_set_str(r20269, "5", 10, MPFR_RNDN);
        mpfr_init(r20270);
        mpfr_init(r20271);
        mpfr_init(r20272);
        mpfr_init(r20273);
        mpfr_init(r20274);
        mpfr_init(r20275);
        mpfr_init(r20276);
        mpfr_init(r20277);
        mpfr_init(r20278);
        mpfr_init(r20279);
}

double f_dm(double re, double im) {
        mpfr_set_d(r20256, re, MPFR_RNDN);
        mpfr_cos(r20257, r20256, MPFR_RNDN);
        ;
        mpfr_mul(r20259, r20257, r20258, MPFR_RNDN);
        ;
        mpfr_neg(r20261, r20260, MPFR_RNDN);
        mpfr_set_d(r20262, im, MPFR_RNDN);
        mpfr_mul(r20263, r20262, r20262, MPFR_RNDN); mpfr_mul(r20263, r20263, r20262, MPFR_RNDN);
        mpfr_mul(r20264, r20261, r20263, MPFR_RNDN);
        mpfr_neg(r20265, r20262, MPFR_RNDN);
        mpfr_add(r20266, r20264, r20265, MPFR_RNDN);
        ;
        mpfr_neg(r20268, r20267, MPFR_RNDN);
        ;
        mpfr_pow(r20270, r20262, r20269, MPFR_RNDN);
        mpfr_mul(r20271, r20268, r20270, MPFR_RNDN);
        mpfr_add(r20272, r20266, r20271, MPFR_RNDN);
        mpfr_exp(r20273, r20262, MPFR_RNDN);
        mpfr_sqrt(r20274, r20273, MPFR_RNDN);
        mpfr_exp(r20275, r20265, MPFR_RNDN);
        mpfr_sqrt(r20276, r20275, MPFR_RNDN);
        mpfr_add(r20277, r20274, r20276, MPFR_RNDN);
        mpfr_mul(r20278, r20272, r20277, MPFR_RNDN);
        mpfr_mul(r20279, r20259, r20278, MPFR_RNDN);
        return mpfr_get_d(r20279, MPFR_RNDN);
}