#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 r9034 = x;
        float r9035 = y;
        float r9036 = z;
        float r9037 = log(r9036);
        float r9038 = r9035 * r9037;
        float r9039 = t;
        float r9040 = 1.0f;
        float r9041 = r9039 - r9040;
        float r9042 = a;
        float r9043 = log(r9042);
        float r9044 = r9041 * r9043;
        float r9045 = r9038 + r9044;
        float r9046 = b;
        float r9047 = r9045 - r9046;
        float r9048 = exp(r9047);
        float r9049 = r9034 * r9048;
        float r9050 = r9049 / r9035;
        return r9050;
}

double f_id(double x, double y, double z, double t, double a, double b) {
        double r9051 = x;
        double r9052 = y;
        double r9053 = z;
        double r9054 = log(r9053);
        double r9055 = r9052 * r9054;
        double r9056 = t;
        double r9057 = 1.0;
        double r9058 = r9056 - r9057;
        double r9059 = a;
        double r9060 = log(r9059);
        double r9061 = r9058 * r9060;
        double r9062 = r9055 + r9061;
        double r9063 = b;
        double r9064 = r9062 - r9063;
        double r9065 = exp(r9064);
        double r9066 = r9051 * r9065;
        double r9067 = r9066 / r9052;
        return r9067;
}


double f_of(float x, float y, float z, float t, float a, float b) {
        float r9068 = x;
        float r9069 = a;
        float r9070 = log(r9069);
        float r9071 = t;
        float r9072 = 1.0f;
        float r9073 = r9071 - r9072;
        float r9074 = r9070 * r9073;
        float r9075 = z;
        float r9076 = log(r9075);
        float r9077 = y;
        float r9078 = r9076 * r9077;
        float r9079 = r9074 + r9078;
        float r9080 = b;
        float r9081 = r9079 - r9080;
        float r9082 = exp(r9081);
        float r9083 = cbrt(r9082);
        float r9084 = r9083 * r9083;
        float r9085 = r9083 * r9084;
        float r9086 = r9068 * r9085;
        float r9087 = r9086 / r9077;
        return r9087;
}

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

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 r9108, r9109, r9110, r9111, r9112, r9113, r9114, r9115, r9116, r9117, r9118, r9119, r9120, r9121, r9122, r9123, r9124;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9108);
        mpfr_init(r9109);
        mpfr_init(r9110);
        mpfr_init(r9111);
        mpfr_init(r9112);
        mpfr_init(r9113);
        mpfr_init_set_str(r9114, "1.0", 10, MPFR_RNDN);
        mpfr_init(r9115);
        mpfr_init(r9116);
        mpfr_init(r9117);
        mpfr_init(r9118);
        mpfr_init(r9119);
        mpfr_init(r9120);
        mpfr_init(r9121);
        mpfr_init(r9122);
        mpfr_init(r9123);
        mpfr_init(r9124);
}

double f_im(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r9108, x, MPFR_RNDN);
        mpfr_set_d(r9109, y, MPFR_RNDN);
        mpfr_set_d(r9110, z, MPFR_RNDN);
        mpfr_log(r9111, r9110, MPFR_RNDN);
        mpfr_mul(r9112, r9109, r9111, MPFR_RNDN);
        mpfr_set_d(r9113, t, MPFR_RNDN);
        ;
        mpfr_sub(r9115, r9113, r9114, MPFR_RNDN);
        mpfr_set_d(r9116, a, MPFR_RNDN);
        mpfr_log(r9117, r9116, MPFR_RNDN);
        mpfr_mul(r9118, r9115, r9117, MPFR_RNDN);
        mpfr_add(r9119, r9112, r9118, MPFR_RNDN);
        mpfr_set_d(r9120, b, MPFR_RNDN);
        mpfr_sub(r9121, r9119, r9120, MPFR_RNDN);
        mpfr_exp(r9122, r9121, MPFR_RNDN);
        mpfr_mul(r9123, r9108, r9122, MPFR_RNDN);
        mpfr_div(r9124, r9123, r9109, MPFR_RNDN);
        return mpfr_get_d(r9124, MPFR_RNDN);
}

static mpfr_t r9125, r9126, r9127, r9128, r9129, r9130, r9131, r9132, r9133, r9134, r9135, r9136, r9137, r9138, r9139, r9140, r9141, r9142, r9143, r9144;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9125);
        mpfr_init(r9126);
        mpfr_init(r9127);
        mpfr_init(r9128);
        mpfr_init_set_str(r9129, "1.0", 10, MPFR_RNDN);
        mpfr_init(r9130);
        mpfr_init(r9131);
        mpfr_init(r9132);
        mpfr_init(r9133);
        mpfr_init(r9134);
        mpfr_init(r9135);
        mpfr_init(r9136);
        mpfr_init(r9137);
        mpfr_init(r9138);
        mpfr_init(r9139);
        mpfr_init(r9140);
        mpfr_init(r9141);
        mpfr_init(r9142);
        mpfr_init(r9143);
        mpfr_init(r9144);
}

double f_fm(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r9125, x, MPFR_RNDN);
        mpfr_set_d(r9126, a, MPFR_RNDN);
        mpfr_log(r9127, r9126, MPFR_RNDN);
        mpfr_set_d(r9128, t, MPFR_RNDN);
        ;
        mpfr_sub(r9130, r9128, r9129, MPFR_RNDN);
        mpfr_mul(r9131, r9127, r9130, MPFR_RNDN);
        mpfr_set_d(r9132, z, MPFR_RNDN);
        mpfr_log(r9133, r9132, MPFR_RNDN);
        mpfr_set_d(r9134, y, MPFR_RNDN);
        mpfr_mul(r9135, r9133, r9134, MPFR_RNDN);
        mpfr_add(r9136, r9131, r9135, MPFR_RNDN);
        mpfr_set_d(r9137, b, MPFR_RNDN);
        mpfr_sub(r9138, r9136, r9137, MPFR_RNDN);
        mpfr_exp(r9139, r9138, MPFR_RNDN);
        mpfr_cbrt(r9140, r9139, MPFR_RNDN);
        mpfr_mul(r9141, r9140, r9140, MPFR_RNDN);
        mpfr_mul(r9142, r9140, r9141, MPFR_RNDN);
        mpfr_mul(r9143, r9125, r9142, MPFR_RNDN);
        mpfr_div(r9144, r9143, r9134, MPFR_RNDN);
        return mpfr_get_d(r9144, MPFR_RNDN);
}

static mpfr_t r9145, r9146, r9147, r9148, r9149, r9150, r9151, r9152, r9153, r9154, r9155, r9156, r9157, r9158, r9159, r9160, r9161, r9162, r9163, r9164;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9145);
        mpfr_init(r9146);
        mpfr_init(r9147);
        mpfr_init(r9148);
        mpfr_init_set_str(r9149, "1.0", 10, MPFR_RNDN);
        mpfr_init(r9150);
        mpfr_init(r9151);
        mpfr_init(r9152);
        mpfr_init(r9153);
        mpfr_init(r9154);
        mpfr_init(r9155);
        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);
}

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