#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 r20726 = x; float r20727 = y; float r20728 = r20726 + r20727; float r20729 = log(r20728); float r20730 = z; float r20731 = log(r20730); float r20732 = r20729 + r20731; float r20733 = t; float r20734 = r20732 - r20733; float r20735 = a; float r20736 = 0.5; float r20737 = r20735 - r20736; float r20738 = log(r20733); float r20739 = r20737 * r20738; float r20740 = r20734 + r20739; return r20740; } double f_id(double x, double y, double z, double t, double a) { double r20741 = x; double r20742 = y; double r20743 = r20741 + r20742; double r20744 = log(r20743); double r20745 = z; double r20746 = log(r20745); double r20747 = r20744 + r20746; double r20748 = t; double r20749 = r20747 - r20748; double r20750 = a; double r20751 = 0.5; double r20752 = r20750 - r20751; double r20753 = log(r20748); double r20754 = r20752 * r20753; double r20755 = r20749 + r20754; return r20755; } double f_of(float x, float y, float z, float t, float a) { float r20756 = x; float r20757 = y; float r20758 = r20756 + r20757; float r20759 = log(r20758); float r20760 = 3; float r20761 = pow(r20759, r20760); float r20762 = z; float r20763 = log(r20762); float r20764 = pow(r20763, r20760); float r20765 = r20761 + r20764; float r20766 = r20757 + r20756; float r20767 = log(r20766); float r20768 = r20763 - r20767; float r20769 = r20768 * r20763; float r20770 = r20767 * r20767; float r20771 = r20769 + r20770; float r20772 = r20765 / r20771; float r20773 = t; float r20774 = r20772 - r20773; float r20775 = a; float r20776 = 0.5; float r20777 = r20775 - r20776; float r20778 = log(r20773); float r20779 = r20777 * r20778; float r20780 = r20774 + r20779; return r20780; } double f_od(double x, double y, double z, double t, double a) { double r20781 = x; double r20782 = y; double r20783 = r20781 + r20782; double r20784 = log(r20783); double r20785 = 3; double r20786 = pow(r20784, r20785); double r20787 = z; double r20788 = log(r20787); double r20789 = pow(r20788, r20785); double r20790 = r20786 + r20789; double r20791 = r20782 + r20781; double r20792 = log(r20791); double r20793 = r20788 - r20792; double r20794 = r20793 * r20788; double r20795 = r20792 * r20792; double r20796 = r20794 + r20795; double r20797 = r20790 / r20796; double r20798 = t; double r20799 = r20797 - r20798; double r20800 = a; double r20801 = 0.5; double r20802 = r20800 - r20801; double r20803 = log(r20798); double r20804 = r20802 * r20803; double r20805 = r20799 + r20804; return r20805; } 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 r20806, r20807, r20808, r20809, r20810, r20811, r20812, r20813, r20814, r20815, r20816, r20817, r20818, r20819, r20820; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r20806); mpfr_init(r20807); mpfr_init(r20808); mpfr_init(r20809); mpfr_init(r20810); mpfr_init(r20811); mpfr_init(r20812); mpfr_init(r20813); mpfr_init(r20814); mpfr_init(r20815); mpfr_init_set_str(r20816, "0.5", 10, MPFR_RNDN); mpfr_init(r20817); mpfr_init(r20818); mpfr_init(r20819); mpfr_init(r20820); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r20806, x, MPFR_RNDN); mpfr_set_d(r20807, y, MPFR_RNDN); mpfr_add(r20808, r20806, r20807, MPFR_RNDN); mpfr_log(r20809, r20808, MPFR_RNDN); mpfr_set_d(r20810, z, MPFR_RNDN); mpfr_log(r20811, r20810, MPFR_RNDN); mpfr_add(r20812, r20809, r20811, MPFR_RNDN); mpfr_set_d(r20813, t, MPFR_RNDN); mpfr_sub(r20814, r20812, r20813, MPFR_RNDN); mpfr_set_d(r20815, a, MPFR_RNDN); ; mpfr_sub(r20817, r20815, r20816, MPFR_RNDN); mpfr_log(r20818, r20813, MPFR_RNDN); mpfr_mul(r20819, r20817, r20818, MPFR_RNDN); mpfr_add(r20820, r20814, r20819, MPFR_RNDN); return mpfr_get_d(r20820, MPFR_RNDN); } static mpfr_t r20821, r20822, r20823, r20824, r20825, r20826, r20827, r20828, r20829, r20830, r20831, r20832, r20833, r20834, r20835, r20836, r20837, r20838, r20839, r20840, r20841, r20842, r20843, r20844, r20845; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r20821); mpfr_init(r20822); mpfr_init(r20823); mpfr_init(r20824); mpfr_init_set_str(r20825, "3", 10, MPFR_RNDN); mpfr_init(r20826); mpfr_init(r20827); mpfr_init(r20828); mpfr_init(r20829); mpfr_init(r20830); mpfr_init(r20831); mpfr_init(r20832); mpfr_init(r20833); mpfr_init(r20834); mpfr_init(r20835); mpfr_init(r20836); mpfr_init(r20837); mpfr_init(r20838); mpfr_init(r20839); mpfr_init(r20840); mpfr_init_set_str(r20841, "0.5", 10, MPFR_RNDN); mpfr_init(r20842); mpfr_init(r20843); mpfr_init(r20844); mpfr_init(r20845); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r20821, x, MPFR_RNDN); mpfr_set_d(r20822, y, MPFR_RNDN); mpfr_add(r20823, r20821, r20822, MPFR_RNDN); mpfr_log(r20824, r20823, MPFR_RNDN); ; mpfr_pow(r20826, r20824, r20825, MPFR_RNDN); mpfr_set_d(r20827, z, MPFR_RNDN); mpfr_log(r20828, r20827, MPFR_RNDN); mpfr_pow(r20829, r20828, r20825, MPFR_RNDN); mpfr_add(r20830, r20826, r20829, MPFR_RNDN); mpfr_add(r20831, r20822, r20821, MPFR_RNDN); mpfr_log(r20832, r20831, MPFR_RNDN); mpfr_sub(r20833, r20828, r20832, MPFR_RNDN); mpfr_mul(r20834, r20833, r20828, MPFR_RNDN); mpfr_mul(r20835, r20832, r20832, MPFR_RNDN); mpfr_add(r20836, r20834, r20835, MPFR_RNDN); mpfr_div(r20837, r20830, r20836, MPFR_RNDN); mpfr_set_d(r20838, t, MPFR_RNDN); mpfr_sub(r20839, r20837, r20838, MPFR_RNDN); mpfr_set_d(r20840, a, MPFR_RNDN); ; mpfr_sub(r20842, r20840, r20841, MPFR_RNDN); mpfr_log(r20843, r20838, MPFR_RNDN); mpfr_mul(r20844, r20842, r20843, MPFR_RNDN); mpfr_add(r20845, r20839, r20844, MPFR_RNDN); return mpfr_get_d(r20845, MPFR_RNDN); } static mpfr_t r20846, r20847, r20848, r20849, r20850, r20851, r20852, r20853, r20854, r20855, r20856, r20857, r20858, r20859, r20860, r20861, r20862, r20863, r20864, r20865, r20866, r20867, r20868, r20869, r20870; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r20846); mpfr_init(r20847); mpfr_init(r20848); mpfr_init(r20849); mpfr_init_set_str(r20850, "3", 10, MPFR_RNDN); mpfr_init(r20851); mpfr_init(r20852); mpfr_init(r20853); mpfr_init(r20854); mpfr_init(r20855); mpfr_init(r20856); mpfr_init(r20857); mpfr_init(r20858); mpfr_init(r20859); mpfr_init(r20860); mpfr_init(r20861); mpfr_init(r20862); mpfr_init(r20863); mpfr_init(r20864); mpfr_init(r20865); mpfr_init_set_str(r20866, "0.5", 10, MPFR_RNDN); mpfr_init(r20867); mpfr_init(r20868); mpfr_init(r20869); mpfr_init(r20870); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r20846, x, MPFR_RNDN); mpfr_set_d(r20847, y, MPFR_RNDN); mpfr_add(r20848, r20846, r20847, MPFR_RNDN); mpfr_log(r20849, r20848, MPFR_RNDN); ; mpfr_pow(r20851, r20849, r20850, MPFR_RNDN); mpfr_set_d(r20852, z, MPFR_RNDN); mpfr_log(r20853, r20852, MPFR_RNDN); mpfr_pow(r20854, r20853, r20850, MPFR_RNDN); mpfr_add(r20855, r20851, r20854, MPFR_RNDN); mpfr_add(r20856, r20847, r20846, MPFR_RNDN); mpfr_log(r20857, r20856, MPFR_RNDN); mpfr_sub(r20858, r20853, r20857, MPFR_RNDN); mpfr_mul(r20859, r20858, r20853, MPFR_RNDN); mpfr_mul(r20860, r20857, r20857, MPFR_RNDN); mpfr_add(r20861, r20859, r20860, MPFR_RNDN); mpfr_div(r20862, r20855, r20861, MPFR_RNDN); mpfr_set_d(r20863, t, MPFR_RNDN); mpfr_sub(r20864, r20862, r20863, MPFR_RNDN); mpfr_set_d(r20865, a, MPFR_RNDN); ; mpfr_sub(r20867, r20865, r20866, MPFR_RNDN); mpfr_log(r20868, r20863, MPFR_RNDN); mpfr_mul(r20869, r20867, r20868, MPFR_RNDN); mpfr_add(r20870, r20864, r20869, MPFR_RNDN); return mpfr_get_d(r20870, MPFR_RNDN); }