#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r27030 = 2; float r27031 = sqrt(r27030); float r27032 = 4; float r27033 = r27031 / r27032; float r27034 = 1; float r27035 = 3; float r27036 = v; float r27037 = r27036 * r27036; float r27038 = r27035 * r27037; float r27039 = r27034 - r27038; float r27040 = sqrt(r27039); float r27041 = r27033 * r27040; float r27042 = r27034 - r27037; float r27043 = r27041 * r27042; return r27043; } double f_id(double v) { double r27044 = 2; double r27045 = sqrt(r27044); double r27046 = 4; double r27047 = r27045 / r27046; double r27048 = 1; double r27049 = 3; double r27050 = v; double r27051 = r27050 * r27050; double r27052 = r27049 * r27051; double r27053 = r27048 - r27052; double r27054 = sqrt(r27053); double r27055 = r27047 * r27054; double r27056 = r27048 - r27051; double r27057 = r27055 * r27056; return r27057; } double f_of(float v) { float r27058 = 2; float r27059 = sqrt(r27058); float r27060 = 4; float r27061 = r27059 / r27060; float r27062 = 1; float r27063 = 3; float r27064 = v; float r27065 = r27064 * r27064; float r27066 = r27063 * r27065; float r27067 = exp(r27066); float r27068 = log(r27067); float r27069 = r27062 - r27068; float r27070 = sqrt(r27069); float r27071 = r27061 * r27070; float r27072 = r27062 - r27065; float r27073 = r27071 * r27072; return r27073; } double f_od(double v) { double r27074 = 2; double r27075 = sqrt(r27074); double r27076 = 4; double r27077 = r27075 / r27076; double r27078 = 1; double r27079 = 3; double r27080 = v; double r27081 = r27080 * r27080; double r27082 = r27079 * r27081; double r27083 = exp(r27082); double r27084 = log(r27083); double r27085 = r27078 - r27084; double r27086 = sqrt(r27085); double r27087 = r27077 * r27086; double r27088 = r27078 - r27081; double r27089 = r27087 * r27088; return r27089; } 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 r27090, r27091, r27092, r27093, r27094, r27095, r27096, r27097, r27098, r27099, r27100, r27101, r27102, r27103; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r27090, "2", 10, MPFR_RNDN); mpfr_init(r27091); mpfr_init_set_str(r27092, "4", 10, MPFR_RNDN); mpfr_init(r27093); mpfr_init_set_str(r27094, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27095, "3", 10, MPFR_RNDN); mpfr_init(r27096); mpfr_init(r27097); mpfr_init(r27098); mpfr_init(r27099); mpfr_init(r27100); mpfr_init(r27101); mpfr_init(r27102); mpfr_init(r27103); } double f_im(double v) { ; mpfr_sqrt(r27091, r27090, MPFR_RNDN); ; mpfr_div(r27093, r27091, r27092, MPFR_RNDN); ; ; mpfr_set_d(r27096, v, MPFR_RNDN); mpfr_mul(r27097, r27096, r27096, MPFR_RNDN); mpfr_mul(r27098, r27095, r27097, MPFR_RNDN); mpfr_sub(r27099, r27094, r27098, MPFR_RNDN); mpfr_sqrt(r27100, r27099, MPFR_RNDN); mpfr_mul(r27101, r27093, r27100, MPFR_RNDN); mpfr_sub(r27102, r27094, r27097, MPFR_RNDN); mpfr_mul(r27103, r27101, r27102, MPFR_RNDN); return mpfr_get_d(r27103, MPFR_RNDN); } static mpfr_t r27104, r27105, r27106, r27107, r27108, r27109, r27110, r27111, r27112, r27113, r27114, r27115, r27116, r27117, r27118, r27119; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27104, "2", 10, MPFR_RNDN); mpfr_init(r27105); mpfr_init_set_str(r27106, "4", 10, MPFR_RNDN); mpfr_init(r27107); mpfr_init_set_str(r27108, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27109, "3", 10, MPFR_RNDN); mpfr_init(r27110); mpfr_init(r27111); mpfr_init(r27112); mpfr_init(r27113); mpfr_init(r27114); mpfr_init(r27115); mpfr_init(r27116); mpfr_init(r27117); mpfr_init(r27118); mpfr_init(r27119); } double f_fm(double v) { ; mpfr_sqrt(r27105, r27104, MPFR_RNDN); ; mpfr_div(r27107, r27105, r27106, MPFR_RNDN); ; ; mpfr_set_d(r27110, v, MPFR_RNDN); mpfr_mul(r27111, r27110, r27110, MPFR_RNDN); mpfr_mul(r27112, r27109, r27111, MPFR_RNDN); mpfr_exp(r27113, r27112, MPFR_RNDN); mpfr_log(r27114, r27113, MPFR_RNDN); mpfr_sub(r27115, r27108, r27114, MPFR_RNDN); mpfr_sqrt(r27116, r27115, MPFR_RNDN); mpfr_mul(r27117, r27107, r27116, MPFR_RNDN); mpfr_sub(r27118, r27108, r27111, MPFR_RNDN); mpfr_mul(r27119, r27117, r27118, MPFR_RNDN); return mpfr_get_d(r27119, MPFR_RNDN); } static mpfr_t r27120, r27121, r27122, r27123, r27124, r27125, r27126, r27127, r27128, r27129, r27130, r27131, r27132, r27133, r27134, r27135; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27120, "2", 10, MPFR_RNDN); mpfr_init(r27121); mpfr_init_set_str(r27122, "4", 10, MPFR_RNDN); mpfr_init(r27123); mpfr_init_set_str(r27124, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27125, "3", 10, MPFR_RNDN); mpfr_init(r27126); mpfr_init(r27127); mpfr_init(r27128); mpfr_init(r27129); mpfr_init(r27130); mpfr_init(r27131); mpfr_init(r27132); mpfr_init(r27133); mpfr_init(r27134); mpfr_init(r27135); } double f_dm(double v) { ; mpfr_sqrt(r27121, r27120, MPFR_RNDN); ; mpfr_div(r27123, r27121, r27122, MPFR_RNDN); ; ; mpfr_set_d(r27126, v, MPFR_RNDN); mpfr_mul(r27127, r27126, r27126, MPFR_RNDN); mpfr_mul(r27128, r27125, r27127, MPFR_RNDN); mpfr_exp(r27129, r27128, MPFR_RNDN); mpfr_log(r27130, r27129, MPFR_RNDN); mpfr_sub(r27131, r27124, r27130, MPFR_RNDN); mpfr_sqrt(r27132, r27131, MPFR_RNDN); mpfr_mul(r27133, r27123, r27132, MPFR_RNDN); mpfr_sub(r27134, r27124, r27127, MPFR_RNDN); mpfr_mul(r27135, r27133, r27134, MPFR_RNDN); return mpfr_get_d(r27135, MPFR_RNDN); }