#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 r7975 = 0.5f;
        float r7976 = re;
        float r7977 = sin(r7976);
        float r7978 = r7975 * r7977;
        float r7979 = 0.0f;
        float r7980 = im;
        float r7981 = r7979 - r7980;
        float r7982 = exp(r7981);
        float r7983 = exp(r7980);
        float r7984 = r7982 + r7983;
        float r7985 = r7978 * r7984;
        return r7985;
}

double f_id(double re, double im) {
        double r7986 = 0.5;
        double r7987 = re;
        double r7988 = sin(r7987);
        double r7989 = r7986 * r7988;
        double r7990 = 0.0;
        double r7991 = im;
        double r7992 = r7990 - r7991;
        double r7993 = exp(r7992);
        double r7994 = exp(r7991);
        double r7995 = r7993 + r7994;
        double r7996 = r7989 * r7995;
        return r7996;
}


double f_of(float re, float im) {
        float r7997 = 0.5f;
        float r7998 = re;
        float r7999 = sin(r7998);
        float r8000 = r7997 * r7999;
        float r8001 = im;
        float r8002 = exp(r8001);
        float r8003 = r8000 / r8002;
        float r8004 = fma(r8000, r8002, r8003);
        return r8004;
}

double f_od(double re, double im) {
        double r8005 = 0.5;
        double r8006 = re;
        double r8007 = sin(r8006);
        double r8008 = r8005 * r8007;
        double r8009 = im;
        double r8010 = exp(r8009);
        double r8011 = r8008 / r8010;
        double r8012 = fma(r8008, r8010, r8011);
        return r8012;
}

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 r8013, r8014, r8015, r8016, r8017, r8018, r8019, r8020, r8021, r8022, r8023;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8013, "0.5", 10, MPFR_RNDN);
        mpfr_init(r8014);
        mpfr_init(r8015);
        mpfr_init(r8016);
        mpfr_init_set_str(r8017, "0", 10, MPFR_RNDN);
        mpfr_init(r8018);
        mpfr_init(r8019);
        mpfr_init(r8020);
        mpfr_init(r8021);
        mpfr_init(r8022);
        mpfr_init(r8023);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r8014, re, MPFR_RNDN);
        mpfr_sin(r8015, r8014, MPFR_RNDN);
        mpfr_mul(r8016, r8013, r8015, MPFR_RNDN);
        ;
        mpfr_set_d(r8018, im, MPFR_RNDN);
        mpfr_sub(r8019, r8017, r8018, MPFR_RNDN);
        mpfr_exp(r8020, r8019, MPFR_RNDN);
        mpfr_exp(r8021, r8018, MPFR_RNDN);
        mpfr_add(r8022, r8020, r8021, MPFR_RNDN);
        mpfr_mul(r8023, r8016, r8022, MPFR_RNDN);
        return mpfr_get_d(r8023, MPFR_RNDN);
}

static mpfr_t r8024, r8025, r8026, r8027, r8028, r8029, r8030, r8031;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8024, "0.5", 10, MPFR_RNDN);
        mpfr_init(r8025);
        mpfr_init(r8026);
        mpfr_init(r8027);
        mpfr_init(r8028);
        mpfr_init(r8029);
        mpfr_init(r8030);
        mpfr_init(r8031);
}

double f_fm(double re, double im) {
        ;
        mpfr_set_d(r8025, re, MPFR_RNDN);
        mpfr_sin(r8026, r8025, MPFR_RNDN);
        mpfr_mul(r8027, r8024, r8026, MPFR_RNDN);
        mpfr_set_d(r8028, im, MPFR_RNDN);
        mpfr_exp(r8029, r8028, MPFR_RNDN);
        mpfr_div(r8030, r8027, r8029, MPFR_RNDN);
        mpfr_fma(r8031, r8027, r8029, r8030, MPFR_RNDN);
        return mpfr_get_d(r8031, MPFR_RNDN);
}

static mpfr_t r8032, r8033, r8034, r8035, r8036, r8037, r8038, r8039;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8032, "0.5", 10, MPFR_RNDN);
        mpfr_init(r8033);
        mpfr_init(r8034);
        mpfr_init(r8035);
        mpfr_init(r8036);
        mpfr_init(r8037);
        mpfr_init(r8038);
        mpfr_init(r8039);
}

double f_dm(double re, double im) {
        ;
        mpfr_set_d(r8033, re, MPFR_RNDN);
        mpfr_sin(r8034, r8033, MPFR_RNDN);
        mpfr_mul(r8035, r8032, r8034, MPFR_RNDN);
        mpfr_set_d(r8036, im, MPFR_RNDN);
        mpfr_exp(r8037, r8036, MPFR_RNDN);
        mpfr_div(r8038, r8035, r8037, MPFR_RNDN);
        mpfr_fma(r8039, r8035, r8037, r8038, MPFR_RNDN);
        return mpfr_get_d(r8039, MPFR_RNDN);
}