#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 r27079 = a;
        float r27080 = exp(r27079);
        float r27081 = b;
        float r27082 = exp(r27081);
        float r27083 = r27080 + r27082;
        float r27084 = r27080 / r27083;
        return r27084;
}

double f_id(double a, double b) {
        double r27085 = a;
        double r27086 = exp(r27085);
        double r27087 = b;
        double r27088 = exp(r27087);
        double r27089 = r27086 + r27088;
        double r27090 = r27086 / r27089;
        return r27090;
}


double f_of(float a, float b) {
        float r27091 = a;
        float r27092 = exp(r27091);
        float r27093 = b;
        float r27094 = exp(r27093);
        float r27095 = r27092 + r27094;
        float r27096 = log(r27095);
        float r27097 = r27091 - r27096;
        float r27098 = exp(r27097);
        return r27098;
}

double f_od(double a, double b) {
        double r27099 = a;
        double r27100 = exp(r27099);
        double r27101 = b;
        double r27102 = exp(r27101);
        double r27103 = r27100 + r27102;
        double r27104 = log(r27103);
        double r27105 = r27099 - r27104;
        double r27106 = exp(r27105);
        return r27106;
}

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 r27107, r27108, r27109, r27110, r27111, r27112;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r27107);
        mpfr_init(r27108);
        mpfr_init(r27109);
        mpfr_init(r27110);
        mpfr_init(r27111);
        mpfr_init(r27112);
}

double f_im(double a, double b) {
        mpfr_set_d(r27107, a, MPFR_RNDN);
        mpfr_exp(r27108, r27107, MPFR_RNDN);
        mpfr_set_d(r27109, b, MPFR_RNDN);
        mpfr_exp(r27110, r27109, MPFR_RNDN);
        mpfr_add(r27111, r27108, r27110, MPFR_RNDN);
        mpfr_div(r27112, r27108, r27111, MPFR_RNDN);
        return mpfr_get_d(r27112, MPFR_RNDN);
}

static mpfr_t r27113, r27114, r27115, r27116, r27117, r27118, r27119, r27120;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27113);
        mpfr_init(r27114);
        mpfr_init(r27115);
        mpfr_init(r27116);
        mpfr_init(r27117);
        mpfr_init(r27118);
        mpfr_init(r27119);
        mpfr_init(r27120);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27113, a, MPFR_RNDN);
        mpfr_exp(r27114, r27113, MPFR_RNDN);
        mpfr_set_d(r27115, b, MPFR_RNDN);
        mpfr_exp(r27116, r27115, MPFR_RNDN);
        mpfr_add(r27117, r27114, r27116, MPFR_RNDN);
        mpfr_log(r27118, r27117, MPFR_RNDN);
        mpfr_sub(r27119, r27113, r27118, MPFR_RNDN);
        mpfr_exp(r27120, r27119, MPFR_RNDN);
        return mpfr_get_d(r27120, MPFR_RNDN);
}

static mpfr_t r27121, r27122, r27123, r27124, r27125, r27126, r27127, r27128;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r27121);
        mpfr_init(r27122);
        mpfr_init(r27123);
        mpfr_init(r27124);
        mpfr_init(r27125);
        mpfr_init(r27126);
        mpfr_init(r27127);
        mpfr_init(r27128);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27121, a, MPFR_RNDN);
        mpfr_exp(r27122, r27121, MPFR_RNDN);
        mpfr_set_d(r27123, b, MPFR_RNDN);
        mpfr_exp(r27124, r27123, MPFR_RNDN);
        mpfr_add(r27125, r27122, r27124, MPFR_RNDN);
        mpfr_log(r27126, r27125, MPFR_RNDN);
        mpfr_sub(r27127, r27121, r27126, MPFR_RNDN);
        mpfr_exp(r27128, r27127, MPFR_RNDN);
        return mpfr_get_d(r27128, MPFR_RNDN);
}