#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "math.sin on complex, real part";

double f_if(float re, float im) {
        float r20859 = 0.5;
        float r20860 = re;
        float r20861 = sin(r20860);
        float r20862 = r20859 * r20861;
        float r20863 = 0;
        float r20864 = im;
        float r20865 = r20863 - r20864;
        float r20866 = exp(r20865);
        float r20867 = exp(r20864);
        float r20868 = r20866 + r20867;
        float r20869 = r20862 * r20868;
        return r20869;
}

double f_id(double re, double im) {
        double r20870 = 0.5;
        double r20871 = re;
        double r20872 = sin(r20871);
        double r20873 = r20870 * r20872;
        double r20874 = 0;
        double r20875 = im;
        double r20876 = r20874 - r20875;
        double r20877 = exp(r20876);
        double r20878 = exp(r20875);
        double r20879 = r20877 + r20878;
        double r20880 = r20873 * r20879;
        return r20880;
}


double f_of(float re, float im) {
        float r20881 = re;
        float r20882 = sin(r20881);
        float r20883 = im;
        float r20884 = exp(r20883);
        float r20885 = r20882 / r20884;
        float r20886 = 0.5;
        float r20887 = r20885 * r20886;
        float r20888 = r20886 * r20882;
        float r20889 = r20888 * r20884;
        float r20890 = r20887 + r20889;
        return r20890;
}

double f_od(double re, double im) {
        double r20891 = re;
        double r20892 = sin(r20891);
        double r20893 = im;
        double r20894 = exp(r20893);
        double r20895 = r20892 / r20894;
        double r20896 = 0.5;
        double r20897 = r20895 * r20896;
        double r20898 = r20896 * r20892;
        double r20899 = r20898 * r20894;
        double r20900 = r20897 + r20899;
        return r20900;
}

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 r20901, r20902, r20903, r20904, r20905, r20906, r20907, r20908, r20909, r20910, r20911;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r20901, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20902);
        mpfr_init(r20903);
        mpfr_init(r20904);
        mpfr_init_set_str(r20905, "0", 10, MPFR_RNDN);
        mpfr_init(r20906);
        mpfr_init(r20907);
        mpfr_init(r20908);
        mpfr_init(r20909);
        mpfr_init(r20910);
        mpfr_init(r20911);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r20902, re, MPFR_RNDN);
        mpfr_sin(r20903, r20902, MPFR_RNDN);
        mpfr_mul(r20904, r20901, r20903, MPFR_RNDN);
        ;
        mpfr_set_d(r20906, im, MPFR_RNDN);
        mpfr_sub(r20907, r20905, r20906, MPFR_RNDN);
        mpfr_exp(r20908, r20907, MPFR_RNDN);
        mpfr_exp(r20909, r20906, MPFR_RNDN);
        mpfr_add(r20910, r20908, r20909, MPFR_RNDN);
        mpfr_mul(r20911, r20904, r20910, MPFR_RNDN);
        return mpfr_get_d(r20911, MPFR_RNDN);
}

static mpfr_t r20912, r20913, r20914, r20915, r20916, r20917, r20918, r20919, r20920, r20921;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20912);
        mpfr_init(r20913);
        mpfr_init(r20914);
        mpfr_init(r20915);
        mpfr_init(r20916);
        mpfr_init_set_str(r20917, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20918);
        mpfr_init(r20919);
        mpfr_init(r20920);
        mpfr_init(r20921);
}

double f_fm(double re, double im) {
        mpfr_set_d(r20912, re, MPFR_RNDN);
        mpfr_sin(r20913, r20912, MPFR_RNDN);
        mpfr_set_d(r20914, im, MPFR_RNDN);
        mpfr_exp(r20915, r20914, MPFR_RNDN);
        mpfr_div(r20916, r20913, r20915, MPFR_RNDN);
        ;
        mpfr_mul(r20918, r20916, r20917, MPFR_RNDN);
        mpfr_mul(r20919, r20917, r20913, MPFR_RNDN);
        mpfr_mul(r20920, r20919, r20915, MPFR_RNDN);
        mpfr_add(r20921, r20918, r20920, MPFR_RNDN);
        return mpfr_get_d(r20921, MPFR_RNDN);
}

static mpfr_t r20922, r20923, r20924, r20925, r20926, r20927, r20928, r20929, r20930, r20931;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r20922);
        mpfr_init(r20923);
        mpfr_init(r20924);
        mpfr_init(r20925);
        mpfr_init(r20926);
        mpfr_init_set_str(r20927, "0.5", 10, MPFR_RNDN);
        mpfr_init(r20928);
        mpfr_init(r20929);
        mpfr_init(r20930);
        mpfr_init(r20931);
}

double f_dm(double re, double im) {
        mpfr_set_d(r20922, re, MPFR_RNDN);
        mpfr_sin(r20923, r20922, MPFR_RNDN);
        mpfr_set_d(r20924, im, MPFR_RNDN);
        mpfr_exp(r20925, r20924, MPFR_RNDN);
        mpfr_div(r20926, r20923, r20925, MPFR_RNDN);
        ;
        mpfr_mul(r20928, r20926, r20927, MPFR_RNDN);
        mpfr_mul(r20929, r20927, r20923, MPFR_RNDN);
        mpfr_mul(r20930, r20929, r20925, MPFR_RNDN);
        mpfr_add(r20931, r20928, r20930, MPFR_RNDN);
        return mpfr_get_d(r20931, MPFR_RNDN);
}