#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 r15918 = x; float r15919 = y; float r15920 = r15918 * r15919; float r15921 = z; float r15922 = r15920 + r15921; float r15923 = r15922 * r15919; float r15924 = 27464.763671875f; float r15925 = r15923 + r15924; float r15926 = r15925 * r15919; float r15927 = 230661.515625f; float r15928 = r15926 + r15927; float r15929 = r15928 * r15919; float r15930 = t; float r15931 = r15929 + r15930; float r15932 = a; float r15933 = r15919 + r15932; float r15934 = r15933 * r15919; float r15935 = b; float r15936 = r15934 + r15935; float r15937 = r15936 * r15919; float r15938 = c; float r15939 = r15937 + r15938; float r15940 = r15939 * r15919; float r15941 = i; float r15942 = r15940 + r15941; float r15943 = r15931 / r15942; return r15943; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r15944 = x; double r15945 = y; double r15946 = r15944 * r15945; double r15947 = z; double r15948 = r15946 + r15947; double r15949 = r15948 * r15945; double r15950 = 27464.763671875; double r15951 = r15949 + r15950; double r15952 = r15951 * r15945; double r15953 = 230661.515625; double r15954 = r15952 + r15953; double r15955 = r15954 * r15945; double r15956 = t; double r15957 = r15955 + r15956; double r15958 = a; double r15959 = r15945 + r15958; double r15960 = r15959 * r15945; double r15961 = b; double r15962 = r15960 + r15961; double r15963 = r15962 * r15945; double r15964 = c; double r15965 = r15963 + r15964; double r15966 = r15965 * r15945; double r15967 = i; double r15968 = r15966 + r15967; double r15969 = r15957 / r15968; return r15969; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r15970 = y; float r15971 = r15970 * r15970; float r15972 = z; float r15973 = r15971 * r15972; float r15974 = 3.0f; float r15975 = pow(r15970, r15974); float r15976 = x; float r15977 = r15975 * r15976; float r15978 = r15973 + r15977; float r15979 = 230661.515625f; float r15980 = r15978 + r15979; float r15981 = r15980 * r15970; float r15982 = t; float r15983 = r15981 + r15982; float r15984 = a; float r15985 = r15970 + r15984; float r15986 = r15985 * r15970; float r15987 = b; float r15988 = r15986 + r15987; float r15989 = r15988 * r15970; float r15990 = 1.0f; float r15991 = pow(r15989, r15990); float r15992 = c; float r15993 = r15991 + r15992; float r15994 = r15993 * r15970; float r15995 = i; float r15996 = r15994 + r15995; float r15997 = r15983 / r15996; return r15997; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r15998 = y; double r15999 = r15998 * r15998; double r16000 = z; double r16001 = r15999 * r16000; double r16002 = 3.0; double r16003 = pow(r15998, r16002); double r16004 = x; double r16005 = r16003 * r16004; double r16006 = r16001 + r16005; double r16007 = 230661.515625; double r16008 = r16006 + r16007; double r16009 = r16008 * r15998; double r16010 = t; double r16011 = r16009 + r16010; double r16012 = a; double r16013 = r15998 + r16012; double r16014 = r16013 * r15998; double r16015 = b; double r16016 = r16014 + r16015; double r16017 = r16016 * r15998; double r16018 = 1.0; double r16019 = pow(r16017, r16018); double r16020 = c; double r16021 = r16019 + r16020; double r16022 = r16021 * r15998; double r16023 = i; double r16024 = r16022 + r16023; double r16025 = r16011 / r16024; return r16025; } 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 r16026, r16027, r16028, r16029, r16030, r16031, r16032, r16033, r16034, r16035, r16036, r16037, r16038, r16039, r16040, r16041, r16042, r16043, r16044, r16045, r16046, r16047, r16048, r16049, r16050, r16051; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16026); mpfr_init(r16027); mpfr_init(r16028); mpfr_init(r16029); mpfr_init(r16030); mpfr_init(r16031); mpfr_init_set_str(r16032, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r16033); mpfr_init(r16034); mpfr_init_set_str(r16035, "230661.510616", 10, MPFR_RNDN); mpfr_init(r16036); mpfr_init(r16037); mpfr_init(r16038); mpfr_init(r16039); mpfr_init(r16040); mpfr_init(r16041); mpfr_init(r16042); mpfr_init(r16043); mpfr_init(r16044); mpfr_init(r16045); mpfr_init(r16046); mpfr_init(r16047); mpfr_init(r16048); mpfr_init(r16049); mpfr_init(r16050); mpfr_init(r16051); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r16026, x, MPFR_RNDN); mpfr_set_d(r16027, y, MPFR_RNDN); mpfr_mul(r16028, r16026, r16027, MPFR_RNDN); mpfr_set_d(r16029, z, MPFR_RNDN); mpfr_add(r16030, r16028, r16029, MPFR_RNDN); mpfr_mul(r16031, r16030, r16027, MPFR_RNDN); ; mpfr_add(r16033, r16031, r16032, MPFR_RNDN); mpfr_mul(r16034, r16033, r16027, MPFR_RNDN); ; mpfr_add(r16036, r16034, r16035, MPFR_RNDN); mpfr_mul(r16037, r16036, r16027, MPFR_RNDN); mpfr_set_d(r16038, t, MPFR_RNDN); mpfr_add(r16039, r16037, r16038, MPFR_RNDN); mpfr_set_d(r16040, a, MPFR_RNDN); mpfr_add(r16041, r16027, r16040, MPFR_RNDN); mpfr_mul(r16042, r16041, r16027, MPFR_RNDN); mpfr_set_d(r16043, b, MPFR_RNDN); mpfr_add(r16044, r16042, r16043, MPFR_RNDN); mpfr_mul(r16045, r16044, r16027, MPFR_RNDN); mpfr_set_d(r16046, c, MPFR_RNDN); mpfr_add(r16047, r16045, r16046, MPFR_RNDN); mpfr_mul(r16048, r16047, r16027, MPFR_RNDN); mpfr_set_d(r16049, i, MPFR_RNDN); mpfr_add(r16050, r16048, r16049, MPFR_RNDN); mpfr_div(r16051, r16039, r16050, MPFR_RNDN); return mpfr_get_d(r16051, MPFR_RNDN); } static mpfr_t r16052, r16053, r16054, r16055, r16056, r16057, r16058, r16059, r16060, r16061, r16062, r16063, r16064, r16065, r16066, r16067, r16068, r16069, r16070, r16071, r16072, r16073, r16074, r16075, r16076, r16077, r16078, r16079; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16052); mpfr_init(r16053); mpfr_init(r16054); mpfr_init(r16055); mpfr_init_set_str(r16056, "3", 10, MPFR_RNDN); mpfr_init(r16057); mpfr_init(r16058); mpfr_init(r16059); mpfr_init(r16060); mpfr_init_set_str(r16061, "230661.510616", 10, MPFR_RNDN); mpfr_init(r16062); mpfr_init(r16063); mpfr_init(r16064); mpfr_init(r16065); mpfr_init(r16066); mpfr_init(r16067); mpfr_init(r16068); mpfr_init(r16069); mpfr_init(r16070); mpfr_init(r16071); mpfr_init_set_str(r16072, "1", 10, MPFR_RNDN); mpfr_init(r16073); mpfr_init(r16074); mpfr_init(r16075); mpfr_init(r16076); mpfr_init(r16077); mpfr_init(r16078); mpfr_init(r16079); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r16052, y, MPFR_RNDN); mpfr_sqr(r16053, r16052, MPFR_RNDN); mpfr_set_d(r16054, z, MPFR_RNDN); mpfr_mul(r16055, r16053, r16054, MPFR_RNDN); ; mpfr_pow(r16057, r16052, r16056, MPFR_RNDN); mpfr_set_d(r16058, x, MPFR_RNDN); mpfr_mul(r16059, r16057, r16058, MPFR_RNDN); mpfr_add(r16060, r16055, r16059, MPFR_RNDN); ; mpfr_add(r16062, r16060, r16061, MPFR_RNDN); mpfr_mul(r16063, r16062, r16052, MPFR_RNDN); mpfr_set_d(r16064, t, MPFR_RNDN); mpfr_add(r16065, r16063, r16064, MPFR_RNDN); mpfr_set_d(r16066, a, MPFR_RNDN); mpfr_add(r16067, r16052, r16066, MPFR_RNDN); mpfr_mul(r16068, r16067, r16052, MPFR_RNDN); mpfr_set_d(r16069, b, MPFR_RNDN); mpfr_add(r16070, r16068, r16069, MPFR_RNDN); mpfr_mul(r16071, r16070, r16052, MPFR_RNDN); ; mpfr_pow(r16073, r16071, r16072, MPFR_RNDN); mpfr_set_d(r16074, c, MPFR_RNDN); mpfr_add(r16075, r16073, r16074, MPFR_RNDN); mpfr_mul(r16076, r16075, r16052, MPFR_RNDN); mpfr_set_d(r16077, i, MPFR_RNDN); mpfr_add(r16078, r16076, r16077, MPFR_RNDN); mpfr_div(r16079, r16065, r16078, MPFR_RNDN); return mpfr_get_d(r16079, MPFR_RNDN); } static mpfr_t r16080, r16081, r16082, r16083, r16084, r16085, r16086, r16087, r16088, r16089, r16090, r16091, r16092, r16093, r16094, r16095, r16096, r16097, r16098, r16099, r16100, r16101, r16102, r16103, r16104, r16105, r16106, r16107; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16080); mpfr_init(r16081); mpfr_init(r16082); mpfr_init(r16083); mpfr_init_set_str(r16084, "3", 10, MPFR_RNDN); mpfr_init(r16085); mpfr_init(r16086); mpfr_init(r16087); mpfr_init(r16088); mpfr_init_set_str(r16089, "230661.510616", 10, MPFR_RNDN); mpfr_init(r16090); mpfr_init(r16091); mpfr_init(r16092); mpfr_init(r16093); mpfr_init(r16094); mpfr_init(r16095); mpfr_init(r16096); mpfr_init(r16097); mpfr_init(r16098); mpfr_init(r16099); mpfr_init_set_str(r16100, "1", 10, MPFR_RNDN); mpfr_init(r16101); mpfr_init(r16102); mpfr_init(r16103); mpfr_init(r16104); mpfr_init(r16105); mpfr_init(r16106); mpfr_init(r16107); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r16080, y, MPFR_RNDN); mpfr_sqr(r16081, r16080, MPFR_RNDN); mpfr_set_d(r16082, z, MPFR_RNDN); mpfr_mul(r16083, r16081, r16082, MPFR_RNDN); ; mpfr_pow(r16085, r16080, r16084, MPFR_RNDN); mpfr_set_d(r16086, x, MPFR_RNDN); mpfr_mul(r16087, r16085, r16086, MPFR_RNDN); mpfr_add(r16088, r16083, r16087, MPFR_RNDN); ; mpfr_add(r16090, r16088, r16089, MPFR_RNDN); mpfr_mul(r16091, r16090, r16080, MPFR_RNDN); mpfr_set_d(r16092, t, MPFR_RNDN); mpfr_add(r16093, r16091, r16092, MPFR_RNDN); mpfr_set_d(r16094, a, MPFR_RNDN); mpfr_add(r16095, r16080, r16094, MPFR_RNDN); mpfr_mul(r16096, r16095, r16080, MPFR_RNDN); mpfr_set_d(r16097, b, MPFR_RNDN); mpfr_add(r16098, r16096, r16097, MPFR_RNDN); mpfr_mul(r16099, r16098, r16080, MPFR_RNDN); ; mpfr_pow(r16101, r16099, r16100, MPFR_RNDN); mpfr_set_d(r16102, c, MPFR_RNDN); mpfr_add(r16103, r16101, r16102, MPFR_RNDN); mpfr_mul(r16104, r16103, r16080, MPFR_RNDN); mpfr_set_d(r16105, i, MPFR_RNDN); mpfr_add(r16106, r16104, r16105, MPFR_RNDN); mpfr_div(r16107, r16093, r16106, MPFR_RNDN); return mpfr_get_d(r16107, MPFR_RNDN); }