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

char *name = "Bulmash initializePoisson";

double f_if(float NdChar, float Ec, float Vef, float EDonor, float mu, float KbT, float NaChar, float Ev, float EAccept) {
        float r27664 = NdChar;
        float r27665 = 1;
        float r27666 = Ec;
        float r27667 = Vef;
        float r27668 = r27666 - r27667;
        float r27669 = EDonor;
        float r27670 = r27668 - r27669;
        float r27671 = mu;
        float r27672 = r27670 - r27671;
        float r27673 = -r27672;
        float r27674 = KbT;
        float r27675 = r27673 / r27674;
        float r27676 = exp(r27675);
        float r27677 = r27665 + r27676;
        float r27678 = r27664 / r27677;
        float r27679 = NaChar;
        float r27680 = Ev;
        float r27681 = r27680 + r27667;
        float r27682 = EAccept;
        float r27683 = r27681 + r27682;
        float r27684 = -r27671;
        float r27685 = r27683 + r27684;
        float r27686 = r27685 / r27674;
        float r27687 = exp(r27686);
        float r27688 = r27665 + r27687;
        float r27689 = r27679 / r27688;
        float r27690 = r27678 + r27689;
        return r27690;
}

double f_id(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        double r27691 = NdChar;
        double r27692 = 1;
        double r27693 = Ec;
        double r27694 = Vef;
        double r27695 = r27693 - r27694;
        double r27696 = EDonor;
        double r27697 = r27695 - r27696;
        double r27698 = mu;
        double r27699 = r27697 - r27698;
        double r27700 = -r27699;
        double r27701 = KbT;
        double r27702 = r27700 / r27701;
        double r27703 = exp(r27702);
        double r27704 = r27692 + r27703;
        double r27705 = r27691 / r27704;
        double r27706 = NaChar;
        double r27707 = Ev;
        double r27708 = r27707 + r27694;
        double r27709 = EAccept;
        double r27710 = r27708 + r27709;
        double r27711 = -r27698;
        double r27712 = r27710 + r27711;
        double r27713 = r27712 / r27701;
        double r27714 = exp(r27713);
        double r27715 = r27692 + r27714;
        double r27716 = r27706 / r27715;
        double r27717 = r27705 + r27716;
        return r27717;
}


double f_of(float NdChar, float Ec, float Vef, float EDonor, float mu, float KbT, float NaChar, float Ev, float EAccept) {
        float r27718 = NdChar;
        float r27719 = 1;
        float r27720 = Vef;
        float r27721 = Ec;
        float r27722 = r27720 - r27721;
        float r27723 = mu;
        float r27724 = EDonor;
        float r27725 = r27723 + r27724;
        float r27726 = r27722 + r27725;
        float r27727 = KbT;
        float r27728 = r27726 / r27727;
        float r27729 = exp(r27728);
        float r27730 = log1p(r27729);
        float r27731 = expm1(r27730);
        float r27732 = r27719 + r27731;
        float r27733 = r27718 / r27732;
        float r27734 = NaChar;
        float r27735 = Ev;
        float r27736 = r27735 + r27720;
        float r27737 = EAccept;
        float r27738 = r27736 + r27737;
        float r27739 = -r27723;
        float r27740 = r27738 + r27739;
        float r27741 = r27740 / r27727;
        float r27742 = exp(r27741);
        float r27743 = r27719 + r27742;
        float r27744 = r27734 / r27743;
        float r27745 = r27733 + r27744;
        return r27745;
}

double f_od(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        double r27746 = NdChar;
        double r27747 = 1;
        double r27748 = Vef;
        double r27749 = Ec;
        double r27750 = r27748 - r27749;
        double r27751 = mu;
        double r27752 = EDonor;
        double r27753 = r27751 + r27752;
        double r27754 = r27750 + r27753;
        double r27755 = KbT;
        double r27756 = r27754 / r27755;
        double r27757 = exp(r27756);
        double r27758 = log1p(r27757);
        double r27759 = expm1(r27758);
        double r27760 = r27747 + r27759;
        double r27761 = r27746 / r27760;
        double r27762 = NaChar;
        double r27763 = Ev;
        double r27764 = r27763 + r27748;
        double r27765 = EAccept;
        double r27766 = r27764 + r27765;
        double r27767 = -r27751;
        double r27768 = r27766 + r27767;
        double r27769 = r27768 / r27755;
        double r27770 = exp(r27769);
        double r27771 = r27747 + r27770;
        double r27772 = r27762 / r27771;
        double r27773 = r27761 + r27772;
        return r27773;
}

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 r27774, r27775, r27776, r27777, r27778, r27779, r27780, r27781, r27782, r27783, r27784, r27785, r27786, r27787, r27788, r27789, r27790, r27791, r27792, r27793, r27794, r27795, r27796, r27797, r27798, r27799, r27800;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r27774);
        mpfr_init_set_str(r27775, "1", 10, MPFR_RNDN);
        mpfr_init(r27776);
        mpfr_init(r27777);
        mpfr_init(r27778);
        mpfr_init(r27779);
        mpfr_init(r27780);
        mpfr_init(r27781);
        mpfr_init(r27782);
        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);
        mpfr_init(r27792);
        mpfr_init(r27793);
        mpfr_init(r27794);
        mpfr_init(r27795);
        mpfr_init(r27796);
        mpfr_init(r27797);
        mpfr_init(r27798);
        mpfr_init(r27799);
        mpfr_init(r27800);
}

