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