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

char *name = "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2";

double f_if(float x, float y, float z, float t, float a, float b) {
        float r9075 = x;
        float r9076 = y;
        float r9077 = z;
        float r9078 = log(r9077);
        float r9079 = r9076 * r9078;
        float r9080 = t;
        float r9081 = 1.0f;
        float r9082 = r9080 - r9081;
        float r9083 = a;
        float r9084 = log(r9083);
        float r9085 = r9082 * r9084;
        float r9086 = r9079 + r9085;
        float r9087 = b;
        float r9088 = r9086 - r9087;
        float r9089 = exp(r9088);
        float r9090 = r9075 * r9089;
        float r9091 = r9090 / r9076;
        return r9091;
}

double f_id(double x, double y, double z, double t, double a, double b) {
        double r9092 = x;
        double r9093 = y;
        double r9094 = z;
        double r9095 = log(r9094);
        double r9096 = r9093 * r9095;
        double r9097 = t;
        double r9098 = 1.0;
        double r9099 = r9097 - r9098;
        double r9100 = a;
        double r9101 = log(r9100);
        double r9102 = r9099 * r9101;
        double r9103 = r9096 + r9102;
        double r9104 = b;
        double r9105 = r9103 - r9104;
        double r9106 = exp(r9105);
        double r9107 = r9092 * r9106;
        double r9108 = r9107 / r9093;
        return r9108;
}


double f_of(float x, float y, float z, float t, float a, float b) {
        float r9109 = x;
        float r9110 = a;
        float r9111 = log(r9110);
        float r9112 = t;
        float r9113 = 1.0f;
        float r9114 = r9112 - r9113;
        float r9115 = r9111 * r9114;
        float r9116 = z;
        float r9117 = log(r9116);
        float r9118 = y;
        float r9119 = r9117 * r9118;
        float r9120 = r9115 + r9119;
        float r9121 = b;
        float r9122 = r9120 - r9121;
        float r9123 = exp(r9122);
        float r9124 = cbrt(r9123);
        float r9125 = r9124 * r9124;
        float r9126 = r9124 * r9125;
        float r9127 = r9109 * r9126;
        float r9128 = r9127 / r9118;
        return r9128;
}

double f_od(double x, double y, double z, double t, double a, double b) {
        double r9129 = x;
        double r9130 = a;
        double r9131 = log(r9130);
        double r9132 = t;
        double r9133 = 1.0;
        double r9134 = r9132 - r9133;
        double r9135 = r9131 * r9134;
        double r9136 = z;
        double r9137 = log(r9136);
        double r9138 = y;
        double r9139 = r9137 * r9138;
        double r9140 = r9135 + r9139;
        double r9141 = b;
        double r9142 = r9140 - r9141;
        double r9143 = exp(r9142);
        double r9144 = cbrt(r9143);
        double r9145 = r9144 * r9144;
        double r9146 = r9144 * r9145;
        double r9147 = r9129 * r9146;
        double r9148 = r9147 / r9138;
        return r9148;
}

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 r9149, r9150, r9151, r9152, r9153, r9154, r9155, r9156, r9157, r9158, r9159, r9160, r9161, r9162, r9163, r9164, r9165;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9149);
        mpfr_init(r9150);
        mpfr_init(r9151);
        mpfr_init(r9152);
        mpfr_init(r9153);
        mpfr_init(r9154);
        mpfr_init_set_str(r9155, "1.0", 10, MPFR_RNDN);
        mpfr_init(r9156);
        mpfr_init(r9157);
        mpfr_init(r9158);
        mpfr_init(r9159);
        mpfr_init(r9160);
        mpfr_init(r9161);
        mpfr_init(r9162);
        mpfr_init(r9163);
        mpfr_init(r9164);
        mpfr_init(r9165);
}

double f_im(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r9149, x, MPFR_RNDN);
        mpfr_set_d(r9150, y, MPFR_RNDN);
        mpfr_set_d(r9151, z, MPFR_RNDN);
        mpfr_log(r9152, r9151, MPFR_RNDN);
        mpfr_mul(r9153, r9150, r9152, MPFR_RNDN);
        mpfr_set_d(r9154, t, MPFR_RNDN);
        ;
        mpfr_sub(r9156, r9154, r9155, MPFR_RNDN);
        mpfr_set_d(r9157, a, MPFR_RNDN);
        mpfr_log(r9158, r9157, MPFR_RNDN);
        mpfr_mul(r9159, r9156, r9158, MPFR_RNDN);
        mpfr_add(r9160, r9153, r9159, MPFR_RNDN);
        mpfr_set_d(r9161, b, MPFR_RNDN);
        mpfr_sub(r9162, r9160, r9161, MPFR_RNDN);
        mpfr_exp(r9163, r9162, MPFR_RNDN);
        mpfr_mul(r9164, r9149, r9163, MPFR_RNDN);
        mpfr_div(r9165, r9164, r9150, MPFR_RNDN);
        return mpfr_get_d(r9165, MPFR_RNDN);
}

