#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 r26740 = d1;
        float r26741 = d2;
        float r26742 = r26740 * r26741;
        float r26743 = d3;
        float r26744 = r26740 * r26743;
        float r26745 = r26742 - r26744;
        float r26746 = d4;
        float r26747 = r26746 * r26740;
        float r26748 = r26745 + r26747;
        float r26749 = r26740 * r26740;
        float r26750 = r26748 - r26749;
        return r26750;
}

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


double f_of(float d1, float d2, float d3, float d4) {
        float r26762 = d4;
        float r26763 = d3;
        float r26764 = r26762 - r26763;
        float r26765 = d2;
        float r26766 = d1;
        float r26767 = r26765 - r26766;
        float r26768 = r26764 + r26767;
        float r26769 = r26768 * r26766;
        return r26769;
}

double f_od(double d1, double d2, double d3, double d4) {
        double r26770 = d4;
        double r26771 = d3;
        double r26772 = r26770 - r26771;
        double r26773 = d2;
        double r26774 = d1;
        double r26775 = r26773 - r26774;
        double r26776 = r26772 + r26775;
        double r26777 = r26776 * r26774;
        return r26777;
}

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 r26778, r26779, r26780, r26781, r26782, r26783, r26784, r26785, r26786, r26787, r26788;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r26778);
        mpfr_init(r26779);
        mpfr_init(r26780);
        mpfr_init(r26781);
        mpfr_init(r26782);
        mpfr_init(r26783);
        mpfr_init(r26784);
        mpfr_init(r26785);
        mpfr_init(r26786);
        mpfr_init(r26787);
        mpfr_init(r26788);
}

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

static mpfr_t r26789, r26790, r26791, r26792, r26793, r26794, r26795, r26796;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r26789);
        mpfr_init(r26790);
        mpfr_init(r26791);
        mpfr_init(r26792);
        mpfr_init(r26793);
        mpfr_init(r26794);
        mpfr_init(r26795);
        mpfr_init(r26796);
}

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

static mpfr_t r26797, r26798, r26799, r26800, r26801, r26802, r26803, r26804;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r26797);
        mpfr_init(r26798);
        mpfr_init(r26799);
        mpfr_init(r26800);
        mpfr_init(r26801);
        mpfr_init(r26802);
        mpfr_init(r26803);
        mpfr_init(r26804);
}

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