#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 r19745 = 0.5f;
        float r19746 = re;
        float r19747 = sin(r19746);
        float r19748 = r19745 * r19747;
        float r19749 = im;
        float r19750 = -r19749;
        float r19751 = exp(r19750);
        float r19752 = exp(r19749);
        float r19753 = r19751 - r19752;
        float r19754 = r19748 * r19753;
        return r19754;
}

double f_id(double re, double im) {
        double r19755 = 0.5;
        double r19756 = re;
        double r19757 = sin(r19756);
        double r19758 = r19755 * r19757;
        double r19759 = im;
        double r19760 = -r19759;
        double r19761 = exp(r19760);
        double r19762 = exp(r19759);
        double r19763 = r19761 - r19762;
        double r19764 = r19758 * r19763;
        return r19764;
}


double f_of(float re, float im) {
        float r19765 = re;
        float r19766 = sin(r19765);
        float r19767 = 0.5f;
        float r19768 = r19766 * r19767;
        float r19769 = -r19768;
        float r19770 = im;
        float r19771 = exp(r19770);
        float r19772 = sqrt(r19771);
        float r19773 = -r19770;
        float r19774 = exp(r19773);
        float r19775 = sqrt(r19774);
        float r19776 = r19772 + r19775;
        float r19777 = r19769 * r19776;
        float r19778 = 0.0005208333333333333f;
        float r19779 = 5.0f;
        float r19780 = pow(r19770, r19779);
        float r19781 = r19778 * r19780;
        float r19782 = 0.041666666666666664f;
        float r19783 = r19770 * (r19770 * r19770);
        float r19784 = r19782 * r19783;
        float r19785 = r19781 + r19784;
        float r19786 = r19785 + r19770;
        float r19787 = r19777 * r19786;
        return r19787;
}

double f_od(double re, double im) {
        double r19788 = re;
        double r19789 = sin(r19788);
        double r19790 = 0.5;
        double r19791 = r19789 * r19790;
        double r19792 = -r19791;
        double r19793 = im;
        double r19794 = exp(r19793);
        double r19795 = sqrt(r19794);
        double r19796 = -r19793;
        double r19797 = exp(r19796);
        double r19798 = sqrt(r19797);
        double r19799 = r19795 + r19798;
        double r19800 = r19792 * r19799;
        double r19801 = 0.0005208333333333333;
        double r19802 = 5.0;
        double r19803 = pow(r19793, r19802);
        double r19804 = r19801 * r19803;
        double r19805 = 0.041666666666666664;
        double r19806 = r19793 * (r19793 * r19793);
        double r19807 = r19805 * r19806;
        double r19808 = r19804 + r19807;
        double r19809 = r19808 + r19793;
        double r19810 = r19800 * r19809;
        return r19810;
}

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 r19811, r19812, r19813, r19814, r19815, r19816, r19817, r19818, r19819, r19820;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19811, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19812);
        mpfr_init(r19813);
        mpfr_init(r19814);
        mpfr_init(r19815);
        mpfr_init(r19816);
        mpfr_init(r19817);
        mpfr_init(r19818);
        mpfr_init(r19819);
        mpfr_init(r19820);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r19812, re, MPFR_RNDN);
        mpfr_sin(r19813, r19812, MPFR_RNDN);
        mpfr_mul(r19814, r19811, r19813, MPFR_RNDN);
        mpfr_set_d(r19815, im, MPFR_RNDN);
        mpfr_neg(r19816, r19815, MPFR_RNDN);
        mpfr_exp(r19817, r19816, MPFR_RNDN);
        mpfr_exp(r19818, r19815, MPFR_RNDN);
        mpfr_sub(r19819, r19817, r19818, MPFR_RNDN);
        mpfr_mul(r19820, r19814, r19819, MPFR_RNDN);
        return mpfr_get_d(r19820, MPFR_RNDN);
}

static mpfr_t r19821, r19822, r19823, r19824, r19825, r19826, r19827, r19828, r19829, r19830, r19831, r19832, r19833, r19834, r19835, r19836, r19837, r19838, r19839, r19840, r19841, r19842, r19843;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19821);
        mpfr_init(r19822);
        mpfr_init_set_str(r19823, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19824);
        mpfr_init(r19825);
        mpfr_init(r19826);
        mpfr_init(r19827);
        mpfr_init(r19828);
        mpfr_init(r19829);
        mpfr_init(r19830);
        mpfr_init(r19831);
        mpfr_init(r19832);
        mpfr_init(r19833);
        mpfr_init_set_str(r19834, "1/1920", 10, MPFR_RNDN);
        mpfr_init_set_str(r19835, "5", 10, MPFR_RNDN);
        mpfr_init(r19836);
        mpfr_init(r19837);
        mpfr_init_set_str(r19838, "1/24", 10, MPFR_RNDN);
        mpfr_init(r19839);
        mpfr_init(r19840);
        mpfr_init(r19841);
        mpfr_init(r19842);
        mpfr_init(r19843);
}

