#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 r16373 = x; float r16374 = y; float r16375 = r16373 * r16374; float r16376 = z; float r16377 = r16375 + r16376; float r16378 = r16377 * r16374; float r16379 = 27464.7644705f; float r16380 = r16378 + r16379; float r16381 = r16380 * r16374; float r16382 = 230661.510616f; float r16383 = r16381 + r16382; float r16384 = r16383 * r16374; float r16385 = t; float r16386 = r16384 + r16385; float r16387 = a; float r16388 = r16374 + r16387; float r16389 = r16388 * r16374; float r16390 = b; float r16391 = r16389 + r16390; float r16392 = r16391 * r16374; float r16393 = c; float r16394 = r16392 + r16393; float r16395 = r16394 * r16374; float r16396 = i; float r16397 = r16395 + r16396; float r16398 = r16386 / r16397; return r16398; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r16399 = x; double r16400 = y; double r16401 = r16399 * r16400; double r16402 = z; double r16403 = r16401 + r16402; double r16404 = r16403 * r16400; double r16405 = 27464.7644705; double r16406 = r16404 + r16405; double r16407 = r16406 * r16400; double r16408 = 230661.510616; double r16409 = r16407 + r16408; double r16410 = r16409 * r16400; double r16411 = t; double r16412 = r16410 + r16411; double r16413 = a; double r16414 = r16400 + r16413; double r16415 = r16414 * r16400; double r16416 = b; double r16417 = r16415 + r16416; double r16418 = r16417 * r16400; double r16419 = c; double r16420 = r16418 + r16419; double r16421 = r16420 * r16400; double r16422 = i; double r16423 = r16421 + r16422; double r16424 = r16412 / r16423; return r16424; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r16425 = 1.0f; float r16426 = y; float r16427 = r16426 * r16426; float r16428 = a; float r16429 = r16426 + r16428; float r16430 = b; float r16431 = fma(r16429, r16426, r16430); float r16432 = c; float r16433 = i; float r16434 = fma(r16426, r16432, r16433); float r16435 = fma(r16427, r16431, r16434); float r16436 = x; float r16437 = z; float r16438 = fma(r16426, r16436, r16437); float r16439 = 27464.7644705f; float r16440 = 230661.510616f; float r16441 = fma(r16426, r16439, r16440); float r16442 = fma(r16438, r16427, r16441); float r16443 = t; float r16444 = fma(r16442, r16426, r16443); float r16445 = r16435 / r16444; float r16446 = r16425 / r16445; return r16446; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r16447 = 1.0; double r16448 = y; double r16449 = r16448 * r16448; double r16450 = a; double r16451 = r16448 + r16450; double r16452 = b; double r16453 = fma(r16451, r16448, r16452); double r16454 = c; double r16455 = i; double r16456 = fma(r16448, r16454, r16455); double r16457 = fma(r16449, r16453, r16456); double r16458 = x; double r16459 = z; double r16460 = fma(r16448, r16458, r16459); double r16461 = 27464.7644705; double r16462 = 230661.510616; double r16463 = fma(r16448, r16461, r16462); double r16464 = fma(r16460, r16449, r16463); double r16465 = t; double r16466 = fma(r16464, r16448, r16465); double r16467 = r16457 / r16466; double r16468 = r16447 / r16467; return r16468; } 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 r16469, r16470, r16471, r16472, r16473, r16474, r16475, r16476, r16477, r16478, r16479, r16480, r16481, r16482, r16483, r16484, r16485, r16486, r16487, r16488, r16489, r16490, r16491, r16492, r16493, r16494; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16469); mpfr_init(r16470); mpfr_init(r16471); mpfr_init(r16472); mpfr_init(r16473); mpfr_init(r16474); mpfr_init_set_str(r16475, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r16476); mpfr_init(r16477); mpfr_init_set_str(r16478, "230661.510616", 10, MPFR_RNDN); mpfr_init(r16479); mpfr_init(r16480); mpfr_init(r16481); mpfr_init(r16482); mpfr_init(r16483); mpfr_init(r16484); mpfr_init(r16485); mpfr_init(r16486); mpfr_init(r16487); mpfr_init(r16488); mpfr_init(r16489); mpfr_init(r16490); mpfr_init(r16491); mpfr_init(r16492); mpfr_init(r16493); mpfr_init(r16494); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r16469, x, MPFR_RNDN); mpfr_set_d(r16470, y, MPFR_RNDN); mpfr_mul(r16471, r16469, r16470, MPFR_RNDN); mpfr_set_d(r16472, z, MPFR_RNDN); mpfr_add(r16473, r16471, r16472, MPFR_RNDN); mpfr_mul(r16474, r16473, r16470, MPFR_RNDN); ; mpfr_add(r16476, r16474, r16475, MPFR_RNDN); mpfr_mul(r16477, r16476, r16470, MPFR_RNDN); ; mpfr_add(r16479, r16477, r16478, MPFR_RNDN); mpfr_mul(r16480, r16479, r16470, MPFR_RNDN); mpfr_set_d(r16481, t, MPFR_RNDN); mpfr_add(r16482, r16480, r16481, MPFR_RNDN); mpfr_set_d(r16483, a, MPFR_RNDN); mpfr_add(r16484, r16470, r16483, MPFR_RNDN); mpfr_mul(r16485, r16484, r16470, MPFR_RNDN); mpfr_set_d(r16486, b, MPFR_RNDN); mpfr_add(r16487, r16485, r16486, MPFR_RNDN); mpfr_mul(r16488, r16487, r16470, MPFR_RNDN); mpfr_set_d(r16489, c, MPFR_RNDN); mpfr_add(r16490, r16488, r16489, MPFR_RNDN); mpfr_mul(r16491, r16490, r16470, MPFR_RNDN); mpfr_set_d(r16492, i, MPFR_RNDN); mpfr_add(r16493, r16491, r16492, MPFR_RNDN); mpfr_div(r16494, r16482, r16493, MPFR_RNDN); return mpfr_get_d(r16494, MPFR_RNDN); } static mpfr_t r16495, r16496, r16497, r16498, r16499, r16500, r16501, r16502, r16503, r16504, r16505, r16506, r16507, r16508, r16509, r16510, r16511, r16512, r16513, r16514, r16515, r16516; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r16495, "1", 10, MPFR_RNDN); mpfr_init(r16496); mpfr_init(r16497); mpfr_init(r16498); mpfr_init(r16499); mpfr_init(r16500); mpfr_init(r16501); mpfr_init(r16502); mpfr_init(r16503); mpfr_init(r16504); mpfr_init(r16505); mpfr_init(r16506); mpfr_init(r16507); mpfr_init(r16508); mpfr_init_set_str(r16509, "27464.7644705", 10, MPFR_RNDN); mpfr_init_set_str(r16510, "230661.510616", 10, MPFR_RNDN); mpfr_init(r16511); mpfr_init(r16512); mpfr_init(r16513); mpfr_init(r16514); mpfr_init(r16515); mpfr_init(r16516); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { ; mpfr_set_d(r16496, y, MPFR_RNDN); mpfr_sqr(r16497, r16496, MPFR_RNDN); mpfr_set_d(r16498, a, MPFR_RNDN); mpfr_add(r16499, r16496, r16498, MPFR_RNDN); mpfr_set_d(r16500, b, MPFR_RNDN); mpfr_fma(r16501, r16499, r16496, r16500, MPFR_RNDN); mpfr_set_d(r16502, c, MPFR_RNDN); mpfr_set_d(r16503, i, MPFR_RNDN); mpfr_fma(r16504, r16496, r16502, r16503, MPFR_RNDN); mpfr_fma(r16505, r16497, r16501, r16504, MPFR_RNDN); mpfr_set_d(r16506, x, MPFR_RNDN); mpfr_set_d(r16507, z, MPFR_RNDN); mpfr_fma(r16508, r16496, r16506, r16507, MPFR_RNDN); ; ; mpfr_fma(r16511, r16496, r16509, r16510, MPFR_RNDN); mpfr_fma(r16512, r16508, r16497, r16511, MPFR_RNDN); mpfr_set_d(r16513, t, MPFR_RNDN); mpfr_fma(r16514, r16512, r16496, r16513, MPFR_RNDN); mpfr_div(r16515, r16505, r16514, MPFR_RNDN); mpfr_div(r16516, r16495, r16515, MPFR_RNDN); return mpfr_get_d(r16516, MPFR_RNDN); } static mpfr_t r16517, r16518, r16519, r16520, r16521, r16522, r16523, r16524, r16525, r16526, r16527, r16528, r16529, r16530, r16531, r16532, r16533, r16534, r16535, r16536, r16537, r16538; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r16517, "1", 10, MPFR_RNDN); mpfr_init(r16518); mpfr_init(r16519); mpfr_init(r16520); mpfr_init(r16521); mpfr_init(r16522); mpfr_init(r16523); mpfr_init(r16524); mpfr_init(r16525); mpfr_init(r16526); mpfr_init(r16527); mpfr_init(r16528); mpfr_init(r16529); mpfr_init(r16530); mpfr_init_set_str(r16531, "27464.7644705", 10, MPFR_RNDN); mpfr_init_set_str(r16532, "230661.510616", 10, MPFR_RNDN); mpfr_init(r16533); mpfr_init(r16534); mpfr_init(r16535); mpfr_init(r16536); mpfr_init(r16537); mpfr_init(r16538); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { ; mpfr_set_d(r16518, y, MPFR_RNDN); mpfr_sqr(r16519, r16518, MPFR_RNDN); mpfr_set_d(r16520, a, MPFR_RNDN); mpfr_add(r16521, r16518, r16520, MPFR_RNDN); mpfr_set_d(r16522, b, MPFR_RNDN); mpfr_fma(r16523, r16521, r16518, r16522, MPFR_RNDN); mpfr_set_d(r16524, c, MPFR_RNDN); mpfr_set_d(r16525, i, MPFR_RNDN); mpfr_fma(r16526, r16518, r16524, r16525, MPFR_RNDN); mpfr_fma(r16527, r16519, r16523, r16526, MPFR_RNDN); mpfr_set_d(r16528, x, MPFR_RNDN); mpfr_set_d(r16529, z, MPFR_RNDN); mpfr_fma(r16530, r16518, r16528, r16529, MPFR_RNDN); ; ; mpfr_fma(r16533, r16518, r16531, r16532, MPFR_RNDN); mpfr_fma(r16534, r16530, r16519, r16533, MPFR_RNDN); mpfr_set_d(r16535, t, MPFR_RNDN); mpfr_fma(r16536, r16534, r16518, r16535, MPFR_RNDN); mpfr_div(r16537, r16527, r16536, MPFR_RNDN); mpfr_div(r16538, r16517, r16537, MPFR_RNDN); return mpfr_get_d(r16538, MPFR_RNDN); }