#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 r24890 = NdChar;
        float r24891 = 1;
        float r24892 = Ec;
        float r24893 = Vef;
        float r24894 = r24892 - r24893;
        float r24895 = EDonor;
        float r24896 = r24894 - r24895;
        float r24897 = mu;
        float r24898 = r24896 - r24897;
        float r24899 = -r24898;
        float r24900 = KbT;
        float r24901 = r24899 / r24900;
        float r24902 = exp(r24901);
        float r24903 = r24891 + r24902;
        float r24904 = r24890 / r24903;
        float r24905 = NaChar;
        float r24906 = Ev;
        float r24907 = r24906 + r24893;
        float r24908 = EAccept;
        float r24909 = r24907 + r24908;
        float r24910 = -r24897;
        float r24911 = r24909 + r24910;
        float r24912 = r24911 / r24900;
        float r24913 = exp(r24912);
        float r24914 = r24891 + r24913;
        float r24915 = r24905 / r24914;
        float r24916 = r24904 + r24915;
        return r24916;
}

double f_id(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        double r24917 = NdChar;
        double r24918 = 1;
        double r24919 = Ec;
        double r24920 = Vef;
        double r24921 = r24919 - r24920;
        double r24922 = EDonor;
        double r24923 = r24921 - r24922;
        double r24924 = mu;
        double r24925 = r24923 - r24924;
        double r24926 = -r24925;
        double r24927 = KbT;
        double r24928 = r24926 / r24927;
        double r24929 = exp(r24928);
        double r24930 = r24918 + r24929;
        double r24931 = r24917 / r24930;
        double r24932 = NaChar;
        double r24933 = Ev;
        double r24934 = r24933 + r24920;
        double r24935 = EAccept;
        double r24936 = r24934 + r24935;
        double r24937 = -r24924;
        double r24938 = r24936 + r24937;
        double r24939 = r24938 / r24927;
        double r24940 = exp(r24939);
        double r24941 = r24918 + r24940;
        double r24942 = r24932 / r24941;
        double r24943 = r24931 + r24942;
        return r24943;
}


double f_of(float NdChar, float Ec, float Vef, float EDonor, float mu, float KbT, float NaChar, float Ev, float EAccept) {
        float r24944 = NdChar;
        float r24945 = 1;
        float r24946 = Ec;
        float r24947 = Vef;
        float r24948 = r24946 - r24947;
        float r24949 = EDonor;
        float r24950 = r24948 - r24949;
        float r24951 = mu;
        float r24952 = r24950 - r24951;
        float r24953 = -r24952;
        float r24954 = KbT;
        float r24955 = r24953 / r24954;
        float r24956 = exp(r24955);
        float r24957 = r24945 + r24956;
        float r24958 = r24944 / r24957;
        float r24959 = Ev;
        float r24960 = r24959 + r24947;
        float r24961 = EAccept;
        float r24962 = r24961 - r24951;
        float r24963 = r24960 + r24962;
        float r24964 = r24963 / r24954;
        float r24965 = exp(r24964);
        float r24966 = r24945 + r24965;
        float r24967 = NaChar;
        float r24968 = r24966 / r24967;
        float r24969 = r24945 / r24968;
        float r24970 = r24958 + r24969;
        return r24970;
}

double f_od(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        double r24971 = NdChar;
        double r24972 = 1;
        double r24973 = Ec;
        double r24974 = Vef;
        double r24975 = r24973 - r24974;
        double r24976 = EDonor;
        double r24977 = r24975 - r24976;
        double r24978 = mu;
        double r24979 = r24977 - r24978;
        double r24980 = -r24979;
        double r24981 = KbT;
        double r24982 = r24980 / r24981;
        double r24983 = exp(r24982);
        double r24984 = r24972 + r24983;
        double r24985 = r24971 / r24984;
        double r24986 = Ev;
        double r24987 = r24986 + r24974;
        double r24988 = EAccept;
        double r24989 = r24988 - r24978;
        double r24990 = r24987 + r24989;
        double r24991 = r24990 / r24981;
        double r24992 = exp(r24991);
        double r24993 = r24972 + r24992;
        double r24994 = NaChar;
        double r24995 = r24993 / r24994;
        double r24996 = r24972 / r24995;
        double r24997 = r24985 + r24996;
        return r24997;
}

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 r24998, r24999, r25000, r25001, r25002, r25003, r25004, r25005, r25006, r25007, r25008, r25009, r25010, r25011, r25012, r25013, r25014, r25015, r25016, r25017, r25018, r25019, r25020, r25021, r25022, r25023, r25024;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(848);
        mpfr_init(r24998);
        mpfr_init_set_str(r24999, "1", 10, MPFR_RNDN);
        mpfr_init(r25000);
        mpfr_init(r25001);
        mpfr_init(r25002);
        mpfr_init(r25003);
        mpfr_init(r25004);
        mpfr_init(r25005);
        mpfr_init(r25006);
        mpfr_init(r25007);
        mpfr_init(r25008);
        mpfr_init(r25009);
        mpfr_init(r25010);
        mpfr_init(r25011);
        mpfr_init(r25012);
        mpfr_init(r25013);
        mpfr_init(r25014);
        mpfr_init(r25015);
        mpfr_init(r25016);
        mpfr_init(r25017);
        mpfr_init(r25018);
        mpfr_init(r25019);
        mpfr_init(r25020);
        mpfr_init(r25021);
        mpfr_init(r25022);
        mpfr_init(r25023);
        mpfr_init(r25024);
}

