#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 r22786 = x; float r22787 = y; float r22788 = r22786 + r22787; float r22789 = log(r22788); float r22790 = z; float r22791 = log(r22790); float r22792 = r22789 + r22791; float r22793 = t; float r22794 = r22792 - r22793; float r22795 = a; float r22796 = 0.5; float r22797 = r22795 - r22796; float r22798 = log(r22793); float r22799 = r22797 * r22798; float r22800 = r22794 + r22799; return r22800; } double f_id(double x, double y, double z, double t, double a) { double r22801 = x; double r22802 = y; double r22803 = r22801 + r22802; double r22804 = log(r22803); double r22805 = z; double r22806 = log(r22805); double r22807 = r22804 + r22806; double r22808 = t; double r22809 = r22807 - r22808; double r22810 = a; double r22811 = 0.5; double r22812 = r22810 - r22811; double r22813 = log(r22808); double r22814 = r22812 * r22813; double r22815 = r22809 + r22814; return r22815; } double f_of(float x, float y, float z, float t, float a) { float r22816 = x; float r22817 = y; float r22818 = r22816 + r22817; float r22819 = log(r22818); float r22820 = a; float r22821 = 0.5; float r22822 = r22820 - r22821; float r22823 = t; float r22824 = log(r22823); float r22825 = z; float r22826 = log(r22825); float r22827 = fma(r22822, r22824, r22826); float r22828 = r22819 + r22827; float r22829 = r22828 - r22823; float r22830 = 0; float r22831 = r22829 + r22830; return r22831; } double f_od(double x, double y, double z, double t, double a) { double r22832 = x; double r22833 = y; double r22834 = r22832 + r22833; double r22835 = log(r22834); double r22836 = a; double r22837 = 0.5; double r22838 = r22836 - r22837; double r22839 = t; double r22840 = log(r22839); double r22841 = z; double r22842 = log(r22841); double r22843 = fma(r22838, r22840, r22842); double r22844 = r22835 + r22843; double r22845 = r22844 - r22839; double r22846 = 0; double r22847 = r22845 + r22846; return r22847; } 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 r22848, r22849, r22850, r22851, r22852, r22853, r22854, r22855, r22856, r22857, r22858, r22859, r22860, r22861, r22862; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22848); mpfr_init(r22849); mpfr_init(r22850); mpfr_init(r22851); mpfr_init(r22852); mpfr_init(r22853); mpfr_init(r22854); mpfr_init(r22855); mpfr_init(r22856); mpfr_init(r22857); mpfr_init_set_str(r22858, "0.5", 10, MPFR_RNDN); mpfr_init(r22859); mpfr_init(r22860); mpfr_init(r22861); mpfr_init(r22862); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r22848, x, MPFR_RNDN); mpfr_set_d(r22849, y, MPFR_RNDN); mpfr_add(r22850, r22848, r22849, MPFR_RNDN); mpfr_log(r22851, r22850, MPFR_RNDN); mpfr_set_d(r22852, z, MPFR_RNDN); mpfr_log(r22853, r22852, MPFR_RNDN); mpfr_add(r22854, r22851, r22853, MPFR_RNDN); mpfr_set_d(r22855, t, MPFR_RNDN); mpfr_sub(r22856, r22854, r22855, MPFR_RNDN); mpfr_set_d(r22857, a, MPFR_RNDN); ; mpfr_sub(r22859, r22857, r22858, MPFR_RNDN); mpfr_log(r22860, r22855, MPFR_RNDN); mpfr_mul(r22861, r22859, r22860, MPFR_RNDN); mpfr_add(r22862, r22856, r22861, MPFR_RNDN); return mpfr_get_d(r22862, MPFR_RNDN); } static mpfr_t r22863, r22864, r22865, r22866, r22867, r22868, r22869, r22870, r22871, r22872, r22873, r22874, r22875, r22876, r22877, r22878; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22863); mpfr_init(r22864); mpfr_init(r22865); mpfr_init(r22866); mpfr_init(r22867); mpfr_init_set_str(r22868, "0.5", 10, MPFR_RNDN); mpfr_init(r22869); mpfr_init(r22870); mpfr_init(r22871); mpfr_init(r22872); mpfr_init(r22873); mpfr_init(r22874); mpfr_init(r22875); mpfr_init(r22876); mpfr_init_set_str(r22877, "0", 10, MPFR_RNDN); mpfr_init(r22878); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r22863, x, MPFR_RNDN); mpfr_set_d(r22864, y, MPFR_RNDN); mpfr_add(r22865, r22863, r22864, MPFR_RNDN); mpfr_log(r22866, r22865, MPFR_RNDN); mpfr_set_d(r22867, a, MPFR_RNDN); ; mpfr_sub(r22869, r22867, r22868, MPFR_RNDN); mpfr_set_d(r22870, t, MPFR_RNDN); mpfr_log(r22871, r22870, MPFR_RNDN); mpfr_set_d(r22872, z, MPFR_RNDN); mpfr_log(r22873, r22872, MPFR_RNDN); mpfr_fma(r22874, r22869, r22871, r22873, MPFR_RNDN); mpfr_add(r22875, r22866, r22874, MPFR_RNDN); mpfr_sub(r22876, r22875, r22870, MPFR_RNDN); ; mpfr_add(r22878, r22876, r22877, MPFR_RNDN); return mpfr_get_d(r22878, MPFR_RNDN); } static mpfr_t r22879, r22880, r22881, r22882, r22883, r22884, r22885, r22886, r22887, r22888, r22889, r22890, r22891, r22892, r22893, r22894; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22879); mpfr_init(r22880); mpfr_init(r22881); mpfr_init(r22882); mpfr_init(r22883); mpfr_init_set_str(r22884, "0.5", 10, MPFR_RNDN); mpfr_init(r22885); mpfr_init(r22886); mpfr_init(r22887); mpfr_init(r22888); mpfr_init(r22889); mpfr_init(r22890); mpfr_init(r22891); mpfr_init(r22892); mpfr_init_set_str(r22893, "0", 10, MPFR_RNDN); mpfr_init(r22894); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r22879, x, MPFR_RNDN); mpfr_set_d(r22880, y, MPFR_RNDN); mpfr_add(r22881, r22879, r22880, MPFR_RNDN); mpfr_log(r22882, r22881, MPFR_RNDN); mpfr_set_d(r22883, a, MPFR_RNDN); ; mpfr_sub(r22885, r22883, r22884, MPFR_RNDN); mpfr_set_d(r22886, t, MPFR_RNDN); mpfr_log(r22887, r22886, MPFR_RNDN); mpfr_set_d(r22888, z, MPFR_RNDN); mpfr_log(r22889, r22888, MPFR_RNDN); mpfr_fma(r22890, r22885, r22887, r22889, MPFR_RNDN); mpfr_add(r22891, r22882, r22890, MPFR_RNDN); mpfr_sub(r22892, r22891, r22886, MPFR_RNDN); ; mpfr_add(r22894, r22892, r22893, MPFR_RNDN); return mpfr_get_d(r22894, MPFR_RNDN); }