#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r27372 = 2; float r27373 = sqrt(r27372); float r27374 = 4; float r27375 = r27373 / r27374; float r27376 = 1; float r27377 = 3; float r27378 = v; float r27379 = r27378 * r27378; float r27380 = r27377 * r27379; float r27381 = r27376 - r27380; float r27382 = sqrt(r27381); float r27383 = r27375 * r27382; float r27384 = r27376 - r27379; float r27385 = r27383 * r27384; return r27385; } double f_id(double v) { double r27386 = 2; double r27387 = sqrt(r27386); double r27388 = 4; double r27389 = r27387 / r27388; double r27390 = 1; double r27391 = 3; double r27392 = v; double r27393 = r27392 * r27392; double r27394 = r27391 * r27393; double r27395 = r27390 - r27394; double r27396 = sqrt(r27395); double r27397 = r27389 * r27396; double r27398 = r27390 - r27393; double r27399 = r27397 * r27398; return r27399; } double f_of(float v) { float r27400 = 2; float r27401 = sqrt(r27400); float r27402 = 4; float r27403 = r27401 / r27402; float r27404 = log(r27403); float r27405 = 1; float r27406 = 3; float r27407 = v; float r27408 = r27407 * r27407; float r27409 = r27406 * r27408; float r27410 = r27405 - r27409; float r27411 = sqrt(r27410); float r27412 = log(r27411); float r27413 = r27404 + r27412; float r27414 = exp(r27413); float r27415 = r27405 - r27408; float r27416 = r27414 * r27415; return r27416; } double f_od(double v) { double r27417 = 2; double r27418 = sqrt(r27417); double r27419 = 4; double r27420 = r27418 / r27419; double r27421 = log(r27420); double r27422 = 1; double r27423 = 3; double r27424 = v; double r27425 = r27424 * r27424; double r27426 = r27423 * r27425; double r27427 = r27422 - r27426; double r27428 = sqrt(r27427); double r27429 = log(r27428); double r27430 = r27421 + r27429; double r27431 = exp(r27430); double r27432 = r27422 - r27425; double r27433 = r27431 * r27432; return r27433; } 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 r27434, r27435, r27436, r27437, r27438, r27439, r27440, r27441, r27442, r27443, r27444, r27445, r27446, r27447; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r27434, "2", 10, MPFR_RNDN); mpfr_init(r27435); mpfr_init_set_str(r27436, "4", 10, MPFR_RNDN); mpfr_init(r27437); mpfr_init_set_str(r27438, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27439, "3", 10, MPFR_RNDN); mpfr_init(r27440); mpfr_init(r27441); mpfr_init(r27442); mpfr_init(r27443); mpfr_init(r27444); mpfr_init(r27445); mpfr_init(r27446); mpfr_init(r27447); } double f_im(double v) { ; mpfr_sqrt(r27435, r27434, MPFR_RNDN); ; mpfr_div(r27437, r27435, r27436, MPFR_RNDN); ; ; mpfr_set_d(r27440, v, MPFR_RNDN); mpfr_mul(r27441, r27440, r27440, MPFR_RNDN); mpfr_mul(r27442, r27439, r27441, MPFR_RNDN); mpfr_sub(r27443, r27438, r27442, MPFR_RNDN); mpfr_sqrt(r27444, r27443, MPFR_RNDN); mpfr_mul(r27445, r27437, r27444, MPFR_RNDN); mpfr_sub(r27446, r27438, r27441, MPFR_RNDN); mpfr_mul(r27447, r27445, r27446, MPFR_RNDN); return mpfr_get_d(r27447, MPFR_RNDN); } static mpfr_t r27448, r27449, r27450, r27451, r27452, r27453, r27454, r27455, r27456, r27457, r27458, r27459, r27460, r27461, r27462, r27463, r27464; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27448, "2", 10, MPFR_RNDN); mpfr_init(r27449); mpfr_init_set_str(r27450, "4", 10, MPFR_RNDN); mpfr_init(r27451); mpfr_init(r27452); mpfr_init_set_str(r27453, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27454, "3", 10, MPFR_RNDN); mpfr_init(r27455); mpfr_init(r27456); mpfr_init(r27457); mpfr_init(r27458); mpfr_init(r27459); mpfr_init(r27460); mpfr_init(r27461); mpfr_init(r27462); mpfr_init(r27463); mpfr_init(r27464); } double f_fm(double v) { ; mpfr_sqrt(r27449, r27448, MPFR_RNDN); ; mpfr_div(r27451, r27449, r27450, MPFR_RNDN); mpfr_log(r27452, r27451, MPFR_RNDN); ; ; mpfr_set_d(r27455, v, MPFR_RNDN); mpfr_mul(r27456, r27455, r27455, MPFR_RNDN); mpfr_mul(r27457, r27454, r27456, MPFR_RNDN); mpfr_sub(r27458, r27453, r27457, MPFR_RNDN); mpfr_sqrt(r27459, r27458, MPFR_RNDN); mpfr_log(r27460, r27459, MPFR_RNDN); mpfr_add(r27461, r27452, r27460, MPFR_RNDN); mpfr_exp(r27462, r27461, MPFR_RNDN); mpfr_sub(r27463, r27453, r27456, MPFR_RNDN); mpfr_mul(r27464, r27462, r27463, MPFR_RNDN); return mpfr_get_d(r27464, MPFR_RNDN); } static mpfr_t r27465, r27466, r27467, r27468, r27469, r27470, r27471, r27472, r27473, r27474, r27475, r27476, r27477, r27478, r27479, r27480, r27481; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27465, "2", 10, MPFR_RNDN); mpfr_init(r27466); mpfr_init_set_str(r27467, "4", 10, MPFR_RNDN); mpfr_init(r27468); mpfr_init(r27469); mpfr_init_set_str(r27470, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27471, "3", 10, MPFR_RNDN); mpfr_init(r27472); mpfr_init(r27473); mpfr_init(r27474); mpfr_init(r27475); mpfr_init(r27476); mpfr_init(r27477); mpfr_init(r27478); mpfr_init(r27479); mpfr_init(r27480); mpfr_init(r27481); } double f_dm(double v) { ; mpfr_sqrt(r27466, r27465, MPFR_RNDN); ; mpfr_div(r27468, r27466, r27467, MPFR_RNDN); mpfr_log(r27469, r27468, MPFR_RNDN); ; ; mpfr_set_d(r27472, v, MPFR_RNDN); mpfr_mul(r27473, r27472, r27472, MPFR_RNDN); mpfr_mul(r27474, r27471, r27473, MPFR_RNDN); mpfr_sub(r27475, r27470, r27474, MPFR_RNDN); mpfr_sqrt(r27476, r27475, MPFR_RNDN); mpfr_log(r27477, r27476, MPFR_RNDN); mpfr_add(r27478, r27469, r27477, MPFR_RNDN); mpfr_exp(r27479, r27478, MPFR_RNDN); mpfr_sub(r27480, r27470, r27473, MPFR_RNDN); mpfr_mul(r27481, r27479, r27480, MPFR_RNDN); return mpfr_get_d(r27481, MPFR_RNDN); }