#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 r21746 = NdChar;
        float r21747 = 1.0f;
        float r21748 = Ec;
        float r21749 = Vef;
        float r21750 = r21748 - r21749;
        float r21751 = EDonor;
        float r21752 = r21750 - r21751;
        float r21753 = mu;
        float r21754 = r21752 - r21753;
        float r21755 = -r21754;
        float r21756 = KbT;
        float r21757 = r21755 / r21756;
        float r21758 = exp(r21757);
        float r21759 = r21747 + r21758;
        float r21760 = r21746 / r21759;
        float r21761 = NaChar;
        float r21762 = Ev;
        float r21763 = r21762 + r21749;
        float r21764 = EAccept;
        float r21765 = r21763 + r21764;
        float r21766 = -r21753;
        float r21767 = r21765 + r21766;
        float r21768 = r21767 / r21756;
        float r21769 = exp(r21768);
        float r21770 = r21747 + r21769;
        float r21771 = r21761 / r21770;
        float r21772 = r21760 + r21771;
        return r21772;
}

double f_id(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        double r21773 = NdChar;
        double r21774 = 1.0;
        double r21775 = Ec;
        double r21776 = Vef;
        double r21777 = r21775 - r21776;
        double r21778 = EDonor;
        double r21779 = r21777 - r21778;
        double r21780 = mu;
        double r21781 = r21779 - r21780;
        double r21782 = -r21781;
        double r21783 = KbT;
        double r21784 = r21782 / r21783;
        double r21785 = exp(r21784);
        double r21786 = r21774 + r21785;
        double r21787 = r21773 / r21786;
        double r21788 = NaChar;
        double r21789 = Ev;
        double r21790 = r21789 + r21776;
        double r21791 = EAccept;
        double r21792 = r21790 + r21791;
        double r21793 = -r21780;
        double r21794 = r21792 + r21793;
        double r21795 = r21794 / r21783;
        double r21796 = exp(r21795);
        double r21797 = r21774 + r21796;
        double r21798 = r21788 / r21797;
        double r21799 = r21787 + r21798;
        return r21799;
}


double f_of(float NdChar, float Ec, float Vef, float EDonor, float mu, float KbT, float NaChar, float Ev, float EAccept) {
        float r21800 = NdChar;
        float r21801 = 1.0f;
        float r21802 = Ec;
        float r21803 = Vef;
        float r21804 = r21802 - r21803;
        float r21805 = EDonor;
        float r21806 = r21804 - r21805;
        float r21807 = mu;
        float r21808 = r21806 - r21807;
        float r21809 = -r21808;
        float r21810 = KbT;
        float r21811 = r21809 / r21810;
        float r21812 = exp(r21811);
        float r21813 = r21801 + r21812;
        float r21814 = r21800 / r21813;
        float r21815 = NaChar;
        float r21816 = Ev;
        float r21817 = r21816 + r21803;
        float r21818 = EAccept;
        float r21819 = r21817 + r21818;
        float r21820 = -r21807;
        float r21821 = r21819 + r21820;
        float r21822 = r21821 / r21810;
        float r21823 = exp(r21822);
        float r21824 = r21801 + r21823;
        float r21825 = r21815 / r21824;
        float r21826 = r21814 + r21825;
        return r21826;
}

double f_od(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        double r21827 = NdChar;
        double r21828 = 1.0;
        double r21829 = Ec;
        double r21830 = Vef;
        double r21831 = r21829 - r21830;
        double r21832 = EDonor;
        double r21833 = r21831 - r21832;
        double r21834 = mu;
        double r21835 = r21833 - r21834;
        double r21836 = -r21835;
        double r21837 = KbT;
        double r21838 = r21836 / r21837;
        double r21839 = exp(r21838);
        double r21840 = r21828 + r21839;
        double r21841 = r21827 / r21840;
        double r21842 = NaChar;
        double r21843 = Ev;
        double r21844 = r21843 + r21830;
        double r21845 = EAccept;
        double r21846 = r21844 + r21845;
        double r21847 = -r21834;
        double r21848 = r21846 + r21847;
        double r21849 = r21848 / r21837;
        double r21850 = exp(r21849);
        double r21851 = r21828 + r21850;
        double r21852 = r21842 / r21851;
        double r21853 = r21841 + r21852;
        return r21853;
}

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 r21854, r21855, r21856, r21857, r21858, r21859, r21860, r21861, r21862, r21863, r21864, r21865, r21866, r21867, r21868, r21869, r21870, r21871, r21872, r21873, r21874, r21875, r21876, r21877, r21878, r21879, r21880;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r21854);
        mpfr_init_set_str(r21855, "1", 10, MPFR_RNDN);
        mpfr_init(r21856);
        mpfr_init(r21857);
        mpfr_init(r21858);
        mpfr_init(r21859);
        mpfr_init(r21860);
        mpfr_init(r21861);
        mpfr_init(r21862);
        mpfr_init(r21863);
        mpfr_init(r21864);
        mpfr_init(r21865);
        mpfr_init(r21866);
        mpfr_init(r21867);
        mpfr_init(r21868);
        mpfr_init(r21869);
        mpfr_init(r21870);
        mpfr_init(r21871);
        mpfr_init(r21872);
        mpfr_init(r21873);
        mpfr_init(r21874);
        mpfr_init(r21875);
        mpfr_init(r21876);
        mpfr_init(r21877);
        mpfr_init(r21878);
        mpfr_init(r21879);
        mpfr_init(r21880);
}

