#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (23)"; double f_if(float F, float B, float x) { float r21475 = x; float r21476 = 1; float r21477 = B; float r21478 = tan(r21477); float r21479 = r21476 / r21478; float r21480 = r21475 * r21479; float r21481 = -r21480; float r21482 = F; float r21483 = sin(r21477); float r21484 = r21482 / r21483; float r21485 = r21482 * r21482; float r21486 = 2; float r21487 = r21485 + r21486; float r21488 = r21486 * r21475; float r21489 = r21487 + r21488; float r21490 = r21476 / r21486; float r21491 = -r21490; float r21492 = pow(r21489, r21491); float r21493 = r21484 * r21492; float r21494 = r21481 + r21493; return r21494; } double f_id(double F, double B, double x) { double r21495 = x; double r21496 = 1; double r21497 = B; double r21498 = tan(r21497); double r21499 = r21496 / r21498; double r21500 = r21495 * r21499; double r21501 = -r21500; double r21502 = F; double r21503 = sin(r21497); double r21504 = r21502 / r21503; double r21505 = r21502 * r21502; double r21506 = 2; double r21507 = r21505 + r21506; double r21508 = r21506 * r21495; double r21509 = r21507 + r21508; double r21510 = r21496 / r21506; double r21511 = -r21510; double r21512 = pow(r21509, r21511); double r21513 = r21504 * r21512; double r21514 = r21501 + r21513; return r21514; } double f_of(float F, float B, float x) { float r21515 = x; float r21516 = 2; float r21517 = F; float r21518 = fma(r21517, r21517, r21516); float r21519 = fma(r21515, r21516, r21518); float r21520 = sqrt(r21519); float r21521 = 1; float r21522 = r21521 / r21516; float r21523 = -r21522; float r21524 = pow(r21520, r21523); float r21525 = r21524 * r21524; float r21526 = B; float r21527 = sin(r21526); float r21528 = r21517 / r21527; float r21529 = -r21515; float r21530 = tan(r21526); float r21531 = r21529 / r21530; float r21532 = fma(r21525, r21528, r21531); return r21532; } double f_od(double F, double B, double x) { double r21533 = x; double r21534 = 2; double r21535 = F; double r21536 = fma(r21535, r21535, r21534); double r21537 = fma(r21533, r21534, r21536); double r21538 = sqrt(r21537); double r21539 = 1; double r21540 = r21539 / r21534; double r21541 = -r21540; double r21542 = pow(r21538, r21541); double r21543 = r21542 * r21542; double r21544 = B; double r21545 = sin(r21544); double r21546 = r21535 / r21545; double r21547 = -r21533; double r21548 = tan(r21544); double r21549 = r21547 / r21548; double r21550 = fma(r21543, r21546, r21549); return r21550; } 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 r21551, r21552, r21553, r21554, r21555, r21556, r21557, r21558, r21559, r21560, r21561, r21562, r21563, r21564, r21565, r21566, r21567, r21568, r21569, r21570; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r21551); mpfr_init_set_str(r21552, "1", 10, MPFR_RNDN); mpfr_init(r21553); mpfr_init(r21554); mpfr_init(r21555); mpfr_init(r21556); mpfr_init(r21557); mpfr_init(r21558); mpfr_init(r21559); mpfr_init(r21560); mpfr_init(r21561); mpfr_init_set_str(r21562, "2", 10, MPFR_RNDN); mpfr_init(r21563); mpfr_init(r21564); mpfr_init(r21565); mpfr_init(r21566); mpfr_init(r21567); mpfr_init(r21568); mpfr_init(r21569); mpfr_init(r21570); } double f_im(double F, double B, double x) { mpfr_set_d(r21551, x, MPFR_RNDN); ; mpfr_set_d(r21553, B, MPFR_RNDN); mpfr_tan(r21554, r21553, MPFR_RNDN); mpfr_div(r21555, r21552, r21554, MPFR_RNDN); mpfr_mul(r21556, r21551, r21555, MPFR_RNDN); mpfr_neg(r21557, r21556, MPFR_RNDN); mpfr_set_d(r21558, F, MPFR_RNDN); mpfr_sin(r21559, r21553, MPFR_RNDN); mpfr_div(r21560, r21558, r21559, MPFR_RNDN); mpfr_mul(r21561, r21558, r21558, MPFR_RNDN); ; mpfr_add(r21563, r21561, r21562, MPFR_RNDN); mpfr_mul(r21564, r21562, r21551, MPFR_RNDN); mpfr_add(r21565, r21563, r21564, MPFR_RNDN); mpfr_div(r21566, r21552, r21562, MPFR_RNDN); mpfr_neg(r21567, r21566, MPFR_RNDN); mpfr_pow(r21568, r21565, r21567, MPFR_RNDN); mpfr_mul(r21569, r21560, r21568, MPFR_RNDN); mpfr_add(r21570, r21557, r21569, MPFR_RNDN); return mpfr_get_d(r21570, MPFR_RNDN); } static mpfr_t r21571, r21572, r21573, r21574, r21575, r21576, r21577, r21578, r21579, r21580, r21581, r21582, r21583, r21584, r21585, r21586, r21587, r21588; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r21571); mpfr_init_set_str(r21572, "2", 10, MPFR_RNDN); mpfr_init(r21573); mpfr_init(r21574); mpfr_init(r21575); mpfr_init(r21576); mpfr_init_set_str(r21577, "1", 10, MPFR_RNDN); mpfr_init(r21578); mpfr_init(r21579); mpfr_init(r21580); mpfr_init(r21581); mpfr_init(r21582); mpfr_init(r21583); mpfr_init(r21584); mpfr_init(r21585); mpfr_init(r21586); mpfr_init(r21587); mpfr_init(r21588); } double f_fm(double F, double B, double x) { mpfr_set_d(r21571, x, MPFR_RNDN); ; mpfr_set_d(r21573, F, MPFR_RNDN); mpfr_fma(r21574, r21573, r21573, r21572, MPFR_RNDN); mpfr_fma(r21575, r21571, r21572, r21574, MPFR_RNDN); mpfr_sqrt(r21576, r21575, MPFR_RNDN); ; mpfr_div(r21578, r21577, r21572, MPFR_RNDN); mpfr_neg(r21579, r21578, MPFR_RNDN); mpfr_pow(r21580, r21576, r21579, MPFR_RNDN); mpfr_mul(r21581, r21580, r21580, MPFR_RNDN); mpfr_set_d(r21582, B, MPFR_RNDN); mpfr_sin(r21583, r21582, MPFR_RNDN); mpfr_div(r21584, r21573, r21583, MPFR_RNDN); mpfr_neg(r21585, r21571, MPFR_RNDN); mpfr_tan(r21586, r21582, MPFR_RNDN); mpfr_div(r21587, r21585, r21586, MPFR_RNDN); mpfr_fma(r21588, r21581, r21584, r21587, MPFR_RNDN); return mpfr_get_d(r21588, MPFR_RNDN); } static mpfr_t r21589, r21590, r21591, r21592, r21593, r21594, r21595, r21596, r21597, r21598, r21599, r21600, r21601, r21602, r21603, r21604, r21605, r21606; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r21589); mpfr_init_set_str(r21590, "2", 10, MPFR_RNDN); mpfr_init(r21591); mpfr_init(r21592); mpfr_init(r21593); mpfr_init(r21594); mpfr_init_set_str(r21595, "1", 10, MPFR_RNDN); mpfr_init(r21596); mpfr_init(r21597); mpfr_init(r21598); mpfr_init(r21599); mpfr_init(r21600); mpfr_init(r21601); mpfr_init(r21602); mpfr_init(r21603); mpfr_init(r21604); mpfr_init(r21605); mpfr_init(r21606); } double f_dm(double F, double B, double x) { mpfr_set_d(r21589, x, MPFR_RNDN); ; mpfr_set_d(r21591, F, MPFR_RNDN); mpfr_fma(r21592, r21591, r21591, r21590, MPFR_RNDN); mpfr_fma(r21593, r21589, r21590, r21592, MPFR_RNDN); mpfr_sqrt(r21594, r21593, MPFR_RNDN); ; mpfr_div(r21596, r21595, r21590, MPFR_RNDN); mpfr_neg(r21597, r21596, MPFR_RNDN); mpfr_pow(r21598, r21594, r21597, MPFR_RNDN); mpfr_mul(r21599, r21598, r21598, MPFR_RNDN); mpfr_set_d(r21600, B, MPFR_RNDN); mpfr_sin(r21601, r21600, MPFR_RNDN); mpfr_div(r21602, r21591, r21601, MPFR_RNDN); mpfr_neg(r21603, r21589, MPFR_RNDN); mpfr_tan(r21604, r21600, MPFR_RNDN); mpfr_div(r21605, r21603, r21604, MPFR_RNDN); mpfr_fma(r21606, r21599, r21602, r21605, MPFR_RNDN); return mpfr_get_d(r21606, MPFR_RNDN); }