#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 r19895 = 0.5f;
        float r19896 = re;
        float r19897 = cos(r19896);
        float r19898 = r19895 * r19897;
        float r19899 = 0.0f;
        float r19900 = im;
        float r19901 = r19899 - r19900;
        float r19902 = exp(r19901);
        float r19903 = exp(r19900);
        float r19904 = r19902 - r19903;
        float r19905 = r19898 * r19904;
        return r19905;
}

double f_id(double re, double im) {
        double r19906 = 0.5;
        double r19907 = re;
        double r19908 = cos(r19907);
        double r19909 = r19906 * r19908;
        double r19910 = 0.0;
        double r19911 = im;
        double r19912 = r19910 - r19911;
        double r19913 = exp(r19912);
        double r19914 = exp(r19911);
        double r19915 = r19913 - r19914;
        double r19916 = r19909 * r19915;
        return r19916;
}


double f_of(float re, float im) {
        float r19917 = re;
        float r19918 = cos(r19917);
        float r19919 = 0.5f;
        float r19920 = r19918 * r19919;
        float r19921 = 0.0416666679084301f;
        float r19922 = -r19921;
        float r19923 = im;
        float r19924 = r19923 * (r19923 * r19923);
        float r19925 = r19922 * r19924;
        float r19926 = -r19923;
        float r19927 = r19925 + r19926;
        float r19928 = 0.0005208333604969084f;
        float r19929 = -r19928;
        float r19930 = 5.0f;
        float r19931 = pow(r19923, r19930);
        float r19932 = r19929 * r19931;
        float r19933 = r19927 + r19932;
        float r19934 = exp(r19923);
        float r19935 = sqrt(r19934);
        float r19936 = exp(r19926);
        float r19937 = sqrt(r19936);
        float r19938 = r19935 + r19937;
        float r19939 = r19933 * r19938;
        float r19940 = r19920 * r19939;
        return r19940;
}

double f_od(double re, double im) {
        double r19941 = re;
        double r19942 = cos(r19941);
        double r19943 = 0.5;
        double r19944 = r19942 * r19943;
        double r19945 = 0.0416666679084301;
        double r19946 = -r19945;
        double r19947 = im;
        double r19948 = r19947 * (r19947 * r19947);
        double r19949 = r19946 * r19948;
        double r19950 = -r19947;
        double r19951 = r19949 + r19950;
        double r19952 = 0.0005208333604969084;
        double r19953 = -r19952;
        double r19954 = 5.0;
        double r19955 = pow(r19947, r19954);
        double r19956 = r19953 * r19955;
        double r19957 = r19951 + r19956;
        double r19958 = exp(r19947);
        double r19959 = sqrt(r19958);
        double r19960 = exp(r19950);
        double r19961 = sqrt(r19960);
        double r19962 = r19959 + r19961;
        double r19963 = r19957 * r19962;
        double r19964 = r19944 * r19963;
        return r19964;
}

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 r19965, r19966, r19967, r19968, r19969, r19970, r19971, r19972, r19973, r19974, r19975;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19965, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19966);
        mpfr_init(r19967);
        mpfr_init(r19968);
        mpfr_init_set_str(r19969, "0", 10, MPFR_RNDN);
        mpfr_init(r19970);
        mpfr_init(r19971);
        mpfr_init(r19972);
        mpfr_init(r19973);
        mpfr_init(r19974);
        mpfr_init(r19975);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r19966, re, MPFR_RNDN);
        mpfr_cos(r19967, r19966, MPFR_RNDN);
        mpfr_mul(r19968, r19965, r19967, MPFR_RNDN);
        ;
        mpfr_set_d(r19970, im, MPFR_RNDN);
        mpfr_sub(r19971, r19969, r19970, MPFR_RNDN);
        mpfr_exp(r19972, r19971, MPFR_RNDN);
        mpfr_exp(r19973, r19970, MPFR_RNDN);
        mpfr_sub(r19974, r19972, r19973, MPFR_RNDN);
        mpfr_mul(r19975, r19968, r19974, MPFR_RNDN);
        return mpfr_get_d(r19975, MPFR_RNDN);
}

static mpfr_t r19976, r19977, r19978, r19979, r19980, r19981, r19982, r19983, r19984, r19985, r19986, r19987, r19988, r19989, r19990, r19991, r19992, r19993, r19994, r19995, r19996, r19997, r19998, r19999;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19976);
        mpfr_init(r19977);
        mpfr_init_set_str(r19978, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19979);
        mpfr_init_set_str(r19980, "1/24", 10, MPFR_RNDN);
        mpfr_init(r19981);
        mpfr_init(r19982);
        mpfr_init(r19983);
        mpfr_init(r19984);
        mpfr_init(r19985);
        mpfr_init(r19986);
        mpfr_init_set_str(r19987, "1/1920", 10, MPFR_RNDN);
        mpfr_init(r19988);
        mpfr_init_set_str(r19989, "5", 10, MPFR_RNDN);
        mpfr_init(r19990);
        mpfr_init(r19991);
        mpfr_init(r19992);
        mpfr_init(r19993);
        mpfr_init(r19994);
        mpfr_init(r19995);
        mpfr_init(r19996);
        mpfr_init(r19997);
        mpfr_init(r19998);
        mpfr_init(r19999);
}

double f_fm(double re, double im) {
        mpfr_set_d(r19976, re, MPFR_RNDN);
        mpfr_cos(r19977, r19976, MPFR_RNDN);
        ;
        mpfr_mul(r19979, r19977, r19978, MPFR_RNDN);
        ;
        mpfr_neg(r19981, r19980, MPFR_RNDN);
        mpfr_set_d(r19982, im, MPFR_RNDN);
        mpfr_mul(r19983, r19982, r19982, MPFR_RNDN); mpfr_mul(r19983, r19983, r19982, MPFR_RNDN);
        mpfr_mul(r19984, r19981, r19983, MPFR_RNDN);
        mpfr_neg(r19985, r19982, MPFR_RNDN);
        mpfr_add(r19986, r19984, r19985, MPFR_RNDN);
        ;
        mpfr_neg(r19988, r19987, MPFR_RNDN);
        ;
        mpfr_pow(r19990, r19982, r19989, MPFR_RNDN);
        mpfr_mul(r19991, r19988, r19990, MPFR_RNDN);
        mpfr_add(r19992, r19986, r19991, MPFR_RNDN);
        mpfr_exp(r19993, r19982, MPFR_RNDN);
        mpfr_sqrt(r19994, r19993, MPFR_RNDN);
        mpfr_exp(r19995, r19985, MPFR_RNDN);
        mpfr_sqrt(r19996, r19995, MPFR_RNDN);
        mpfr_add(r19997, r19994, r19996, MPFR_RNDN);
        mpfr_mul(r19998, r19992, r19997, MPFR_RNDN);
        mpfr_mul(r19999, r19979, r19998, MPFR_RNDN);
        return mpfr_get_d(r19999, MPFR_RNDN);
}

static mpfr_t r20000, r20001, r20002, r20003, r20004, r20005, r20006, r20007, r20008, r20009, r20010, r20011, r20012, r20013, r20014, r20015, r20016, r20017, r20018, r20019, r20020, r20021, r20022, r20023;

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

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