#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 r19915 = 0.5f;
        float r19916 = re;
        float r19917 = cos(r19916);
        float r19918 = r19915 * r19917;
        float r19919 = 0.0f;
        float r19920 = im;
        float r19921 = r19919 - r19920;
        float r19922 = exp(r19921);
        float r19923 = exp(r19920);
        float r19924 = r19922 - r19923;
        float r19925 = r19918 * r19924;
        return r19925;
}

double f_id(double re, double im) {
        double r19926 = 0.5;
        double r19927 = re;
        double r19928 = cos(r19927);
        double r19929 = r19926 * r19928;
        double r19930 = 0.0;
        double r19931 = im;
        double r19932 = r19930 - r19931;
        double r19933 = exp(r19932);
        double r19934 = exp(r19931);
        double r19935 = r19933 - r19934;
        double r19936 = r19929 * r19935;
        return r19936;
}


double f_of(float re, float im) {
        float r19937 = re;
        float r19938 = cos(r19937);
        float r19939 = 0.5f;
        float r19940 = r19938 * r19939;
        float r19941 = 0.0416666679084301f;
        float r19942 = -r19941;
        float r19943 = im;
        float r19944 = r19943 * (r19943 * r19943);
        float r19945 = r19942 * r19944;
        float r19946 = -r19943;
        float r19947 = r19945 + r19946;
        float r19948 = 0.0005208333604969084f;
        float r19949 = -r19948;
        float r19950 = 5.0f;
        float r19951 = pow(r19943, r19950);
        float r19952 = r19949 * r19951;
        float r19953 = r19947 + r19952;
        float r19954 = exp(r19943);
        float r19955 = sqrt(r19954);
        float r19956 = exp(r19946);
        float r19957 = sqrt(r19956);
        float r19958 = r19955 + r19957;
        float r19959 = r19953 * r19958;
        float r19960 = r19940 * r19959;
        return r19960;
}

double f_od(double re, double im) {
        double r19961 = re;
        double r19962 = cos(r19961);
        double r19963 = 0.5;
        double r19964 = r19962 * r19963;
        double r19965 = 0.0416666679084301;
        double r19966 = -r19965;
        double r19967 = im;
        double r19968 = r19967 * (r19967 * r19967);
        double r19969 = r19966 * r19968;
        double r19970 = -r19967;
        double r19971 = r19969 + r19970;
        double r19972 = 0.0005208333604969084;
        double r19973 = -r19972;
        double r19974 = 5.0;
        double r19975 = pow(r19967, r19974);
        double r19976 = r19973 * r19975;
        double r19977 = r19971 + r19976;
        double r19978 = exp(r19967);
        double r19979 = sqrt(r19978);
        double r19980 = exp(r19970);
        double r19981 = sqrt(r19980);
        double r19982 = r19979 + r19981;
        double r19983 = r19977 * r19982;
        double r19984 = r19964 * r19983;
        return r19984;
}

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 r19985, r19986, r19987, r19988, r19989, r19990, r19991, r19992, r19993, r19994, r19995;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19985, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19986);
        mpfr_init(r19987);
        mpfr_init(r19988);
        mpfr_init_set_str(r19989, "0", 10, MPFR_RNDN);
        mpfr_init(r19990);
        mpfr_init(r19991);
        mpfr_init(r19992);
        mpfr_init(r19993);
        mpfr_init(r19994);
        mpfr_init(r19995);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r19986, re, MPFR_RNDN);
        mpfr_cos(r19987, r19986, MPFR_RNDN);
        mpfr_mul(r19988, r19985, r19987, MPFR_RNDN);
        ;
        mpfr_set_d(r19990, im, MPFR_RNDN);
        mpfr_sub(r19991, r19989, r19990, MPFR_RNDN);
        mpfr_exp(r19992, r19991, MPFR_RNDN);
        mpfr_exp(r19993, r19990, MPFR_RNDN);
        mpfr_sub(r19994, r19992, r19993, MPFR_RNDN);
        mpfr_mul(r19995, r19988, r19994, MPFR_RNDN);
        return mpfr_get_d(r19995, MPFR_RNDN);
}

static mpfr_t r19996, r19997, r19998, r19999, r20000, r20001, r20002, r20003, r20004, r20005, r20006, r20007, r20008, r20009, r20010, r20011, r20012, r20013, r20014, r20015, r20016, r20017, r20018, r20019;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19996);
        mpfr_init(r19997);
        mpfr_init_set_str(r19998, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19999);
        mpfr_init_set_str(r20000, "1/24", 10, MPFR_RNDN);
        mpfr_init(r20001);
        mpfr_init(r20002);
        mpfr_init(r20003);
        mpfr_init(r20004);
        mpfr_init(r20005);
        mpfr_init(r20006);
        mpfr_init_set_str(r20007, "1/1920", 10, MPFR_RNDN);
        mpfr_init(r20008);
        mpfr_init_set_str(r20009, "5", 10, MPFR_RNDN);
        mpfr_init(r20010);
        mpfr_init(r20011);
        mpfr_init(r20012);
        mpfr_init(r20013);
        mpfr_init(r20014);
        mpfr_init(r20015);
        mpfr_init(r20016);
        mpfr_init(r20017);
        mpfr_init(r20018);
        mpfr_init(r20019);
}

