#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 r28991 = d1;
        float r28992 = d2;
        float r28993 = r28991 * r28992;
        float r28994 = d3;
        float r28995 = r28991 * r28994;
        float r28996 = r28993 - r28995;
        float r28997 = d4;
        float r28998 = r28997 * r28991;
        float r28999 = r28996 + r28998;
        float r29000 = r28991 * r28991;
        float r29001 = r28999 - r29000;
        return r29001;
}

double f_id(double d1, double d2, double d3, double d4) {
        double r29002 = d1;
        double r29003 = d2;
        double r29004 = r29002 * r29003;
        double r29005 = d3;
        double r29006 = r29002 * r29005;
        double r29007 = r29004 - r29006;
        double r29008 = d4;
        double r29009 = r29008 * r29002;
        double r29010 = r29007 + r29009;
        double r29011 = r29002 * r29002;
        double r29012 = r29010 - r29011;
        return r29012;
}


double f_of(float d1, float d2, float d3, float d4) {
        float r29013 = d4;
        float r29014 = d3;
        float r29015 = r29013 - r29014;
        float r29016 = d2;
        float r29017 = d1;
        float r29018 = r29016 - r29017;
        float r29019 = r29015 + r29018;
        float r29020 = r29019 * r29017;
        return r29020;
}

double f_od(double d1, double d2, double d3, double d4) {
        double r29021 = d4;
        double r29022 = d3;
        double r29023 = r29021 - r29022;
        double r29024 = d2;
        double r29025 = d1;
        double r29026 = r29024 - r29025;
        double r29027 = r29023 + r29026;
        double r29028 = r29027 * r29025;
        return r29028;
}

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 r29029, r29030, r29031, r29032, r29033, r29034, r29035, r29036, r29037, r29038, r29039;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r29029);
        mpfr_init(r29030);
        mpfr_init(r29031);
        mpfr_init(r29032);
        mpfr_init(r29033);
        mpfr_init(r29034);
        mpfr_init(r29035);
        mpfr_init(r29036);
        mpfr_init(r29037);
        mpfr_init(r29038);
        mpfr_init(r29039);
}

double f_im(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r29029, d1, MPFR_RNDN);
        mpfr_set_d(r29030, d2, MPFR_RNDN);
        mpfr_mul(r29031, r29029, r29030, MPFR_RNDN);
        mpfr_set_d(r29032, d3, MPFR_RNDN);
        mpfr_mul(r29033, r29029, r29032, MPFR_RNDN);
        mpfr_sub(r29034, r29031, r29033, MPFR_RNDN);
        mpfr_set_d(r29035, d4, MPFR_RNDN);
        mpfr_mul(r29036, r29035, r29029, MPFR_RNDN);
        mpfr_add(r29037, r29034, r29036, MPFR_RNDN);
        mpfr_mul(r29038, r29029, r29029, MPFR_RNDN);
        mpfr_sub(r29039, r29037, r29038, MPFR_RNDN);
        return mpfr_get_d(r29039, MPFR_RNDN);
}

static mpfr_t r29040, r29041, r29042, r29043, r29044, r29045, r29046, r29047;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r29040);
        mpfr_init(r29041);
        mpfr_init(r29042);
        mpfr_init(r29043);
        mpfr_init(r29044);
        mpfr_init(r29045);
        mpfr_init(r29046);
        mpfr_init(r29047);
}

double f_fm(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r29040, d4, MPFR_RNDN);
        mpfr_set_d(r29041, d3, MPFR_RNDN);
        mpfr_sub(r29042, r29040, r29041, MPFR_RNDN);
        mpfr_set_d(r29043, d2, MPFR_RNDN);
        mpfr_set_d(r29044, d1, MPFR_RNDN);
        mpfr_sub(r29045, r29043, r29044, MPFR_RNDN);
        mpfr_add(r29046, r29042, r29045, MPFR_RNDN);
        mpfr_mul(r29047, r29046, r29044, MPFR_RNDN);
        return mpfr_get_d(r29047, MPFR_RNDN);
}

static mpfr_t r29048, r29049, r29050, r29051, r29052, r29053, r29054, r29055;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r29048);
        mpfr_init(r29049);
        mpfr_init(r29050);
        mpfr_init(r29051);
        mpfr_init(r29052);
        mpfr_init(r29053);
        mpfr_init(r29054);
        mpfr_init(r29055);
}

double f_dm(double d1, double d2, double d3, double d4) {
        mpfr_set_d(r29048, d4, MPFR_RNDN);
        mpfr_set_d(r29049, d3, MPFR_RNDN);
        mpfr_sub(r29050, r29048, r29049, MPFR_RNDN);
        mpfr_set_d(r29051, d2, MPFR_RNDN);
        mpfr_set_d(r29052, d1, MPFR_RNDN);
        mpfr_sub(r29053, r29051, r29052, MPFR_RNDN);
        mpfr_add(r29054, r29050, r29053, MPFR_RNDN);
        mpfr_mul(r29055, r29054, r29052, MPFR_RNDN);
        return mpfr_get_d(r29055, MPFR_RNDN);
}