#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r27302 = 2; float r27303 = sqrt(r27302); float r27304 = 4; float r27305 = r27303 / r27304; float r27306 = 1; float r27307 = 3; float r27308 = v; float r27309 = r27308 * r27308; float r27310 = r27307 * r27309; float r27311 = r27306 - r27310; float r27312 = sqrt(r27311); float r27313 = r27305 * r27312; float r27314 = r27306 - r27309; float r27315 = r27313 * r27314; return r27315; } double f_id(double v) { double r27316 = 2; double r27317 = sqrt(r27316); double r27318 = 4; double r27319 = r27317 / r27318; double r27320 = 1; double r27321 = 3; double r27322 = v; double r27323 = r27322 * r27322; double r27324 = r27321 * r27323; double r27325 = r27320 - r27324; double r27326 = sqrt(r27325); double r27327 = r27319 * r27326; double r27328 = r27320 - r27323; double r27329 = r27327 * r27328; return r27329; } double f_of(float v) { float r27330 = 2; float r27331 = sqrt(r27330); float r27332 = 4; float r27333 = r27331 / r27332; float r27334 = 1; float r27335 = 3; float r27336 = v; float r27337 = r27336 * r27336; float r27338 = r27335 * r27337; float r27339 = exp(r27338); float r27340 = log(r27339); float r27341 = r27334 - r27340; float r27342 = sqrt(r27341); float r27343 = r27333 * r27342; float r27344 = r27334 - r27337; float r27345 = r27343 * r27344; return r27345; } double f_od(double v) { double r27346 = 2; double r27347 = sqrt(r27346); double r27348 = 4; double r27349 = r27347 / r27348; double r27350 = 1; double r27351 = 3; double r27352 = v; double r27353 = r27352 * r27352; double r27354 = r27351 * r27353; double r27355 = exp(r27354); double r27356 = log(r27355); double r27357 = r27350 - r27356; double r27358 = sqrt(r27357); double r27359 = r27349 * r27358; double r27360 = r27350 - r27353; double r27361 = r27359 * r27360; return r27361; } 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 r27362, r27363, r27364, r27365, r27366, r27367, r27368, r27369, r27370, r27371, r27372, r27373, r27374, r27375; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r27362, "2", 10, MPFR_RNDN); mpfr_init(r27363); mpfr_init_set_str(r27364, "4", 10, MPFR_RNDN); mpfr_init(r27365); mpfr_init_set_str(r27366, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27367, "3", 10, MPFR_RNDN); mpfr_init(r27368); mpfr_init(r27369); mpfr_init(r27370); mpfr_init(r27371); mpfr_init(r27372); mpfr_init(r27373); mpfr_init(r27374); mpfr_init(r27375); } double f_im(double v) { ; mpfr_sqrt(r27363, r27362, MPFR_RNDN); ; mpfr_div(r27365, r27363, r27364, MPFR_RNDN); ; ; mpfr_set_d(r27368, v, MPFR_RNDN); mpfr_mul(r27369, r27368, r27368, MPFR_RNDN); mpfr_mul(r27370, r27367, r27369, MPFR_RNDN); mpfr_sub(r27371, r27366, r27370, MPFR_RNDN); mpfr_sqrt(r27372, r27371, MPFR_RNDN); mpfr_mul(r27373, r27365, r27372, MPFR_RNDN); mpfr_sub(r27374, r27366, r27369, MPFR_RNDN); mpfr_mul(r27375, r27373, r27374, MPFR_RNDN); return mpfr_get_d(r27375, MPFR_RNDN); } static mpfr_t r27376, r27377, r27378, r27379, r27380, r27381, r27382, r27383, r27384, r27385, r27386, r27387, r27388, r27389, r27390, r27391; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27376, "2", 10, MPFR_RNDN); mpfr_init(r27377); mpfr_init_set_str(r27378, "4", 10, MPFR_RNDN); mpfr_init(r27379); mpfr_init_set_str(r27380, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27381, "3", 10, MPFR_RNDN); mpfr_init(r27382); mpfr_init(r27383); mpfr_init(r27384); mpfr_init(r27385); mpfr_init(r27386); mpfr_init(r27387); mpfr_init(r27388); mpfr_init(r27389); mpfr_init(r27390); mpfr_init(r27391); } double f_fm(double v) { ; mpfr_sqrt(r27377, r27376, MPFR_RNDN); ; mpfr_div(r27379, r27377, r27378, MPFR_RNDN); ; ; mpfr_set_d(r27382, v, MPFR_RNDN); mpfr_mul(r27383, r27382, r27382, MPFR_RNDN); mpfr_mul(r27384, r27381, r27383, MPFR_RNDN); mpfr_exp(r27385, r27384, MPFR_RNDN); mpfr_log(r27386, r27385, MPFR_RNDN); mpfr_sub(r27387, r27380, r27386, MPFR_RNDN); mpfr_sqrt(r27388, r27387, MPFR_RNDN); mpfr_mul(r27389, r27379, r27388, MPFR_RNDN); mpfr_sub(r27390, r27380, r27383, MPFR_RNDN); mpfr_mul(r27391, r27389, r27390, MPFR_RNDN); return mpfr_get_d(r27391, MPFR_RNDN); } static mpfr_t r27392, r27393, r27394, r27395, r27396, r27397, r27398, r27399, r27400, r27401, r27402, r27403, r27404, r27405, r27406, r27407; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27392, "2", 10, MPFR_RNDN); mpfr_init(r27393); mpfr_init_set_str(r27394, "4", 10, MPFR_RNDN); mpfr_init(r27395); mpfr_init_set_str(r27396, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27397, "3", 10, MPFR_RNDN); mpfr_init(r27398); mpfr_init(r27399); mpfr_init(r27400); mpfr_init(r27401); mpfr_init(r27402); mpfr_init(r27403); mpfr_init(r27404); mpfr_init(r27405); mpfr_init(r27406); mpfr_init(r27407); } double f_dm(double v) { ; mpfr_sqrt(r27393, r27392, MPFR_RNDN); ; mpfr_div(r27395, r27393, r27394, MPFR_RNDN); ; ; mpfr_set_d(r27398, v, MPFR_RNDN); mpfr_mul(r27399, r27398, r27398, MPFR_RNDN); mpfr_mul(r27400, r27397, r27399, MPFR_RNDN); mpfr_exp(r27401, r27400, MPFR_RNDN); mpfr_log(r27402, r27401, MPFR_RNDN); mpfr_sub(r27403, r27396, r27402, MPFR_RNDN); mpfr_sqrt(r27404, r27403, MPFR_RNDN); mpfr_mul(r27405, r27395, r27404, MPFR_RNDN); mpfr_sub(r27406, r27396, r27399, MPFR_RNDN); mpfr_mul(r27407, r27405, r27406, MPFR_RNDN); return mpfr_get_d(r27407, MPFR_RNDN); }