#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 r27654 = x;
        float r27655 = y;
        float r27656 = z;
        float r27657 = log(r27656);
        float r27658 = r27655 * r27657;
        float r27659 = t;
        float r27660 = 1.0;
        float r27661 = r27659 - r27660;
        float r27662 = a;
        float r27663 = log(r27662);
        float r27664 = r27661 * r27663;
        float r27665 = r27658 + r27664;
        float r27666 = b;
        float r27667 = r27665 - r27666;
        float r27668 = exp(r27667);
        float r27669 = r27654 * r27668;
        float r27670 = r27669 / r27655;
        return r27670;
}

double f_id(double x, double y, double z, double t, double a, double b) {
        double r27671 = x;
        double r27672 = y;
        double r27673 = z;
        double r27674 = log(r27673);
        double r27675 = r27672 * r27674;
        double r27676 = t;
        double r27677 = 1.0;
        double r27678 = r27676 - r27677;
        double r27679 = a;
        double r27680 = log(r27679);
        double r27681 = r27678 * r27680;
        double r27682 = r27675 + r27681;
        double r27683 = b;
        double r27684 = r27682 - r27683;
        double r27685 = exp(r27684);
        double r27686 = r27671 * r27685;
        double r27687 = r27686 / r27672;
        return r27687;
}


double f_of(float x, float y, float z, float t, float a, float b) {
        float r27688 = t;
        float r27689 = -1486774482200.535;
        bool r27690 = r27688 <= r27689;
        float r27691 = x;
        float r27692 = a;
        float r27693 = 1.0;
        float r27694 = r27688 - r27693;
        float r27695 = pow(r27692, r27694);
        float r27696 = b;
        float r27697 = 1;
        float r27698 = r27696 + r27697;
        float r27699 = y;
        float r27700 = z;
        float r27701 = log(r27700);
        float r27702 = r27699 * r27701;
        float r27703 = r27698 - r27702;
        float r27704 = r27695 / r27703;
        float r27705 = r27691 * r27704;
        float r27706 = r27705 / r27699;
        float r27707 = 1.324949139524979e-102;
        bool r27708 = r27688 <= r27707;
        float r27709 = r27691 / r27699;
        float r27710 = r27701 * r27699;
        float r27711 = r27696 - r27710;
        float r27712 = exp(r27711);
        float r27713 = r27695 / r27712;
        float r27714 = r27709 * r27713;
        float r27715 = r27708 ? r27714 : r27706;
        float r27716 = r27690 ? r27706 : r27715;
        return r27716;
}

double f_od(double x, double y, double z, double t, double a, double b) {
        double r27717 = t;
        double r27718 = -1486774482200.535;
        bool r27719 = r27717 <= r27718;
        double r27720 = x;
        double r27721 = a;
        double r27722 = 1.0;
        double r27723 = r27717 - r27722;
        double r27724 = pow(r27721, r27723);
        double r27725 = b;
        double r27726 = 1;
        double r27727 = r27725 + r27726;
        double r27728 = y;
        double r27729 = z;
        double r27730 = log(r27729);
        double r27731 = r27728 * r27730;
        double r27732 = r27727 - r27731;
        double r27733 = r27724 / r27732;
        double r27734 = r27720 * r27733;
        double r27735 = r27734 / r27728;
        double r27736 = 1.324949139524979e-102;
        bool r27737 = r27717 <= r27736;
        double r27738 = r27720 / r27728;
        double r27739 = r27730 * r27728;
        double r27740 = r27725 - r27739;
        double r27741 = exp(r27740);
        double r27742 = r27724 / r27741;
        double r27743 = r27738 * r27742;
        double r27744 = r27737 ? r27743 : r27735;
        double r27745 = r27719 ? r27735 : r27744;
        return r27745;
}

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 r27746, r27747, r27748, r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756, r27757, r27758, r27759, r27760, r27761, r27762;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r27746);
        mpfr_init(r27747);
        mpfr_init(r27748);
        mpfr_init(r27749);
        mpfr_init(r27750);
        mpfr_init(r27751);
        mpfr_init_set_str(r27752, "1.0", 10, MPFR_RNDN);
        mpfr_init(r27753);
        mpfr_init(r27754);
        mpfr_init(r27755);
        mpfr_init(r27756);
        mpfr_init(r27757);
        mpfr_init(r27758);
        mpfr_init(r27759);
        mpfr_init(r27760);
        mpfr_init(r27761);
        mpfr_init(r27762);
}

