#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 r18925 = 0.5f;
        float r18926 = re;
        float r18927 = sin(r18926);
        float r18928 = r18925 * r18927;
        float r18929 = im;
        float r18930 = -r18929;
        float r18931 = exp(r18930);
        float r18932 = exp(r18929);
        float r18933 = r18931 - r18932;
        float r18934 = r18928 * r18933;
        return r18934;
}

double f_id(double re, double im) {
        double r18935 = 0.5;
        double r18936 = re;
        double r18937 = sin(r18936);
        double r18938 = r18935 * r18937;
        double r18939 = im;
        double r18940 = -r18939;
        double r18941 = exp(r18940);
        double r18942 = exp(r18939);
        double r18943 = r18941 - r18942;
        double r18944 = r18938 * r18943;
        return r18944;
}


double f_of(float re, float im) {
        float r18945 = im;
        float r18946 = r18945 * (r18945 * r18945);
        float r18947 = 0.3333333333333333f;
        float r18948 = 5.0f;
        float r18949 = pow(r18945, r18948);
        float r18950 = 0.016666666666666666f;
        float r18951 = 2.0f;
        float r18952 = r18945 * r18951;
        float r18953 = fma(r18949, r18950, r18952);
        float r18954 = fma(r18946, r18947, r18953);
        float r18955 = re;
        float r18956 = sin(r18955);
        float r18957 = 0.5f;
        float r18958 = -r18957;
        float r18959 = r18956 * r18958;
        float r18960 = r18954 * r18959;
        return r18960;
}

double f_od(double re, double im) {
        double r18961 = im;
        double r18962 = r18961 * (r18961 * r18961);
        double r18963 = 0.3333333333333333;
        double r18964 = 5.0;
        double r18965 = pow(r18961, r18964);
        double r18966 = 0.016666666666666666;
        double r18967 = 2.0;
        double r18968 = r18961 * r18967;
        double r18969 = fma(r18965, r18966, r18968);
        double r18970 = fma(r18962, r18963, r18969);
        double r18971 = re;
        double r18972 = sin(r18971);
        double r18973 = 0.5;
        double r18974 = -r18973;
        double r18975 = r18972 * r18974;
        double r18976 = r18970 * r18975;
        return r18976;
}

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 r18977, r18978, r18979, r18980, r18981, r18982, r18983, r18984, r18985, r18986;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18977, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18978);
        mpfr_init(r18979);
        mpfr_init(r18980);
        mpfr_init(r18981);
        mpfr_init(r18982);
        mpfr_init(r18983);
        mpfr_init(r18984);
        mpfr_init(r18985);
        mpfr_init(r18986);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r18978, re, MPFR_RNDN);
        mpfr_sin(r18979, r18978, MPFR_RNDN);
        mpfr_mul(r18980, r18977, r18979, MPFR_RNDN);
        mpfr_set_d(r18981, im, MPFR_RNDN);
        mpfr_neg(r18982, r18981, MPFR_RNDN);
        mpfr_exp(r18983, r18982, MPFR_RNDN);
        mpfr_exp(r18984, r18981, MPFR_RNDN);
        mpfr_sub(r18985, r18983, r18984, MPFR_RNDN);
        mpfr_mul(r18986, r18980, r18985, MPFR_RNDN);
        return mpfr_get_d(r18986, MPFR_RNDN);
}

static mpfr_t r18987, r18988, r18989, r18990, r18991, r18992, r18993, r18994, r18995, r18996, r18997, r18998, r18999, r19000, r19001, r19002;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18987);
        mpfr_init(r18988);
        mpfr_init_set_str(r18989, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r18990, "5", 10, MPFR_RNDN);
        mpfr_init(r18991);
        mpfr_init_set_str(r18992, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r18993, "2", 10, MPFR_RNDN);
        mpfr_init(r18994);
        mpfr_init(r18995);
        mpfr_init(r18996);
        mpfr_init(r18997);
        mpfr_init(r18998);
        mpfr_init_set_str(r18999, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19000);
        mpfr_init(r19001);
        mpfr_init(r19002);
}

double f_fm(double re, double im) {
        mpfr_set_d(r18987, im, MPFR_RNDN);
        mpfr_mul(r18988, r18987, r18987, MPFR_RNDN); mpfr_mul(r18988, r18988, r18987, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r18991, r18987, r18990, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18994, r18987, r18993, MPFR_RNDN);
        mpfr_fma(r18995, r18991, r18992, r18994, MPFR_RNDN);
        mpfr_fma(r18996, r18988, r18989, r18995, MPFR_RNDN);
        mpfr_set_d(r18997, re, MPFR_RNDN);
        mpfr_sin(r18998, r18997, MPFR_RNDN);
        ;
        mpfr_neg(r19000, r18999, MPFR_RNDN);
        mpfr_mul(r19001, r18998, r19000, MPFR_RNDN);
        mpfr_mul(r19002, r18996, r19001, MPFR_RNDN);
        return mpfr_get_d(r19002, MPFR_RNDN);
}

static mpfr_t r19003, r19004, r19005, r19006, r19007, r19008, r19009, r19010, r19011, r19012, r19013, r19014, r19015, r19016, r19017, r19018;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19003);
        mpfr_init(r19004);
        mpfr_init_set_str(r19005, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r19006, "5", 10, MPFR_RNDN);
        mpfr_init(r19007);
        mpfr_init_set_str(r19008, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r19009, "2", 10, MPFR_RNDN);
        mpfr_init(r19010);
        mpfr_init(r19011);
        mpfr_init(r19012);
        mpfr_init(r19013);
        mpfr_init(r19014);
        mpfr_init_set_str(r19015, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19016);
        mpfr_init(r19017);
        mpfr_init(r19018);
}

double f_dm(double re, double im) {
        mpfr_set_d(r19003, im, MPFR_RNDN);
        mpfr_mul(r19004, r19003, r19003, MPFR_RNDN); mpfr_mul(r19004, r19004, r19003, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r19007, r19003, r19006, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r19010, r19003, r19009, MPFR_RNDN);
        mpfr_fma(r19011, r19007, r19008, r19010, MPFR_RNDN);
        mpfr_fma(r19012, r19004, r19005, r19011, MPFR_RNDN);
        mpfr_set_d(r19013, re, MPFR_RNDN);
        mpfr_sin(r19014, r19013, MPFR_RNDN);
        ;
        mpfr_neg(r19016, r19015, MPFR_RNDN);
        mpfr_mul(r19017, r19014, r19016, MPFR_RNDN);
        mpfr_mul(r19018, r19012, r19017, MPFR_RNDN);
        return mpfr_get_d(r19018, MPFR_RNDN);
}