double f_im(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        mpfr_set_d(r27774, NdChar, MPFR_RNDN);
        ;
        mpfr_set_d(r27776, Ec, MPFR_RNDN);
        mpfr_set_d(r27777, Vef, MPFR_RNDN);
        mpfr_sub(r27778, r27776, r27777, MPFR_RNDN);
        mpfr_set_d(r27779, EDonor, MPFR_RNDN);
        mpfr_sub(r27780, r27778, r27779, MPFR_RNDN);
        mpfr_set_d(r27781, mu, MPFR_RNDN);
        mpfr_sub(r27782, r27780, r27781, MPFR_RNDN);
        mpfr_neg(r27783, r27782, MPFR_RNDN);
        mpfr_set_d(r27784, KbT, MPFR_RNDN);
        mpfr_div(r27785, r27783, r27784, MPFR_RNDN);
        mpfr_exp(r27786, r27785, MPFR_RNDN);
        mpfr_add(r27787, r27775, r27786, MPFR_RNDN);
        mpfr_div(r27788, r27774, r27787, MPFR_RNDN);
        mpfr_set_d(r27789, NaChar, MPFR_RNDN);
        mpfr_set_d(r27790, Ev, MPFR_RNDN);
        mpfr_add(r27791, r27790, r27777, MPFR_RNDN);
        mpfr_set_d(r27792, EAccept, MPFR_RNDN);
        mpfr_add(r27793, r27791, r27792, MPFR_RNDN);
        mpfr_neg(r27794, r27781, MPFR_RNDN);
        mpfr_add(r27795, r27793, r27794, MPFR_RNDN);
        mpfr_div(r27796, r27795, r27784, MPFR_RNDN);
        mpfr_exp(r27797, r27796, MPFR_RNDN);
        mpfr_add(r27798, r27775, r27797, MPFR_RNDN);
        mpfr_div(r27799, r27789, r27798, MPFR_RNDN);
        mpfr_add(r27800, r27788, r27799, MPFR_RNDN);
        return mpfr_get_d(r27800, MPFR_RNDN);
}

static mpfr_t r27801, r27802, r27803, r27804, r27805, r27806, r27807, r27808, r27809, r27810, r27811, r27812, r27813, r27814, r27815, r27816, r27817, r27818, r27819, r27820, r27821, r27822, r27823, r27824, r27825, r27826, r27827, r27828;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27801);
        mpfr_init_set_str(r27802, "1", 10, MPFR_RNDN);
        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(r27811);
        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);
        mpfr_init(r27821);
        mpfr_init(r27822);
        mpfr_init(r27823);
        mpfr_init(r27824);
        mpfr_init(r27825);
        mpfr_init(r27826);
        mpfr_init(r27827);
        mpfr_init(r27828);
}

