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

char *name = "Quotient of sum of exps";

double f_if(float a, float b) {
        float r27325 = a;
        float r27326 = exp(r27325);
        float r27327 = b;
        float r27328 = exp(r27327);
        float r27329 = r27326 + r27328;
        float r27330 = r27326 / r27329;
        return r27330;
}

double f_id(double a, double b) {
        double r27331 = a;
        double r27332 = exp(r27331);
        double r27333 = b;
        double r27334 = exp(r27333);
        double r27335 = r27332 + r27334;
        double r27336 = r27332 / r27335;
        return r27336;
}


double f_of(float a, float b) {
        float r27337 = a;
        float r27338 = exp(r27337);
        float r27339 = b;
        float r27340 = exp(r27339);
        float r27341 = r27338 + r27340;
        float r27342 = log(r27341);
        float r27343 = r27337 - r27342;
        float r27344 = exp(r27343);
        return r27344;
}

double f_od(double a, double b) {
        double r27345 = a;
        double r27346 = exp(r27345);
        double r27347 = b;
        double r27348 = exp(r27347);
        double r27349 = r27346 + r27348;
        double r27350 = log(r27349);
        double r27351 = r27345 - r27350;
        double r27352 = exp(r27351);
        return r27352;
}

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 r27353, r27354, r27355, r27356, r27357, r27358;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r27353);
        mpfr_init(r27354);
        mpfr_init(r27355);
        mpfr_init(r27356);
        mpfr_init(r27357);
        mpfr_init(r27358);
}

double f_im(double a, double b) {
        mpfr_set_d(r27353, a, MPFR_RNDN);
        mpfr_exp(r27354, r27353, MPFR_RNDN);
        mpfr_set_d(r27355, b, MPFR_RNDN);
        mpfr_exp(r27356, r27355, MPFR_RNDN);
        mpfr_add(r27357, r27354, r27356, MPFR_RNDN);
        mpfr_div(r27358, r27354, r27357, MPFR_RNDN);
        return mpfr_get_d(r27358, MPFR_RNDN);
}

static mpfr_t r27359, r27360, r27361, r27362, r27363, r27364, r27365, r27366;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27359);
        mpfr_init(r27360);
        mpfr_init(r27361);
        mpfr_init(r27362);
        mpfr_init(r27363);
        mpfr_init(r27364);
        mpfr_init(r27365);
        mpfr_init(r27366);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27359, a, MPFR_RNDN);
        mpfr_exp(r27360, r27359, MPFR_RNDN);
        mpfr_set_d(r27361, b, MPFR_RNDN);
        mpfr_exp(r27362, r27361, MPFR_RNDN);
        mpfr_add(r27363, r27360, r27362, MPFR_RNDN);
        mpfr_log(r27364, r27363, MPFR_RNDN);
        mpfr_sub(r27365, r27359, r27364, MPFR_RNDN);
        mpfr_exp(r27366, r27365, MPFR_RNDN);
        return mpfr_get_d(r27366, MPFR_RNDN);
}

static mpfr_t r27367, r27368, r27369, r27370, r27371, r27372, r27373, r27374;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27367);
        mpfr_init(r27368);
        mpfr_init(r27369);
        mpfr_init(r27370);
        mpfr_init(r27371);
        mpfr_init(r27372);
        mpfr_init(r27373);
        mpfr_init(r27374);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27367, a, MPFR_RNDN);
        mpfr_exp(r27368, r27367, MPFR_RNDN);
        mpfr_set_d(r27369, b, MPFR_RNDN);
        mpfr_exp(r27370, r27369, MPFR_RNDN);
        mpfr_add(r27371, r27368, r27370, MPFR_RNDN);
        mpfr_log(r27372, r27371, MPFR_RNDN);
        mpfr_sub(r27373, r27367, r27372, MPFR_RNDN);
        mpfr_exp(r27374, r27373, MPFR_RNDN);
        return mpfr_get_d(r27374, MPFR_RNDN);
}