#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 r24788 = 0.5;
        float r24789 = re;
        float r24790 = cos(r24789);
        float r24791 = r24788 * r24790;
        float r24792 = 0;
        float r24793 = im;
        float r24794 = r24792 - r24793;
        float r24795 = exp(r24794);
        float r24796 = exp(r24793);
        float r24797 = r24795 - r24796;
        float r24798 = r24791 * r24797;
        return r24798;
}

double f_id(double re, double im) {
        double r24799 = 0.5;
        double r24800 = re;
        double r24801 = cos(r24800);
        double r24802 = r24799 * r24801;
        double r24803 = 0;
        double r24804 = im;
        double r24805 = r24803 - r24804;
        double r24806 = exp(r24805);
        double r24807 = exp(r24804);
        double r24808 = r24806 - r24807;
        double r24809 = r24802 * r24808;
        return r24809;
}


double f_of(float re, float im) {
        float r24810 = im;
        float r24811 = 1/3;
        float r24812 = r24811 * r24810;
        float r24813 = 2;
        float r24814 = fma(r24812, r24810, r24813);
        float r24815 = 5;
        float r24816 = pow(r24810, r24815);
        float r24817 = 1/60;
        float r24818 = r24816 * r24817;
        float r24819 = fma(r24810, r24814, r24818);
        float r24820 = re;
        float r24821 = cos(r24820);
        float r24822 = 0.5;
        float r24823 = -r24822;
        float r24824 = r24821 * r24823;
        float r24825 = r24819 * r24824;
        return r24825;
}

double f_od(double re, double im) {
        double r24826 = im;
        double r24827 = 1/3;
        double r24828 = r24827 * r24826;
        double r24829 = 2;
        double r24830 = fma(r24828, r24826, r24829);
        double r24831 = 5;
        double r24832 = pow(r24826, r24831);
        double r24833 = 1/60;
        double r24834 = r24832 * r24833;
        double r24835 = fma(r24826, r24830, r24834);
        double r24836 = re;
        double r24837 = cos(r24836);
        double r24838 = 0.5;
        double r24839 = -r24838;
        double r24840 = r24837 * r24839;
        double r24841 = r24835 * r24840;
        return r24841;
}

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 r24842, r24843, r24844, r24845, r24846, r24847, r24848, r24849, r24850, r24851, r24852;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r24842, "0.5", 10, MPFR_RNDN);
        mpfr_init(r24843);
        mpfr_init(r24844);
        mpfr_init(r24845);
        mpfr_init_set_str(r24846, "0", 10, MPFR_RNDN);
        mpfr_init(r24847);
        mpfr_init(r24848);
        mpfr_init(r24849);
        mpfr_init(r24850);
        mpfr_init(r24851);
        mpfr_init(r24852);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r24843, re, MPFR_RNDN);
        mpfr_cos(r24844, r24843, MPFR_RNDN);
        mpfr_mul(r24845, r24842, r24844, MPFR_RNDN);
        ;
        mpfr_set_d(r24847, im, MPFR_RNDN);
        mpfr_sub(r24848, r24846, r24847, MPFR_RNDN);
        mpfr_exp(r24849, r24848, MPFR_RNDN);
        mpfr_exp(r24850, r24847, MPFR_RNDN);
        mpfr_sub(r24851, r24849, r24850, MPFR_RNDN);
        mpfr_mul(r24852, r24845, r24851, MPFR_RNDN);
        return mpfr_get_d(r24852, MPFR_RNDN);
}

static mpfr_t r24853, r24854, r24855, r24856, r24857, r24858, r24859, r24860, r24861, r24862, r24863, r24864, r24865, r24866, r24867, r24868;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r24853);
        mpfr_init_set_str(r24854, "1/3", 10, MPFR_RNDN);
        mpfr_init(r24855);
        mpfr_init_set_str(r24856, "2", 10, MPFR_RNDN);
        mpfr_init(r24857);
        mpfr_init_set_str(r24858, "5", 10, MPFR_RNDN);
        mpfr_init(r24859);
        mpfr_init_set_str(r24860, "1/60", 10, MPFR_RNDN);
        mpfr_init(r24861);
        mpfr_init(r24862);
        mpfr_init(r24863);
        mpfr_init(r24864);
        mpfr_init_set_str(r24865, "0.5", 10, MPFR_RNDN);
        mpfr_init(r24866);
        mpfr_init(r24867);
        mpfr_init(r24868);
}

double f_fm(double re, double im) {
        mpfr_set_d(r24853, im, MPFR_RNDN);
        ;
        mpfr_mul(r24855, r24854, r24853, MPFR_RNDN);
        ;
        mpfr_fma(r24857, r24855, r24853, r24856, MPFR_RNDN);
        ;
        mpfr_pow(r24859, r24853, r24858, MPFR_RNDN);
        ;
        mpfr_mul(r24861, r24859, r24860, MPFR_RNDN);
        mpfr_fma(r24862, r24853, r24857, r24861, MPFR_RNDN);
        mpfr_set_d(r24863, re, MPFR_RNDN);
        mpfr_cos(r24864, r24863, MPFR_RNDN);
        ;
        mpfr_neg(r24866, r24865, MPFR_RNDN);
        mpfr_mul(r24867, r24864, r24866, MPFR_RNDN);
        mpfr_mul(r24868, r24862, r24867, MPFR_RNDN);
        return mpfr_get_d(r24868, MPFR_RNDN);
}

static mpfr_t r24869, r24870, r24871, r24872, r24873, r24874, r24875, r24876, r24877, r24878, r24879, r24880, r24881, r24882, r24883, r24884;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r24869);
        mpfr_init_set_str(r24870, "1/3", 10, MPFR_RNDN);
        mpfr_init(r24871);
        mpfr_init_set_str(r24872, "2", 10, MPFR_RNDN);
        mpfr_init(r24873);
        mpfr_init_set_str(r24874, "5", 10, MPFR_RNDN);
        mpfr_init(r24875);
        mpfr_init_set_str(r24876, "1/60", 10, MPFR_RNDN);
        mpfr_init(r24877);
        mpfr_init(r24878);
        mpfr_init(r24879);
        mpfr_init(r24880);
        mpfr_init_set_str(r24881, "0.5", 10, MPFR_RNDN);
        mpfr_init(r24882);
        mpfr_init(r24883);
        mpfr_init(r24884);
}

double f_dm(double re, double im) {
        mpfr_set_d(r24869, im, MPFR_RNDN);
        ;
        mpfr_mul(r24871, r24870, r24869, MPFR_RNDN);
        ;
        mpfr_fma(r24873, r24871, r24869, r24872, MPFR_RNDN);
        ;
        mpfr_pow(r24875, r24869, r24874, MPFR_RNDN);
        ;
        mpfr_mul(r24877, r24875, r24876, MPFR_RNDN);
        mpfr_fma(r24878, r24869, r24873, r24877, MPFR_RNDN);
        mpfr_set_d(r24879, re, MPFR_RNDN);
        mpfr_cos(r24880, r24879, MPFR_RNDN);
        ;
        mpfr_neg(r24882, r24881, MPFR_RNDN);
        mpfr_mul(r24883, r24880, r24882, MPFR_RNDN);
        mpfr_mul(r24884, r24878, r24883, MPFR_RNDN);
        return mpfr_get_d(r24884, MPFR_RNDN);
}