#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r24335 = 4; float r24336 = 3; float r24337 = atan2(1.0, 0.0); float r24338 = r24336 * r24337; float r24339 = 1; float r24340 = v; float r24341 = r24340 * r24340; float r24342 = r24339 - r24341; float r24343 = r24338 * r24342; float r24344 = 2; float r24345 = 6; float r24346 = r24345 * r24341; float r24347 = r24344 - r24346; float r24348 = sqrt(r24347); float r24349 = r24343 * r24348; float r24350 = r24335 / r24349; return r24350; } double f_id(double v) { double r24351 = 4; double r24352 = 3; double r24353 = atan2(1.0, 0.0); double r24354 = r24352 * r24353; double r24355 = 1; double r24356 = v; double r24357 = r24356 * r24356; double r24358 = r24355 - r24357; double r24359 = r24354 * r24358; double r24360 = 2; double r24361 = 6; double r24362 = r24361 * r24357; double r24363 = r24360 - r24362; double r24364 = sqrt(r24363); double r24365 = r24359 * r24364; double r24366 = r24351 / r24365; return r24366; } double f_of(float v) { float r24367 = 4; float r24368 = atan2(1.0, 0.0); float r24369 = 3; float r24370 = r24368 * r24369; float r24371 = r24367 / r24370; float r24372 = 1; float r24373 = v; float r24374 = r24373 * r24373; float r24375 = r24372 - r24374; float r24376 = r24371 / r24375; float r24377 = 6; float r24378 = -r24377; float r24379 = 2; float r24380 = fma(r24378, r24374, r24379); float r24381 = sqrt(r24380); float r24382 = r24376 / r24381; float r24383 = pow(r24382, r24369); float r24384 = cbrt(r24383); return r24384; } double f_od(double v) { double r24385 = 4; double r24386 = atan2(1.0, 0.0); double r24387 = 3; double r24388 = r24386 * r24387; double r24389 = r24385 / r24388; double r24390 = 1; double r24391 = v; double r24392 = r24391 * r24391; double r24393 = r24390 - r24392; double r24394 = r24389 / r24393; double r24395 = 6; double r24396 = -r24395; double r24397 = 2; double r24398 = fma(r24396, r24392, r24397); double r24399 = sqrt(r24398); double r24400 = r24394 / r24399; double r24401 = pow(r24400, r24387); double r24402 = cbrt(r24401); return r24402; } 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 r24403, r24404, r24405, r24406, r24407, r24408, r24409, r24410, r24411, r24412, r24413, r24414, r24415, r24416, r24417, r24418; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r24403, "4", 10, MPFR_RNDN); mpfr_init_set_str(r24404, "3", 10, MPFR_RNDN); mpfr_init(r24405); mpfr_init(r24406); mpfr_init_set_str(r24407, "1", 10, MPFR_RNDN); mpfr_init(r24408); mpfr_init(r24409); mpfr_init(r24410); mpfr_init(r24411); mpfr_init_set_str(r24412, "2", 10, MPFR_RNDN); mpfr_init_set_str(r24413, "6", 10, MPFR_RNDN); mpfr_init(r24414); mpfr_init(r24415); mpfr_init(r24416); mpfr_init(r24417); mpfr_init(r24418); } double f_im(double v) { ; ; mpfr_const_pi(r24405, MPFR_RNDN); mpfr_mul(r24406, r24404, r24405, MPFR_RNDN); ; mpfr_set_d(r24408, v, MPFR_RNDN); mpfr_mul(r24409, r24408, r24408, MPFR_RNDN); mpfr_sub(r24410, r24407, r24409, MPFR_RNDN); mpfr_mul(r24411, r24406, r24410, MPFR_RNDN); ; ; mpfr_mul(r24414, r24413, r24409, MPFR_RNDN); mpfr_sub(r24415, r24412, r24414, MPFR_RNDN); mpfr_sqrt(r24416, r24415, MPFR_RNDN); mpfr_mul(r24417, r24411, r24416, MPFR_RNDN); mpfr_div(r24418, r24403, r24417, MPFR_RNDN); return mpfr_get_d(r24418, MPFR_RNDN); } static mpfr_t r24419, r24420, r24421, r24422, r24423, r24424, r24425, r24426, r24427, r24428, r24429, r24430, r24431, r24432, r24433, r24434, r24435, r24436; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24419, "4", 10, MPFR_RNDN); mpfr_init(r24420); mpfr_init_set_str(r24421, "3", 10, MPFR_RNDN); mpfr_init(r24422); mpfr_init(r24423); mpfr_init_set_str(r24424, "1", 10, MPFR_RNDN); mpfr_init(r24425); mpfr_init(r24426); mpfr_init(r24427); mpfr_init(r24428); mpfr_init_set_str(r24429, "6", 10, MPFR_RNDN); mpfr_init(r24430); mpfr_init_set_str(r24431, "2", 10, MPFR_RNDN); mpfr_init(r24432); mpfr_init(r24433); mpfr_init(r24434); mpfr_init(r24435); mpfr_init(r24436); } double f_fm(double v) { ; mpfr_const_pi(r24420, MPFR_RNDN); ; mpfr_mul(r24422, r24420, r24421, MPFR_RNDN); mpfr_div(r24423, r24419, r24422, MPFR_RNDN); ; mpfr_set_d(r24425, v, MPFR_RNDN); mpfr_mul(r24426, r24425, r24425, MPFR_RNDN); mpfr_sub(r24427, r24424, r24426, MPFR_RNDN); mpfr_div(r24428, r24423, r24427, MPFR_RNDN); ; mpfr_neg(r24430, r24429, MPFR_RNDN); ; mpfr_fma(r24432, r24430, r24426, r24431, MPFR_RNDN); mpfr_sqrt(r24433, r24432, MPFR_RNDN); mpfr_div(r24434, r24428, r24433, MPFR_RNDN); mpfr_pow(r24435, r24434, r24421, MPFR_RNDN); mpfr_cbrt(r24436, r24435, MPFR_RNDN); return mpfr_get_d(r24436, MPFR_RNDN); } static mpfr_t r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446, r24447, r24448, r24449, r24450, r24451, r24452, r24453, r24454; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24437, "4", 10, MPFR_RNDN); mpfr_init(r24438); mpfr_init_set_str(r24439, "3", 10, MPFR_RNDN); mpfr_init(r24440); mpfr_init(r24441); mpfr_init_set_str(r24442, "1", 10, MPFR_RNDN); mpfr_init(r24443); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init_set_str(r24447, "6", 10, MPFR_RNDN); mpfr_init(r24448); mpfr_init_set_str(r24449, "2", 10, MPFR_RNDN); mpfr_init(r24450); mpfr_init(r24451); mpfr_init(r24452); mpfr_init(r24453); mpfr_init(r24454); } double f_dm(double v) { ; mpfr_const_pi(r24438, MPFR_RNDN); ; mpfr_mul(r24440, r24438, r24439, MPFR_RNDN); mpfr_div(r24441, r24437, r24440, MPFR_RNDN); ; mpfr_set_d(r24443, v, MPFR_RNDN); mpfr_mul(r24444, r24443, r24443, MPFR_RNDN); mpfr_sub(r24445, r24442, r24444, MPFR_RNDN); mpfr_div(r24446, r24441, r24445, MPFR_RNDN); ; mpfr_neg(r24448, r24447, MPFR_RNDN); ; mpfr_fma(r24450, r24448, r24444, r24449, MPFR_RNDN); mpfr_sqrt(r24451, r24450, MPFR_RNDN); mpfr_div(r24452, r24446, r24451, MPFR_RNDN); mpfr_pow(r24453, r24452, r24439, MPFR_RNDN); mpfr_cbrt(r24454, r24453, MPFR_RNDN); return mpfr_get_d(r24454, MPFR_RNDN); }