#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r26993 = 4; float r26994 = 3; float r26995 = atan2(1.0, 0.0); float r26996 = r26994 * r26995; float r26997 = 1; float r26998 = v; float r26999 = r26998 * r26998; float r27000 = r26997 - r26999; float r27001 = r26996 * r27000; float r27002 = 2; float r27003 = 6; float r27004 = r27003 * r26999; float r27005 = r27002 - r27004; float r27006 = sqrt(r27005); float r27007 = r27001 * r27006; float r27008 = r26993 / r27007; return r27008; } double f_id(double v) { double r27009 = 4; double r27010 = 3; double r27011 = atan2(1.0, 0.0); double r27012 = r27010 * r27011; double r27013 = 1; double r27014 = v; double r27015 = r27014 * r27014; double r27016 = r27013 - r27015; double r27017 = r27012 * r27016; double r27018 = 2; double r27019 = 6; double r27020 = r27019 * r27015; double r27021 = r27018 - r27020; double r27022 = sqrt(r27021); double r27023 = r27017 * r27022; double r27024 = r27009 / r27023; return r27024; } double f_of(float v) { float r27025 = 4; float r27026 = 3; float r27027 = atan2(1.0, 0.0); float r27028 = r27026 * r27027; float r27029 = 1; float r27030 = v; float r27031 = r27030 * r27030; float r27032 = r27029 - r27031; float r27033 = r27028 * r27032; float r27034 = r27025 / r27033; float r27035 = 2; float r27036 = 6; float r27037 = r27036 * r27031; float r27038 = r27035 - r27037; float r27039 = sqrt(r27038); float r27040 = r27034 / r27039; return r27040; } double f_od(double v) { double r27041 = 4; double r27042 = 3; double r27043 = atan2(1.0, 0.0); double r27044 = r27042 * r27043; double r27045 = 1; double r27046 = v; double r27047 = r27046 * r27046; double r27048 = r27045 - r27047; double r27049 = r27044 * r27048; double r27050 = r27041 / r27049; double r27051 = 2; double r27052 = 6; double r27053 = r27052 * r27047; double r27054 = r27051 - r27053; double r27055 = sqrt(r27054); double r27056 = r27050 / r27055; return r27056; } 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 r27057, r27058, r27059, r27060, r27061, r27062, r27063, r27064, r27065, r27066, r27067, r27068, r27069, r27070, r27071, r27072; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r27057, "4", 10, MPFR_RNDN); mpfr_init_set_str(r27058, "3", 10, MPFR_RNDN); mpfr_init(r27059); mpfr_init(r27060); mpfr_init_set_str(r27061, "1", 10, MPFR_RNDN); mpfr_init(r27062); mpfr_init(r27063); mpfr_init(r27064); mpfr_init(r27065); mpfr_init_set_str(r27066, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27067, "6", 10, MPFR_RNDN); mpfr_init(r27068); mpfr_init(r27069); mpfr_init(r27070); mpfr_init(r27071); mpfr_init(r27072); } double f_im(double v) { ; ; mpfr_const_pi(r27059, MPFR_RNDN); mpfr_mul(r27060, r27058, r27059, MPFR_RNDN); ; mpfr_set_d(r27062, v, MPFR_RNDN); mpfr_mul(r27063, r27062, r27062, MPFR_RNDN); mpfr_sub(r27064, r27061, r27063, MPFR_RNDN); mpfr_mul(r27065, r27060, r27064, MPFR_RNDN); ; ; mpfr_mul(r27068, r27067, r27063, MPFR_RNDN); mpfr_sub(r27069, r27066, r27068, MPFR_RNDN); mpfr_sqrt(r27070, r27069, MPFR_RNDN); mpfr_mul(r27071, r27065, r27070, MPFR_RNDN); mpfr_div(r27072, r27057, r27071, MPFR_RNDN); return mpfr_get_d(r27072, MPFR_RNDN); } static mpfr_t r27073, r27074, r27075, r27076, r27077, r27078, r27079, r27080, r27081, r27082, r27083, r27084, r27085, r27086, r27087, r27088; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27073, "4", 10, MPFR_RNDN); mpfr_init_set_str(r27074, "3", 10, MPFR_RNDN); mpfr_init(r27075); mpfr_init(r27076); mpfr_init_set_str(r27077, "1", 10, MPFR_RNDN); mpfr_init(r27078); mpfr_init(r27079); mpfr_init(r27080); mpfr_init(r27081); mpfr_init(r27082); mpfr_init_set_str(r27083, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27084, "6", 10, MPFR_RNDN); mpfr_init(r27085); mpfr_init(r27086); mpfr_init(r27087); mpfr_init(r27088); } double f_fm(double v) { ; ; mpfr_const_pi(r27075, MPFR_RNDN); mpfr_mul(r27076, r27074, r27075, MPFR_RNDN); ; mpfr_set_d(r27078, v, MPFR_RNDN); mpfr_mul(r27079, r27078, r27078, MPFR_RNDN); mpfr_sub(r27080, r27077, r27079, MPFR_RNDN); mpfr_mul(r27081, r27076, r27080, MPFR_RNDN); mpfr_div(r27082, r27073, r27081, MPFR_RNDN); ; ; mpfr_mul(r27085, r27084, r27079, MPFR_RNDN); mpfr_sub(r27086, r27083, r27085, MPFR_RNDN); mpfr_sqrt(r27087, r27086, MPFR_RNDN); mpfr_div(r27088, r27082, r27087, MPFR_RNDN); return mpfr_get_d(r27088, MPFR_RNDN); } static mpfr_t r27089, r27090, r27091, r27092, r27093, r27094, r27095, r27096, r27097, r27098, r27099, r27100, r27101, r27102, r27103, r27104; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27089, "4", 10, MPFR_RNDN); mpfr_init_set_str(r27090, "3", 10, MPFR_RNDN); mpfr_init(r27091); mpfr_init(r27092); mpfr_init_set_str(r27093, "1", 10, MPFR_RNDN); mpfr_init(r27094); mpfr_init(r27095); mpfr_init(r27096); mpfr_init(r27097); mpfr_init(r27098); mpfr_init_set_str(r27099, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27100, "6", 10, MPFR_RNDN); mpfr_init(r27101); mpfr_init(r27102); mpfr_init(r27103); mpfr_init(r27104); } double f_dm(double v) { ; ; mpfr_const_pi(r27091, MPFR_RNDN); mpfr_mul(r27092, r27090, r27091, MPFR_RNDN); ; mpfr_set_d(r27094, v, MPFR_RNDN); mpfr_mul(r27095, r27094, r27094, MPFR_RNDN); mpfr_sub(r27096, r27093, r27095, MPFR_RNDN); mpfr_mul(r27097, r27092, r27096, MPFR_RNDN); mpfr_div(r27098, r27089, r27097, MPFR_RNDN); ; ; mpfr_mul(r27101, r27100, r27095, MPFR_RNDN); mpfr_sub(r27102, r27099, r27101, MPFR_RNDN); mpfr_sqrt(r27103, r27102, MPFR_RNDN); mpfr_div(r27104, r27098, r27103, MPFR_RNDN); return mpfr_get_d(r27104, MPFR_RNDN); }