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

char *name = "FastMath test3";

double f_if(float d1, float d2, float d3) {
        float r26001 = d1;
        float r26002 = 3;
        float r26003 = r26001 * r26002;
        float r26004 = d2;
        float r26005 = r26001 * r26004;
        float r26006 = r26003 + r26005;
        float r26007 = d3;
        float r26008 = r26001 * r26007;
        float r26009 = r26006 + r26008;
        return r26009;
}

double f_id(double d1, double d2, double d3) {
        double r26010 = d1;
        double r26011 = 3;
        double r26012 = r26010 * r26011;
        double r26013 = d2;
        double r26014 = r26010 * r26013;
        double r26015 = r26012 + r26014;
        double r26016 = d3;
        double r26017 = r26010 * r26016;
        double r26018 = r26015 + r26017;
        return r26018;
}


double f_of(float d1, float d2, float d3) {
        float r26019 = d1;
        float r26020 = 3;
        float r26021 = r26019 * r26020;
        float r26022 = d2;
        float r26023 = r26019 * r26022;
        float r26024 = r26021 + r26023;
        float r26025 = d3;
        float r26026 = r26019 * r26025;
        float r26027 = r26024 + r26026;
        return r26027;
}

double f_od(double d1, double d2, double d3) {
        double r26028 = d1;
        double r26029 = 3;
        double r26030 = r26028 * r26029;
        double r26031 = d2;
        double r26032 = r26028 * r26031;
        double r26033 = r26030 + r26032;
        double r26034 = d3;
        double r26035 = r26028 * r26034;
        double r26036 = r26033 + r26035;
        return r26036;
}

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 r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26037);
        mpfr_init_set_str(r26038, "3", 10, MPFR_RNDN);
        mpfr_init(r26039);
        mpfr_init(r26040);
        mpfr_init(r26041);
        mpfr_init(r26042);
        mpfr_init(r26043);
        mpfr_init(r26044);
        mpfr_init(r26045);
}

double f_im(double d1, double d2, double d3) {
        mpfr_set_d(r26037, d1, MPFR_RNDN);
        ;
        mpfr_mul(r26039, r26037, r26038, MPFR_RNDN);
        mpfr_set_d(r26040, d2, MPFR_RNDN);
        mpfr_mul(r26041, r26037, r26040, MPFR_RNDN);
        mpfr_add(r26042, r26039, r26041, MPFR_RNDN);
        mpfr_set_d(r26043, d3, MPFR_RNDN);
        mpfr_mul(r26044, r26037, r26043, MPFR_RNDN);
        mpfr_add(r26045, r26042, r26044, MPFR_RNDN);
        return mpfr_get_d(r26045, MPFR_RNDN);
}

static mpfr_t r26046, r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26046);
        mpfr_init_set_str(r26047, "3", 10, MPFR_RNDN);
        mpfr_init(r26048);
        mpfr_init(r26049);
        mpfr_init(r26050);
        mpfr_init(r26051);
        mpfr_init(r26052);
        mpfr_init(r26053);
        mpfr_init(r26054);
}

double f_fm(double d1, double d2, double d3) {
        mpfr_set_d(r26046, d1, MPFR_RNDN);
        ;
        mpfr_mul(r26048, r26046, r26047, MPFR_RNDN);
        mpfr_set_d(r26049, d2, MPFR_RNDN);
        mpfr_mul(r26050, r26046, r26049, MPFR_RNDN);
        mpfr_add(r26051, r26048, r26050, MPFR_RNDN);
        mpfr_set_d(r26052, d3, MPFR_RNDN);
        mpfr_mul(r26053, r26046, r26052, MPFR_RNDN);
        mpfr_add(r26054, r26051, r26053, MPFR_RNDN);
        return mpfr_get_d(r26054, MPFR_RNDN);
}

static mpfr_t r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26055);
        mpfr_init_set_str(r26056, "3", 10, MPFR_RNDN);
        mpfr_init(r26057);
        mpfr_init(r26058);
        mpfr_init(r26059);
        mpfr_init(r26060);
        mpfr_init(r26061);
        mpfr_init(r26062);
        mpfr_init(r26063);
}

double f_dm(double d1, double d2, double d3) {
        mpfr_set_d(r26055, d1, MPFR_RNDN);
        ;
        mpfr_mul(r26057, r26055, r26056, MPFR_RNDN);
        mpfr_set_d(r26058, d2, MPFR_RNDN);
        mpfr_mul(r26059, r26055, r26058, MPFR_RNDN);
        mpfr_add(r26060, r26057, r26059, MPFR_RNDN);
        mpfr_set_d(r26061, d3, MPFR_RNDN);
        mpfr_mul(r26062, r26055, r26061, MPFR_RNDN);
        mpfr_add(r26063, r26060, r26062, MPFR_RNDN);
        return mpfr_get_d(r26063, MPFR_RNDN);
}