#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 r20210 = d1;
        float r20211 = d2;
        float r20212 = r20210 * r20211;
        float r20213 = d3;
        float r20214 = r20210 * r20213;
        float r20215 = r20212 - r20214;
        float r20216 = d4;
        float r20217 = r20216 * r20210;
        float r20218 = r20215 + r20217;
        float r20219 = r20210 * r20210;
        float r20220 = r20218 - r20219;
        return r20220;
}

double f_id(double d1, double d2, double d3, double d4) {
        double r20221 = d1;
        double r20222 = d2;
        double r20223 = r20221 * r20222;
        double r20224 = d3;
        double r20225 = r20221 * r20224;
        double r20226 = r20223 - r20225;
        double r20227 = d4;
        double r20228 = r20227 * r20221;
        double r20229 = r20226 + r20228;
        double r20230 = r20221 * r20221;
        double r20231 = r20229 - r20230;
        return r20231;
}


double f_of(float d1, float d2, float d3, float d4) {
        float r20232 = d1;
        float r20233 = d4;
        float r20234 = d2;
        float r20235 = r20233 + r20234;
        float r20236 = d3;
        float r20237 = r20236 + r20232;
        float r20238 = r20235 - r20237;
        float r20239 = r20232 * r20238;
        return r20239;
}

double f_od(double d1, double d2, double d3, double d4) {
        double r20240 = d1;
        double r20241 = d4;
        double r20242 = d2;
        double r20243 = r20241 + r20242;
        double r20244 = d3;
        double r20245 = r20244 + r20240;
        double r20246 = r20243 - r20245;
        double r20247 = r20240 * r20246;
        return r20247;
}

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 r20248, r20249, r20250, r20251, r20252, r20253, r20254, r20255, r20256, r20257, r20258;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r20248);
        mpfr_init(r20249);
        mpfr_init(r20250);
        mpfr_init(r20251);
        mpfr_init(r20252);
        mpfr_init(r20253);
        mpfr_init(r20254);
        mpfr_init(r20255);
        mpfr_init(r20256);
        mpfr_init(r20257);
        mpfr_init(r20258);
}

double f_im(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r20248, d1, MPFR_RNDN);
        mpfr_set_d(r20249, d2, MPFR_RNDN);
        mpfr_mul(r20250, r20248, r20249, MPFR_RNDN);
        mpfr_set_d(r20251, d3, MPFR_RNDN);
        mpfr_mul(r20252, r20248, r20251, MPFR_RNDN);
        mpfr_sub(r20253, r20250, r20252, MPFR_RNDN);
        mpfr_set_d(r20254, d4, MPFR_RNDN);
        mpfr_mul(r20255, r20254, r20248, MPFR_RNDN);
        mpfr_add(r20256, r20253, r20255, MPFR_RNDN);
        mpfr_mul(r20257, r20248, r20248, MPFR_RNDN);
        mpfr_sub(r20258, r20256, r20257, MPFR_RNDN);
        return mpfr_get_d(r20258, MPFR_RNDN);
}

static mpfr_t r20259, r20260, r20261, r20262, r20263, r20264, r20265, r20266;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20259);
        mpfr_init(r20260);
        mpfr_init(r20261);
        mpfr_init(r20262);
        mpfr_init(r20263);
        mpfr_init(r20264);
        mpfr_init(r20265);
        mpfr_init(r20266);
}

double f_fm(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r20259, d1, MPFR_RNDN);
        mpfr_set_d(r20260, d4, MPFR_RNDN);
        mpfr_set_d(r20261, d2, MPFR_RNDN);
        mpfr_add(r20262, r20260, r20261, MPFR_RNDN);
        mpfr_set_d(r20263, d3, MPFR_RNDN);
        mpfr_add(r20264, r20263, r20259, MPFR_RNDN);
        mpfr_sub(r20265, r20262, r20264, MPFR_RNDN);
        mpfr_mul(r20266, r20259, r20265, MPFR_RNDN);
        return mpfr_get_d(r20266, MPFR_RNDN);
}

static mpfr_t r20267, r20268, r20269, r20270, r20271, r20272, r20273, r20274;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20267);
        mpfr_init(r20268);
        mpfr_init(r20269);
        mpfr_init(r20270);
        mpfr_init(r20271);
        mpfr_init(r20272);
        mpfr_init(r20273);
        mpfr_init(r20274);
}

double f_dm(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r20267, d1, MPFR_RNDN);
        mpfr_set_d(r20268, d4, MPFR_RNDN);
        mpfr_set_d(r20269, d2, MPFR_RNDN);
        mpfr_add(r20270, r20268, r20269, MPFR_RNDN);
        mpfr_set_d(r20271, d3, MPFR_RNDN);
        mpfr_add(r20272, r20271, r20267, MPFR_RNDN);
        mpfr_sub(r20273, r20270, r20272, MPFR_RNDN);
        mpfr_mul(r20274, r20267, r20273, MPFR_RNDN);
        return mpfr_get_d(r20274, MPFR_RNDN);
}