#include #include #include #include #include char *name = "Migdal et al, Equation (51)"; double f_if(float k, float n) { float r20005 = k; float r20006 = sqrt(r20005); float r20007 = 1.0/r20006; float r20008 = 2.0f; float r20009 = atan2(1.0, 0.0); float r20010 = r20008 * r20009; float r20011 = n; float r20012 = r20010 * r20011; float r20013 = 1.0f; float r20014 = r20013 - r20005; float r20015 = r20014 / r20008; float r20016 = pow(r20012, r20015); float r20017 = r20007 * r20016; return r20017; } double f_id(double k, double n) { double r20018 = k; double r20019 = sqrt(r20018); double r20020 = 1.0/r20019; double r20021 = 2.0; double r20022 = atan2(1.0, 0.0); double r20023 = r20021 * r20022; double r20024 = n; double r20025 = r20023 * r20024; double r20026 = 1.0; double r20027 = r20026 - r20018; double r20028 = r20027 / r20021; double r20029 = pow(r20025, r20028); double r20030 = r20020 * r20029; return r20030; } double f_of(float k, float n) { float r20031 = k; float r20032 = sqrt(r20031); float r20033 = 1.0/r20032; float r20034 = 2.0f; float r20035 = atan2(1.0, 0.0); float r20036 = r20034 * r20035; float r20037 = n; float r20038 = r20036 * r20037; float r20039 = 1.0f; float r20040 = r20039 / r20034; float r20041 = pow(r20038, r20040); float r20042 = r20031 / r20034; float r20043 = pow(r20038, r20042); float r20044 = r20041 / r20043; float r20045 = r20033 * r20044; return r20045; } double f_od(double k, double n) { double r20046 = k; double r20047 = sqrt(r20046); double r20048 = 1.0/r20047; double r20049 = 2.0; double r20050 = atan2(1.0, 0.0); double r20051 = r20049 * r20050; double r20052 = n; double r20053 = r20051 * r20052; double r20054 = 1.0; double r20055 = r20054 / r20049; double r20056 = pow(r20053, r20055); double r20057 = r20046 / r20049; double r20058 = pow(r20053, r20057); double r20059 = r20056 / r20058; double r20060 = r20048 * r20059; return r20060; } 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 r20061, r20062, r20063, r20064, r20065, r20066, r20067, r20068, r20069, r20070, r20071, r20072, r20073; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r20061); mpfr_init(r20062); mpfr_init(r20063); mpfr_init_set_str(r20064, "2", 10, MPFR_RNDN); mpfr_init(r20065); mpfr_init(r20066); mpfr_init(r20067); mpfr_init(r20068); mpfr_init_set_str(r20069, "1", 10, MPFR_RNDN); mpfr_init(r20070); mpfr_init(r20071); mpfr_init(r20072); mpfr_init(r20073); } double f_im(double k, double n) { mpfr_set_d(r20061, k, MPFR_RNDN); mpfr_sqrt(r20062, r20061, MPFR_RNDN); mpfr_ui_div(r20063, 1, r20062, MPFR_RNDN); ; mpfr_const_pi(r20065, MPFR_RNDN); mpfr_mul(r20066, r20064, r20065, MPFR_RNDN); mpfr_set_d(r20067, n, MPFR_RNDN); mpfr_mul(r20068, r20066, r20067, MPFR_RNDN); ; mpfr_sub(r20070, r20069, r20061, MPFR_RNDN); mpfr_div(r20071, r20070, r20064, MPFR_RNDN); mpfr_pow(r20072, r20068, r20071, MPFR_RNDN); mpfr_mul(r20073, r20063, r20072, MPFR_RNDN); return mpfr_get_d(r20073, MPFR_RNDN); } static mpfr_t r20074, r20075, r20076, r20077, r20078, r20079, r20080, r20081, r20082, r20083, r20084, r20085, r20086, r20087, r20088; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r20074); mpfr_init(r20075); mpfr_init(r20076); mpfr_init_set_str(r20077, "2", 10, MPFR_RNDN); mpfr_init(r20078); mpfr_init(r20079); mpfr_init(r20080); mpfr_init(r20081); mpfr_init_set_str(r20082, "1", 10, MPFR_RNDN); mpfr_init(r20083); mpfr_init(r20084); mpfr_init(r20085); mpfr_init(r20086); mpfr_init(r20087); mpfr_init(r20088); } double f_fm(double k, double n) { mpfr_set_d(r20074, k, MPFR_RNDN); mpfr_sqrt(r20075, r20074, MPFR_RNDN); mpfr_ui_div(r20076, 1, r20075, MPFR_RNDN); ; mpfr_const_pi(r20078, MPFR_RNDN); mpfr_mul(r20079, r20077, r20078, MPFR_RNDN); mpfr_set_d(r20080, n, MPFR_RNDN); mpfr_mul(r20081, r20079, r20080, MPFR_RNDN); ; mpfr_div(r20083, r20082, r20077, MPFR_RNDN); mpfr_pow(r20084, r20081, r20083, MPFR_RNDN); mpfr_div(r20085, r20074, r20077, MPFR_RNDN); mpfr_pow(r20086, r20081, r20085, MPFR_RNDN); mpfr_div(r20087, r20084, r20086, MPFR_RNDN); mpfr_mul(r20088, r20076, r20087, MPFR_RNDN); return mpfr_get_d(r20088, MPFR_RNDN); } static mpfr_t r20089, r20090, r20091, r20092, r20093, r20094, r20095, r20096, r20097, r20098, r20099, r20100, r20101, r20102, r20103; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r20089); mpfr_init(r20090); mpfr_init(r20091); mpfr_init_set_str(r20092, "2", 10, MPFR_RNDN); mpfr_init(r20093); mpfr_init(r20094); mpfr_init(r20095); mpfr_init(r20096); mpfr_init_set_str(r20097, "1", 10, MPFR_RNDN); mpfr_init(r20098); mpfr_init(r20099); mpfr_init(r20100); mpfr_init(r20101); mpfr_init(r20102); mpfr_init(r20103); } double f_dm(double k, double n) { mpfr_set_d(r20089, k, MPFR_RNDN); mpfr_sqrt(r20090, r20089, MPFR_RNDN); mpfr_ui_div(r20091, 1, r20090, MPFR_RNDN); ; mpfr_const_pi(r20093, MPFR_RNDN); mpfr_mul(r20094, r20092, r20093, MPFR_RNDN); mpfr_set_d(r20095, n, MPFR_RNDN); mpfr_mul(r20096, r20094, r20095, MPFR_RNDN); ; mpfr_div(r20098, r20097, r20092, MPFR_RNDN); mpfr_pow(r20099, r20096, r20098, MPFR_RNDN); mpfr_div(r20100, r20089, r20092, MPFR_RNDN); mpfr_pow(r20101, r20096, r20100, MPFR_RNDN); mpfr_div(r20102, r20099, r20101, MPFR_RNDN); mpfr_mul(r20103, r20091, r20102, MPFR_RNDN); return mpfr_get_d(r20103, MPFR_RNDN); }