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

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

double f_if(float re, float im) {
        float r8963 = 0.5;
        float r8964 = re;
        float r8965 = cos(r8964);
        float r8966 = r8963 * r8965;
        float r8967 = im;
        float r8968 = -r8967;
        float r8969 = exp(r8968);
        float r8970 = exp(r8967);
        float r8971 = r8969 + r8970;
        float r8972 = r8966 * r8971;
        return r8972;
}

double f_id(double re, double im) {
        double r8973 = 0.5;
        double r8974 = re;
        double r8975 = cos(r8974);
        double r8976 = r8973 * r8975;
        double r8977 = im;
        double r8978 = -r8977;
        double r8979 = exp(r8978);
        double r8980 = exp(r8977);
        double r8981 = r8979 + r8980;
        double r8982 = r8976 * r8981;
        return r8982;
}


double f_of(float re, float im) {
        float r8983 = 0.5;
        float r8984 = re;
        float r8985 = cos(r8984);
        float r8986 = r8983 * r8985;
        float r8987 = im;
        float r8988 = -r8987;
        float r8989 = exp(r8988);
        float r8990 = exp(r8987);
        float r8991 = r8989 + r8990;
        float r8992 = r8986 * r8991;
        return r8992;
}

double f_od(double re, double im) {
        double r8993 = 0.5;
        double r8994 = re;
        double r8995 = cos(r8994);
        double r8996 = r8993 * r8995;
        double r8997 = im;
        double r8998 = -r8997;
        double r8999 = exp(r8998);
        double r9000 = exp(r8997);
        double r9001 = r8999 + r9000;
        double r9002 = r8996 * r9001;
        return r9002;
}

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 r9003, r9004, r9005, r9006, r9007, r9008, r9009, r9010, r9011, r9012;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9003, "0.5", 10, MPFR_RNDN);
        mpfr_init(r9004);
        mpfr_init(r9005);
        mpfr_init(r9006);
        mpfr_init(r9007);
        mpfr_init(r9008);
        mpfr_init(r9009);
        mpfr_init(r9010);
        mpfr_init(r9011);
        mpfr_init(r9012);
}

double f_im(double re, double im) {
        ;
        mpfr_set_d(r9004, re, MPFR_RNDN);
        mpfr_cos(r9005, r9004, MPFR_RNDN);
        mpfr_mul(r9006, r9003, r9005, MPFR_RNDN);
        mpfr_set_d(r9007, im, MPFR_RNDN);
        mpfr_neg(r9008, r9007, MPFR_RNDN);
        mpfr_exp(r9009, r9008, MPFR_RNDN);
        mpfr_exp(r9010, r9007, MPFR_RNDN);
        mpfr_add(r9011, r9009, r9010, MPFR_RNDN);
        mpfr_mul(r9012, r9006, r9011, MPFR_RNDN);
        return mpfr_get_d(r9012, MPFR_RNDN);
}

static mpfr_t r9013, r9014, r9015, r9016, r9017, r9018, r9019, r9020, r9021, r9022;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9013, "0.5", 10, MPFR_RNDN);
        mpfr_init(r9014);
        mpfr_init(r9015);
        mpfr_init(r9016);
        mpfr_init(r9017);
        mpfr_init(r9018);
        mpfr_init(r9019);
        mpfr_init(r9020);
        mpfr_init(r9021);
        mpfr_init(r9022);
}

double f_fm(double re, double im) {
        ;
        mpfr_set_d(r9014, re, MPFR_RNDN);
        mpfr_cos(r9015, r9014, MPFR_RNDN);
        mpfr_mul(r9016, r9013, r9015, MPFR_RNDN);
        mpfr_set_d(r9017, im, MPFR_RNDN);
        mpfr_neg(r9018, r9017, MPFR_RNDN);
        mpfr_exp(r9019, r9018, MPFR_RNDN);
        mpfr_exp(r9020, r9017, MPFR_RNDN);
        mpfr_add(r9021, r9019, r9020, MPFR_RNDN);
        mpfr_mul(r9022, r9016, r9021, MPFR_RNDN);
        return mpfr_get_d(r9022, MPFR_RNDN);
}

static mpfr_t r9023, r9024, r9025, r9026, r9027, r9028, r9029, r9030, r9031, r9032;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9023, "0.5", 10, MPFR_RNDN);
        mpfr_init(r9024);
        mpfr_init(r9025);
        mpfr_init(r9026);
        mpfr_init(r9027);
        mpfr_init(r9028);
        mpfr_init(r9029);
        mpfr_init(r9030);
        mpfr_init(r9031);
        mpfr_init(r9032);
}

double f_dm(double re, double im) {
        ;
        mpfr_set_d(r9024, re, MPFR_RNDN);
        mpfr_cos(r9025, r9024, MPFR_RNDN);
        mpfr_mul(r9026, r9023, r9025, MPFR_RNDN);
        mpfr_set_d(r9027, im, MPFR_RNDN);
        mpfr_neg(r9028, r9027, MPFR_RNDN);
        mpfr_exp(r9029, r9028, MPFR_RNDN);
        mpfr_exp(r9030, r9027, MPFR_RNDN);
        mpfr_add(r9031, r9029, r9030, MPFR_RNDN);
        mpfr_mul(r9032, r9026, r9031, MPFR_RNDN);
        return mpfr_get_d(r9032, MPFR_RNDN);
}

