#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 r11731 = d1;
        float r11732 = 10.0f;
        float r11733 = r11731 * r11732;
        float r11734 = d2;
        float r11735 = r11731 * r11734;
        float r11736 = r11733 + r11735;
        float r11737 = 20.0f;
        float r11738 = r11731 * r11737;
        float r11739 = r11736 + r11738;
        return r11739;
}

double f_id(double d1, double d2) {
        double r11740 = d1;
        double r11741 = 10.0;
        double r11742 = r11740 * r11741;
        double r11743 = d2;
        double r11744 = r11740 * r11743;
        double r11745 = r11742 + r11744;
        double r11746 = 20.0;
        double r11747 = r11740 * r11746;
        double r11748 = r11745 + r11747;
        return r11748;
}


double f_of(float d1, float d2) {
        float r11749 = d2;
        float r11750 = 20.0f;
        float r11751 = 10.0f;
        float r11752 = r11750 + r11751;
        float r11753 = r11749 + r11752;
        float r11754 = d1;
        float r11755 = r11753 * r11754;
        return r11755;
}

double f_od(double d1, double d2) {
        double r11756 = d2;
        double r11757 = 20.0;
        double r11758 = 10.0;
        double r11759 = r11757 + r11758;
        double r11760 = r11756 + r11759;
        double r11761 = d1;
        double r11762 = r11760 * r11761;
        return r11762;
}

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 r11763, r11764, r11765, r11766, r11767, r11768, r11769, r11770, r11771;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r11763);
        mpfr_init_set_str(r11764, "10", 10, MPFR_RNDN);
        mpfr_init(r11765);
        mpfr_init(r11766);
        mpfr_init(r11767);
        mpfr_init(r11768);
        mpfr_init_set_str(r11769, "20", 10, MPFR_RNDN);
        mpfr_init(r11770);
        mpfr_init(r11771);
}

double f_im(double d1, double d2) {
        mpfr_set_d(r11763, d1, MPFR_RNDN);
        ;
        mpfr_mul(r11765, r11763, r11764, MPFR_RNDN);
        mpfr_set_d(r11766, d2, MPFR_RNDN);
        mpfr_mul(r11767, r11763, r11766, MPFR_RNDN);
        mpfr_add(r11768, r11765, r11767, MPFR_RNDN);
        ;
        mpfr_mul(r11770, r11763, r11769, MPFR_RNDN);
        mpfr_add(r11771, r11768, r11770, MPFR_RNDN);
        return mpfr_get_d(r11771, MPFR_RNDN);
}

static mpfr_t r11772, r11773, r11774, r11775, r11776, r11777, r11778;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r11772);
        mpfr_init_set_str(r11773, "20", 10, MPFR_RNDN);
        mpfr_init_set_str(r11774, "10", 10, MPFR_RNDN);
        mpfr_init(r11775);
        mpfr_init(r11776);
        mpfr_init(r11777);
        mpfr_init(r11778);
}

double f_fm(double d1, double d2) {
        mpfr_set_d(r11772, d2, MPFR_RNDN);
        ;
        ;
        mpfr_add(r11775, r11773, r11774, MPFR_RNDN);
        mpfr_add(r11776, r11772, r11775, MPFR_RNDN);
        mpfr_set_d(r11777, d1, MPFR_RNDN);
        mpfr_mul(r11778, r11776, r11777, MPFR_RNDN);
        return mpfr_get_d(r11778, MPFR_RNDN);
}

static mpfr_t r11779, r11780, r11781, r11782, r11783, r11784, r11785;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r11779);
        mpfr_init_set_str(r11780, "20", 10, MPFR_RNDN);
        mpfr_init_set_str(r11781, "10", 10, MPFR_RNDN);
        mpfr_init(r11782);
        mpfr_init(r11783);
        mpfr_init(r11784);
        mpfr_init(r11785);
}

double f_dm(double d1, double d2) {
        mpfr_set_d(r11779, d2, MPFR_RNDN);
        ;
        ;
        mpfr_add(r11782, r11780, r11781, MPFR_RNDN);
        mpfr_add(r11783, r11779, r11782, MPFR_RNDN);
        mpfr_set_d(r11784, d1, MPFR_RNDN);
        mpfr_mul(r11785, r11783, r11784, MPFR_RNDN);
        return mpfr_get_d(r11785, MPFR_RNDN);
}