double f_im(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        mpfr_set_d(r21854, NdChar, MPFR_RNDN);
        ;
        mpfr_set_d(r21856, Ec, MPFR_RNDN);
        mpfr_set_d(r21857, Vef, MPFR_RNDN);
        mpfr_sub(r21858, r21856, r21857, MPFR_RNDN);
        mpfr_set_d(r21859, EDonor, MPFR_RNDN);
        mpfr_sub(r21860, r21858, r21859, MPFR_RNDN);
        mpfr_set_d(r21861, mu, MPFR_RNDN);
        mpfr_sub(r21862, r21860, r21861, MPFR_RNDN);
        mpfr_neg(r21863, r21862, MPFR_RNDN);
        mpfr_set_d(r21864, KbT, MPFR_RNDN);
        mpfr_div(r21865, r21863, r21864, MPFR_RNDN);
        mpfr_exp(r21866, r21865, MPFR_RNDN);
        mpfr_add(r21867, r21855, r21866, MPFR_RNDN);
        mpfr_div(r21868, r21854, r21867, MPFR_RNDN);
        mpfr_set_d(r21869, NaChar, MPFR_RNDN);
        mpfr_set_d(r21870, Ev, MPFR_RNDN);
        mpfr_add(r21871, r21870, r21857, MPFR_RNDN);
        mpfr_set_d(r21872, EAccept, MPFR_RNDN);
        mpfr_add(r21873, r21871, r21872, MPFR_RNDN);
        mpfr_neg(r21874, r21861, MPFR_RNDN);
        mpfr_add(r21875, r21873, r21874, MPFR_RNDN);
        mpfr_div(r21876, r21875, r21864, MPFR_RNDN);
        mpfr_exp(r21877, r21876, MPFR_RNDN);
        mpfr_add(r21878, r21855, r21877, MPFR_RNDN);
        mpfr_div(r21879, r21869, r21878, MPFR_RNDN);
        mpfr_add(r21880, r21868, r21879, MPFR_RNDN);
        return mpfr_get_d(r21880, MPFR_RNDN);
}

static mpfr_t r21881, r21882, r21883, r21884, r21885, r21886, r21887, r21888, r21889, r21890, r21891, r21892, r21893, r21894, r21895, r21896, r21897, r21898, r21899, r21900, r21901, r21902, r21903, r21904, r21905, r21906, r21907;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r21881);
        mpfr_init_set_str(r21882, "1", 10, MPFR_RNDN);
        mpfr_init(r21883);
        mpfr_init(r21884);
        mpfr_init(r21885);
        mpfr_init(r21886);
        mpfr_init(r21887);
        mpfr_init(r21888);
        mpfr_init(r21889);
        mpfr_init(r21890);
        mpfr_init(r21891);
        mpfr_init(r21892);
        mpfr_init(r21893);
        mpfr_init(r21894);
        mpfr_init(r21895);
        mpfr_init(r21896);
        mpfr_init(r21897);
        mpfr_init(r21898);
        mpfr_init(r21899);
        mpfr_init(r21900);
        mpfr_init(r21901);
        mpfr_init(r21902);
        mpfr_init(r21903);
        mpfr_init(r21904);
        mpfr_init(r21905);
        mpfr_init(r21906);
        mpfr_init(r21907);
}

