#include #include #include #include #include char *name = "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2"; double f_if(float x, float y, float z, float t, float a) { float r22762 = x; float r22763 = y; float r22764 = r22762 + r22763; float r22765 = log(r22764); float r22766 = z; float r22767 = log(r22766); float r22768 = r22765 + r22767; float r22769 = t; float r22770 = r22768 - r22769; float r22771 = a; float r22772 = 0.5; float r22773 = r22771 - r22772; float r22774 = log(r22769); float r22775 = r22773 * r22774; float r22776 = r22770 + r22775; return r22776; } double f_id(double x, double y, double z, double t, double a) { double r22777 = x; double r22778 = y; double r22779 = r22777 + r22778; double r22780 = log(r22779); double r22781 = z; double r22782 = log(r22781); double r22783 = r22780 + r22782; double r22784 = t; double r22785 = r22783 - r22784; double r22786 = a; double r22787 = 0.5; double r22788 = r22786 - r22787; double r22789 = log(r22784); double r22790 = r22788 * r22789; double r22791 = r22785 + r22790; return r22791; } double f_of(float x, float y, float z, float t, float a) { float r22792 = x; float r22793 = y; float r22794 = r22792 + r22793; float r22795 = log(r22794); float r22796 = z; float r22797 = log(r22796); float r22798 = t; float r22799 = r22797 - r22798; float r22800 = a; float r22801 = 0.5; float r22802 = r22800 - r22801; float r22803 = log(r22798); float r22804 = r22802 * r22803; float r22805 = r22799 + r22804; float r22806 = r22795 + r22805; return r22806; } double f_od(double x, double y, double z, double t, double a) { double r22807 = x; double r22808 = y; double r22809 = r22807 + r22808; double r22810 = log(r22809); double r22811 = z; double r22812 = log(r22811); double r22813 = t; double r22814 = r22812 - r22813; double r22815 = a; double r22816 = 0.5; double r22817 = r22815 - r22816; double r22818 = log(r22813); double r22819 = r22817 * r22818; double r22820 = r22814 + r22819; double r22821 = r22810 + r22820; return r22821; } 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 r22822, r22823, r22824, r22825, r22826, r22827, r22828, r22829, r22830, r22831, r22832, r22833, r22834, r22835, r22836; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r22822); mpfr_init(r22823); mpfr_init(r22824); mpfr_init(r22825); mpfr_init(r22826); mpfr_init(r22827); mpfr_init(r22828); mpfr_init(r22829); mpfr_init(r22830); mpfr_init(r22831); mpfr_init_set_str(r22832, "0.5", 10, MPFR_RNDN); mpfr_init(r22833); mpfr_init(r22834); mpfr_init(r22835); mpfr_init(r22836); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r22822, x, MPFR_RNDN); mpfr_set_d(r22823, y, MPFR_RNDN); mpfr_add(r22824, r22822, r22823, MPFR_RNDN); mpfr_log(r22825, r22824, MPFR_RNDN); mpfr_set_d(r22826, z, MPFR_RNDN); mpfr_log(r22827, r22826, MPFR_RNDN); mpfr_add(r22828, r22825, r22827, MPFR_RNDN); mpfr_set_d(r22829, t, MPFR_RNDN); mpfr_sub(r22830, r22828, r22829, MPFR_RNDN); mpfr_set_d(r22831, a, MPFR_RNDN); ; mpfr_sub(r22833, r22831, r22832, MPFR_RNDN); mpfr_log(r22834, r22829, MPFR_RNDN); mpfr_mul(r22835, r22833, r22834, MPFR_RNDN); mpfr_add(r22836, r22830, r22835, MPFR_RNDN); return mpfr_get_d(r22836, MPFR_RNDN); } static mpfr_t r22837, r22838, r22839, r22840, r22841, r22842, r22843, r22844, r22845, r22846, r22847, r22848, r22849, r22850, r22851; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r22837); mpfr_init(r22838); mpfr_init(r22839); mpfr_init(r22840); mpfr_init(r22841); mpfr_init(r22842); mpfr_init(r22843); mpfr_init(r22844); mpfr_init(r22845); mpfr_init_set_str(r22846, "0.5", 10, MPFR_RNDN); mpfr_init(r22847); mpfr_init(r22848); mpfr_init(r22849); mpfr_init(r22850); mpfr_init(r22851); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r22837, x, MPFR_RNDN); mpfr_set_d(r22838, y, MPFR_RNDN); mpfr_add(r22839, r22837, r22838, MPFR_RNDN); mpfr_log(r22840, r22839, MPFR_RNDN); mpfr_set_d(r22841, z, MPFR_RNDN); mpfr_log(r22842, r22841, MPFR_RNDN); mpfr_set_d(r22843, t, MPFR_RNDN); mpfr_sub(r22844, r22842, r22843, MPFR_RNDN); mpfr_set_d(r22845, a, MPFR_RNDN); ; mpfr_sub(r22847, r22845, r22846, MPFR_RNDN); mpfr_log(r22848, r22843, MPFR_RNDN); mpfr_mul(r22849, r22847, r22848, MPFR_RNDN); mpfr_add(r22850, r22844, r22849, MPFR_RNDN); mpfr_add(r22851, r22840, r22850, MPFR_RNDN); return mpfr_get_d(r22851, MPFR_RNDN); } static mpfr_t r22852, r22853, r22854, r22855, r22856, r22857, r22858, r22859, r22860, r22861, r22862, r22863, r22864, r22865, r22866; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r22852); mpfr_init(r22853); mpfr_init(r22854); mpfr_init(r22855); mpfr_init(r22856); mpfr_init(r22857); mpfr_init(r22858); mpfr_init(r22859); mpfr_init(r22860); mpfr_init_set_str(r22861, "0.5", 10, MPFR_RNDN); mpfr_init(r22862); mpfr_init(r22863); mpfr_init(r22864); mpfr_init(r22865); mpfr_init(r22866); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r22852, x, MPFR_RNDN); mpfr_set_d(r22853, y, MPFR_RNDN); mpfr_add(r22854, r22852, r22853, MPFR_RNDN); mpfr_log(r22855, r22854, MPFR_RNDN); mpfr_set_d(r22856, z, MPFR_RNDN); mpfr_log(r22857, r22856, MPFR_RNDN); mpfr_set_d(r22858, t, MPFR_RNDN); mpfr_sub(r22859, r22857, r22858, MPFR_RNDN); mpfr_set_d(r22860, a, MPFR_RNDN); ; mpfr_sub(r22862, r22860, r22861, MPFR_RNDN); mpfr_log(r22863, r22858, MPFR_RNDN); mpfr_mul(r22864, r22862, r22863, MPFR_RNDN); mpfr_add(r22865, r22859, r22864, MPFR_RNDN); mpfr_add(r22866, r22855, r22865, MPFR_RNDN); return mpfr_get_d(r22866, MPFR_RNDN); }