#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 r27913 = a;
        float r27914 = log(r27913);
        float r27915 = b;
        float r27916 = log(r27915);
        float r27917 = r27914 + r27916;
        float r27918 = exp(r27917);
        return r27918;
}

double f_id(double a, double b) {
        double r27919 = a;
        double r27920 = log(r27919);
        double r27921 = b;
        double r27922 = log(r27921);
        double r27923 = r27920 + r27922;
        double r27924 = exp(r27923);
        return r27924;
}


double f_of(float a, float b) {
        float r27925 = a;
        float r27926 = b;
        float r27927 = r27925 * r27926;
        return r27927;
}

double f_od(double a, double b) {
        double r27928 = a;
        double r27929 = b;
        double r27930 = r27928 * r27929;
        return r27930;
}

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 r27931, r27932, r27933, r27934, r27935, r27936;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r27931);
        mpfr_init(r27932);
        mpfr_init(r27933);
        mpfr_init(r27934);
        mpfr_init(r27935);
        mpfr_init(r27936);
}

double f_im(double a, double b) {
        mpfr_set_d(r27931, a, MPFR_RNDN);
        mpfr_log(r27932, r27931, MPFR_RNDN);
        mpfr_set_d(r27933, b, MPFR_RNDN);
        mpfr_log(r27934, r27933, MPFR_RNDN);
        mpfr_add(r27935, r27932, r27934, MPFR_RNDN);
        mpfr_exp(r27936, r27935, MPFR_RNDN);
        return mpfr_get_d(r27936, MPFR_RNDN);
}

static mpfr_t r27937, r27938, r27939;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27937);
        mpfr_init(r27938);
        mpfr_init(r27939);
}

double f_fm(double a, double b) {
        mpfr_set_d(r27937, a, MPFR_RNDN);
        mpfr_set_d(r27938, b, MPFR_RNDN);
        mpfr_mul(r27939, r27937, r27938, MPFR_RNDN);
        return mpfr_get_d(r27939, MPFR_RNDN);
}

static mpfr_t r27940, r27941, r27942;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27940);
        mpfr_init(r27941);
        mpfr_init(r27942);
}

double f_dm(double a, double b) {
        mpfr_set_d(r27940, a, MPFR_RNDN);
        mpfr_set_d(r27941, b, MPFR_RNDN);
        mpfr_mul(r27942, r27940, r27941, MPFR_RNDN);
        return mpfr_get_d(r27942, MPFR_RNDN);
}