double f_im(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        mpfr_set_d(r24998, NdChar, MPFR_RNDN);
        ;
        mpfr_set_d(r25000, Ec, MPFR_RNDN);
        mpfr_set_d(r25001, Vef, MPFR_RNDN);
        mpfr_sub(r25002, r25000, r25001, MPFR_RNDN);
        mpfr_set_d(r25003, EDonor, MPFR_RNDN);
        mpfr_sub(r25004, r25002, r25003, MPFR_RNDN);
        mpfr_set_d(r25005, mu, MPFR_RNDN);
        mpfr_sub(r25006, r25004, r25005, MPFR_RNDN);
        mpfr_neg(r25007, r25006, MPFR_RNDN);
        mpfr_set_d(r25008, KbT, MPFR_RNDN);
        mpfr_div(r25009, r25007, r25008, MPFR_RNDN);
        mpfr_exp(r25010, r25009, MPFR_RNDN);
        mpfr_add(r25011, r24999, r25010, MPFR_RNDN);
        mpfr_div(r25012, r24998, r25011, MPFR_RNDN);
        mpfr_set_d(r25013, NaChar, MPFR_RNDN);
        mpfr_set_d(r25014, Ev, MPFR_RNDN);
        mpfr_add(r25015, r25014, r25001, MPFR_RNDN);
        mpfr_set_d(r25016, EAccept, MPFR_RNDN);
        mpfr_add(r25017, r25015, r25016, MPFR_RNDN);
        mpfr_neg(r25018, r25005, MPFR_RNDN);
        mpfr_add(r25019, r25017, r25018, MPFR_RNDN);
        mpfr_div(r25020, r25019, r25008, MPFR_RNDN);
        mpfr_exp(r25021, r25020, MPFR_RNDN);
        mpfr_add(r25022, r24999, r25021, MPFR_RNDN);
        mpfr_div(r25023, r25013, r25022, MPFR_RNDN);
        mpfr_add(r25024, r25012, r25023, MPFR_RNDN);
        return mpfr_get_d(r25024, MPFR_RNDN);
}

static mpfr_t r25025, r25026, r25027, r25028, r25029, r25030, r25031, r25032, r25033, r25034, r25035, r25036, r25037, r25038, r25039, r25040, r25041, r25042, r25043, r25044, r25045, r25046, r25047, r25048, r25049, r25050, r25051;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(848);
        mpfr_init(r25025);
        mpfr_init_set_str(r25026, "1", 10, MPFR_RNDN);
        mpfr_init(r25027);
        mpfr_init(r25028);
        mpfr_init(r25029);
        mpfr_init(r25030);
        mpfr_init(r25031);
        mpfr_init(r25032);
        mpfr_init(r25033);
        mpfr_init(r25034);
        mpfr_init(r25035);
        mpfr_init(r25036);
        mpfr_init(r25037);
        mpfr_init(r25038);
        mpfr_init(r25039);
        mpfr_init(r25040);
        mpfr_init(r25041);
        mpfr_init(r25042);
        mpfr_init(r25043);
        mpfr_init(r25044);
        mpfr_init(r25045);
        mpfr_init(r25046);
        mpfr_init(r25047);
        mpfr_init(r25048);
        mpfr_init(r25049);
        mpfr_init(r25050);
        mpfr_init(r25051);
}

