#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 r27772 = a;
        float r27773 = log(r27772);
        float r27774 = b;
        float r27775 = log(r27774);
        float r27776 = r27773 + r27775;
        float r27777 = exp(r27776);
        return r27777;
}

double f_id(double a, double b) {
        double r27778 = a;
        double r27779 = log(r27778);
        double r27780 = b;
        double r27781 = log(r27780);
        double r27782 = r27779 + r27781;
        double r27783 = exp(r27782);
        return r27783;
}


double f_of(float a, float b) {
        float r27784 = a;
        float r27785 = b;
        float r27786 = r27784 * r27785;
        return r27786;
}

double f_od(double a, double b) {
        double r27787 = a;
        double r27788 = b;
        double r27789 = r27787 * r27788;
        return r27789;
}

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 r27790, r27791, r27792, r27793, r27794, r27795;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r27790);
        mpfr_init(r27791);
        mpfr_init(r27792);
        mpfr_init(r27793);
        mpfr_init(r27794);
        mpfr_init(r27795);
}

double f_im(double a, double b) {
        mpfr_set_d(r27790, a, MPFR_RNDN);
        mpfr_log(r27791, r27790, MPFR_RNDN);
        mpfr_set_d(r27792, b, MPFR_RNDN);
        mpfr_log(r27793, r27792, MPFR_RNDN);
        mpfr_add(r27794, r27791, r27793, MPFR_RNDN);
        mpfr_exp(r27795, r27794, MPFR_RNDN);
        return mpfr_get_d(r27795, MPFR_RNDN);
}

static mpfr_t r27796, r27797, r27798;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27796);
        mpfr_init(r27797);
        mpfr_init(r27798);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27796, a, MPFR_RNDN);
        mpfr_set_d(r27797, b, MPFR_RNDN);
        mpfr_mul(r27798, r27796, r27797, MPFR_RNDN);
        return mpfr_get_d(r27798, MPFR_RNDN);
}

static mpfr_t r27799, r27800, r27801;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27799);
        mpfr_init(r27800);
        mpfr_init(r27801);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27799, a, MPFR_RNDN);
        mpfr_set_d(r27800, b, MPFR_RNDN);
        mpfr_mul(r27801, r27799, r27800, MPFR_RNDN);
        return mpfr_get_d(r27801, MPFR_RNDN);
}