#include #include #include #include #include char *name = "Migdal et al, Equation (51)"; double f_if(float k, float n) { float r20053 = k; float r20054 = sqrt(r20053); float r20055 = 1.0/r20054; float r20056 = 2.0f; float r20057 = atan2(1.0, 0.0); float r20058 = r20056 * r20057; float r20059 = n; float r20060 = r20058 * r20059; float r20061 = 1.0f; float r20062 = r20061 - r20053; float r20063 = r20062 / r20056; float r20064 = pow(r20060, r20063); float r20065 = r20055 * r20064; return r20065; } double f_id(double k, double n) { double r20066 = k; double r20067 = sqrt(r20066); double r20068 = 1.0/r20067; double r20069 = 2.0; double r20070 = atan2(1.0, 0.0); double r20071 = r20069 * r20070; double r20072 = n; double r20073 = r20071 * r20072; double r20074 = 1.0; double r20075 = r20074 - r20066; double r20076 = r20075 / r20069; double r20077 = pow(r20073, r20076); double r20078 = r20068 * r20077; return r20078; } double f_of(float k, float n) { float r20079 = k; float r20080 = sqrt(r20079); float r20081 = 1.0/r20080; float r20082 = 2.0f; float r20083 = atan2(1.0, 0.0); float r20084 = r20082 * r20083; float r20085 = n; float r20086 = r20084 * r20085; float r20087 = 1.0f; float r20088 = r20087 / r20082; float r20089 = pow(r20086, r20088); float r20090 = r20079 / r20082; float r20091 = pow(r20086, r20090); float r20092 = r20089 / r20091; float r20093 = r20081 * r20092; return r20093; } double f_od(double k, double n) { double r20094 = k; double r20095 = sqrt(r20094); double r20096 = 1.0/r20095; double r20097 = 2.0; double r20098 = atan2(1.0, 0.0); double r20099 = r20097 * r20098; double r20100 = n; double r20101 = r20099 * r20100; double r20102 = 1.0; double r20103 = r20102 / r20097; double r20104 = pow(r20101, r20103); double r20105 = r20094 / r20097; double r20106 = pow(r20101, r20105); double r20107 = r20104 / r20106; double r20108 = r20096 * r20107; return r20108; } 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 r20109, r20110, r20111, r20112, r20113, r20114, r20115, r20116, r20117, r20118, r20119, r20120, r20121; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r20109); mpfr_init(r20110); mpfr_init(r20111); mpfr_init_set_str(r20112, "2", 10, MPFR_RNDN); mpfr_init(r20113); mpfr_init(r20114); mpfr_init(r20115); mpfr_init(r20116); mpfr_init_set_str(r20117, "1", 10, MPFR_RNDN); mpfr_init(r20118); mpfr_init(r20119); mpfr_init(r20120); mpfr_init(r20121); } double f_im(double k, double n) { mpfr_set_d(r20109, k, MPFR_RNDN); mpfr_sqrt(r20110, r20109, MPFR_RNDN); mpfr_ui_div(r20111, 1, r20110, MPFR_RNDN); ; mpfr_const_pi(r20113, MPFR_RNDN); mpfr_mul(r20114, r20112, r20113, MPFR_RNDN); mpfr_set_d(r20115, n, MPFR_RNDN); mpfr_mul(r20116, r20114, r20115, MPFR_RNDN); ; mpfr_sub(r20118, r20117, r20109, MPFR_RNDN); mpfr_div(r20119, r20118, r20112, MPFR_RNDN); mpfr_pow(r20120, r20116, r20119, MPFR_RNDN); mpfr_mul(r20121, r20111, r20120, MPFR_RNDN); return mpfr_get_d(r20121, MPFR_RNDN); } static mpfr_t r20122, r20123, r20124, r20125, r20126, r20127, r20128, r20129, r20130, r20131, r20132, r20133, r20134, r20135, r20136; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r20122); mpfr_init(r20123); mpfr_init(r20124); mpfr_init_set_str(r20125, "2", 10, MPFR_RNDN); mpfr_init(r20126); mpfr_init(r20127); mpfr_init(r20128); mpfr_init(r20129); mpfr_init_set_str(r20130, "1", 10, MPFR_RNDN); mpfr_init(r20131); mpfr_init(r20132); mpfr_init(r20133); mpfr_init(r20134); mpfr_init(r20135); mpfr_init(r20136); } double f_fm(double k, double n) { mpfr_set_d(r20122, k, MPFR_RNDN); mpfr_sqrt(r20123, r20122, MPFR_RNDN); mpfr_ui_div(r20124, 1, r20123, MPFR_RNDN); ; mpfr_const_pi(r20126, MPFR_RNDN); mpfr_mul(r20127, r20125, r20126, MPFR_RNDN); mpfr_set_d(r20128, n, MPFR_RNDN); mpfr_mul(r20129, r20127, r20128, MPFR_RNDN); ; mpfr_div(r20131, r20130, r20125, MPFR_RNDN); mpfr_pow(r20132, r20129, r20131, MPFR_RNDN); mpfr_div(r20133, r20122, r20125, MPFR_RNDN); mpfr_pow(r20134, r20129, r20133, MPFR_RNDN); mpfr_div(r20135, r20132, r20134, MPFR_RNDN); mpfr_mul(r20136, r20124, r20135, MPFR_RNDN); return mpfr_get_d(r20136, MPFR_RNDN); } static mpfr_t r20137, r20138, r20139, r20140, r20141, r20142, r20143, r20144, r20145, r20146, r20147, r20148, r20149, r20150, r20151; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r20137); mpfr_init(r20138); mpfr_init(r20139); mpfr_init_set_str(r20140, "2", 10, MPFR_RNDN); mpfr_init(r20141); mpfr_init(r20142); mpfr_init(r20143); mpfr_init(r20144); mpfr_init_set_str(r20145, "1", 10, MPFR_RNDN); mpfr_init(r20146); mpfr_init(r20147); mpfr_init(r20148); mpfr_init(r20149); mpfr_init(r20150); mpfr_init(r20151); } double f_dm(double k, double n) { mpfr_set_d(r20137, k, MPFR_RNDN); mpfr_sqrt(r20138, r20137, MPFR_RNDN); mpfr_ui_div(r20139, 1, r20138, MPFR_RNDN); ; mpfr_const_pi(r20141, MPFR_RNDN); mpfr_mul(r20142, r20140, r20141, MPFR_RNDN); mpfr_set_d(r20143, n, MPFR_RNDN); mpfr_mul(r20144, r20142, r20143, MPFR_RNDN); ; mpfr_div(r20146, r20145, r20140, MPFR_RNDN); mpfr_pow(r20147, r20144, r20146, MPFR_RNDN); mpfr_div(r20148, r20137, r20140, MPFR_RNDN); mpfr_pow(r20149, r20144, r20148, MPFR_RNDN); mpfr_div(r20150, r20147, r20149, MPFR_RNDN); mpfr_mul(r20151, r20139, r20150, MPFR_RNDN); return mpfr_get_d(r20151, MPFR_RNDN); }