#include #include #include #include #include char *name = "cos(2*x)/(cos^2(x)*sin^2(x))"; double f_if(float x, float cos, float sin) { float r22495 = 2; float r22496 = x; float r22497 = r22495 * r22496; float r22498 = cos(r22497); float r22499 = cos; float r22500 = pow(r22499, r22495); float r22501 = sin; float r22502 = pow(r22501, r22495); float r22503 = r22496 * r22502; float r22504 = r22503 * r22496; float r22505 = r22500 * r22504; float r22506 = r22498 / r22505; return r22506; } double f_id(double x, double cos, double sin) { double r22507 = 2; double r22508 = x; double r22509 = r22507 * r22508; double r22510 = cos(r22509); double r22511 = cos; double r22512 = pow(r22511, r22507); double r22513 = sin; double r22514 = pow(r22513, r22507); double r22515 = r22508 * r22514; double r22516 = r22515 * r22508; double r22517 = r22512 * r22516; double r22518 = r22510 / r22517; return r22518; } double f_of(float x, float cos, float sin) { float r22519 = 2; float r22520 = x; float r22521 = r22519 * r22520; float r22522 = cos(r22521); float r22523 = cos; float r22524 = pow(r22523, r22519); float r22525 = sin; float r22526 = pow(r22525, r22519); float r22527 = r22520 * r22526; float r22528 = r22527 * r22520; float r22529 = r22524 * r22528; float r22530 = r22522 / r22529; float r22531 = +inf.0; bool r22532 = r22530 <= r22531; float r22533 = r22520 * r22525; float r22534 = r22523 * r22533; float r22535 = fabs(r22534); float r22536 = pow(r22535, r22519); float r22537 = r22522 / r22536; float r22538 = r22525 * r22523; float r22539 = r22538 * r22520; float r22540 = fabs(r22539); float r22541 = fabs(r22540); float r22542 = r22541 * r22541; float r22543 = r22522 / r22542; float r22544 = r22532 ? r22537 : r22543; return r22544; } double f_od(double x, double cos, double sin) { double r22545 = 2; double r22546 = x; double r22547 = r22545 * r22546; double r22548 = cos(r22547); double r22549 = cos; double r22550 = pow(r22549, r22545); double r22551 = sin; double r22552 = pow(r22551, r22545); double r22553 = r22546 * r22552; double r22554 = r22553 * r22546; double r22555 = r22550 * r22554; double r22556 = r22548 / r22555; double r22557 = +inf.0; bool r22558 = r22556 <= r22557; double r22559 = r22546 * r22551; double r22560 = r22549 * r22559; double r22561 = fabs(r22560); double r22562 = pow(r22561, r22545); double r22563 = r22548 / r22562; double r22564 = r22551 * r22549; double r22565 = r22564 * r22546; double r22566 = fabs(r22565); double r22567 = fabs(r22566); double r22568 = r22567 * r22567; double r22569 = r22548 / r22568; double r22570 = r22558 ? r22563 : r22569; return r22570; } 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 r22571, r22572, r22573, r22574, r22575, r22576, r22577, r22578, r22579, r22580, r22581, r22582; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22571, "2", 10, MPFR_RNDN); mpfr_init(r22572); mpfr_init(r22573); mpfr_init(r22574); mpfr_init(r22575); mpfr_init(r22576); mpfr_init(r22577); mpfr_init(r22578); mpfr_init(r22579); mpfr_init(r22580); mpfr_init(r22581); mpfr_init(r22582); } double f_im(double x, double cos, double sin) { ; mpfr_set_d(r22572, x, MPFR_RNDN); mpfr_mul(r22573, r22571, r22572, MPFR_RNDN); mpfr_cos(r22574, r22573, MPFR_RNDN); mpfr_set_d(r22575, cos, MPFR_RNDN); mpfr_pow(r22576, r22575, r22571, MPFR_RNDN); mpfr_set_d(r22577, sin, MPFR_RNDN); mpfr_pow(r22578, r22577, r22571, MPFR_RNDN); mpfr_mul(r22579, r22572, r22578, MPFR_RNDN); mpfr_mul(r22580, r22579, r22572, MPFR_RNDN); mpfr_mul(r22581, r22576, r22580, MPFR_RNDN); mpfr_div(r22582, r22574, r22581, MPFR_RNDN); return mpfr_get_d(r22582, MPFR_RNDN); } static mpfr_t r22583, r22584, r22585, r22586, r22587, r22588, r22589, r22590, r22591, r22592, r22593, r22594, r22595, r22596, r22597, r22598, r22599, r22600, r22601, r22602, r22603, r22604, r22605, r22606, r22607, r22608; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22583, "2", 10, MPFR_RNDN); mpfr_init(r22584); mpfr_init(r22585); mpfr_init(r22586); mpfr_init(r22587); mpfr_init(r22588); mpfr_init(r22589); mpfr_init(r22590); mpfr_init(r22591); mpfr_init(r22592); mpfr_init(r22593); mpfr_init(r22594); mpfr_init_set_str(r22595, "+inf.0", 10, MPFR_RNDN); mpfr_init(r22596); mpfr_init(r22597); mpfr_init(r22598); mpfr_init(r22599); mpfr_init(r22600); mpfr_init(r22601); mpfr_init(r22602); mpfr_init(r22603); mpfr_init(r22604); mpfr_init(r22605); mpfr_init(r22606); mpfr_init(r22607); mpfr_init(r22608); } double f_fm(double x, double cos, double sin) { ; mpfr_set_d(r22584, x, MPFR_RNDN); mpfr_mul(r22585, r22583, r22584, MPFR_RNDN); mpfr_cos(r22586, r22585, MPFR_RNDN); mpfr_set_d(r22587, cos, MPFR_RNDN); mpfr_pow(r22588, r22587, r22583, MPFR_RNDN); mpfr_set_d(r22589, sin, MPFR_RNDN); mpfr_pow(r22590, r22589, r22583, MPFR_RNDN); mpfr_mul(r22591, r22584, r22590, MPFR_RNDN); mpfr_mul(r22592, r22591, r22584, MPFR_RNDN); mpfr_mul(r22593, r22588, r22592, MPFR_RNDN); mpfr_div(r22594, r22586, r22593, MPFR_RNDN); ; mpfr_set_si(r22596, mpfr_cmp(r22594, r22595) <= 0, MPFR_RNDN); mpfr_mul(r22597, r22584, r22589, MPFR_RNDN); mpfr_mul(r22598, r22587, r22597, MPFR_RNDN); mpfr_abs(r22599, r22598, MPFR_RNDN); mpfr_pow(r22600, r22599, r22583, MPFR_RNDN); mpfr_div(r22601, r22586, r22600, MPFR_RNDN); mpfr_mul(r22602, r22589, r22587, MPFR_RNDN); mpfr_mul(r22603, r22602, r22584, MPFR_RNDN); mpfr_abs(r22604, r22603, MPFR_RNDN); mpfr_abs(r22605, r22604, MPFR_RNDN); mpfr_mul(r22606, r22605, r22605, MPFR_RNDN); mpfr_div(r22607, r22586, r22606, MPFR_RNDN); if (mpfr_get_si(r22596, MPFR_RNDN)) { mpfr_set(r22608, r22601, MPFR_RNDN); } else { mpfr_set(r22608, r22607, MPFR_RNDN); }; return mpfr_get_d(r22608, MPFR_RNDN); } static mpfr_t r22609, r22610, r22611, r22612, r22613, r22614, r22615, r22616, r22617, r22618, r22619, r22620, r22621, r22622, r22623, r22624, r22625, r22626, r22627, r22628, r22629, r22630, r22631, r22632, r22633, r22634; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22609, "2", 10, MPFR_RNDN); mpfr_init(r22610); mpfr_init(r22611); mpfr_init(r22612); mpfr_init(r22613); mpfr_init(r22614); mpfr_init(r22615); mpfr_init(r22616); mpfr_init(r22617); mpfr_init(r22618); mpfr_init(r22619); mpfr_init(r22620); mpfr_init_set_str(r22621, "+inf.0", 10, MPFR_RNDN); mpfr_init(r22622); mpfr_init(r22623); mpfr_init(r22624); mpfr_init(r22625); mpfr_init(r22626); mpfr_init(r22627); mpfr_init(r22628); mpfr_init(r22629); mpfr_init(r22630); mpfr_init(r22631); mpfr_init(r22632); mpfr_init(r22633); mpfr_init(r22634); } double f_dm(double x, double cos, double sin) { ; mpfr_set_d(r22610, x, MPFR_RNDN); mpfr_mul(r22611, r22609, r22610, MPFR_RNDN); mpfr_cos(r22612, r22611, MPFR_RNDN); mpfr_set_d(r22613, cos, MPFR_RNDN); mpfr_pow(r22614, r22613, r22609, MPFR_RNDN); mpfr_set_d(r22615, sin, MPFR_RNDN); mpfr_pow(r22616, r22615, r22609, MPFR_RNDN); mpfr_mul(r22617, r22610, r22616, MPFR_RNDN); mpfr_mul(r22618, r22617, r22610, MPFR_RNDN); mpfr_mul(r22619, r22614, r22618, MPFR_RNDN); mpfr_div(r22620, r22612, r22619, MPFR_RNDN); ; mpfr_set_si(r22622, mpfr_cmp(r22620, r22621) <= 0, MPFR_RNDN); mpfr_mul(r22623, r22610, r22615, MPFR_RNDN); mpfr_mul(r22624, r22613, r22623, MPFR_RNDN); mpfr_abs(r22625, r22624, MPFR_RNDN); mpfr_pow(r22626, r22625, r22609, MPFR_RNDN); mpfr_div(r22627, r22612, r22626, MPFR_RNDN); mpfr_mul(r22628, r22615, r22613, MPFR_RNDN); mpfr_mul(r22629, r22628, r22610, MPFR_RNDN); mpfr_abs(r22630, r22629, MPFR_RNDN); mpfr_abs(r22631, r22630, MPFR_RNDN); mpfr_mul(r22632, r22631, r22631, MPFR_RNDN); mpfr_div(r22633, r22612, r22632, MPFR_RNDN); if (mpfr_get_si(r22622, MPFR_RNDN)) { mpfr_set(r22634, r22627, MPFR_RNDN); } else { mpfr_set(r22634, r22633, MPFR_RNDN); }; return mpfr_get_d(r22634, MPFR_RNDN); }