double f_fm(double re, double im) {
        mpfr_set_d(r19821, re, MPFR_RNDN);
        mpfr_sin(r19822, r19821, MPFR_RNDN);
        ;
        mpfr_mul(r19824, r19822, r19823, MPFR_RNDN);
        mpfr_neg(r19825, r19824, MPFR_RNDN);
        mpfr_set_d(r19826, im, MPFR_RNDN);
        mpfr_exp(r19827, r19826, MPFR_RNDN);
        mpfr_sqrt(r19828, r19827, MPFR_RNDN);
        mpfr_neg(r19829, r19826, MPFR_RNDN);
        mpfr_exp(r19830, r19829, MPFR_RNDN);
        mpfr_sqrt(r19831, r19830, MPFR_RNDN);
        mpfr_add(r19832, r19828, r19831, MPFR_RNDN);
        mpfr_mul(r19833, r19825, r19832, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r19836, r19826, r19835, MPFR_RNDN);
        mpfr_mul(r19837, r19834, r19836, MPFR_RNDN);
        ;
        mpfr_mul(r19839, r19826, r19826, MPFR_RNDN); mpfr_mul(r19839, r19839, r19826, MPFR_RNDN);
        mpfr_mul(r19840, r19838, r19839, MPFR_RNDN);
        mpfr_add(r19841, r19837, r19840, MPFR_RNDN);
        mpfr_add(r19842, r19841, r19826, MPFR_RNDN);
        mpfr_mul(r19843, r19833, r19842, MPFR_RNDN);
        return mpfr_get_d(r19843, MPFR_RNDN);
}

static mpfr_t r19844, r19845, r19846, r19847, r19848, r19849, r19850, r19851, r19852, r19853, r19854, r19855, r19856, r19857, r19858, r19859, r19860, r19861, r19862, r19863, r19864, r19865, r19866;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19844);
        mpfr_init(r19845);
        mpfr_init_set_str(r19846, "0.5", 10, MPFR_RNDN);
        mpfr_init(r19847);
        mpfr_init(r19848);
        mpfr_init(r19849);
        mpfr_init(r19850);
        mpfr_init(r19851);
        mpfr_init(r19852);
        mpfr_init(r19853);
        mpfr_init(r19854);
        mpfr_init(r19855);
        mpfr_init(r19856);
        mpfr_init_set_str(r19857, "1/1920", 10, MPFR_RNDN);
        mpfr_init_set_str(r19858, "5", 10, MPFR_RNDN);
        mpfr_init(r19859);
        mpfr_init(r19860);
        mpfr_init_set_str(r19861, "1/24", 10, MPFR_RNDN);
        mpfr_init(r19862);
        mpfr_init(r19863);
        mpfr_init(r19864);
        mpfr_init(r19865);
        mpfr_init(r19866);
}

double f_dm(double re, double im) {
        mpfr_set_d(r19844, re, MPFR_RNDN);
        mpfr_sin(r19845, r19844, MPFR_RNDN);
        ;
        mpfr_mul(r19847, r19845, r19846, MPFR_RNDN);
        mpfr_neg(r19848, r19847, MPFR_RNDN);
        mpfr_set_d(r19849, im, MPFR_RNDN);
        mpfr_exp(r19850, r19849, MPFR_RNDN);
        mpfr_sqrt(r19851, r19850, MPFR_RNDN);
        mpfr_neg(r19852, r19849, MPFR_RNDN);
        mpfr_exp(r19853, r19852, MPFR_RNDN);
        mpfr_sqrt(r19854, r19853, MPFR_RNDN);
        mpfr_add(r19855, r19851, r19854, MPFR_RNDN);
        mpfr_mul(r19856, r19848, r19855, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r19859, r19849, r19858, MPFR_RNDN);
        mpfr_mul(r19860, r19857, r19859, MPFR_RNDN);
        ;
        mpfr_mul(r19862, r19849, r19849, MPFR_RNDN); mpfr_mul(r19862, r19862, r19849, MPFR_RNDN);
        mpfr_mul(r19863, r19861, r19862, MPFR_RNDN);
        mpfr_add(r19864, r19860, r19863, MPFR_RNDN);
        mpfr_add(r19865, r19864, r19849, MPFR_RNDN);
        mpfr_mul(r19866, r19856, r19865, MPFR_RNDN);
        return mpfr_get_d(r19866, MPFR_RNDN);
}