#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 r25454 = a;
        float r25455 = log(r25454);
        float r25456 = b;
        float r25457 = log(r25456);
        float r25458 = r25455 + r25457;
        float r25459 = exp(r25458);
        return r25459;
}

double f_id(double a, double b) {
        double r25460 = a;
        double r25461 = log(r25460);
        double r25462 = b;
        double r25463 = log(r25462);
        double r25464 = r25461 + r25463;
        double r25465 = exp(r25464);
        return r25465;
}


double f_of(float a, float b) {
        float r25466 = a;
        float r25467 = b;
        float r25468 = r25466 * r25467;
        return r25468;
}

double f_od(double a, double b) {
        double r25469 = a;
        double r25470 = b;
        double r25471 = r25469 * r25470;
        return r25471;
}

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 r25472, r25473, r25474, r25475, r25476, r25477;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r25472);
        mpfr_init(r25473);
        mpfr_init(r25474);
        mpfr_init(r25475);
        mpfr_init(r25476);
        mpfr_init(r25477);
}

double f_im(double a, double b) {
        mpfr_set_d(r25472, a, MPFR_RNDN);
        mpfr_log(r25473, r25472, MPFR_RNDN);
        mpfr_set_d(r25474, b, MPFR_RNDN);
        mpfr_log(r25475, r25474, MPFR_RNDN);
        mpfr_add(r25476, r25473, r25475, MPFR_RNDN);
        mpfr_exp(r25477, r25476, MPFR_RNDN);
        return mpfr_get_d(r25477, MPFR_RNDN);
}

static mpfr_t r25478, r25479, r25480;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r25478);
        mpfr_init(r25479);
        mpfr_init(r25480);
}

double f_fm(double a, double b) {
        mpfr_set_d(r25478, a, MPFR_RNDN);
        mpfr_set_d(r25479, b, MPFR_RNDN);
        mpfr_mul(r25480, r25478, r25479, MPFR_RNDN);
        return mpfr_get_d(r25480, MPFR_RNDN);
}

static mpfr_t r25481, r25482, r25483;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r25481);
        mpfr_init(r25482);
        mpfr_init(r25483);
}

double f_dm(double a, double b) {
        mpfr_set_d(r25481, a, MPFR_RNDN);
        mpfr_set_d(r25482, b, MPFR_RNDN);
        mpfr_mul(r25483, r25481, r25482, MPFR_RNDN);
        return mpfr_get_d(r25483, MPFR_RNDN);
}