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

char *name = "math.exp on complex, imaginary part";

double f_if(float re, float im) {
        float r9139 = re;
        float r9140 = exp(r9139);
        float r9141 = im;
        float r9142 = sin(r9141);
        float r9143 = r9140 * r9142;
        return r9143;
}

double f_id(double re, double im) {
        double r9144 = re;
        double r9145 = exp(r9144);
        double r9146 = im;
        double r9147 = sin(r9146);
        double r9148 = r9145 * r9147;
        return r9148;
}


double f_of(float re, float im) {
        float r9149 = re;
        float r9150 = exp(r9149);
        float r9151 = sqrt(r9150);
        float r9152 = im;
        float r9153 = sin(r9152);
        float r9154 = r9153 * r9151;
        float r9155 = r9151 * r9154;
        return r9155;
}

double f_od(double re, double im) {
        double r9156 = re;
        double r9157 = exp(r9156);
        double r9158 = sqrt(r9157);
        double r9159 = im;
        double r9160 = sin(r9159);
        double r9161 = r9160 * r9158;
        double r9162 = r9158 * r9161;
        return r9162;
}

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 r9163, r9164, r9165, r9166, r9167;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r9163);
        mpfr_init(r9164);
        mpfr_init(r9165);
        mpfr_init(r9166);
        mpfr_init(r9167);
}

double f_im(double re, double im) {
        mpfr_set_d(r9163, re, MPFR_RNDN);
        mpfr_exp(r9164, r9163, MPFR_RNDN);
        mpfr_set_d(r9165, im, MPFR_RNDN);
        mpfr_sin(r9166, r9165, MPFR_RNDN);
        mpfr_mul(r9167, r9164, r9166, MPFR_RNDN);
        return mpfr_get_d(r9167, MPFR_RNDN);
}

static mpfr_t r9168, r9169, r9170, r9171, r9172, r9173, r9174;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r9168);
        mpfr_init(r9169);
        mpfr_init(r9170);
        mpfr_init(r9171);
        mpfr_init(r9172);
        mpfr_init(r9173);
        mpfr_init(r9174);
}

double f_fm(double re, double im) {
        mpfr_set_d(r9168, re, MPFR_RNDN);
        mpfr_exp(r9169, r9168, MPFR_RNDN);
        mpfr_sqrt(r9170, r9169, MPFR_RNDN);
        mpfr_set_d(r9171, im, MPFR_RNDN);
        mpfr_sin(r9172, r9171, MPFR_RNDN);
        mpfr_mul(r9173, r9172, r9170, MPFR_RNDN);
        mpfr_mul(r9174, r9170, r9173, MPFR_RNDN);
        return mpfr_get_d(r9174, MPFR_RNDN);
}

static mpfr_t r9175, r9176, r9177, r9178, r9179, r9180, r9181;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r9175);
        mpfr_init(r9176);
        mpfr_init(r9177);
        mpfr_init(r9178);
        mpfr_init(r9179);
        mpfr_init(r9180);
        mpfr_init(r9181);
}

double f_dm(double re, double im) {
        mpfr_set_d(r9175, re, MPFR_RNDN);
        mpfr_exp(r9176, r9175, MPFR_RNDN);
        mpfr_sqrt(r9177, r9176, MPFR_RNDN);
        mpfr_set_d(r9178, im, MPFR_RNDN);
        mpfr_sin(r9179, r9178, MPFR_RNDN);
        mpfr_mul(r9180, r9179, r9177, MPFR_RNDN);
        mpfr_mul(r9181, r9177, r9180, MPFR_RNDN);
        return mpfr_get_d(r9181, MPFR_RNDN);
}