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

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

double f_if(float m, float v) {
        float r8050 = m;
        float r8051 = 1;
        float r8052 = r8051 - r8050;
        float r8053 = r8050 * r8052;
        float r8054 = v;
        float r8055 = r8053 / r8054;
        float r8056 = r8055 - r8051;
        float r8057 = r8056 * r8050;
        return r8057;
}

double f_id(double m, double v) {
        double r8058 = m;
        double r8059 = 1;
        double r8060 = r8059 - r8058;
        double r8061 = r8058 * r8060;
        double r8062 = v;
        double r8063 = r8061 / r8062;
        double r8064 = r8063 - r8059;
        double r8065 = r8064 * r8058;
        return r8065;
}


double f_of(float m, float v) {
        float r8066 = m;
        float r8067 = r8066 * r8066;
        float r8068 = r8066 - r8067;
        float r8069 = v;
        float r8070 = r8066 / r8069;
        float r8071 = r8068 * r8070;
        float r8072 = r8071 - r8066;
        return r8072;
}

double f_od(double m, double v) {
        double r8073 = m;
        double r8074 = r8073 * r8073;
        double r8075 = r8073 - r8074;
        double r8076 = v;
        double r8077 = r8073 / r8076;
        double r8078 = r8075 * r8077;
        double r8079 = r8078 - r8073;
        return r8079;
}

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 r8080, r8081, r8082, r8083, r8084, r8085, r8086, r8087;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r8080);
        mpfr_init_set_str(r8081, "1", 10, MPFR_RNDN);
        mpfr_init(r8082);
        mpfr_init(r8083);
        mpfr_init(r8084);
        mpfr_init(r8085);
        mpfr_init(r8086);
        mpfr_init(r8087);
}

double f_im(double m, double v) {
        mpfr_set_d(r8080, m, MPFR_RNDN);
        ;
        mpfr_sub(r8082, r8081, r8080, MPFR_RNDN);
        mpfr_mul(r8083, r8080, r8082, MPFR_RNDN);
        mpfr_set_d(r8084, v, MPFR_RNDN);
        mpfr_div(r8085, r8083, r8084, MPFR_RNDN);
        mpfr_sub(r8086, r8085, r8081, MPFR_RNDN);
        mpfr_mul(r8087, r8086, r8080, MPFR_RNDN);
        return mpfr_get_d(r8087, MPFR_RNDN);
}

static mpfr_t r8088, r8089, r8090, r8091, r8092, r8093, r8094;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r8088);
        mpfr_init(r8089);
        mpfr_init(r8090);
        mpfr_init(r8091);
        mpfr_init(r8092);
        mpfr_init(r8093);
        mpfr_init(r8094);
}

double f_fm(double m, double v) {
        mpfr_set_d(r8088, m, MPFR_RNDN);
        mpfr_mul(r8089, r8088, r8088, MPFR_RNDN);
        mpfr_sub(r8090, r8088, r8089, MPFR_RNDN);
        mpfr_set_d(r8091, v, MPFR_RNDN);
        mpfr_div(r8092, r8088, r8091, MPFR_RNDN);
        mpfr_mul(r8093, r8090, r8092, MPFR_RNDN);
        mpfr_sub(r8094, r8093, r8088, MPFR_RNDN);
        return mpfr_get_d(r8094, MPFR_RNDN);
}

static mpfr_t r8095, r8096, r8097, r8098, r8099, r8100, r8101;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r8095);
        mpfr_init(r8096);
        mpfr_init(r8097);
        mpfr_init(r8098);
        mpfr_init(r8099);
        mpfr_init(r8100);
        mpfr_init(r8101);
}

double f_dm(double m, double v) {
        mpfr_set_d(r8095, m, MPFR_RNDN);
        mpfr_mul(r8096, r8095, r8095, MPFR_RNDN);
        mpfr_sub(r8097, r8095, r8096, MPFR_RNDN);
        mpfr_set_d(r8098, v, MPFR_RNDN);
        mpfr_div(r8099, r8095, r8098, MPFR_RNDN);
        mpfr_mul(r8100, r8097, r8099, MPFR_RNDN);
        mpfr_sub(r8101, r8100, r8095, MPFR_RNDN);
        return mpfr_get_d(r8101, MPFR_RNDN);
}