#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r27411 = 2; float r27412 = sqrt(r27411); float r27413 = 4; float r27414 = r27412 / r27413; float r27415 = 1; float r27416 = 3; float r27417 = v; float r27418 = r27417 * r27417; float r27419 = r27416 * r27418; float r27420 = r27415 - r27419; float r27421 = sqrt(r27420); float r27422 = r27414 * r27421; float r27423 = r27415 - r27418; float r27424 = r27422 * r27423; return r27424; } double f_id(double v) { double r27425 = 2; double r27426 = sqrt(r27425); double r27427 = 4; double r27428 = r27426 / r27427; double r27429 = 1; double r27430 = 3; double r27431 = v; double r27432 = r27431 * r27431; double r27433 = r27430 * r27432; double r27434 = r27429 - r27433; double r27435 = sqrt(r27434); double r27436 = r27428 * r27435; double r27437 = r27429 - r27432; double r27438 = r27436 * r27437; return r27438; } double f_of(float v) { float r27439 = 2; float r27440 = sqrt(r27439); float r27441 = 4; float r27442 = r27440 / r27441; float r27443 = 1; float r27444 = 3; float r27445 = v; float r27446 = r27445 * r27445; float r27447 = r27444 * r27446; float r27448 = r27443 - r27447; float r27449 = sqrt(r27448); float r27450 = r27442 * r27449; float r27451 = expm1(r27450); float r27452 = log1p(r27451); float r27453 = r27443 - r27446; float r27454 = r27452 * r27453; return r27454; } double f_od(double v) { double r27455 = 2; double r27456 = sqrt(r27455); double r27457 = 4; double r27458 = r27456 / r27457; double r27459 = 1; double r27460 = 3; double r27461 = v; double r27462 = r27461 * r27461; double r27463 = r27460 * r27462; double r27464 = r27459 - r27463; double r27465 = sqrt(r27464); double r27466 = r27458 * r27465; double r27467 = expm1(r27466); double r27468 = log1p(r27467); double r27469 = r27459 - r27462; double r27470 = r27468 * r27469; return r27470; } 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 r27471, r27472, r27473, r27474, r27475, r27476, r27477, r27478, r27479, r27480, r27481, r27482, r27483, r27484; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27471, "2", 10, MPFR_RNDN); mpfr_init(r27472); mpfr_init_set_str(r27473, "4", 10, MPFR_RNDN); mpfr_init(r27474); mpfr_init_set_str(r27475, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27476, "3", 10, MPFR_RNDN); mpfr_init(r27477); mpfr_init(r27478); mpfr_init(r27479); mpfr_init(r27480); mpfr_init(r27481); mpfr_init(r27482); mpfr_init(r27483); mpfr_init(r27484); } double f_im(double v) { ; mpfr_sqrt(r27472, r27471, MPFR_RNDN); ; mpfr_div(r27474, r27472, r27473, MPFR_RNDN); ; ; mpfr_set_d(r27477, v, MPFR_RNDN); mpfr_mul(r27478, r27477, r27477, MPFR_RNDN); mpfr_mul(r27479, r27476, r27478, MPFR_RNDN); mpfr_sub(r27480, r27475, r27479, MPFR_RNDN); mpfr_sqrt(r27481, r27480, MPFR_RNDN); mpfr_mul(r27482, r27474, r27481, MPFR_RNDN); mpfr_sub(r27483, r27475, r27478, MPFR_RNDN); mpfr_mul(r27484, r27482, r27483, MPFR_RNDN); return mpfr_get_d(r27484, MPFR_RNDN); } static mpfr_t r27485, r27486, r27487, r27488, r27489, r27490, r27491, r27492, r27493, r27494, r27495, r27496, r27497, r27498, r27499, r27500; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27485, "2", 10, MPFR_RNDN); mpfr_init(r27486); mpfr_init_set_str(r27487, "4", 10, MPFR_RNDN); mpfr_init(r27488); mpfr_init_set_str(r27489, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27490, "3", 10, MPFR_RNDN); mpfr_init(r27491); mpfr_init(r27492); mpfr_init(r27493); mpfr_init(r27494); mpfr_init(r27495); mpfr_init(r27496); mpfr_init(r27497); mpfr_init(r27498); mpfr_init(r27499); mpfr_init(r27500); } double f_fm(double v) { ; mpfr_sqrt(r27486, r27485, MPFR_RNDN); ; mpfr_div(r27488, r27486, r27487, MPFR_RNDN); ; ; mpfr_set_d(r27491, v, MPFR_RNDN); mpfr_mul(r27492, r27491, r27491, MPFR_RNDN); mpfr_mul(r27493, r27490, r27492, MPFR_RNDN); mpfr_sub(r27494, r27489, r27493, MPFR_RNDN); mpfr_sqrt(r27495, r27494, MPFR_RNDN); mpfr_mul(r27496, r27488, r27495, MPFR_RNDN); mpfr_expm1(r27497, r27496, MPFR_RNDN); mpfr_log1p(r27498, r27497, MPFR_RNDN); mpfr_sub(r27499, r27489, r27492, MPFR_RNDN); mpfr_mul(r27500, r27498, r27499, MPFR_RNDN); return mpfr_get_d(r27500, MPFR_RNDN); } static mpfr_t r27501, r27502, r27503, r27504, r27505, r27506, r27507, r27508, r27509, r27510, r27511, r27512, r27513, r27514, r27515, r27516; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27501, "2", 10, MPFR_RNDN); mpfr_init(r27502); mpfr_init_set_str(r27503, "4", 10, MPFR_RNDN); mpfr_init(r27504); mpfr_init_set_str(r27505, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27506, "3", 10, MPFR_RNDN); mpfr_init(r27507); mpfr_init(r27508); mpfr_init(r27509); mpfr_init(r27510); mpfr_init(r27511); mpfr_init(r27512); mpfr_init(r27513); mpfr_init(r27514); mpfr_init(r27515); mpfr_init(r27516); } double f_dm(double v) { ; mpfr_sqrt(r27502, r27501, MPFR_RNDN); ; mpfr_div(r27504, r27502, r27503, MPFR_RNDN); ; ; mpfr_set_d(r27507, v, MPFR_RNDN); mpfr_mul(r27508, r27507, r27507, MPFR_RNDN); mpfr_mul(r27509, r27506, r27508, MPFR_RNDN); mpfr_sub(r27510, r27505, r27509, MPFR_RNDN); mpfr_sqrt(r27511, r27510, MPFR_RNDN); mpfr_mul(r27512, r27504, r27511, MPFR_RNDN); mpfr_expm1(r27513, r27512, MPFR_RNDN); mpfr_log1p(r27514, r27513, MPFR_RNDN); mpfr_sub(r27515, r27505, r27508, MPFR_RNDN); mpfr_mul(r27516, r27514, r27515, MPFR_RNDN); return mpfr_get_d(r27516, MPFR_RNDN); }