#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 r23564 = d1;
        float r23565 = 10.0f;
        float r23566 = r23564 * r23565;
        float r23567 = d2;
        float r23568 = r23564 * r23567;
        float r23569 = r23566 + r23568;
        float r23570 = 20.0f;
        float r23571 = r23564 * r23570;
        float r23572 = r23569 + r23571;
        return r23572;
}

double f_id(double d1, double d2) {
        double r23573 = d1;
        double r23574 = 10.0;
        double r23575 = r23573 * r23574;
        double r23576 = d2;
        double r23577 = r23573 * r23576;
        double r23578 = r23575 + r23577;
        double r23579 = 20.0;
        double r23580 = r23573 * r23579;
        double r23581 = r23578 + r23580;
        return r23581;
}


double f_of(float d1, float d2) {
        float r23582 = d2;
        float r23583 = 10.0f;
        float r23584 = r23582 + r23583;
        float r23585 = 20.0f;
        float r23586 = r23584 + r23585;
        float r23587 = d1;
        float r23588 = r23586 * r23587;
        return r23588;
}

double f_od(double d1, double d2) {
        double r23589 = d2;
        double r23590 = 10.0;
        double r23591 = r23589 + r23590;
        double r23592 = 20.0;
        double r23593 = r23591 + r23592;
        double r23594 = d1;
        double r23595 = r23593 * r23594;
        return r23595;
}

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 r23596, r23597, r23598, r23599, r23600, r23601, r23602, r23603, r23604;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r23596);
        mpfr_init_set_str(r23597, "10", 10, MPFR_RNDN);
        mpfr_init(r23598);
        mpfr_init(r23599);
        mpfr_init(r23600);
        mpfr_init(r23601);
        mpfr_init_set_str(r23602, "20", 10, MPFR_RNDN);
        mpfr_init(r23603);
        mpfr_init(r23604);
}

double f_im(double d1, double d2) {
        mpfr_set_d(r23596, d1, MPFR_RNDN);
        ;
        mpfr_mul(r23598, r23596, r23597, MPFR_RNDN);
        mpfr_set_d(r23599, d2, MPFR_RNDN);
        mpfr_mul(r23600, r23596, r23599, MPFR_RNDN);
        mpfr_add(r23601, r23598, r23600, MPFR_RNDN);
        ;
        mpfr_mul(r23603, r23596, r23602, MPFR_RNDN);
        mpfr_add(r23604, r23601, r23603, MPFR_RNDN);
        return mpfr_get_d(r23604, MPFR_RNDN);
}

static mpfr_t r23605, r23606, r23607, r23608, r23609, r23610, r23611;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r23605);
        mpfr_init_set_str(r23606, "10", 10, MPFR_RNDN);
        mpfr_init(r23607);
        mpfr_init_set_str(r23608, "20", 10, MPFR_RNDN);
        mpfr_init(r23609);
        mpfr_init(r23610);
        mpfr_init(r23611);
}

double f_fm(double d1, double d2) {
        mpfr_set_d(r23605, d2, MPFR_RNDN);
        ;
        mpfr_add(r23607, r23605, r23606, MPFR_RNDN);
        ;
        mpfr_add(r23609, r23607, r23608, MPFR_RNDN);
        mpfr_set_d(r23610, d1, MPFR_RNDN);
        mpfr_mul(r23611, r23609, r23610, MPFR_RNDN);
        return mpfr_get_d(r23611, MPFR_RNDN);
}

static mpfr_t r23612, r23613, r23614, r23615, r23616, r23617, r23618;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r23612);
        mpfr_init_set_str(r23613, "10", 10, MPFR_RNDN);
        mpfr_init(r23614);
        mpfr_init_set_str(r23615, "20", 10, MPFR_RNDN);
        mpfr_init(r23616);
        mpfr_init(r23617);
        mpfr_init(r23618);
}

double f_dm(double d1, double d2) {
        mpfr_set_d(r23612, d2, MPFR_RNDN);
        ;
        mpfr_add(r23614, r23612, r23613, MPFR_RNDN);
        ;
        mpfr_add(r23616, r23614, r23615, MPFR_RNDN);
        mpfr_set_d(r23617, d1, MPFR_RNDN);
        mpfr_mul(r23618, r23616, r23617, MPFR_RNDN);
        return mpfr_get_d(r23618, MPFR_RNDN);
}