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

char *name = "Exp of sum of logs";

double f_if(float a, float b) {
        float r27736 = a;
        float r27737 = log(r27736);
        float r27738 = b;
        float r27739 = log(r27738);
        float r27740 = r27737 + r27739;
        float r27741 = exp(r27740);
        return r27741;
}

double f_id(double a, double b) {
        double r27742 = a;
        double r27743 = log(r27742);
        double r27744 = b;
        double r27745 = log(r27744);
        double r27746 = r27743 + r27745;
        double r27747 = exp(r27746);
        return r27747;
}


double f_of(float a, float b) {
        float r27748 = a;
        float r27749 = b;
        float r27750 = r27748 * r27749;
        return r27750;
}

double f_od(double a, double b) {
        double r27751 = a;
        double r27752 = b;
        double r27753 = r27751 * r27752;
        return r27753;
}

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 r27754, r27755, r27756, r27757, r27758, r27759;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r27754);
        mpfr_init(r27755);
        mpfr_init(r27756);
        mpfr_init(r27757);
        mpfr_init(r27758);
        mpfr_init(r27759);
}

double f_im(double a, double b) {
        mpfr_set_d(r27754, a, MPFR_RNDN);
        mpfr_log(r27755, r27754, MPFR_RNDN);
        mpfr_set_d(r27756, b, MPFR_RNDN);
        mpfr_log(r27757, r27756, MPFR_RNDN);
        mpfr_add(r27758, r27755, r27757, MPFR_RNDN);
        mpfr_exp(r27759, r27758, MPFR_RNDN);
        return mpfr_get_d(r27759, MPFR_RNDN);
}

static mpfr_t r27760, r27761, r27762;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27760);
        mpfr_init(r27761);
        mpfr_init(r27762);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27760, a, MPFR_RNDN);
        mpfr_set_d(r27761, b, MPFR_RNDN);
        mpfr_mul(r27762, r27760, r27761, MPFR_RNDN);
        return mpfr_get_d(r27762, MPFR_RNDN);
}

static mpfr_t r27763, r27764, r27765;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27763);
        mpfr_init(r27764);
        mpfr_init(r27765);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27763, a, MPFR_RNDN);
        mpfr_set_d(r27764, b, MPFR_RNDN);
        mpfr_mul(r27765, r27763, r27764, MPFR_RNDN);
        return mpfr_get_d(r27765, MPFR_RNDN);
}