#include #include #include #include #include char *name = "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i) { float r25216 = x; float r25217 = y; float r25218 = r25216 * r25217; float r25219 = z; float r25220 = r25218 + r25219; float r25221 = r25220 * r25217; float r25222 = 27464.7644705; float r25223 = r25221 + r25222; float r25224 = r25223 * r25217; float r25225 = 230661.510616; float r25226 = r25224 + r25225; float r25227 = r25226 * r25217; float r25228 = t; float r25229 = r25227 + r25228; float r25230 = a; float r25231 = r25217 + r25230; float r25232 = r25231 * r25217; float r25233 = b; float r25234 = r25232 + r25233; float r25235 = r25234 * r25217; float r25236 = c; float r25237 = r25235 + r25236; float r25238 = r25237 * r25217; float r25239 = i; float r25240 = r25238 + r25239; float r25241 = r25229 / r25240; return r25241; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r25242 = x; double r25243 = y; double r25244 = r25242 * r25243; double r25245 = z; double r25246 = r25244 + r25245; double r25247 = r25246 * r25243; double r25248 = 27464.7644705; double r25249 = r25247 + r25248; double r25250 = r25249 * r25243; double r25251 = 230661.510616; double r25252 = r25250 + r25251; double r25253 = r25252 * r25243; double r25254 = t; double r25255 = r25253 + r25254; double r25256 = a; double r25257 = r25243 + r25256; double r25258 = r25257 * r25243; double r25259 = b; double r25260 = r25258 + r25259; double r25261 = r25260 * r25243; double r25262 = c; double r25263 = r25261 + r25262; double r25264 = r25263 * r25243; double r25265 = i; double r25266 = r25264 + r25265; double r25267 = r25255 / r25266; return r25267; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r25268 = x; float r25269 = y; float r25270 = r25268 * r25269; float r25271 = z; float r25272 = r25270 + r25271; float r25273 = r25272 * r25269; float r25274 = 27464.7644705; float r25275 = r25273 + r25274; float r25276 = r25275 * r25269; float r25277 = 230661.510616; float r25278 = r25276 + r25277; float r25279 = r25278 * r25269; float r25280 = t; float r25281 = r25279 + r25280; float r25282 = 1; float r25283 = a; float r25284 = r25269 + r25283; float r25285 = r25284 * r25269; float r25286 = b; float r25287 = r25285 + r25286; float r25288 = r25287 * r25269; float r25289 = c; float r25290 = r25288 + r25289; float r25291 = r25290 * r25269; float r25292 = i; float r25293 = r25291 + r25292; float r25294 = r25282 / r25293; float r25295 = r25281 * r25294; return r25295; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r25296 = x; double r25297 = y; double r25298 = r25296 * r25297; double r25299 = z; double r25300 = r25298 + r25299; double r25301 = r25300 * r25297; double r25302 = 27464.7644705; double r25303 = r25301 + r25302; double r25304 = r25303 * r25297; double r25305 = 230661.510616; double r25306 = r25304 + r25305; double r25307 = r25306 * r25297; double r25308 = t; double r25309 = r25307 + r25308; double r25310 = 1; double r25311 = a; double r25312 = r25297 + r25311; double r25313 = r25312 * r25297; double r25314 = b; double r25315 = r25313 + r25314; double r25316 = r25315 * r25297; double r25317 = c; double r25318 = r25316 + r25317; double r25319 = r25318 * r25297; double r25320 = i; double r25321 = r25319 + r25320; double r25322 = r25310 / r25321; double r25323 = r25309 * r25322; return r25323; } 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 r25324, r25325, r25326, r25327, r25328, r25329, r25330, r25331, r25332, r25333, r25334, r25335, r25336, r25337, r25338, r25339, r25340, r25341, r25342, r25343, r25344, r25345, r25346, r25347, r25348, r25349; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25324); mpfr_init(r25325); mpfr_init(r25326); mpfr_init(r25327); mpfr_init(r25328); mpfr_init(r25329); mpfr_init_set_str(r25330, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r25331); mpfr_init(r25332); mpfr_init_set_str(r25333, "230661.510616", 10, MPFR_RNDN); mpfr_init(r25334); mpfr_init(r25335); mpfr_init(r25336); mpfr_init(r25337); mpfr_init(r25338); mpfr_init(r25339); mpfr_init(r25340); mpfr_init(r25341); mpfr_init(r25342); mpfr_init(r25343); mpfr_init(r25344); mpfr_init(r25345); mpfr_init(r25346); mpfr_init(r25347); mpfr_init(r25348); mpfr_init(r25349); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r25324, x, MPFR_RNDN); mpfr_set_d(r25325, y, MPFR_RNDN); mpfr_mul(r25326, r25324, r25325, MPFR_RNDN); mpfr_set_d(r25327, z, MPFR_RNDN); mpfr_add(r25328, r25326, r25327, MPFR_RNDN); mpfr_mul(r25329, r25328, r25325, MPFR_RNDN); ; mpfr_add(r25331, r25329, r25330, MPFR_RNDN); mpfr_mul(r25332, r25331, r25325, MPFR_RNDN); ; mpfr_add(r25334, r25332, r25333, MPFR_RNDN); mpfr_mul(r25335, r25334, r25325, MPFR_RNDN); mpfr_set_d(r25336, t, MPFR_RNDN); mpfr_add(r25337, r25335, r25336, MPFR_RNDN); mpfr_set_d(r25338, a, MPFR_RNDN); mpfr_add(r25339, r25325, r25338, MPFR_RNDN); mpfr_mul(r25340, r25339, r25325, MPFR_RNDN); mpfr_set_d(r25341, b, MPFR_RNDN); mpfr_add(r25342, r25340, r25341, MPFR_RNDN); mpfr_mul(r25343, r25342, r25325, MPFR_RNDN); mpfr_set_d(r25344, c, MPFR_RNDN); mpfr_add(r25345, r25343, r25344, MPFR_RNDN); mpfr_mul(r25346, r25345, r25325, MPFR_RNDN); mpfr_set_d(r25347, i, MPFR_RNDN); mpfr_add(r25348, r25346, r25347, MPFR_RNDN); mpfr_div(r25349, r25337, r25348, MPFR_RNDN); return mpfr_get_d(r25349, MPFR_RNDN); } static mpfr_t r25350, r25351, r25352, r25353, r25354, r25355, r25356, r25357, r25358, r25359, r25360, r25361, r25362, r25363, r25364, r25365, r25366, r25367, r25368, r25369, r25370, r25371, r25372, r25373, r25374, r25375, r25376, r25377; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25350); mpfr_init(r25351); mpfr_init(r25352); mpfr_init(r25353); mpfr_init(r25354); mpfr_init(r25355); mpfr_init_set_str(r25356, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r25357); mpfr_init(r25358); mpfr_init_set_str(r25359, "230661.510616", 10, MPFR_RNDN); mpfr_init(r25360); mpfr_init(r25361); mpfr_init(r25362); mpfr_init(r25363); mpfr_init_set_str(r25364, "1", 10, MPFR_RNDN); mpfr_init(r25365); mpfr_init(r25366); mpfr_init(r25367); mpfr_init(r25368); mpfr_init(r25369); mpfr_init(r25370); mpfr_init(r25371); mpfr_init(r25372); mpfr_init(r25373); mpfr_init(r25374); mpfr_init(r25375); mpfr_init(r25376); mpfr_init(r25377); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r25350, x, MPFR_RNDN); mpfr_set_d(r25351, y, MPFR_RNDN); mpfr_mul(r25352, r25350, r25351, MPFR_RNDN); mpfr_set_d(r25353, z, MPFR_RNDN); mpfr_add(r25354, r25352, r25353, MPFR_RNDN); mpfr_mul(r25355, r25354, r25351, MPFR_RNDN); ; mpfr_add(r25357, r25355, r25356, MPFR_RNDN); mpfr_mul(r25358, r25357, r25351, MPFR_RNDN); ; mpfr_add(r25360, r25358, r25359, MPFR_RNDN); mpfr_mul(r25361, r25360, r25351, MPFR_RNDN); mpfr_set_d(r25362, t, MPFR_RNDN); mpfr_add(r25363, r25361, r25362, MPFR_RNDN); ; mpfr_set_d(r25365, a, MPFR_RNDN); mpfr_add(r25366, r25351, r25365, MPFR_RNDN); mpfr_mul(r25367, r25366, r25351, MPFR_RNDN); mpfr_set_d(r25368, b, MPFR_RNDN); mpfr_add(r25369, r25367, r25368, MPFR_RNDN); mpfr_mul(r25370, r25369, r25351, MPFR_RNDN); mpfr_set_d(r25371, c, MPFR_RNDN); mpfr_add(r25372, r25370, r25371, MPFR_RNDN); mpfr_mul(r25373, r25372, r25351, MPFR_RNDN); mpfr_set_d(r25374, i, MPFR_RNDN); mpfr_add(r25375, r25373, r25374, MPFR_RNDN); mpfr_div(r25376, r25364, r25375, MPFR_RNDN); mpfr_mul(r25377, r25363, r25376, MPFR_RNDN); return mpfr_get_d(r25377, MPFR_RNDN); } static mpfr_t r25378, r25379, r25380, r25381, r25382, r25383, r25384, r25385, r25386, r25387, r25388, r25389, r25390, r25391, r25392, r25393, r25394, r25395, r25396, r25397, r25398, r25399, r25400, r25401, r25402, r25403, r25404, r25405; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25378); mpfr_init(r25379); mpfr_init(r25380); mpfr_init(r25381); mpfr_init(r25382); mpfr_init(r25383); mpfr_init_set_str(r25384, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r25385); mpfr_init(r25386); mpfr_init_set_str(r25387, "230661.510616", 10, MPFR_RNDN); mpfr_init(r25388); mpfr_init(r25389); mpfr_init(r25390); mpfr_init(r25391); mpfr_init_set_str(r25392, "1", 10, MPFR_RNDN); mpfr_init(r25393); mpfr_init(r25394); mpfr_init(r25395); mpfr_init(r25396); mpfr_init(r25397); mpfr_init(r25398); mpfr_init(r25399); mpfr_init(r25400); mpfr_init(r25401); mpfr_init(r25402); mpfr_init(r25403); mpfr_init(r25404); mpfr_init(r25405); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r25378, x, MPFR_RNDN); mpfr_set_d(r25379, y, MPFR_RNDN); mpfr_mul(r25380, r25378, r25379, MPFR_RNDN); mpfr_set_d(r25381, z, MPFR_RNDN); mpfr_add(r25382, r25380, r25381, MPFR_RNDN); mpfr_mul(r25383, r25382, r25379, MPFR_RNDN); ; mpfr_add(r25385, r25383, r25384, MPFR_RNDN); mpfr_mul(r25386, r25385, r25379, MPFR_RNDN); ; mpfr_add(r25388, r25386, r25387, MPFR_RNDN); mpfr_mul(r25389, r25388, r25379, MPFR_RNDN); mpfr_set_d(r25390, t, MPFR_RNDN); mpfr_add(r25391, r25389, r25390, MPFR_RNDN); ; mpfr_set_d(r25393, a, MPFR_RNDN); mpfr_add(r25394, r25379, r25393, MPFR_RNDN); mpfr_mul(r25395, r25394, r25379, MPFR_RNDN); mpfr_set_d(r25396, b, MPFR_RNDN); mpfr_add(r25397, r25395, r25396, MPFR_RNDN); mpfr_mul(r25398, r25397, r25379, MPFR_RNDN); mpfr_set_d(r25399, c, MPFR_RNDN); mpfr_add(r25400, r25398, r25399, MPFR_RNDN); mpfr_mul(r25401, r25400, r25379, MPFR_RNDN); mpfr_set_d(r25402, i, MPFR_RNDN); mpfr_add(r25403, r25401, r25402, MPFR_RNDN); mpfr_div(r25404, r25392, r25403, MPFR_RNDN); mpfr_mul(r25405, r25391, r25404, MPFR_RNDN); return mpfr_get_d(r25405, MPFR_RNDN); }