#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 r8469 = re;
        float r8470 = exp(r8469);
        float r8471 = im;
        float r8472 = sin(r8471);
        float r8473 = r8470 * r8472;
        return r8473;
}

double f_id(double re, double im) {
        double r8474 = re;
        double r8475 = exp(r8474);
        double r8476 = im;
        double r8477 = sin(r8476);
        double r8478 = r8475 * r8477;
        return r8478;
}


double f_of(float re, float im) {
        float r8479 = im;
        float r8480 = sin(r8479);
        float r8481 = re;
        float r8482 = exp(r8481);
        float r8483 = r8480 * r8482;
        return r8483;
}

double f_od(double re, double im) {
        double r8484 = im;
        double r8485 = sin(r8484);
        double r8486 = re;
        double r8487 = exp(r8486);
        double r8488 = r8485 * r8487;
        return r8488;
}

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 r8489, r8490, r8491, r8492, r8493;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r8489);
        mpfr_init(r8490);
        mpfr_init(r8491);
        mpfr_init(r8492);
        mpfr_init(r8493);
}

double f_im(double re, double im) {
        mpfr_set_d(r8489, re, MPFR_RNDN);
        mpfr_exp(r8490, r8489, MPFR_RNDN);
        mpfr_set_d(r8491, im, MPFR_RNDN);
        mpfr_sin(r8492, r8491, MPFR_RNDN);
        mpfr_mul(r8493, r8490, r8492, MPFR_RNDN);
        return mpfr_get_d(r8493, MPFR_RNDN);
}

static mpfr_t r8494, r8495, r8496, r8497, r8498;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r8494);
        mpfr_init(r8495);
        mpfr_init(r8496);
        mpfr_init(r8497);
        mpfr_init(r8498);
}

double f_fm(double re, double im) {
        mpfr_set_d(r8494, im, MPFR_RNDN);
        mpfr_sin(r8495, r8494, MPFR_RNDN);
        mpfr_set_d(r8496, re, MPFR_RNDN);
        mpfr_exp(r8497, r8496, MPFR_RNDN);
        mpfr_mul(r8498, r8495, r8497, MPFR_RNDN);
        return mpfr_get_d(r8498, MPFR_RNDN);
}

static mpfr_t r8499, r8500, r8501, r8502, r8503;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r8499);
        mpfr_init(r8500);
        mpfr_init(r8501);
        mpfr_init(r8502);
        mpfr_init(r8503);
}

double f_dm(double re, double im) {
        mpfr_set_d(r8499, im, MPFR_RNDN);
        mpfr_sin(r8500, r8499, MPFR_RNDN);
        mpfr_set_d(r8501, re, MPFR_RNDN);
        mpfr_exp(r8502, r8501, MPFR_RNDN);
        mpfr_mul(r8503, r8500, r8502, MPFR_RNDN);
        return mpfr_get_d(r8503, MPFR_RNDN);
}