double f_im(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r27746, x, MPFR_RNDN);
        mpfr_set_d(r27747, y, MPFR_RNDN);
        mpfr_set_d(r27748, z, MPFR_RNDN);
        mpfr_log(r27749, r27748, MPFR_RNDN);
        mpfr_mul(r27750, r27747, r27749, MPFR_RNDN);
        mpfr_set_d(r27751, t, MPFR_RNDN);
        ;
        mpfr_sub(r27753, r27751, r27752, MPFR_RNDN);
        mpfr_set_d(r27754, a, MPFR_RNDN);
        mpfr_log(r27755, r27754, MPFR_RNDN);
        mpfr_mul(r27756, r27753, r27755, MPFR_RNDN);
        mpfr_add(r27757, r27750, r27756, MPFR_RNDN);
        mpfr_set_d(r27758, b, MPFR_RNDN);
        mpfr_sub(r27759, r27757, r27758, MPFR_RNDN);
        mpfr_exp(r27760, r27759, MPFR_RNDN);
        mpfr_mul(r27761, r27746, r27760, MPFR_RNDN);
        mpfr_div(r27762, r27761, r27747, MPFR_RNDN);
        return mpfr_get_d(r27762, MPFR_RNDN);
}

static mpfr_t r27763, r27764, r27765, r27766, r27767, r27768, r27769, r27770, r27771, r27772, r27773, r27774, r27775, r27776, r27777, r27778, r27779, r27780, r27781, r27782, r27783, r27784, r27785, r27786, r27787, r27788, r27789, r27790, r27791;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r27763);
        mpfr_init_set_str(r27764, "-1486774482200.535", 10, MPFR_RNDN);
        mpfr_init(r27765);
        mpfr_init(r27766);
        mpfr_init(r27767);
        mpfr_init_set_str(r27768, "1.0", 10, MPFR_RNDN);
        mpfr_init(r27769);
        mpfr_init(r27770);
        mpfr_init(r27771);
        mpfr_init_set_str(r27772, "1", 10, MPFR_RNDN);
        mpfr_init(r27773);
        mpfr_init(r27774);
        mpfr_init(r27775);
        mpfr_init(r27776);
        mpfr_init(r27777);
        mpfr_init(r27778);
        mpfr_init(r27779);
        mpfr_init(r27780);
        mpfr_init(r27781);
        mpfr_init_set_str(r27782, "1.324949139524979e-102", 10, MPFR_RNDN);
        mpfr_init(r27783);
        mpfr_init(r27784);
        mpfr_init(r27785);
        mpfr_init(r27786);
        mpfr_init(r27787);
        mpfr_init(r27788);
        mpfr_init(r27789);
        mpfr_init(r27790);
        mpfr_init(r27791);
}

double f_fm(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r27763, t, MPFR_RNDN);
        ;
        mpfr_set_si(r27765, mpfr_cmp(r27763, r27764) <= 0, MPFR_RNDN);
        mpfr_set_d(r27766, x, MPFR_RNDN);
        mpfr_set_d(r27767, a, MPFR_RNDN);
        ;
        mpfr_sub(r27769, r27763, r27768, MPFR_RNDN);
        mpfr_pow(r27770, r27767, r27769, MPFR_RNDN);
        mpfr_set_d(r27771, b, MPFR_RNDN);
        ;
        mpfr_add(r27773, r27771, r27772, MPFR_RNDN);
        mpfr_set_d(r27774, y, MPFR_RNDN);
        mpfr_set_d(r27775, z, MPFR_RNDN);
        mpfr_log(r27776, r27775, MPFR_RNDN);
        mpfr_mul(r27777, r27774, r27776, MPFR_RNDN);
        mpfr_sub(r27778, r27773, r27777, MPFR_RNDN);
        mpfr_div(r27779, r27770, r27778, MPFR_RNDN);
        mpfr_mul(r27780, r27766, r27779, MPFR_RNDN);
        mpfr_div(r27781, r27780, r27774, MPFR_RNDN);
        ;
        mpfr_set_si(r27783, mpfr_cmp(r27763, r27782) <= 0, MPFR_RNDN);
        mpfr_div(r27784, r27766, r27774, MPFR_RNDN);
        mpfr_mul(r27785, r27776, r27774, MPFR_RNDN);
        mpfr_sub(r27786, r27771, r27785, MPFR_RNDN);
        mpfr_exp(r27787, r27786, MPFR_RNDN);
        mpfr_div(r27788, r27770, r27787, MPFR_RNDN);
        mpfr_mul(r27789, r27784, r27788, MPFR_RNDN);
        if (mpfr_get_si(r27783, MPFR_RNDN)) { mpfr_set(r27790, r27789, MPFR_RNDN); } else { mpfr_set(r27790, r27781, MPFR_RNDN); };
        if (mpfr_get_si(r27765, MPFR_RNDN)) { mpfr_set(r27791, r27781, MPFR_RNDN); } else { mpfr_set(r27791, r27790, MPFR_RNDN); };
        return mpfr_get_d(r27791, MPFR_RNDN);
}

