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

char *name = "b parameter of renormalized beta distribution";

double f_if(float m, float v) {
        float r7786 = m;
        float r7787 = 1;
        float r7788 = r7787 - r7786;
        float r7789 = r7786 * r7788;
        float r7790 = v;
        float r7791 = r7789 / r7790;
        float r7792 = r7791 - r7787;
        float r7793 = r7792 * r7788;
        return r7793;
}

double f_id(double m, double v) {
        double r7794 = m;
        double r7795 = 1;
        double r7796 = r7795 - r7794;
        double r7797 = r7794 * r7796;
        double r7798 = v;
        double r7799 = r7797 / r7798;
        double r7800 = r7799 - r7795;
        double r7801 = r7800 * r7796;
        return r7801;
}


double f_of(float m, float v) {
        float r7802 = m;
        float r7803 = 1;
        float r7804 = r7803 - r7802;
        float r7805 = r7802 * r7804;
        float r7806 = v;
        float r7807 = r7805 / r7806;
        float r7808 = r7807 - r7803;
        float r7809 = r7808 * r7804;
        return r7809;
}

double f_od(double m, double v) {
        double r7810 = m;
        double r7811 = 1;
        double r7812 = r7811 - r7810;
        double r7813 = r7810 * r7812;
        double r7814 = v;
        double r7815 = r7813 / r7814;
        double r7816 = r7815 - r7811;
        double r7817 = r7816 * r7812;
        return r7817;
}

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 r7818, r7819, r7820, r7821, r7822, r7823, r7824, r7825;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r7818);
        mpfr_init_set_str(r7819, "1", 10, MPFR_RNDN);
        mpfr_init(r7820);
        mpfr_init(r7821);
        mpfr_init(r7822);
        mpfr_init(r7823);
        mpfr_init(r7824);
        mpfr_init(r7825);
}

double f_im(double m, double v) {
        mpfr_set_d(r7818, m, MPFR_RNDN);
        ;
        mpfr_sub(r7820, r7819, r7818, MPFR_RNDN);
        mpfr_mul(r7821, r7818, r7820, MPFR_RNDN);
        mpfr_set_d(r7822, v, MPFR_RNDN);
        mpfr_div(r7823, r7821, r7822, MPFR_RNDN);
        mpfr_sub(r7824, r7823, r7819, MPFR_RNDN);
        mpfr_mul(r7825, r7824, r7820, MPFR_RNDN);
        return mpfr_get_d(r7825, MPFR_RNDN);
}

static mpfr_t r7826, r7827, r7828, r7829, r7830, r7831, r7832, r7833;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r7826);
        mpfr_init_set_str(r7827, "1", 10, MPFR_RNDN);
        mpfr_init(r7828);
        mpfr_init(r7829);
        mpfr_init(r7830);
        mpfr_init(r7831);
        mpfr_init(r7832);
        mpfr_init(r7833);
}

double f_fm(double m, double v) {
        mpfr_set_d(r7826, m, MPFR_RNDN);
        ;
        mpfr_sub(r7828, r7827, r7826, MPFR_RNDN);
        mpfr_mul(r7829, r7826, r7828, MPFR_RNDN);
        mpfr_set_d(r7830, v, MPFR_RNDN);
        mpfr_div(r7831, r7829, r7830, MPFR_RNDN);
        mpfr_sub(r7832, r7831, r7827, MPFR_RNDN);
        mpfr_mul(r7833, r7832, r7828, MPFR_RNDN);
        return mpfr_get_d(r7833, MPFR_RNDN);
}

static mpfr_t r7834, r7835, r7836, r7837, r7838, r7839, r7840, r7841;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r7834);
        mpfr_init_set_str(r7835, "1", 10, MPFR_RNDN);
        mpfr_init(r7836);
        mpfr_init(r7837);
        mpfr_init(r7838);
        mpfr_init(r7839);
        mpfr_init(r7840);
        mpfr_init(r7841);
}

double f_dm(double m, double v) {
        mpfr_set_d(r7834, m, MPFR_RNDN);
        ;
        mpfr_sub(r7836, r7835, r7834, MPFR_RNDN);
        mpfr_mul(r7837, r7834, r7836, MPFR_RNDN);
        mpfr_set_d(r7838, v, MPFR_RNDN);
        mpfr_div(r7839, r7837, r7838, MPFR_RNDN);
        mpfr_sub(r7840, r7839, r7835, MPFR_RNDN);
        mpfr_mul(r7841, r7840, r7836, MPFR_RNDN);
        return mpfr_get_d(r7841, MPFR_RNDN);
}