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

char *name = "FastMath test2";

double f_if(float d1, float d2) {
        float r18993 = d1;
        float r18994 = 10.0f;
        float r18995 = r18993 * r18994;
        float r18996 = d2;
        float r18997 = r18993 * r18996;
        float r18998 = r18995 + r18997;
        float r18999 = 20.0f;
        float r19000 = r18993 * r18999;
        float r19001 = r18998 + r19000;
        return r19001;
}

double f_id(double d1, double d2) {
        double r19002 = d1;
        double r19003 = 10.0;
        double r19004 = r19002 * r19003;
        double r19005 = d2;
        double r19006 = r19002 * r19005;
        double r19007 = r19004 + r19006;
        double r19008 = 20.0;
        double r19009 = r19002 * r19008;
        double r19010 = r19007 + r19009;
        return r19010;
}


double f_of(float d1, float d2) {
        float r19011 = d1;
        float r19012 = 20.0f;
        float r19013 = 10.0f;
        float r19014 = r19012 + r19013;
        float r19015 = d2;
        float r19016 = r19011 * r19015;
        float r19017 = fma(r19011, r19014, r19016);
        return r19017;
}

double f_od(double d1, double d2) {
        double r19018 = d1;
        double r19019 = 20.0;
        double r19020 = 10.0;
        double r19021 = r19019 + r19020;
        double r19022 = d2;
        double r19023 = r19018 * r19022;
        double r19024 = fma(r19018, r19021, r19023);
        return r19024;
}

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 r19025, r19026, r19027, r19028, r19029, r19030, r19031, r19032, r19033;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r19025);
        mpfr_init_set_str(r19026, "10", 10, MPFR_RNDN);
        mpfr_init(r19027);
        mpfr_init(r19028);
        mpfr_init(r19029);
        mpfr_init(r19030);
        mpfr_init_set_str(r19031, "20", 10, MPFR_RNDN);
        mpfr_init(r19032);
        mpfr_init(r19033);
}

double f_im(double d1, double d2) {
        mpfr_set_d(r19025, d1, MPFR_RNDN);
        ;
        mpfr_mul(r19027, r19025, r19026, MPFR_RNDN);
        mpfr_set_d(r19028, d2, MPFR_RNDN);
        mpfr_mul(r19029, r19025, r19028, MPFR_RNDN);
        mpfr_add(r19030, r19027, r19029, MPFR_RNDN);
        ;
        mpfr_mul(r19032, r19025, r19031, MPFR_RNDN);
        mpfr_add(r19033, r19030, r19032, MPFR_RNDN);
        return mpfr_get_d(r19033, MPFR_RNDN);
}

static mpfr_t r19034, r19035, r19036, r19037, r19038, r19039, r19040;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19034);
        mpfr_init_set_str(r19035, "20", 10, MPFR_RNDN);
        mpfr_init_set_str(r19036, "10", 10, MPFR_RNDN);
        mpfr_init(r19037);
        mpfr_init(r19038);
        mpfr_init(r19039);
        mpfr_init(r19040);
}

double f_fm(double d1, double d2) {
        mpfr_set_d(r19034, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r19037, r19035, r19036, MPFR_RNDN);
        mpfr_set_d(r19038, d2, MPFR_RNDN);
        mpfr_mul(r19039, r19034, r19038, MPFR_RNDN);
        mpfr_fma(r19040, r19034, r19037, r19039, MPFR_RNDN);
        return mpfr_get_d(r19040, MPFR_RNDN);
}

static mpfr_t r19041, r19042, r19043, r19044, r19045, r19046, r19047;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19041);
        mpfr_init_set_str(r19042, "20", 10, MPFR_RNDN);
        mpfr_init_set_str(r19043, "10", 10, MPFR_RNDN);
        mpfr_init(r19044);
        mpfr_init(r19045);
        mpfr_init(r19046);
        mpfr_init(r19047);
}

double f_dm(double d1, double d2) {
        mpfr_set_d(r19041, d1, MPFR_RNDN);
        ;
        ;
        mpfr_add(r19044, r19042, r19043, MPFR_RNDN);
        mpfr_set_d(r19045, d2, MPFR_RNDN);
        mpfr_mul(r19046, r19041, r19045, MPFR_RNDN);
        mpfr_fma(r19047, r19041, r19044, r19046, MPFR_RNDN);
        return mpfr_get_d(r19047, MPFR_RNDN);
}