#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 r19769 = d1;
        float r19770 = 10.0f;
        float r19771 = r19769 * r19770;
        float r19772 = d2;
        float r19773 = r19769 * r19772;
        float r19774 = r19771 + r19773;
        float r19775 = 20.0f;
        float r19776 = r19769 * r19775;
        float r19777 = r19774 + r19776;
        return r19777;
}

double f_id(double d1, double d2) {
        double r19778 = d1;
        double r19779 = 10.0;
        double r19780 = r19778 * r19779;
        double r19781 = d2;
        double r19782 = r19778 * r19781;
        double r19783 = r19780 + r19782;
        double r19784 = 20.0;
        double r19785 = r19778 * r19784;
        double r19786 = r19783 + r19785;
        return r19786;
}


double f_of(float d1, float d2) {
        float r19787 = d1;
        float r19788 = 10.0f;
        float r19789 = d2;
        float r19790 = r19788 + r19789;
        float r19791 = 20.0f;
        float r19792 = r19790 + r19791;
        float r19793 = r19787 * r19792;
        return r19793;
}

double f_od(double d1, double d2) {
        double r19794 = d1;
        double r19795 = 10.0;
        double r19796 = d2;
        double r19797 = r19795 + r19796;
        double r19798 = 20.0;
        double r19799 = r19797 + r19798;
        double r19800 = r19794 * r19799;
        return r19800;
}

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 r19801, r19802, r19803, r19804, r19805, r19806, r19807, r19808, r19809;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r19801);
        mpfr_init_set_str(r19802, "10", 10, MPFR_RNDN);
        mpfr_init(r19803);
        mpfr_init(r19804);
        mpfr_init(r19805);
        mpfr_init(r19806);
        mpfr_init_set_str(r19807, "20", 10, MPFR_RNDN);
        mpfr_init(r19808);
        mpfr_init(r19809);
}

double f_im(double d1, double d2) {
        mpfr_set_d(r19801, d1, MPFR_RNDN);
        ;
        mpfr_mul(r19803, r19801, r19802, MPFR_RNDN);
        mpfr_set_d(r19804, d2, MPFR_RNDN);
        mpfr_mul(r19805, r19801, r19804, MPFR_RNDN);
        mpfr_add(r19806, r19803, r19805, MPFR_RNDN);
        ;
        mpfr_mul(r19808, r19801, r19807, MPFR_RNDN);
        mpfr_add(r19809, r19806, r19808, MPFR_RNDN);
        return mpfr_get_d(r19809, MPFR_RNDN);
}

static mpfr_t r19810, r19811, r19812, r19813, r19814, r19815, r19816;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19810);
        mpfr_init_set_str(r19811, "10", 10, MPFR_RNDN);
        mpfr_init(r19812);
        mpfr_init(r19813);
        mpfr_init_set_str(r19814, "20", 10, MPFR_RNDN);
        mpfr_init(r19815);
        mpfr_init(r19816);
}

double f_fm(double d1, double d2) {
        mpfr_set_d(r19810, d1, MPFR_RNDN);
        ;
        mpfr_set_d(r19812, d2, MPFR_RNDN);
        mpfr_add(r19813, r19811, r19812, MPFR_RNDN);
        ;
        mpfr_add(r19815, r19813, r19814, MPFR_RNDN);
        mpfr_mul(r19816, r19810, r19815, MPFR_RNDN);
        return mpfr_get_d(r19816, MPFR_RNDN);
}

static mpfr_t r19817, r19818, r19819, r19820, r19821, r19822, r19823;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r19817);
        mpfr_init_set_str(r19818, "10", 10, MPFR_RNDN);
        mpfr_init(r19819);
        mpfr_init(r19820);
        mpfr_init_set_str(r19821, "20", 10, MPFR_RNDN);
        mpfr_init(r19822);
        mpfr_init(r19823);
}

double f_dm(double d1, double d2) {
        mpfr_set_d(r19817, d1, MPFR_RNDN);
        ;
        mpfr_set_d(r19819, d2, MPFR_RNDN);
        mpfr_add(r19820, r19818, r19819, MPFR_RNDN);
        ;
        mpfr_add(r19822, r19820, r19821, MPFR_RNDN);
        mpfr_mul(r19823, r19817, r19822, MPFR_RNDN);
        return mpfr_get_d(r19823, MPFR_RNDN);
}