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

char *name = "FastMath dist4";

double f_if(float d1, float d2, float d3, float d4) {
        float r26743 = d1;
        float r26744 = d2;
        float r26745 = r26743 * r26744;
        float r26746 = d3;
        float r26747 = r26743 * r26746;
        float r26748 = r26745 - r26747;
        float r26749 = d4;
        float r26750 = r26749 * r26743;
        float r26751 = r26748 + r26750;
        float r26752 = r26743 * r26743;
        float r26753 = r26751 - r26752;
        return r26753;
}

double f_id(double d1, double d2, double d3, double d4) {
        double r26754 = d1;
        double r26755 = d2;
        double r26756 = r26754 * r26755;
        double r26757 = d3;
        double r26758 = r26754 * r26757;
        double r26759 = r26756 - r26758;
        double r26760 = d4;
        double r26761 = r26760 * r26754;
        double r26762 = r26759 + r26761;
        double r26763 = r26754 * r26754;
        double r26764 = r26762 - r26763;
        return r26764;
}


double f_of(float d1, float d2, float d3, float d4) {
        float r26765 = d4;
        float r26766 = d3;
        float r26767 = r26765 - r26766;
        float r26768 = d2;
        float r26769 = d1;
        float r26770 = r26768 - r26769;
        float r26771 = r26767 + r26770;
        float r26772 = r26771 * r26769;
        return r26772;
}

double f_od(double d1, double d2, double d3, double d4) {
        double r26773 = d4;
        double r26774 = d3;
        double r26775 = r26773 - r26774;
        double r26776 = d2;
        double r26777 = d1;
        double r26778 = r26776 - r26777;
        double r26779 = r26775 + r26778;
        double r26780 = r26779 * r26777;
        return r26780;
}

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 r26781, r26782, r26783, r26784, r26785, r26786, r26787, r26788, r26789, r26790, r26791;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r26781);
        mpfr_init(r26782);
        mpfr_init(r26783);
        mpfr_init(r26784);
        mpfr_init(r26785);
        mpfr_init(r26786);
        mpfr_init(r26787);
        mpfr_init(r26788);
        mpfr_init(r26789);
        mpfr_init(r26790);
        mpfr_init(r26791);
}

double f_im(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r26781, d1, MPFR_RNDN);
        mpfr_set_d(r26782, d2, MPFR_RNDN);
        mpfr_mul(r26783, r26781, r26782, MPFR_RNDN);
        mpfr_set_d(r26784, d3, MPFR_RNDN);
        mpfr_mul(r26785, r26781, r26784, MPFR_RNDN);
        mpfr_sub(r26786, r26783, r26785, MPFR_RNDN);
        mpfr_set_d(r26787, d4, MPFR_RNDN);
        mpfr_mul(r26788, r26787, r26781, MPFR_RNDN);
        mpfr_add(r26789, r26786, r26788, MPFR_RNDN);
        mpfr_mul(r26790, r26781, r26781, MPFR_RNDN);
        mpfr_sub(r26791, r26789, r26790, MPFR_RNDN);
        return mpfr_get_d(r26791, MPFR_RNDN);
}

static mpfr_t r26792, r26793, r26794, r26795, r26796, r26797, r26798, r26799;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26792);
        mpfr_init(r26793);
        mpfr_init(r26794);
        mpfr_init(r26795);
        mpfr_init(r26796);
        mpfr_init(r26797);
        mpfr_init(r26798);
        mpfr_init(r26799);
}

double f_fm(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r26792, d4, MPFR_RNDN);
        mpfr_set_d(r26793, d3, MPFR_RNDN);
        mpfr_sub(r26794, r26792, r26793, MPFR_RNDN);
        mpfr_set_d(r26795, d2, MPFR_RNDN);
        mpfr_set_d(r26796, d1, MPFR_RNDN);
        mpfr_sub(r26797, r26795, r26796, MPFR_RNDN);
        mpfr_add(r26798, r26794, r26797, MPFR_RNDN);
        mpfr_mul(r26799, r26798, r26796, MPFR_RNDN);
        return mpfr_get_d(r26799, MPFR_RNDN);
}

static mpfr_t r26800, r26801, r26802, r26803, r26804, r26805, r26806, r26807;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r26800);
        mpfr_init(r26801);
        mpfr_init(r26802);
        mpfr_init(r26803);
        mpfr_init(r26804);
        mpfr_init(r26805);
        mpfr_init(r26806);
        mpfr_init(r26807);
}

double f_dm(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r26800, d4, MPFR_RNDN);
        mpfr_set_d(r26801, d3, MPFR_RNDN);
        mpfr_sub(r26802, r26800, r26801, MPFR_RNDN);
        mpfr_set_d(r26803, d2, MPFR_RNDN);
        mpfr_set_d(r26804, d1, MPFR_RNDN);
        mpfr_sub(r26805, r26803, r26804, MPFR_RNDN);
        mpfr_add(r26806, r26802, r26805, MPFR_RNDN);
        mpfr_mul(r26807, r26806, r26804, MPFR_RNDN);
        return mpfr_get_d(r26807, MPFR_RNDN);
}