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

char *name = "Data.Number.Erf:$dmerfcx from erf-2.0.0.0";

double f_if(float x, float y) {
        float r55943 = x;
        float r55944 = y;
        float r55945 = r55944 * r55944;
        float r55946 = exp(r55945);
        float r55947 = r55943 * r55946;
        return r55947;
}

double f_id(double x, double y) {
        double r55948 = x;
        double r55949 = y;
        double r55950 = r55949 * r55949;
        double r55951 = exp(r55950);
        double r55952 = r55948 * r55951;
        return r55952;
}


double f_of(float x, float y) {
        float r55953 = x;
        float r55954 = y;
        float r55955 = r55954 * r55954;
        float r55956 = exp(r55955);
        float r55957 = r55953 * r55956;
        return r55957;
}

double f_od(double x, double y) {
        double r55958 = x;
        double r55959 = y;
        double r55960 = r55959 * r55959;
        double r55961 = exp(r55960);
        double r55962 = r55958 * r55961;
        return r55962;
}

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 r55963, r55964, r55965, r55966, r55967;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r55963);
        mpfr_init(r55964);
        mpfr_init(r55965);
        mpfr_init(r55966);
        mpfr_init(r55967);
}

double f_im(double x, double y) {
        mpfr_set_d(r55963, x, MPFR_RNDN);
        mpfr_set_d(r55964, y, MPFR_RNDN);
        mpfr_mul(r55965, r55964, r55964, MPFR_RNDN);
        mpfr_exp(r55966, r55965, MPFR_RNDN);
        mpfr_mul(r55967, r55963, r55966, MPFR_RNDN);
        return mpfr_get_d(r55967, MPFR_RNDN);
}

static mpfr_t r55968, r55969, r55970, r55971, r55972;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r55968);
        mpfr_init(r55969);
        mpfr_init(r55970);
        mpfr_init(r55971);
        mpfr_init(r55972);
}

double f_fm(double x, double y) {
        mpfr_set_d(r55968, x, MPFR_RNDN);
        mpfr_set_d(r55969, y, MPFR_RNDN);
        mpfr_mul(r55970, r55969, r55969, MPFR_RNDN);
        mpfr_exp(r55971, r55970, MPFR_RNDN);
        mpfr_mul(r55972, r55968, r55971, MPFR_RNDN);
        return mpfr_get_d(r55972, MPFR_RNDN);
}

static mpfr_t r55973, r55974, r55975, r55976, r55977;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r55973);
        mpfr_init(r55974);
        mpfr_init(r55975);
        mpfr_init(r55976);
        mpfr_init(r55977);
}

double f_dm(double x, double y) {
        mpfr_set_d(r55973, x, MPFR_RNDN);
        mpfr_set_d(r55974, y, MPFR_RNDN);
        mpfr_mul(r55975, r55974, r55974, MPFR_RNDN);
        mpfr_exp(r55976, r55975, MPFR_RNDN);
        mpfr_mul(r55977, r55973, r55976, MPFR_RNDN);
        return mpfr_get_d(r55977, MPFR_RNDN);
}