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

char *name = "FastMath dist";

double f_if(float d1, float d2, float d3) {
        float r26640 = d1;
        float r26641 = d2;
        float r26642 = r26640 * r26641;
        float r26643 = d3;
        float r26644 = r26640 * r26643;
        float r26645 = r26642 + r26644;
        return r26645;
}

double f_id(double d1, double d2, double d3) {
        double r26646 = d1;
        double r26647 = d2;
        double r26648 = r26646 * r26647;
        double r26649 = d3;
        double r26650 = r26646 * r26649;
        double r26651 = r26648 + r26650;
        return r26651;
}


double f_of(float d1, float d2, float d3) {
        float r26652 = d1;
        float r26653 = d3;
        float r26654 = d2;
        float r26655 = r26653 + r26654;
        float r26656 = r26652 * r26655;
        return r26656;
}

double f_od(double d1, double d2, double d3) {
        double r26657 = d1;
        double r26658 = d3;
        double r26659 = d2;
        double r26660 = r26658 + r26659;
        double r26661 = r26657 * r26660;
        return r26661;
}

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 r26662, r26663, r26664, r26665, r26666, r26667;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r26662);
        mpfr_init(r26663);
        mpfr_init(r26664);
        mpfr_init(r26665);
        mpfr_init(r26666);
        mpfr_init(r26667);
}

double f_im(double d1, double d2, double d3) {
        mpfr_set_d(r26662, d1, MPFR_RNDN);
        mpfr_set_d(r26663, d2, MPFR_RNDN);
        mpfr_mul(r26664, r26662, r26663, MPFR_RNDN);
        mpfr_set_d(r26665, d3, MPFR_RNDN);
        mpfr_mul(r26666, r26662, r26665, MPFR_RNDN);
        mpfr_add(r26667, r26664, r26666, MPFR_RNDN);
        return mpfr_get_d(r26667, MPFR_RNDN);
}

static mpfr_t r26668, r26669, r26670, r26671, r26672;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r26668);
        mpfr_init(r26669);
        mpfr_init(r26670);
        mpfr_init(r26671);
        mpfr_init(r26672);
}

double f_fm(double d1, double d2, double d3) {
        mpfr_set_d(r26668, d1, MPFR_RNDN);
        mpfr_set_d(r26669, d3, MPFR_RNDN);
        mpfr_set_d(r26670, d2, MPFR_RNDN);
        mpfr_add(r26671, r26669, r26670, MPFR_RNDN);
        mpfr_mul(r26672, r26668, r26671, MPFR_RNDN);
        return mpfr_get_d(r26672, MPFR_RNDN);
}

static mpfr_t r26673, r26674, r26675, r26676, r26677;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r26673);
        mpfr_init(r26674);
        mpfr_init(r26675);
        mpfr_init(r26676);
        mpfr_init(r26677);
}

double f_dm(double d1, double d2, double d3) {
        mpfr_set_d(r26673, d1, MPFR_RNDN);
        mpfr_set_d(r26674, d3, MPFR_RNDN);
        mpfr_set_d(r26675, d2, MPFR_RNDN);
        mpfr_add(r26676, r26674, r26675, MPFR_RNDN);
        mpfr_mul(r26677, r26673, r26676, MPFR_RNDN);
        return mpfr_get_d(r26677, MPFR_RNDN);
}