#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r11274 = 2; float r11275 = sqrt(r11274); float r11276 = 4; float r11277 = r11275 / r11276; float r11278 = 1; float r11279 = 3; float r11280 = v; float r11281 = r11280 * r11280; float r11282 = r11279 * r11281; float r11283 = r11278 - r11282; float r11284 = sqrt(r11283); float r11285 = r11277 * r11284; float r11286 = r11278 - r11281; float r11287 = r11285 * r11286; return r11287; } double f_id(double v) { double r11288 = 2; double r11289 = sqrt(r11288); double r11290 = 4; double r11291 = r11289 / r11290; double r11292 = 1; double r11293 = 3; double r11294 = v; double r11295 = r11294 * r11294; double r11296 = r11293 * r11295; double r11297 = r11292 - r11296; double r11298 = sqrt(r11297); double r11299 = r11291 * r11298; double r11300 = r11292 - r11295; double r11301 = r11299 * r11300; return r11301; } double f_of(float v) { float r11302 = 2; float r11303 = sqrt(r11302); float r11304 = 4; float r11305 = r11303 / r11304; float r11306 = 1; float r11307 = 3; float r11308 = v; float r11309 = r11308 * r11308; float r11310 = r11307 * r11309; float r11311 = exp(r11310); float r11312 = log(r11311); float r11313 = r11306 - r11312; float r11314 = sqrt(r11313); float r11315 = r11305 * r11314; float r11316 = r11306 - r11309; float r11317 = r11315 * r11316; return r11317; } double f_od(double v) { double r11318 = 2; double r11319 = sqrt(r11318); double r11320 = 4; double r11321 = r11319 / r11320; double r11322 = 1; double r11323 = 3; double r11324 = v; double r11325 = r11324 * r11324; double r11326 = r11323 * r11325; double r11327 = exp(r11326); double r11328 = log(r11327); double r11329 = r11322 - r11328; double r11330 = sqrt(r11329); double r11331 = r11321 * r11330; double r11332 = r11322 - r11325; double r11333 = r11331 * r11332; return r11333; } 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 r11334, r11335, r11336, r11337, r11338, r11339, r11340, r11341, r11342, r11343, r11344, r11345, r11346, r11347; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r11334, "2", 10, MPFR_RNDN); mpfr_init(r11335); mpfr_init_set_str(r11336, "4", 10, MPFR_RNDN); mpfr_init(r11337); mpfr_init_set_str(r11338, "1", 10, MPFR_RNDN); mpfr_init_set_str(r11339, "3", 10, MPFR_RNDN); mpfr_init(r11340); mpfr_init(r11341); mpfr_init(r11342); mpfr_init(r11343); mpfr_init(r11344); mpfr_init(r11345); mpfr_init(r11346); mpfr_init(r11347); } double f_im(double v) { ; mpfr_sqrt(r11335, r11334, MPFR_RNDN); ; mpfr_div(r11337, r11335, r11336, MPFR_RNDN); ; ; mpfr_set_d(r11340, v, MPFR_RNDN); mpfr_mul(r11341, r11340, r11340, MPFR_RNDN); mpfr_mul(r11342, r11339, r11341, MPFR_RNDN); mpfr_sub(r11343, r11338, r11342, MPFR_RNDN); mpfr_sqrt(r11344, r11343, MPFR_RNDN); mpfr_mul(r11345, r11337, r11344, MPFR_RNDN); mpfr_sub(r11346, r11338, r11341, MPFR_RNDN); mpfr_mul(r11347, r11345, r11346, MPFR_RNDN); return mpfr_get_d(r11347, MPFR_RNDN); } static mpfr_t r11348, r11349, r11350, r11351, r11352, r11353, r11354, r11355, r11356, r11357, r11358, r11359, r11360, r11361, r11362, r11363; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r11348, "2", 10, MPFR_RNDN); mpfr_init(r11349); mpfr_init_set_str(r11350, "4", 10, MPFR_RNDN); mpfr_init(r11351); mpfr_init_set_str(r11352, "1", 10, MPFR_RNDN); mpfr_init_set_str(r11353, "3", 10, MPFR_RNDN); mpfr_init(r11354); mpfr_init(r11355); mpfr_init(r11356); mpfr_init(r11357); mpfr_init(r11358); mpfr_init(r11359); mpfr_init(r11360); mpfr_init(r11361); mpfr_init(r11362); mpfr_init(r11363); } double f_fm(double v) { ; mpfr_sqrt(r11349, r11348, MPFR_RNDN); ; mpfr_div(r11351, r11349, r11350, MPFR_RNDN); ; ; mpfr_set_d(r11354, v, MPFR_RNDN); mpfr_mul(r11355, r11354, r11354, MPFR_RNDN); mpfr_mul(r11356, r11353, r11355, MPFR_RNDN); mpfr_exp(r11357, r11356, MPFR_RNDN); mpfr_log(r11358, r11357, MPFR_RNDN); mpfr_sub(r11359, r11352, r11358, MPFR_RNDN); mpfr_sqrt(r11360, r11359, MPFR_RNDN); mpfr_mul(r11361, r11351, r11360, MPFR_RNDN); mpfr_sub(r11362, r11352, r11355, MPFR_RNDN); mpfr_mul(r11363, r11361, r11362, MPFR_RNDN); return mpfr_get_d(r11363, MPFR_RNDN); } static mpfr_t r11364, r11365, r11366, r11367, r11368, r11369, r11370, r11371, r11372, r11373, r11374, r11375, r11376, r11377, r11378, r11379; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r11364, "2", 10, MPFR_RNDN); mpfr_init(r11365); mpfr_init_set_str(r11366, "4", 10, MPFR_RNDN); mpfr_init(r11367); mpfr_init_set_str(r11368, "1", 10, MPFR_RNDN); mpfr_init_set_str(r11369, "3", 10, MPFR_RNDN); mpfr_init(r11370); mpfr_init(r11371); mpfr_init(r11372); mpfr_init(r11373); mpfr_init(r11374); mpfr_init(r11375); mpfr_init(r11376); mpfr_init(r11377); mpfr_init(r11378); mpfr_init(r11379); } double f_dm(double v) { ; mpfr_sqrt(r11365, r11364, MPFR_RNDN); ; mpfr_div(r11367, r11365, r11366, MPFR_RNDN); ; ; mpfr_set_d(r11370, v, MPFR_RNDN); mpfr_mul(r11371, r11370, r11370, MPFR_RNDN); mpfr_mul(r11372, r11369, r11371, MPFR_RNDN); mpfr_exp(r11373, r11372, MPFR_RNDN); mpfr_log(r11374, r11373, MPFR_RNDN); mpfr_sub(r11375, r11368, r11374, MPFR_RNDN); mpfr_sqrt(r11376, r11375, MPFR_RNDN); mpfr_mul(r11377, r11367, r11376, MPFR_RNDN); mpfr_sub(r11378, r11368, r11371, MPFR_RNDN); mpfr_mul(r11379, r11377, r11378, MPFR_RNDN); return mpfr_get_d(r11379, MPFR_RNDN); }