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