#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 r18765 = 0.5f;
        float r18766 = re;
        float r18767 = sin(r18766);
        float r18768 = r18765 * r18767;
        float r18769 = im;
        float r18770 = -r18769;
        float r18771 = exp(r18770);
        float r18772 = exp(r18769);
        float r18773 = r18771 - r18772;
        float r18774 = r18768 * r18773;
        return r18774;
}

double f_id(double re, double im) {
        double r18775 = 0.5;
        double r18776 = re;
        double r18777 = sin(r18776);
        double r18778 = r18775 * r18777;
        double r18779 = im;
        double r18780 = -r18779;
        double r18781 = exp(r18780);
        double r18782 = exp(r18779);
        double r18783 = r18781 - r18782;
        double r18784 = r18778 * r18783;
        return r18784;
}


double f_of(float re, float im) {
        float r18785 = im;
        float r18786 = r18785 * (r18785 * r18785);
        float r18787 = 0.3333333432674408f;
        float r18788 = 5.0f;
        float r18789 = pow(r18785, r18788);
        float r18790 = 0.01666666753590107f;
        float r18791 = 2.0f;
        float r18792 = r18785 * r18791;
        float r18793 = fma(r18789, r18790, r18792);
        float r18794 = fma(r18786, r18787, r18793);
        float r18795 = re;
        float r18796 = sin(r18795);
        float r18797 = 0.5f;
        float r18798 = -r18797;
        float r18799 = r18796 * r18798;
        float r18800 = r18794 * r18799;
        return r18800;
}

double f_od(double re, double im) {
        double r18801 = im;
        double r18802 = r18801 * (r18801 * r18801);
        double r18803 = 0.3333333432674408;
        double r18804 = 5.0;
        double r18805 = pow(r18801, r18804);
        double r18806 = 0.01666666753590107;
        double r18807 = 2.0;
        double r18808 = r18801 * r18807;
        double r18809 = fma(r18805, r18806, r18808);
        double r18810 = fma(r18802, r18803, r18809);
        double r18811 = re;
        double r18812 = sin(r18811);
        double r18813 = 0.5;
        double r18814 = -r18813;
        double r18815 = r18812 * r18814;
        double r18816 = r18810 * r18815;
        return r18816;
}

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 r18817, r18818, r18819, r18820, r18821, r18822, r18823, r18824, r18825, r18826;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18817, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18818);
        mpfr_init(r18819);
        mpfr_init(r18820);
        mpfr_init(r18821);
        mpfr_init(r18822);
        mpfr_init(r18823);
        mpfr_init(r18824);
        mpfr_init(r18825);
        mpfr_init(r18826);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r18818, re, MPFR_RNDN);
        mpfr_sin(r18819, r18818, MPFR_RNDN);
        mpfr_mul(r18820, r18817, r18819, MPFR_RNDN);
        mpfr_set_d(r18821, im, MPFR_RNDN);
        mpfr_neg(r18822, r18821, MPFR_RNDN);
        mpfr_exp(r18823, r18822, MPFR_RNDN);
        mpfr_exp(r18824, r18821, MPFR_RNDN);
        mpfr_sub(r18825, r18823, r18824, MPFR_RNDN);
        mpfr_mul(r18826, r18820, r18825, MPFR_RNDN);
        return mpfr_get_d(r18826, MPFR_RNDN);
}

static mpfr_t r18827, r18828, r18829, r18830, r18831, r18832, r18833, r18834, r18835, r18836, r18837, r18838, r18839, r18840, r18841, r18842;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18827);
        mpfr_init(r18828);
        mpfr_init_set_str(r18829, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r18830, "5", 10, MPFR_RNDN);
        mpfr_init(r18831);
        mpfr_init_set_str(r18832, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r18833, "2", 10, MPFR_RNDN);
        mpfr_init(r18834);
        mpfr_init(r18835);
        mpfr_init(r18836);
        mpfr_init(r18837);
        mpfr_init(r18838);
        mpfr_init_set_str(r18839, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18840);
        mpfr_init(r18841);
        mpfr_init(r18842);
}

double f_fm(double re, double im) {
        mpfr_set_d(r18827, im, MPFR_RNDN);
        mpfr_mul(r18828, r18827, r18827, MPFR_RNDN); mpfr_mul(r18828, r18828, r18827, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r18831, r18827, r18830, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18834, r18827, r18833, MPFR_RNDN);
        mpfr_fma(r18835, r18831, r18832, r18834, MPFR_RNDN);
        mpfr_fma(r18836, r18828, r18829, r18835, MPFR_RNDN);
        mpfr_set_d(r18837, re, MPFR_RNDN);
        mpfr_sin(r18838, r18837, MPFR_RNDN);
        ;
        mpfr_neg(r18840, r18839, MPFR_RNDN);
        mpfr_mul(r18841, r18838, r18840, MPFR_RNDN);
        mpfr_mul(r18842, r18836, r18841, MPFR_RNDN);
        return mpfr_get_d(r18842, MPFR_RNDN);
}

static mpfr_t r18843, r18844, r18845, r18846, r18847, r18848, r18849, r18850, r18851, r18852, r18853, r18854, r18855, r18856, r18857, r18858;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18843);
        mpfr_init(r18844);
        mpfr_init_set_str(r18845, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r18846, "5", 10, MPFR_RNDN);
        mpfr_init(r18847);
        mpfr_init_set_str(r18848, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r18849, "2", 10, MPFR_RNDN);
        mpfr_init(r18850);
        mpfr_init(r18851);
        mpfr_init(r18852);
        mpfr_init(r18853);
        mpfr_init(r18854);
        mpfr_init_set_str(r18855, "0.5", 10, MPFR_RNDN);
        mpfr_init(r18856);
        mpfr_init(r18857);
        mpfr_init(r18858);
}

double f_dm(double re, double im) {
        mpfr_set_d(r18843, im, MPFR_RNDN);
        mpfr_mul(r18844, r18843, r18843, MPFR_RNDN); mpfr_mul(r18844, r18844, r18843, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r18847, r18843, r18846, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r18850, r18843, r18849, MPFR_RNDN);
        mpfr_fma(r18851, r18847, r18848, r18850, MPFR_RNDN);
        mpfr_fma(r18852, r18844, r18845, r18851, MPFR_RNDN);
        mpfr_set_d(r18853, re, MPFR_RNDN);
        mpfr_sin(r18854, r18853, MPFR_RNDN);
        ;
        mpfr_neg(r18856, r18855, MPFR_RNDN);
        mpfr_mul(r18857, r18854, r18856, MPFR_RNDN);
        mpfr_mul(r18858, r18852, r18857, MPFR_RNDN);
        return mpfr_get_d(r18858, MPFR_RNDN);
}