static mpfr_t r9166, r9167, r9168, r9169, r9170, r9171, r9172, r9173, r9174, r9175, r9176, r9177, r9178, r9179, r9180, r9181, r9182, r9183, r9184, r9185;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9166);
        mpfr_init(r9167);
        mpfr_init(r9168);
        mpfr_init(r9169);
        mpfr_init_set_str(r9170, "1.0", 10, MPFR_RNDN);
        mpfr_init(r9171);
        mpfr_init(r9172);
        mpfr_init(r9173);
        mpfr_init(r9174);
        mpfr_init(r9175);
        mpfr_init(r9176);
        mpfr_init(r9177);
        mpfr_init(r9178);
        mpfr_init(r9179);
        mpfr_init(r9180);
        mpfr_init(r9181);
        mpfr_init(r9182);
        mpfr_init(r9183);
        mpfr_init(r9184);
        mpfr_init(r9185);
}

double f_fm(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r9166, x, MPFR_RNDN);
        mpfr_set_d(r9167, a, MPFR_RNDN);
        mpfr_log(r9168, r9167, MPFR_RNDN);
        mpfr_set_d(r9169, t, MPFR_RNDN);
        ;
        mpfr_sub(r9171, r9169, r9170, MPFR_RNDN);
        mpfr_mul(r9172, r9168, r9171, MPFR_RNDN);
        mpfr_set_d(r9173, z, MPFR_RNDN);
        mpfr_log(r9174, r9173, MPFR_RNDN);
        mpfr_set_d(r9175, y, MPFR_RNDN);
        mpfr_mul(r9176, r9174, r9175, MPFR_RNDN);
        mpfr_add(r9177, r9172, r9176, MPFR_RNDN);
        mpfr_set_d(r9178, b, MPFR_RNDN);
        mpfr_sub(r9179, r9177, r9178, MPFR_RNDN);
        mpfr_exp(r9180, r9179, MPFR_RNDN);
        mpfr_cbrt(r9181, r9180, MPFR_RNDN);
        mpfr_mul(r9182, r9181, r9181, MPFR_RNDN);
        mpfr_mul(r9183, r9181, r9182, MPFR_RNDN);
        mpfr_mul(r9184, r9166, r9183, MPFR_RNDN);
        mpfr_div(r9185, r9184, r9175, MPFR_RNDN);
        return mpfr_get_d(r9185, MPFR_RNDN);
}

static mpfr_t r9186, r9187, r9188, r9189, r9190, r9191, r9192, r9193, r9194, r9195, r9196, r9197, r9198, r9199, r9200, r9201, r9202, r9203, r9204, r9205;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9186);
        mpfr_init(r9187);
        mpfr_init(r9188);
        mpfr_init(r9189);
        mpfr_init_set_str(r9190, "1.0", 10, MPFR_RNDN);
        mpfr_init(r9191);
        mpfr_init(r9192);
        mpfr_init(r9193);
        mpfr_init(r9194);
        mpfr_init(r9195);
        mpfr_init(r9196);
        mpfr_init(r9197);
        mpfr_init(r9198);
        mpfr_init(r9199);
        mpfr_init(r9200);
        mpfr_init(r9201);
        mpfr_init(r9202);
        mpfr_init(r9203);
        mpfr_init(r9204);
        mpfr_init(r9205);
}

double f_dm(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r9186, x, MPFR_RNDN);
        mpfr_set_d(r9187, a, MPFR_RNDN);
        mpfr_log(r9188, r9187, MPFR_RNDN);
        mpfr_set_d(r9189, t, MPFR_RNDN);
        ;
        mpfr_sub(r9191, r9189, r9190, MPFR_RNDN);
        mpfr_mul(r9192, r9188, r9191, MPFR_RNDN);
        mpfr_set_d(r9193, z, MPFR_RNDN);
        mpfr_log(r9194, r9193, MPFR_RNDN);
        mpfr_set_d(r9195, y, MPFR_RNDN);
        mpfr_mul(r9196, r9194, r9195, MPFR_RNDN);
        mpfr_add(r9197, r9192, r9196, MPFR_RNDN);
        mpfr_set_d(r9198, b, MPFR_RNDN);
        mpfr_sub(r9199, r9197, r9198, MPFR_RNDN);
        mpfr_exp(r9200, r9199, MPFR_RNDN);
        mpfr_cbrt(r9201, r9200, MPFR_RNDN);
        mpfr_mul(r9202, r9201, r9201, MPFR_RNDN);
        mpfr_mul(r9203, r9201, r9202, MPFR_RNDN);
        mpfr_mul(r9204, r9186, r9203, MPFR_RNDN);
        mpfr_div(r9205, r9204, r9195, MPFR_RNDN);
        return mpfr_get_d(r9205, MPFR_RNDN);
}