double f_fm(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        mpfr_set_d(r27801, NdChar, MPFR_RNDN);
        ;
        mpfr_set_d(r27803, Vef, MPFR_RNDN);
        mpfr_set_d(r27804, Ec, MPFR_RNDN);
        mpfr_sub(r27805, r27803, r27804, MPFR_RNDN);
        mpfr_set_d(r27806, mu, MPFR_RNDN);
        mpfr_set_d(r27807, EDonor, MPFR_RNDN);
        mpfr_add(r27808, r27806, r27807, MPFR_RNDN);
        mpfr_add(r27809, r27805, r27808, MPFR_RNDN);
        mpfr_set_d(r27810, KbT, MPFR_RNDN);
        mpfr_div(r27811, r27809, r27810, MPFR_RNDN);
        mpfr_exp(r27812, r27811, MPFR_RNDN);
        mpfr_log1p(r27813, r27812, MPFR_RNDN);
        mpfr_expm1(r27814, r27813, MPFR_RNDN);
        mpfr_add(r27815, r27802, r27814, MPFR_RNDN);
        mpfr_div(r27816, r27801, r27815, MPFR_RNDN);
        mpfr_set_d(r27817, NaChar, MPFR_RNDN);
        mpfr_set_d(r27818, Ev, MPFR_RNDN);
        mpfr_add(r27819, r27818, r27803, MPFR_RNDN);
        mpfr_set_d(r27820, EAccept, MPFR_RNDN);
        mpfr_add(r27821, r27819, r27820, MPFR_RNDN);
        mpfr_neg(r27822, r27806, MPFR_RNDN);
        mpfr_add(r27823, r27821, r27822, MPFR_RNDN);
        mpfr_div(r27824, r27823, r27810, MPFR_RNDN);
        mpfr_exp(r27825, r27824, MPFR_RNDN);
        mpfr_add(r27826, r27802, r27825, MPFR_RNDN);
        mpfr_div(r27827, r27817, r27826, MPFR_RNDN);
        mpfr_add(r27828, r27816, r27827, MPFR_RNDN);
        return mpfr_get_d(r27828, MPFR_RNDN);
}

static mpfr_t r27829, r27830, r27831, r27832, r27833, r27834, r27835, r27836, r27837, r27838, r27839, r27840, r27841, r27842, r27843, r27844, r27845, r27846, r27847, r27848, r27849, r27850, r27851, r27852, r27853, r27854, r27855, r27856;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r27829);
        mpfr_init_set_str(r27830, "1", 10, MPFR_RNDN);
        mpfr_init(r27831);
        mpfr_init(r27832);
        mpfr_init(r27833);
        mpfr_init(r27834);
        mpfr_init(r27835);
        mpfr_init(r27836);
        mpfr_init(r27837);
        mpfr_init(r27838);
        mpfr_init(r27839);
        mpfr_init(r27840);
        mpfr_init(r27841);
        mpfr_init(r27842);
        mpfr_init(r27843);
        mpfr_init(r27844);
        mpfr_init(r27845);
        mpfr_init(r27846);
        mpfr_init(r27847);
        mpfr_init(r27848);
        mpfr_init(r27849);
        mpfr_init(r27850);
        mpfr_init(r27851);
        mpfr_init(r27852);
        mpfr_init(r27853);
        mpfr_init(r27854);
        mpfr_init(r27855);
        mpfr_init(r27856);
}

double f_dm(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        mpfr_set_d(r27829, NdChar, MPFR_RNDN);
        ;
        mpfr_set_d(r27831, Vef, MPFR_RNDN);
        mpfr_set_d(r27832, Ec, MPFR_RNDN);
        mpfr_sub(r27833, r27831, r27832, MPFR_RNDN);
        mpfr_set_d(r27834, mu, MPFR_RNDN);
        mpfr_set_d(r27835, EDonor, MPFR_RNDN);
        mpfr_add(r27836, r27834, r27835, MPFR_RNDN);
        mpfr_add(r27837, r27833, r27836, MPFR_RNDN);
        mpfr_set_d(r27838, KbT, MPFR_RNDN);
        mpfr_div(r27839, r27837, r27838, MPFR_RNDN);
        mpfr_exp(r27840, r27839, MPFR_RNDN);
        mpfr_log1p(r27841, r27840, MPFR_RNDN);
        mpfr_expm1(r27842, r27841, MPFR_RNDN);
        mpfr_add(r27843, r27830, r27842, MPFR_RNDN);
        mpfr_div(r27844, r27829, r27843, MPFR_RNDN);
        mpfr_set_d(r27845, NaChar, MPFR_RNDN);
        mpfr_set_d(r27846, Ev, MPFR_RNDN);
        mpfr_add(r27847, r27846, r27831, MPFR_RNDN);
        mpfr_set_d(r27848, EAccept, MPFR_RNDN);
        mpfr_add(r27849, r27847, r27848, MPFR_RNDN);
        mpfr_neg(r27850, r27834, MPFR_RNDN);
        mpfr_add(r27851, r27849, r27850, MPFR_RNDN);
        mpfr_div(r27852, r27851, r27838, MPFR_RNDN);
        mpfr_exp(r27853, r27852, MPFR_RNDN);
        mpfr_add(r27854, r27830, r27853, MPFR_RNDN);
        mpfr_div(r27855, r27845, r27854, MPFR_RNDN);
        mpfr_add(r27856, r27844, r27855, MPFR_RNDN);
        return mpfr_get_d(r27856, MPFR_RNDN);
}