#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r10397 = 4.0f; float r10398 = 3.0f; float r10399 = atan2(1.0, 0.0); float r10400 = r10398 * r10399; float r10401 = 1.0f; float r10402 = v; float r10403 = r10402 * r10402; float r10404 = r10401 - r10403; float r10405 = r10400 * r10404; float r10406 = 2.0f; float r10407 = 6.0f; float r10408 = r10407 * r10403; float r10409 = r10406 - r10408; float r10410 = sqrt(r10409); float r10411 = r10405 * r10410; float r10412 = r10397 / r10411; return r10412; } double f_id(double v) { double r10413 = 4.0; double r10414 = 3.0; double r10415 = atan2(1.0, 0.0); double r10416 = r10414 * r10415; double r10417 = 1.0; double r10418 = v; double r10419 = r10418 * r10418; double r10420 = r10417 - r10419; double r10421 = r10416 * r10420; double r10422 = 2.0; double r10423 = 6.0; double r10424 = r10423 * r10419; double r10425 = r10422 - r10424; double r10426 = sqrt(r10425); double r10427 = r10421 * r10426; double r10428 = r10413 / r10427; return r10428; } double f_of(float v) { float r10429 = 4.0f; float r10430 = 3.0f; float r10431 = r10429 / r10430; float r10432 = atan2(1.0, 0.0); float r10433 = r10431 / r10432; float r10434 = 1.0f; float r10435 = v; float r10436 = r10435 * r10435; float r10437 = r10434 - r10436; float r10438 = r10433 / r10437; float r10439 = 6.0f; float r10440 = r10439 * r10435; float r10441 = -r10435; float r10442 = 2.0f; float r10443 = fma(r10440, r10441, r10442); float r10444 = sqrt(r10443); float r10445 = r10438 / r10444; return r10445; } double f_od(double v) { double r10446 = 4.0; double r10447 = 3.0; double r10448 = r10446 / r10447; double r10449 = atan2(1.0, 0.0); double r10450 = r10448 / r10449; double r10451 = 1.0; double r10452 = v; double r10453 = r10452 * r10452; double r10454 = r10451 - r10453; double r10455 = r10450 / r10454; double r10456 = 6.0; double r10457 = r10456 * r10452; double r10458 = -r10452; double r10459 = 2.0; double r10460 = fma(r10457, r10458, r10459); double r10461 = sqrt(r10460); double r10462 = r10455 / r10461; return r10462; } 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 r10463, r10464, r10465, r10466, r10467, r10468, r10469, r10470, r10471, r10472, r10473, r10474, r10475, r10476, r10477, r10478; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r10463, "4", 10, MPFR_RNDN); mpfr_init_set_str(r10464, "3", 10, MPFR_RNDN); mpfr_init(r10465); mpfr_init(r10466); mpfr_init_set_str(r10467, "1", 10, MPFR_RNDN); mpfr_init(r10468); mpfr_init(r10469); mpfr_init(r10470); mpfr_init(r10471); mpfr_init_set_str(r10472, "2", 10, MPFR_RNDN); mpfr_init_set_str(r10473, "6", 10, MPFR_RNDN); mpfr_init(r10474); mpfr_init(r10475); mpfr_init(r10476); mpfr_init(r10477); mpfr_init(r10478); } double f_im(double v) { ; ; mpfr_const_pi(r10465, MPFR_RNDN); mpfr_mul(r10466, r10464, r10465, MPFR_RNDN); ; mpfr_set_d(r10468, v, MPFR_RNDN); mpfr_mul(r10469, r10468, r10468, MPFR_RNDN); mpfr_sub(r10470, r10467, r10469, MPFR_RNDN); mpfr_mul(r10471, r10466, r10470, MPFR_RNDN); ; ; mpfr_mul(r10474, r10473, r10469, MPFR_RNDN); mpfr_sub(r10475, r10472, r10474, MPFR_RNDN); mpfr_sqrt(r10476, r10475, MPFR_RNDN); mpfr_mul(r10477, r10471, r10476, MPFR_RNDN); mpfr_div(r10478, r10463, r10477, MPFR_RNDN); return mpfr_get_d(r10478, MPFR_RNDN); } static mpfr_t r10479, r10480, r10481, r10482, r10483, r10484, r10485, r10486, r10487, r10488, r10489, r10490, r10491, r10492, r10493, r10494, r10495; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r10479, "4", 10, MPFR_RNDN); mpfr_init_set_str(r10480, "3", 10, MPFR_RNDN); mpfr_init(r10481); mpfr_init(r10482); mpfr_init(r10483); mpfr_init_set_str(r10484, "1", 10, MPFR_RNDN); mpfr_init(r10485); mpfr_init(r10486); mpfr_init(r10487); mpfr_init(r10488); mpfr_init_set_str(r10489, "6", 10, MPFR_RNDN); mpfr_init(r10490); mpfr_init(r10491); mpfr_init_set_str(r10492, "2", 10, MPFR_RNDN); mpfr_init(r10493); mpfr_init(r10494); mpfr_init(r10495); } double f_fm(double v) { ; ; mpfr_div(r10481, r10479, r10480, MPFR_RNDN); mpfr_const_pi(r10482, MPFR_RNDN); mpfr_div(r10483, r10481, r10482, MPFR_RNDN); ; mpfr_set_d(r10485, v, MPFR_RNDN); mpfr_mul(r10486, r10485, r10485, MPFR_RNDN); mpfr_sub(r10487, r10484, r10486, MPFR_RNDN); mpfr_div(r10488, r10483, r10487, MPFR_RNDN); ; mpfr_mul(r10490, r10489, r10485, MPFR_RNDN); mpfr_neg(r10491, r10485, MPFR_RNDN); ; mpfr_fma(r10493, r10490, r10491, r10492, MPFR_RNDN); mpfr_sqrt(r10494, r10493, MPFR_RNDN); mpfr_div(r10495, r10488, r10494, MPFR_RNDN); return mpfr_get_d(r10495, MPFR_RNDN); } static mpfr_t r10496, r10497, r10498, r10499, r10500, r10501, r10502, r10503, r10504, r10505, r10506, r10507, r10508, r10509, r10510, r10511, r10512; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r10496, "4", 10, MPFR_RNDN); mpfr_init_set_str(r10497, "3", 10, MPFR_RNDN); mpfr_init(r10498); mpfr_init(r10499); mpfr_init(r10500); mpfr_init_set_str(r10501, "1", 10, MPFR_RNDN); mpfr_init(r10502); mpfr_init(r10503); mpfr_init(r10504); mpfr_init(r10505); mpfr_init_set_str(r10506, "6", 10, MPFR_RNDN); mpfr_init(r10507); mpfr_init(r10508); mpfr_init_set_str(r10509, "2", 10, MPFR_RNDN); mpfr_init(r10510); mpfr_init(r10511); mpfr_init(r10512); } double f_dm(double v) { ; ; mpfr_div(r10498, r10496, r10497, MPFR_RNDN); mpfr_const_pi(r10499, MPFR_RNDN); mpfr_div(r10500, r10498, r10499, MPFR_RNDN); ; mpfr_set_d(r10502, v, MPFR_RNDN); mpfr_mul(r10503, r10502, r10502, MPFR_RNDN); mpfr_sub(r10504, r10501, r10503, MPFR_RNDN); mpfr_div(r10505, r10500, r10504, MPFR_RNDN); ; mpfr_mul(r10507, r10506, r10502, MPFR_RNDN); mpfr_neg(r10508, r10502, MPFR_RNDN); ; mpfr_fma(r10510, r10507, r10508, r10509, MPFR_RNDN); mpfr_sqrt(r10511, r10510, MPFR_RNDN); mpfr_div(r10512, r10505, r10511, MPFR_RNDN); return mpfr_get_d(r10512, MPFR_RNDN); }