double f_fm(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        mpfr_set_d(r21881, NdChar, MPFR_RNDN);
        ;
        mpfr_set_d(r21883, Ec, MPFR_RNDN);
        mpfr_set_d(r21884, Vef, MPFR_RNDN);
        mpfr_sub(r21885, r21883, r21884, MPFR_RNDN);
        mpfr_set_d(r21886, EDonor, MPFR_RNDN);
        mpfr_sub(r21887, r21885, r21886, MPFR_RNDN);
        mpfr_set_d(r21888, mu, MPFR_RNDN);
        mpfr_sub(r21889, r21887, r21888, MPFR_RNDN);
        mpfr_neg(r21890, r21889, MPFR_RNDN);
        mpfr_set_d(r21891, KbT, MPFR_RNDN);
        mpfr_div(r21892, r21890, r21891, MPFR_RNDN);
        mpfr_exp(r21893, r21892, MPFR_RNDN);
        mpfr_add(r21894, r21882, r21893, MPFR_RNDN);
        mpfr_div(r21895, r21881, r21894, MPFR_RNDN);
        mpfr_set_d(r21896, NaChar, MPFR_RNDN);
        mpfr_set_d(r21897, Ev, MPFR_RNDN);
        mpfr_add(r21898, r21897, r21884, MPFR_RNDN);
        mpfr_set_d(r21899, EAccept, MPFR_RNDN);
        mpfr_add(r21900, r21898, r21899, MPFR_RNDN);
        mpfr_neg(r21901, r21888, MPFR_RNDN);
        mpfr_add(r21902, r21900, r21901, MPFR_RNDN);
        mpfr_div(r21903, r21902, r21891, MPFR_RNDN);
        mpfr_exp(r21904, r21903, MPFR_RNDN);
        mpfr_add(r21905, r21882, r21904, MPFR_RNDN);
        mpfr_div(r21906, r21896, r21905, MPFR_RNDN);
        mpfr_add(r21907, r21895, r21906, MPFR_RNDN);
        return mpfr_get_d(r21907, MPFR_RNDN);
}

static mpfr_t r21908, r21909, r21910, r21911, r21912, r21913, r21914, r21915, r21916, r21917, r21918, r21919, r21920, r21921, r21922, r21923, r21924, r21925, r21926, r21927, r21928, r21929, r21930, r21931, r21932, r21933, r21934;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r21908);
        mpfr_init_set_str(r21909, "1", 10, MPFR_RNDN);
        mpfr_init(r21910);
        mpfr_init(r21911);
        mpfr_init(r21912);
        mpfr_init(r21913);
        mpfr_init(r21914);
        mpfr_init(r21915);
        mpfr_init(r21916);
        mpfr_init(r21917);
        mpfr_init(r21918);
        mpfr_init(r21919);
        mpfr_init(r21920);
        mpfr_init(r21921);
        mpfr_init(r21922);
        mpfr_init(r21923);
        mpfr_init(r21924);
        mpfr_init(r21925);
        mpfr_init(r21926);
        mpfr_init(r21927);
        mpfr_init(r21928);
        mpfr_init(r21929);
        mpfr_init(r21930);
        mpfr_init(r21931);
        mpfr_init(r21932);
        mpfr_init(r21933);
        mpfr_init(r21934);
}

double f_dm(double NdChar, double Ec, double Vef, double EDonor, double mu, double KbT, double NaChar, double Ev, double EAccept) {
        mpfr_set_d(r21908, NdChar, MPFR_RNDN);
        ;
        mpfr_set_d(r21910, Ec, MPFR_RNDN);
        mpfr_set_d(r21911, Vef, MPFR_RNDN);
        mpfr_sub(r21912, r21910, r21911, MPFR_RNDN);
        mpfr_set_d(r21913, EDonor, MPFR_RNDN);
        mpfr_sub(r21914, r21912, r21913, MPFR_RNDN);
        mpfr_set_d(r21915, mu, MPFR_RNDN);
        mpfr_sub(r21916, r21914, r21915, MPFR_RNDN);
        mpfr_neg(r21917, r21916, MPFR_RNDN);
        mpfr_set_d(r21918, KbT, MPFR_RNDN);
        mpfr_div(r21919, r21917, r21918, MPFR_RNDN);
        mpfr_exp(r21920, r21919, MPFR_RNDN);
        mpfr_add(r21921, r21909, r21920, MPFR_RNDN);
        mpfr_div(r21922, r21908, r21921, MPFR_RNDN);
        mpfr_set_d(r21923, NaChar, MPFR_RNDN);
        mpfr_set_d(r21924, Ev, MPFR_RNDN);
        mpfr_add(r21925, r21924, r21911, MPFR_RNDN);
        mpfr_set_d(r21926, EAccept, MPFR_RNDN);
        mpfr_add(r21927, r21925, r21926, MPFR_RNDN);
        mpfr_neg(r21928, r21915, MPFR_RNDN);
        mpfr_add(r21929, r21927, r21928, MPFR_RNDN);
        mpfr_div(r21930, r21929, r21918, MPFR_RNDN);
        mpfr_exp(r21931, r21930, MPFR_RNDN);
        mpfr_add(r21932, r21909, r21931, MPFR_RNDN);
        mpfr_div(r21933, r21923, r21932, MPFR_RNDN);
        mpfr_add(r21934, r21922, r21933, MPFR_RNDN);
        return mpfr_get_d(r21934, MPFR_RNDN);
}