static mpfr_t r27792, r27793, r27794, r27795, r27796, r27797, r27798, r27799, r27800, r27801, r27802, r27803, r27804, r27805, r27806, r27807, r27808, r27809, r27810, r27811, r27812, r27813, r27814, r27815, r27816, r27817, r27818, r27819, r27820;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r27792);
        mpfr_init_set_str(r27793, "-1486774482200.535", 10, MPFR_RNDN);
        mpfr_init(r27794);
        mpfr_init(r27795);
        mpfr_init(r27796);
        mpfr_init_set_str(r27797, "1.0", 10, MPFR_RNDN);
        mpfr_init(r27798);
        mpfr_init(r27799);
        mpfr_init(r27800);
        mpfr_init_set_str(r27801, "1", 10, MPFR_RNDN);
        mpfr_init(r27802);
        mpfr_init(r27803);
        mpfr_init(r27804);
        mpfr_init(r27805);
        mpfr_init(r27806);
        mpfr_init(r27807);
        mpfr_init(r27808);
        mpfr_init(r27809);
        mpfr_init(r27810);
        mpfr_init_set_str(r27811, "1.324949139524979e-102", 10, MPFR_RNDN);
        mpfr_init(r27812);
        mpfr_init(r27813);
        mpfr_init(r27814);
        mpfr_init(r27815);
        mpfr_init(r27816);
        mpfr_init(r27817);
        mpfr_init(r27818);
        mpfr_init(r27819);
        mpfr_init(r27820);
}

double f_dm(double x, double y, double z, double t, double a, double b) {
        mpfr_set_d(r27792, t, MPFR_RNDN);
        ;
        mpfr_set_si(r27794, mpfr_cmp(r27792, r27793) <= 0, MPFR_RNDN);
        mpfr_set_d(r27795, x, MPFR_RNDN);
        mpfr_set_d(r27796, a, MPFR_RNDN);
        ;
        mpfr_sub(r27798, r27792, r27797, MPFR_RNDN);
        mpfr_pow(r27799, r27796, r27798, MPFR_RNDN);
        mpfr_set_d(r27800, b, MPFR_RNDN);
        ;
        mpfr_add(r27802, r27800, r27801, MPFR_RNDN);
        mpfr_set_d(r27803, y, MPFR_RNDN);
        mpfr_set_d(r27804, z, MPFR_RNDN);
        mpfr_log(r27805, r27804, MPFR_RNDN);
        mpfr_mul(r27806, r27803, r27805, MPFR_RNDN);
        mpfr_sub(r27807, r27802, r27806, MPFR_RNDN);
        mpfr_div(r27808, r27799, r27807, MPFR_RNDN);
        mpfr_mul(r27809, r27795, r27808, MPFR_RNDN);
        mpfr_div(r27810, r27809, r27803, MPFR_RNDN);
        ;
        mpfr_set_si(r27812, mpfr_cmp(r27792, r27811) <= 0, MPFR_RNDN);
        mpfr_div(r27813, r27795, r27803, MPFR_RNDN);
        mpfr_mul(r27814, r27805, r27803, MPFR_RNDN);
        mpfr_sub(r27815, r27800, r27814, MPFR_RNDN);
        mpfr_exp(r27816, r27815, MPFR_RNDN);
        mpfr_div(r27817, r27799, r27816, MPFR_RNDN);
        mpfr_mul(r27818, r27813, r27817, MPFR_RNDN);
        if (mpfr_get_si(r27812, MPFR_RNDN)) { mpfr_set(r27819, r27818, MPFR_RNDN); } else { mpfr_set(r27819, r27810, MPFR_RNDN); };
        if (mpfr_get_si(r27794, MPFR_RNDN)) { mpfr_set(r27820, r27810, MPFR_RNDN); } else { mpfr_set(r27820, r27819, MPFR_RNDN); };
        return mpfr_get_d(r27820, MPFR_RNDN);
}