#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 r22446 = 2; float r22447 = x; float r22448 = r22446 * r22447; float r22449 = cos(r22448); float r22450 = cos; float r22451 = pow(r22450, r22446); float r22452 = sin; float r22453 = pow(r22452, r22446); float r22454 = r22447 * r22453; float r22455 = r22454 * r22447; float r22456 = r22451 * r22455; float r22457 = r22449 / r22456; return r22457; } double f_id(double x, double cos, double sin) { double r22458 = 2; double r22459 = x; double r22460 = r22458 * r22459; double r22461 = cos(r22460); double r22462 = cos; double r22463 = pow(r22462, r22458); double r22464 = sin; double r22465 = pow(r22464, r22458); double r22466 = r22459 * r22465; double r22467 = r22466 * r22459; double r22468 = r22463 * r22467; double r22469 = r22461 / r22468; return r22469; } double f_of(float x, float cos, float sin) { float r22470 = 2; float r22471 = x; float r22472 = r22470 * r22471; float r22473 = cos(r22472); float r22474 = cos; float r22475 = pow(r22474, r22470); float r22476 = sqrt(r22475); float r22477 = sin; float r22478 = r22477 * r22471; float r22479 = r22478 * r22478; float r22480 = fabs(r22474); float r22481 = r22479 * r22480; float r22482 = r22476 * r22481; float r22483 = r22473 / r22482; float r22484 = 1.139874212304426e+303; bool r22485 = r22483 <= r22484; float r22486 = r22471 * r22477; float r22487 = r22474 * r22486; float r22488 = fabs(r22487); float r22489 = pow(r22488, r22470); float r22490 = r22473 / r22489; float r22491 = 1; float r22492 = r22471 * r22474; float r22493 = r22492 * r22477; float r22494 = fabs(r22493); float r22495 = r22494 * r22494; float r22496 = r22491 / r22495; float r22497 = r22473 * r22496; float r22498 = r22485 ? r22490 : r22497; return r22498; } double f_od(double x, double cos, double sin) { double r22499 = 2; double r22500 = x; double r22501 = r22499 * r22500; double r22502 = cos(r22501); double r22503 = cos; double r22504 = pow(r22503, r22499); double r22505 = sqrt(r22504); double r22506 = sin; double r22507 = r22506 * r22500; double r22508 = r22507 * r22507; double r22509 = fabs(r22503); double r22510 = r22508 * r22509; double r22511 = r22505 * r22510; double r22512 = r22502 / r22511; double r22513 = 1.139874212304426e+303; bool r22514 = r22512 <= r22513; double r22515 = r22500 * r22506; double r22516 = r22503 * r22515; double r22517 = fabs(r22516); double r22518 = pow(r22517, r22499); double r22519 = r22502 / r22518; double r22520 = 1; double r22521 = r22500 * r22503; double r22522 = r22521 * r22506; double r22523 = fabs(r22522); double r22524 = r22523 * r22523; double r22525 = r22520 / r22524; double r22526 = r22502 * r22525; double r22527 = r22514 ? r22519 : r22526; return r22527; } 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 r22528, r22529, r22530, r22531, r22532, r22533, r22534, r22535, r22536, r22537, r22538, r22539; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r22528, "2", 10, MPFR_RNDN); mpfr_init(r22529); mpfr_init(r22530); mpfr_init(r22531); mpfr_init(r22532); mpfr_init(r22533); mpfr_init(r22534); mpfr_init(r22535); mpfr_init(r22536); mpfr_init(r22537); mpfr_init(r22538); mpfr_init(r22539); } double f_im(double x, double cos, double sin) { ; mpfr_set_d(r22529, x, MPFR_RNDN); mpfr_mul(r22530, r22528, r22529, MPFR_RNDN); mpfr_cos(r22531, r22530, MPFR_RNDN); mpfr_set_d(r22532, cos, MPFR_RNDN); mpfr_pow(r22533, r22532, r22528, MPFR_RNDN); mpfr_set_d(r22534, sin, MPFR_RNDN); mpfr_pow(r22535, r22534, r22528, MPFR_RNDN); mpfr_mul(r22536, r22529, r22535, MPFR_RNDN); mpfr_mul(r22537, r22536, r22529, MPFR_RNDN); mpfr_mul(r22538, r22533, r22537, MPFR_RNDN); mpfr_div(r22539, r22531, r22538, MPFR_RNDN); return mpfr_get_d(r22539, MPFR_RNDN); } static mpfr_t r22540, r22541, r22542, r22543, r22544, r22545, r22546, r22547, r22548, r22549, r22550, r22551, r22552, r22553, r22554, r22555, r22556, r22557, r22558, r22559, r22560, r22561, r22562, r22563, r22564, r22565, r22566, r22567, r22568; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22540, "2", 10, MPFR_RNDN); mpfr_init(r22541); mpfr_init(r22542); mpfr_init(r22543); mpfr_init(r22544); mpfr_init(r22545); mpfr_init(r22546); mpfr_init(r22547); mpfr_init(r22548); mpfr_init(r22549); mpfr_init(r22550); mpfr_init(r22551); mpfr_init(r22552); mpfr_init(r22553); mpfr_init_set_str(r22554, "1.139874212304426e+303", 10, MPFR_RNDN); mpfr_init(r22555); mpfr_init(r22556); mpfr_init(r22557); mpfr_init(r22558); mpfr_init(r22559); mpfr_init(r22560); mpfr_init_set_str(r22561, "1", 10, MPFR_RNDN); mpfr_init(r22562); mpfr_init(r22563); mpfr_init(r22564); mpfr_init(r22565); mpfr_init(r22566); mpfr_init(r22567); mpfr_init(r22568); } double f_fm(double x, double cos, double sin) { ; mpfr_set_d(r22541, x, MPFR_RNDN); mpfr_mul(r22542, r22540, r22541, MPFR_RNDN); mpfr_cos(r22543, r22542, MPFR_RNDN); mpfr_set_d(r22544, cos, MPFR_RNDN); mpfr_pow(r22545, r22544, r22540, MPFR_RNDN); mpfr_sqrt(r22546, r22545, MPFR_RNDN); mpfr_set_d(r22547, sin, MPFR_RNDN); mpfr_mul(r22548, r22547, r22541, MPFR_RNDN); mpfr_mul(r22549, r22548, r22548, MPFR_RNDN); mpfr_abs(r22550, r22544, MPFR_RNDN); mpfr_mul(r22551, r22549, r22550, MPFR_RNDN); mpfr_mul(r22552, r22546, r22551, MPFR_RNDN); mpfr_div(r22553, r22543, r22552, MPFR_RNDN); ; mpfr_set_si(r22555, mpfr_cmp(r22553, r22554) <= 0, MPFR_RNDN); mpfr_mul(r22556, r22541, r22547, MPFR_RNDN); mpfr_mul(r22557, r22544, r22556, MPFR_RNDN); mpfr_abs(r22558, r22557, MPFR_RNDN); mpfr_pow(r22559, r22558, r22540, MPFR_RNDN); mpfr_div(r22560, r22543, r22559, MPFR_RNDN); ; mpfr_mul(r22562, r22541, r22544, MPFR_RNDN); mpfr_mul(r22563, r22562, r22547, MPFR_RNDN); mpfr_abs(r22564, r22563, MPFR_RNDN); mpfr_mul(r22565, r22564, r22564, MPFR_RNDN); mpfr_div(r22566, r22561, r22565, MPFR_RNDN); mpfr_mul(r22567, r22543, r22566, MPFR_RNDN); if (mpfr_get_si(r22555, MPFR_RNDN)) { mpfr_set(r22568, r22560, MPFR_RNDN); } else { mpfr_set(r22568, r22567, MPFR_RNDN); }; return mpfr_get_d(r22568, MPFR_RNDN); } static mpfr_t r22569, r22570, r22571, r22572, r22573, r22574, r22575, r22576, r22577, r22578, r22579, r22580, r22581, r22582, r22583, r22584, r22585, r22586, r22587, r22588, r22589, r22590, r22591, r22592, r22593, r22594, r22595, r22596, r22597; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22569, "2", 10, MPFR_RNDN); mpfr_init(r22570); mpfr_init(r22571); 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); mpfr_init_set_str(r22583, "1.139874212304426e+303", 10, MPFR_RNDN); mpfr_init(r22584); mpfr_init(r22585); mpfr_init(r22586); mpfr_init(r22587); mpfr_init(r22588); mpfr_init(r22589); mpfr_init_set_str(r22590, "1", 10, MPFR_RNDN); mpfr_init(r22591); mpfr_init(r22592); mpfr_init(r22593); mpfr_init(r22594); mpfr_init(r22595); mpfr_init(r22596); mpfr_init(r22597); } double f_dm(double x, double cos, double sin) { ; mpfr_set_d(r22570, x, MPFR_RNDN); mpfr_mul(r22571, r22569, r22570, MPFR_RNDN); mpfr_cos(r22572, r22571, MPFR_RNDN); mpfr_set_d(r22573, cos, MPFR_RNDN); mpfr_pow(r22574, r22573, r22569, MPFR_RNDN); mpfr_sqrt(r22575, r22574, MPFR_RNDN); mpfr_set_d(r22576, sin, MPFR_RNDN); mpfr_mul(r22577, r22576, r22570, MPFR_RNDN); mpfr_mul(r22578, r22577, r22577, MPFR_RNDN); mpfr_abs(r22579, r22573, MPFR_RNDN); mpfr_mul(r22580, r22578, r22579, MPFR_RNDN); mpfr_mul(r22581, r22575, r22580, MPFR_RNDN); mpfr_div(r22582, r22572, r22581, MPFR_RNDN); ; mpfr_set_si(r22584, mpfr_cmp(r22582, r22583) <= 0, MPFR_RNDN); mpfr_mul(r22585, r22570, r22576, MPFR_RNDN); mpfr_mul(r22586, r22573, r22585, MPFR_RNDN); mpfr_abs(r22587, r22586, MPFR_RNDN); mpfr_pow(r22588, r22587, r22569, MPFR_RNDN); mpfr_div(r22589, r22572, r22588, MPFR_RNDN); ; mpfr_mul(r22591, r22570, r22573, MPFR_RNDN); mpfr_mul(r22592, r22591, r22576, MPFR_RNDN); mpfr_abs(r22593, r22592, MPFR_RNDN); mpfr_mul(r22594, r22593, r22593, MPFR_RNDN); mpfr_div(r22595, r22590, r22594, MPFR_RNDN); mpfr_mul(r22596, r22572, r22595, MPFR_RNDN); if (mpfr_get_si(r22584, MPFR_RNDN)) { mpfr_set(r22597, r22589, MPFR_RNDN); } else { mpfr_set(r22597, r22596, MPFR_RNDN); }; return mpfr_get_d(r22597, MPFR_RNDN); }