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

char *name = "Development.Shake.Progress:message from shake-0.15.5";

double f_if(float x, float y) {
        float r55383 = x;
        float r55384 = 100.0;
        float r55385 = r55383 * r55384;
        float r55386 = y;
        float r55387 = r55383 + r55386;
        float r55388 = r55385 / r55387;
        return r55388;
}

double f_id(double x, double y) {
        double r55389 = x;
        double r55390 = 100.0;
        double r55391 = r55389 * r55390;
        double r55392 = y;
        double r55393 = r55389 + r55392;
        double r55394 = r55391 / r55393;
        return r55394;
}


double f_of(float x, float y) {
        float r55395 = x;
        float r55396 = y;
        float r55397 = r55395 + r55396;
        float r55398 = 100.0;
        float r55399 = r55397 / r55398;
        float r55400 = r55395 / r55399;
        return r55400;
}

double f_od(double x, double y) {
        double r55401 = x;
        double r55402 = y;
        double r55403 = r55401 + r55402;
        double r55404 = 100.0;
        double r55405 = r55403 / r55404;
        double r55406 = r55401 / r55405;
        return r55406;
}

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 r55407, r55408, r55409, r55410, r55411, r55412;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r55407);
        mpfr_init_set_str(r55408, "100.0", 10, MPFR_RNDN);
        mpfr_init(r55409);
        mpfr_init(r55410);
        mpfr_init(r55411);
        mpfr_init(r55412);
}

double f_im(double x, double y) {
        mpfr_set_d(r55407, x, MPFR_RNDN);
        ;
        mpfr_mul(r55409, r55407, r55408, MPFR_RNDN);
        mpfr_set_d(r55410, y, MPFR_RNDN);
        mpfr_add(r55411, r55407, r55410, MPFR_RNDN);
        mpfr_div(r55412, r55409, r55411, MPFR_RNDN);
        return mpfr_get_d(r55412, MPFR_RNDN);
}

static mpfr_t r55413, r55414, r55415, r55416, r55417, r55418;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r55413);
        mpfr_init(r55414);
        mpfr_init(r55415);
        mpfr_init_set_str(r55416, "100.0", 10, MPFR_RNDN);
        mpfr_init(r55417);
        mpfr_init(r55418);
}

double f_fm(double x, double y) {
        mpfr_set_d(r55413, x, MPFR_RNDN);
        mpfr_set_d(r55414, y, MPFR_RNDN);
        mpfr_add(r55415, r55413, r55414, MPFR_RNDN);
        ;
        mpfr_div(r55417, r55415, r55416, MPFR_RNDN);
        mpfr_div(r55418, r55413, r55417, MPFR_RNDN);
        return mpfr_get_d(r55418, MPFR_RNDN);
}

static mpfr_t r55419, r55420, r55421, r55422, r55423, r55424;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r55419);
        mpfr_init(r55420);
        mpfr_init(r55421);
        mpfr_init_set_str(r55422, "100.0", 10, MPFR_RNDN);
        mpfr_init(r55423);
        mpfr_init(r55424);
}

double f_dm(double x, double y) {
        mpfr_set_d(r55419, x, MPFR_RNDN);
        mpfr_set_d(r55420, y, MPFR_RNDN);
        mpfr_add(r55421, r55419, r55420, MPFR_RNDN);
        ;
        mpfr_div(r55423, r55421, r55422, MPFR_RNDN);
        mpfr_div(r55424, r55419, r55423, MPFR_RNDN);
        return mpfr_get_d(r55424, MPFR_RNDN);
}