#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 r28209 = d1;
        float r28210 = d2;
        float r28211 = r28209 * r28210;
        float r28212 = d3;
        float r28213 = r28209 * r28212;
        float r28214 = r28211 - r28213;
        float r28215 = d4;
        float r28216 = r28215 * r28209;
        float r28217 = r28214 + r28216;
        float r28218 = r28209 * r28209;
        float r28219 = r28217 - r28218;
        return r28219;
}

double f_id(double d1, double d2, double d3, double d4) {
        double r28220 = d1;
        double r28221 = d2;
        double r28222 = r28220 * r28221;
        double r28223 = d3;
        double r28224 = r28220 * r28223;
        double r28225 = r28222 - r28224;
        double r28226 = d4;
        double r28227 = r28226 * r28220;
        double r28228 = r28225 + r28227;
        double r28229 = r28220 * r28220;
        double r28230 = r28228 - r28229;
        return r28230;
}


double f_of(float d1, float d2, float d3, float d4) {
        float r28231 = d1;
        float r28232 = d4;
        float r28233 = r28232 - r28231;
        float r28234 = d2;
        float r28235 = d3;
        float r28236 = r28234 - r28235;
        float r28237 = r28233 + r28236;
        float r28238 = r28231 * r28237;
        return r28238;
}

double f_od(double d1, double d2, double d3, double d4) {
        double r28239 = d1;
        double r28240 = d4;
        double r28241 = r28240 - r28239;
        double r28242 = d2;
        double r28243 = d3;
        double r28244 = r28242 - r28243;
        double r28245 = r28241 + r28244;
        double r28246 = r28239 * r28245;
        return r28246;
}

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 r28247, r28248, r28249, r28250, r28251, r28252, r28253, r28254, r28255, r28256, r28257;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r28247);
        mpfr_init(r28248);
        mpfr_init(r28249);
        mpfr_init(r28250);
        mpfr_init(r28251);
        mpfr_init(r28252);
        mpfr_init(r28253);
        mpfr_init(r28254);
        mpfr_init(r28255);
        mpfr_init(r28256);
        mpfr_init(r28257);
}

double f_im(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r28247, d1, MPFR_RNDN);
        mpfr_set_d(r28248, d2, MPFR_RNDN);
        mpfr_mul(r28249, r28247, r28248, MPFR_RNDN);
        mpfr_set_d(r28250, d3, MPFR_RNDN);
        mpfr_mul(r28251, r28247, r28250, MPFR_RNDN);
        mpfr_sub(r28252, r28249, r28251, MPFR_RNDN);
        mpfr_set_d(r28253, d4, MPFR_RNDN);
        mpfr_mul(r28254, r28253, r28247, MPFR_RNDN);
        mpfr_add(r28255, r28252, r28254, MPFR_RNDN);
        mpfr_mul(r28256, r28247, r28247, MPFR_RNDN);
        mpfr_sub(r28257, r28255, r28256, MPFR_RNDN);
        return mpfr_get_d(r28257, MPFR_RNDN);
}

static mpfr_t r28258, r28259, r28260, r28261, r28262, r28263, r28264, r28265;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r28258);
        mpfr_init(r28259);
        mpfr_init(r28260);
        mpfr_init(r28261);
        mpfr_init(r28262);
        mpfr_init(r28263);
        mpfr_init(r28264);
        mpfr_init(r28265);
}

double f_fm(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r28258, d1, MPFR_RNDN);
        mpfr_set_d(r28259, d4, MPFR_RNDN);
        mpfr_sub(r28260, r28259, r28258, MPFR_RNDN);
        mpfr_set_d(r28261, d2, MPFR_RNDN);
        mpfr_set_d(r28262, d3, MPFR_RNDN);
        mpfr_sub(r28263, r28261, r28262, MPFR_RNDN);
        mpfr_add(r28264, r28260, r28263, MPFR_RNDN);
        mpfr_mul(r28265, r28258, r28264, MPFR_RNDN);
        return mpfr_get_d(r28265, MPFR_RNDN);
}

static mpfr_t r28266, r28267, r28268, r28269, r28270, r28271, r28272, r28273;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r28266);
        mpfr_init(r28267);
        mpfr_init(r28268);
        mpfr_init(r28269);
        mpfr_init(r28270);
        mpfr_init(r28271);
        mpfr_init(r28272);
        mpfr_init(r28273);
}

double f_dm(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r28266, d1, MPFR_RNDN);
        mpfr_set_d(r28267, d4, MPFR_RNDN);
        mpfr_sub(r28268, r28267, r28266, MPFR_RNDN);
        mpfr_set_d(r28269, d2, MPFR_RNDN);
        mpfr_set_d(r28270, d3, MPFR_RNDN);
        mpfr_sub(r28271, r28269, r28270, MPFR_RNDN);
        mpfr_add(r28272, r28268, r28271, MPFR_RNDN);
        mpfr_mul(r28273, r28266, r28272, MPFR_RNDN);
        return mpfr_get_d(r28273, MPFR_RNDN);
}