#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r23375 = 2; float r23376 = sqrt(r23375); float r23377 = 4; float r23378 = r23376 / r23377; float r23379 = 1; float r23380 = 3; float r23381 = v; float r23382 = r23381 * r23381; float r23383 = r23380 * r23382; float r23384 = r23379 - r23383; float r23385 = sqrt(r23384); float r23386 = r23378 * r23385; float r23387 = r23379 - r23382; float r23388 = r23386 * r23387; return r23388; } double f_id(double v) { double r23389 = 2; double r23390 = sqrt(r23389); double r23391 = 4; double r23392 = r23390 / r23391; double r23393 = 1; double r23394 = 3; double r23395 = v; double r23396 = r23395 * r23395; double r23397 = r23394 * r23396; double r23398 = r23393 - r23397; double r23399 = sqrt(r23398); double r23400 = r23392 * r23399; double r23401 = r23393 - r23396; double r23402 = r23400 * r23401; return r23402; } double f_of(float v) { float r23403 = 2; float r23404 = sqrt(r23403); float r23405 = 4; float r23406 = r23404 / r23405; float r23407 = log(r23406); float r23408 = 1; float r23409 = 3; float r23410 = v; float r23411 = r23410 * r23410; float r23412 = r23409 * r23411; float r23413 = r23408 - r23412; float r23414 = sqrt(r23413); float r23415 = log(r23414); float r23416 = r23407 + r23415; float r23417 = r23408 - r23411; float r23418 = log(r23417); float r23419 = r23416 + r23418; float r23420 = exp(r23419); return r23420; } double f_od(double v) { double r23421 = 2; double r23422 = sqrt(r23421); double r23423 = 4; double r23424 = r23422 / r23423; double r23425 = log(r23424); double r23426 = 1; double r23427 = 3; double r23428 = v; double r23429 = r23428 * r23428; double r23430 = r23427 * r23429; double r23431 = r23426 - r23430; double r23432 = sqrt(r23431); double r23433 = log(r23432); double r23434 = r23425 + r23433; double r23435 = r23426 - r23429; double r23436 = log(r23435); double r23437 = r23434 + r23436; double r23438 = exp(r23437); return r23438; } 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 r23439, r23440, r23441, r23442, r23443, r23444, r23445, r23446, r23447, r23448, r23449, r23450, r23451, r23452; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r23439, "2", 10, MPFR_RNDN); mpfr_init(r23440); mpfr_init_set_str(r23441, "4", 10, MPFR_RNDN); mpfr_init(r23442); mpfr_init_set_str(r23443, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23444, "3", 10, MPFR_RNDN); mpfr_init(r23445); mpfr_init(r23446); mpfr_init(r23447); mpfr_init(r23448); mpfr_init(r23449); mpfr_init(r23450); mpfr_init(r23451); mpfr_init(r23452); } double f_im(double v) { ; mpfr_sqrt(r23440, r23439, MPFR_RNDN); ; mpfr_div(r23442, r23440, r23441, MPFR_RNDN); ; ; mpfr_set_d(r23445, v, MPFR_RNDN); mpfr_mul(r23446, r23445, r23445, MPFR_RNDN); mpfr_mul(r23447, r23444, r23446, MPFR_RNDN); mpfr_sub(r23448, r23443, r23447, MPFR_RNDN); mpfr_sqrt(r23449, r23448, MPFR_RNDN); mpfr_mul(r23450, r23442, r23449, MPFR_RNDN); mpfr_sub(r23451, r23443, r23446, MPFR_RNDN); mpfr_mul(r23452, r23450, r23451, MPFR_RNDN); return mpfr_get_d(r23452, MPFR_RNDN); } static mpfr_t r23453, r23454, r23455, r23456, r23457, r23458, r23459, r23460, r23461, r23462, r23463, r23464, r23465, r23466, r23467, r23468, r23469, r23470; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23453, "2", 10, MPFR_RNDN); mpfr_init(r23454); mpfr_init_set_str(r23455, "4", 10, MPFR_RNDN); mpfr_init(r23456); mpfr_init(r23457); mpfr_init_set_str(r23458, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23459, "3", 10, MPFR_RNDN); mpfr_init(r23460); mpfr_init(r23461); mpfr_init(r23462); mpfr_init(r23463); mpfr_init(r23464); mpfr_init(r23465); mpfr_init(r23466); mpfr_init(r23467); mpfr_init(r23468); mpfr_init(r23469); mpfr_init(r23470); } double f_fm(double v) { ; mpfr_sqrt(r23454, r23453, MPFR_RNDN); ; mpfr_div(r23456, r23454, r23455, MPFR_RNDN); mpfr_log(r23457, r23456, MPFR_RNDN); ; ; mpfr_set_d(r23460, v, MPFR_RNDN); mpfr_mul(r23461, r23460, r23460, MPFR_RNDN); mpfr_mul(r23462, r23459, r23461, MPFR_RNDN); mpfr_sub(r23463, r23458, r23462, MPFR_RNDN); mpfr_sqrt(r23464, r23463, MPFR_RNDN); mpfr_log(r23465, r23464, MPFR_RNDN); mpfr_add(r23466, r23457, r23465, MPFR_RNDN); mpfr_sub(r23467, r23458, r23461, MPFR_RNDN); mpfr_log(r23468, r23467, MPFR_RNDN); mpfr_add(r23469, r23466, r23468, MPFR_RNDN); mpfr_exp(r23470, r23469, MPFR_RNDN); return mpfr_get_d(r23470, MPFR_RNDN); } static mpfr_t r23471, r23472, r23473, r23474, r23475, r23476, r23477, r23478, r23479, r23480, r23481, r23482, r23483, r23484, r23485, r23486, r23487, r23488; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23471, "2", 10, MPFR_RNDN); mpfr_init(r23472); mpfr_init_set_str(r23473, "4", 10, MPFR_RNDN); mpfr_init(r23474); mpfr_init(r23475); mpfr_init_set_str(r23476, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23477, "3", 10, MPFR_RNDN); mpfr_init(r23478); mpfr_init(r23479); mpfr_init(r23480); mpfr_init(r23481); mpfr_init(r23482); mpfr_init(r23483); mpfr_init(r23484); mpfr_init(r23485); mpfr_init(r23486); mpfr_init(r23487); mpfr_init(r23488); } double f_dm(double v) { ; mpfr_sqrt(r23472, r23471, MPFR_RNDN); ; mpfr_div(r23474, r23472, r23473, MPFR_RNDN); mpfr_log(r23475, r23474, MPFR_RNDN); ; ; mpfr_set_d(r23478, v, MPFR_RNDN); mpfr_mul(r23479, r23478, r23478, MPFR_RNDN); mpfr_mul(r23480, r23477, r23479, MPFR_RNDN); mpfr_sub(r23481, r23476, r23480, MPFR_RNDN); mpfr_sqrt(r23482, r23481, MPFR_RNDN); mpfr_log(r23483, r23482, MPFR_RNDN); mpfr_add(r23484, r23475, r23483, MPFR_RNDN); mpfr_sub(r23485, r23476, r23479, MPFR_RNDN); mpfr_log(r23486, r23485, MPFR_RNDN); mpfr_add(r23487, r23484, r23486, MPFR_RNDN); mpfr_exp(r23488, r23487, MPFR_RNDN); return mpfr_get_d(r23488, MPFR_RNDN); }