#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 r27729 = a;
        float r27730 = log(r27729);
        float r27731 = b;
        float r27732 = log(r27731);
        float r27733 = r27730 + r27732;
        float r27734 = exp(r27733);
        return r27734;
}

double f_id(double a, double b) {
        double r27735 = a;
        double r27736 = log(r27735);
        double r27737 = b;
        double r27738 = log(r27737);
        double r27739 = r27736 + r27738;
        double r27740 = exp(r27739);
        return r27740;
}


double f_of(float a, float b) {
        float r27741 = a;
        float r27742 = b;
        float r27743 = r27741 * r27742;
        return r27743;
}

double f_od(double a, double b) {
        double r27744 = a;
        double r27745 = b;
        double r27746 = r27744 * r27745;
        return r27746;
}

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 r27747, r27748, r27749, r27750, r27751, r27752;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r27747);
        mpfr_init(r27748);
        mpfr_init(r27749);
        mpfr_init(r27750);
        mpfr_init(r27751);
        mpfr_init(r27752);
}

double f_im(double a, double b) {
        mpfr_set_d(r27747, a, MPFR_RNDN);
        mpfr_log(r27748, r27747, MPFR_RNDN);
        mpfr_set_d(r27749, b, MPFR_RNDN);
        mpfr_log(r27750, r27749, MPFR_RNDN);
        mpfr_add(r27751, r27748, r27750, MPFR_RNDN);
        mpfr_exp(r27752, r27751, MPFR_RNDN);
        return mpfr_get_d(r27752, MPFR_RNDN);
}

static mpfr_t r27753, r27754, r27755;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27753);
        mpfr_init(r27754);
        mpfr_init(r27755);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27753, a, MPFR_RNDN);
        mpfr_set_d(r27754, b, MPFR_RNDN);
        mpfr_mul(r27755, r27753, r27754, MPFR_RNDN);
        return mpfr_get_d(r27755, MPFR_RNDN);
}

static mpfr_t r27756, r27757, r27758;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27756);
        mpfr_init(r27757);
        mpfr_init(r27758);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27756, a, MPFR_RNDN);
        mpfr_set_d(r27757, b, MPFR_RNDN);
        mpfr_mul(r27758, r27756, r27757, MPFR_RNDN);
        return mpfr_get_d(r27758, MPFR_RNDN);
}