double f_fm(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        mpfr_set_d(r25025, NdChar, MPFR_RNDN);
        ;
        mpfr_set_d(r25027, Ec, MPFR_RNDN);
        mpfr_set_d(r25028, Vef, MPFR_RNDN);
        mpfr_sub(r25029, r25027, r25028, MPFR_RNDN);
        mpfr_set_d(r25030, EDonor, MPFR_RNDN);
        mpfr_sub(r25031, r25029, r25030, MPFR_RNDN);
        mpfr_set_d(r25032, mu, MPFR_RNDN);
        mpfr_sub(r25033, r25031, r25032, MPFR_RNDN);
        mpfr_neg(r25034, r25033, MPFR_RNDN);
        mpfr_set_d(r25035, KbT, MPFR_RNDN);
        mpfr_div(r25036, r25034, r25035, MPFR_RNDN);
        mpfr_exp(r25037, r25036, MPFR_RNDN);
        mpfr_add(r25038, r25026, r25037, MPFR_RNDN);
        mpfr_div(r25039, r25025, r25038, MPFR_RNDN);
        mpfr_set_d(r25040, Ev, MPFR_RNDN);
        mpfr_add(r25041, r25040, r25028, MPFR_RNDN);
        mpfr_set_d(r25042, EAccept, MPFR_RNDN);
        mpfr_sub(r25043, r25042, r25032, MPFR_RNDN);
        mpfr_add(r25044, r25041, r25043, MPFR_RNDN);
        mpfr_div(r25045, r25044, r25035, MPFR_RNDN);
        mpfr_exp(r25046, r25045, MPFR_RNDN);
        mpfr_add(r25047, r25026, r25046, MPFR_RNDN);
        mpfr_set_d(r25048, NaChar, MPFR_RNDN);
        mpfr_div(r25049, r25047, r25048, MPFR_RNDN);
        mpfr_div(r25050, r25026, r25049, MPFR_RNDN);
        mpfr_add(r25051, r25039, r25050, MPFR_RNDN);
        return mpfr_get_d(r25051, MPFR_RNDN);
}

static mpfr_t r25052, r25053, r25054, r25055, r25056, r25057, r25058, r25059, r25060, r25061, r25062, r25063, r25064, r25065, r25066, r25067, r25068, r25069, r25070, r25071, r25072, r25073, r25074, r25075, r25076, r25077, r25078;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(848);
        mpfr_init(r25052);
        mpfr_init_set_str(r25053, "1", 10, MPFR_RNDN);
        mpfr_init(r25054);
        mpfr_init(r25055);
        mpfr_init(r25056);
        mpfr_init(r25057);
        mpfr_init(r25058);
        mpfr_init(r25059);
        mpfr_init(r25060);
        mpfr_init(r25061);
        mpfr_init(r25062);
        mpfr_init(r25063);
        mpfr_init(r25064);
        mpfr_init(r25065);
        mpfr_init(r25066);
        mpfr_init(r25067);
        mpfr_init(r25068);
        mpfr_init(r25069);
        mpfr_init(r25070);
        mpfr_init(r25071);
        mpfr_init(r25072);
        mpfr_init(r25073);
        mpfr_init(r25074);
        mpfr_init(r25075);
        mpfr_init(r25076);
        mpfr_init(r25077);
        mpfr_init(r25078);
}

double f_dm(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        mpfr_set_d(r25052, NdChar, MPFR_RNDN);
        ;
        mpfr_set_d(r25054, Ec, MPFR_RNDN);
        mpfr_set_d(r25055, Vef, MPFR_RNDN);
        mpfr_sub(r25056, r25054, r25055, MPFR_RNDN);
        mpfr_set_d(r25057, EDonor, MPFR_RNDN);
        mpfr_sub(r25058, r25056, r25057, MPFR_RNDN);
        mpfr_set_d(r25059, mu, MPFR_RNDN);
        mpfr_sub(r25060, r25058, r25059, MPFR_RNDN);
        mpfr_neg(r25061, r25060, MPFR_RNDN);
        mpfr_set_d(r25062, KbT, MPFR_RNDN);
        mpfr_div(r25063, r25061, r25062, MPFR_RNDN);
        mpfr_exp(r25064, r25063, MPFR_RNDN);
        mpfr_add(r25065, r25053, r25064, MPFR_RNDN);
        mpfr_div(r25066, r25052, r25065, MPFR_RNDN);
        mpfr_set_d(r25067, Ev, MPFR_RNDN);
        mpfr_add(r25068, r25067, r25055, MPFR_RNDN);
        mpfr_set_d(r25069, EAccept, MPFR_RNDN);
        mpfr_sub(r25070, r25069, r25059, MPFR_RNDN);
        mpfr_add(r25071, r25068, r25070, MPFR_RNDN);
        mpfr_div(r25072, r25071, r25062, MPFR_RNDN);
        mpfr_exp(r25073, r25072, MPFR_RNDN);
        mpfr_add(r25074, r25053, r25073, MPFR_RNDN);
        mpfr_set_d(r25075, NaChar, MPFR_RNDN);
        mpfr_div(r25076, r25074, r25075, MPFR_RNDN);
        mpfr_div(r25077, r25053, r25076, MPFR_RNDN);
        mpfr_add(r25078, r25066, r25077, MPFR_RNDN);
        return mpfr_get_d(r25078, MPFR_RNDN);
}