double f_fm(double re, double im) {
        mpfr_set_d(r19996, re, MPFR_RNDN);
        mpfr_cos(r19997, r19996, MPFR_RNDN);
        ;
        mpfr_mul(r19999, r19997, r19998, MPFR_RNDN);
        ;
        mpfr_neg(r20001, r20000, MPFR_RNDN);
        mpfr_set_d(r20002, im, MPFR_RNDN);
        mpfr_mul(r20003, r20002, r20002, MPFR_RNDN); mpfr_mul(r20003, r20003, r20002, MPFR_RNDN);
        mpfr_mul(r20004, r20001, r20003, MPFR_RNDN);
        mpfr_neg(r20005, r20002, MPFR_RNDN);
        mpfr_add(r20006, r20004, r20005, MPFR_RNDN);
        ;
        mpfr_neg(r20008, r20007, MPFR_RNDN);
        ;
        mpfr_pow(r20010, r20002, r20009, MPFR_RNDN);
        mpfr_mul(r20011, r20008, r20010, MPFR_RNDN);
        mpfr_add(r20012, r20006, r20011, MPFR_RNDN);
        mpfr_exp(r20013, r20002, MPFR_RNDN);
        mpfr_sqrt(r20014, r20013, MPFR_RNDN);
        mpfr_exp(r20015, r20005, MPFR_RNDN);
        mpfr_sqrt(r20016, r20015, MPFR_RNDN);
        mpfr_add(r20017, r20014, r20016, MPFR_RNDN);
        mpfr_mul(r20018, r20012, r20017, MPFR_RNDN);
        mpfr_mul(r20019, r19999, r20018, MPFR_RNDN);
        return mpfr_get_d(r20019, MPFR_RNDN);
}

static mpfr_t r20020, r20021, r20022, r20023, r20024, r20025, r20026, r20027, r20028, r20029, r20030, r20031, r20032, r20033, r20034, r20035, r20036, r20037, r20038, r20039, r20040, r20041, r20042, r20043;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20020);
        mpfr_init(r20021);
        mpfr_init_set_str(r20022, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20023);
        mpfr_init_set_str(r20024, "1/24", 10, MPFR_RNDN);
        mpfr_init(r20025);
        mpfr_init(r20026);
        mpfr_init(r20027);
        mpfr_init(r20028);
        mpfr_init(r20029);
        mpfr_init(r20030);
        mpfr_init_set_str(r20031, "1/1920", 10, MPFR_RNDN);
        mpfr_init(r20032);
        mpfr_init_set_str(r20033, "5", 10, MPFR_RNDN);
        mpfr_init(r20034);
        mpfr_init(r20035);
        mpfr_init(r20036);
        mpfr_init(r20037);
        mpfr_init(r20038);
        mpfr_init(r20039);
        mpfr_init(r20040);
        mpfr_init(r20041);
        mpfr_init(r20042);
        mpfr_init(r20043);
}

double f_dm(double re, double im) {
        mpfr_set_d(r20020, re, MPFR_RNDN);
        mpfr_cos(r20021, r20020, MPFR_RNDN);
        ;
        mpfr_mul(r20023, r20021, r20022, MPFR_RNDN);
        ;
        mpfr_neg(r20025, r20024, MPFR_RNDN);
        mpfr_set_d(r20026, im, MPFR_RNDN);
        mpfr_mul(r20027, r20026, r20026, MPFR_RNDN); mpfr_mul(r20027, r20027, r20026, MPFR_RNDN);
        mpfr_mul(r20028, r20025, r20027, MPFR_RNDN);
        mpfr_neg(r20029, r20026, MPFR_RNDN);
        mpfr_add(r20030, r20028, r20029, MPFR_RNDN);
        ;
        mpfr_neg(r20032, r20031, MPFR_RNDN);
        ;
        mpfr_pow(r20034, r20026, r20033, MPFR_RNDN);
        mpfr_mul(r20035, r20032, r20034, MPFR_RNDN);
        mpfr_add(r20036, r20030, r20035, MPFR_RNDN);
        mpfr_exp(r20037, r20026, MPFR_RNDN);
        mpfr_sqrt(r20038, r20037, MPFR_RNDN);
        mpfr_exp(r20039, r20029, MPFR_RNDN);
        mpfr_sqrt(r20040, r20039, MPFR_RNDN);
        mpfr_add(r20041, r20038, r20040, MPFR_RNDN);
        mpfr_mul(r20042, r20036, r20041, MPFR_RNDN);
        mpfr_mul(r20043, r20023, r20042, MPFR_RNDN);
        return mpfr_get_d(r20043, MPFR_RNDN);
}