#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r27090 = 2; float r27091 = sqrt(r27090); float r27092 = 4; float r27093 = r27091 / r27092; float r27094 = 1; float r27095 = 3; float r27096 = v; float r27097 = r27096 * r27096; float r27098 = r27095 * r27097; float r27099 = r27094 - r27098; float r27100 = sqrt(r27099); float r27101 = r27093 * r27100; float r27102 = r27094 - r27097; float r27103 = r27101 * r27102; return r27103; } double f_id(double v) { double r27104 = 2; double r27105 = sqrt(r27104); double r27106 = 4; double r27107 = r27105 / r27106; double r27108 = 1; double r27109 = 3; double r27110 = v; double r27111 = r27110 * r27110; double r27112 = r27109 * r27111; double r27113 = r27108 - r27112; double r27114 = sqrt(r27113); double r27115 = r27107 * r27114; double r27116 = r27108 - r27111; double r27117 = r27115 * r27116; return r27117; } double f_of(float v) { float r27118 = 2; float r27119 = sqrt(r27118); float r27120 = 4; float r27121 = r27119 / r27120; float r27122 = 1; float r27123 = 3; float r27124 = v; float r27125 = r27124 * r27124; float r27126 = r27123 * r27125; float r27127 = r27122 - r27126; float r27128 = sqrt(r27127); float r27129 = r27121 * r27128; float r27130 = exp(r27129); float r27131 = log(r27130); float r27132 = r27122 - r27125; float r27133 = r27131 * r27132; return r27133; } double f_od(double v) { double r27134 = 2; double r27135 = sqrt(r27134); double r27136 = 4; double r27137 = r27135 / r27136; double r27138 = 1; double r27139 = 3; double r27140 = v; double r27141 = r27140 * r27140; double r27142 = r27139 * r27141; double r27143 = r27138 - r27142; double r27144 = sqrt(r27143); double r27145 = r27137 * r27144; double r27146 = exp(r27145); double r27147 = log(r27146); double r27148 = r27138 - r27141; double r27149 = r27147 * r27148; return r27149; } 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 r27150, r27151, r27152, r27153, r27154, r27155, r27156, r27157, r27158, r27159, r27160, r27161, r27162, r27163; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27150, "2", 10, MPFR_RNDN); mpfr_init(r27151); mpfr_init_set_str(r27152, "4", 10, MPFR_RNDN); mpfr_init(r27153); mpfr_init_set_str(r27154, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27155, "3", 10, MPFR_RNDN); mpfr_init(r27156); mpfr_init(r27157); mpfr_init(r27158); mpfr_init(r27159); mpfr_init(r27160); mpfr_init(r27161); mpfr_init(r27162); mpfr_init(r27163); } double f_im(double v) { ; mpfr_sqrt(r27151, r27150, MPFR_RNDN); ; mpfr_div(r27153, r27151, r27152, MPFR_RNDN); ; ; mpfr_set_d(r27156, v, MPFR_RNDN); mpfr_mul(r27157, r27156, r27156, MPFR_RNDN); mpfr_mul(r27158, r27155, r27157, MPFR_RNDN); mpfr_sub(r27159, r27154, r27158, MPFR_RNDN); mpfr_sqrt(r27160, r27159, MPFR_RNDN); mpfr_mul(r27161, r27153, r27160, MPFR_RNDN); mpfr_sub(r27162, r27154, r27157, MPFR_RNDN); mpfr_mul(r27163, r27161, r27162, MPFR_RNDN); return mpfr_get_d(r27163, MPFR_RNDN); } static mpfr_t r27164, r27165, r27166, r27167, r27168, r27169, r27170, r27171, r27172, r27173, r27174, r27175, r27176, r27177, r27178, r27179; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27164, "2", 10, MPFR_RNDN); mpfr_init(r27165); mpfr_init_set_str(r27166, "4", 10, MPFR_RNDN); mpfr_init(r27167); mpfr_init_set_str(r27168, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27169, "3", 10, MPFR_RNDN); mpfr_init(r27170); mpfr_init(r27171); mpfr_init(r27172); mpfr_init(r27173); mpfr_init(r27174); mpfr_init(r27175); mpfr_init(r27176); mpfr_init(r27177); mpfr_init(r27178); mpfr_init(r27179); } double f_fm(double v) { ; mpfr_sqrt(r27165, r27164, MPFR_RNDN); ; mpfr_div(r27167, r27165, r27166, MPFR_RNDN); ; ; mpfr_set_d(r27170, v, MPFR_RNDN); mpfr_mul(r27171, r27170, r27170, MPFR_RNDN); mpfr_mul(r27172, r27169, r27171, MPFR_RNDN); mpfr_sub(r27173, r27168, r27172, MPFR_RNDN); mpfr_sqrt(r27174, r27173, MPFR_RNDN); mpfr_mul(r27175, r27167, r27174, MPFR_RNDN); mpfr_exp(r27176, r27175, MPFR_RNDN); mpfr_log(r27177, r27176, MPFR_RNDN); mpfr_sub(r27178, r27168, r27171, MPFR_RNDN); mpfr_mul(r27179, r27177, r27178, MPFR_RNDN); return mpfr_get_d(r27179, MPFR_RNDN); } static mpfr_t r27180, r27181, r27182, r27183, r27184, r27185, r27186, r27187, r27188, r27189, r27190, r27191, r27192, r27193, r27194, r27195; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27180, "2", 10, MPFR_RNDN); mpfr_init(r27181); mpfr_init_set_str(r27182, "4", 10, MPFR_RNDN); mpfr_init(r27183); mpfr_init_set_str(r27184, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27185, "3", 10, MPFR_RNDN); mpfr_init(r27186); mpfr_init(r27187); mpfr_init(r27188); mpfr_init(r27189); mpfr_init(r27190); mpfr_init(r27191); mpfr_init(r27192); mpfr_init(r27193); mpfr_init(r27194); mpfr_init(r27195); } double f_dm(double v) { ; mpfr_sqrt(r27181, r27180, MPFR_RNDN); ; mpfr_div(r27183, r27181, r27182, MPFR_RNDN); ; ; mpfr_set_d(r27186, v, MPFR_RNDN); mpfr_mul(r27187, r27186, r27186, MPFR_RNDN); mpfr_mul(r27188, r27185, r27187, MPFR_RNDN); mpfr_sub(r27189, r27184, r27188, MPFR_RNDN); mpfr_sqrt(r27190, r27189, MPFR_RNDN); mpfr_mul(r27191, r27183, r27190, MPFR_RNDN); mpfr_exp(r27192, r27191, MPFR_RNDN); mpfr_log(r27193, r27192, MPFR_RNDN); mpfr_sub(r27194, r27184, r27187, MPFR_RNDN); mpfr_mul(r27195, r27193, r27194, MPFR_RNDN); return mpfr_get_d(r27195, MPFR_RNDN); }