#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 = sin; float r22475 = cos; float r22476 = r22474 * r22475; float r22477 = r22476 * r22471; float r22478 = fabs(r22477); float r22479 = fabs(r22478); float r22480 = r22479 * r22479; float r22481 = r22473 / r22480; float r22482 = 1.6580614592888172e+277; bool r22483 = r22481 <= r22482; float r22484 = r22471 * r22470; float r22485 = cos(r22484); float r22486 = sqrt(r22485); float r22487 = sqrt(r22473); float r22488 = r22471 * r22474; float r22489 = r22475 * r22488; float r22490 = fabs(r22489); float r22491 = pow(r22490, r22470); float r22492 = r22487 / r22491; float r22493 = r22486 * r22492; float r22494 = r22483 ? r22481 : r22493; return r22494; } double f_od(double x, double cos, double sin) { double r22495 = 2; double r22496 = x; double r22497 = r22495 * r22496; double r22498 = cos(r22497); double r22499 = sin; double r22500 = cos; double r22501 = r22499 * r22500; double r22502 = r22501 * r22496; double r22503 = fabs(r22502); double r22504 = fabs(r22503); double r22505 = r22504 * r22504; double r22506 = r22498 / r22505; double r22507 = 1.6580614592888172e+277; bool r22508 = r22506 <= r22507; double r22509 = r22496 * r22495; double r22510 = cos(r22509); double r22511 = sqrt(r22510); double r22512 = sqrt(r22498); double r22513 = r22496 * r22499; double r22514 = r22500 * r22513; double r22515 = fabs(r22514); double r22516 = pow(r22515, r22495); double r22517 = r22512 / r22516; double r22518 = r22511 * r22517; double r22519 = r22508 ? r22506 : r22518; return r22519; } 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 r22520, r22521, r22522, r22523, r22524, r22525, r22526, r22527, r22528, r22529, r22530, r22531; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r22520, "2", 10, MPFR_RNDN); mpfr_init(r22521); mpfr_init(r22522); mpfr_init(r22523); mpfr_init(r22524); mpfr_init(r22525); mpfr_init(r22526); mpfr_init(r22527); mpfr_init(r22528); mpfr_init(r22529); mpfr_init(r22530); mpfr_init(r22531); } double f_im(double x, double cos, double sin) { ; mpfr_set_d(r22521, x, MPFR_RNDN); mpfr_mul(r22522, r22520, r22521, MPFR_RNDN); mpfr_cos(r22523, r22522, MPFR_RNDN); mpfr_set_d(r22524, cos, MPFR_RNDN); mpfr_pow(r22525, r22524, r22520, MPFR_RNDN); mpfr_set_d(r22526, sin, MPFR_RNDN); mpfr_pow(r22527, r22526, r22520, MPFR_RNDN); mpfr_mul(r22528, r22521, r22527, MPFR_RNDN); mpfr_mul(r22529, r22528, r22521, MPFR_RNDN); mpfr_mul(r22530, r22525, r22529, MPFR_RNDN); mpfr_div(r22531, r22523, r22530, MPFR_RNDN); return mpfr_get_d(r22531, MPFR_RNDN); } static mpfr_t r22532, r22533, r22534, r22535, r22536, r22537, r22538, r22539, r22540, r22541, r22542, r22543, r22544, r22545, r22546, r22547, r22548, r22549, r22550, r22551, r22552, r22553, r22554, r22555, r22556; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22532, "2", 10, MPFR_RNDN); mpfr_init(r22533); mpfr_init(r22534); mpfr_init(r22535); mpfr_init(r22536); mpfr_init(r22537); mpfr_init(r22538); mpfr_init(r22539); mpfr_init(r22540); mpfr_init(r22541); mpfr_init(r22542); mpfr_init(r22543); mpfr_init_set_str(r22544, "1.6580614592888172e+277", 10, MPFR_RNDN); 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(r22554); mpfr_init(r22555); mpfr_init(r22556); } double f_fm(double x, double cos, double sin) { ; mpfr_set_d(r22533, x, MPFR_RNDN); mpfr_mul(r22534, r22532, r22533, MPFR_RNDN); mpfr_cos(r22535, r22534, MPFR_RNDN); mpfr_set_d(r22536, sin, MPFR_RNDN); mpfr_set_d(r22537, cos, MPFR_RNDN); mpfr_mul(r22538, r22536, r22537, MPFR_RNDN); mpfr_mul(r22539, r22538, r22533, MPFR_RNDN); mpfr_abs(r22540, r22539, MPFR_RNDN); mpfr_abs(r22541, r22540, MPFR_RNDN); mpfr_mul(r22542, r22541, r22541, MPFR_RNDN); mpfr_div(r22543, r22535, r22542, MPFR_RNDN); ; mpfr_set_si(r22545, mpfr_cmp(r22543, r22544) <= 0, MPFR_RNDN); mpfr_mul(r22546, r22533, r22532, MPFR_RNDN); mpfr_cos(r22547, r22546, MPFR_RNDN); mpfr_sqrt(r22548, r22547, MPFR_RNDN); mpfr_sqrt(r22549, r22535, MPFR_RNDN); mpfr_mul(r22550, r22533, r22536, MPFR_RNDN); mpfr_mul(r22551, r22537, r22550, MPFR_RNDN); mpfr_abs(r22552, r22551, MPFR_RNDN); mpfr_pow(r22553, r22552, r22532, MPFR_RNDN); mpfr_div(r22554, r22549, r22553, MPFR_RNDN); mpfr_mul(r22555, r22548, r22554, MPFR_RNDN); if (mpfr_get_si(r22545, MPFR_RNDN)) { mpfr_set(r22556, r22543, MPFR_RNDN); } else { mpfr_set(r22556, r22555, MPFR_RNDN); }; return mpfr_get_d(r22556, MPFR_RNDN); } static mpfr_t r22557, r22558, r22559, r22560, r22561, r22562, r22563, r22564, r22565, r22566, r22567, r22568, r22569, r22570, r22571, r22572, r22573, r22574, r22575, r22576, r22577, r22578, r22579, r22580, r22581; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22557, "2", 10, MPFR_RNDN); mpfr_init(r22558); mpfr_init(r22559); mpfr_init(r22560); mpfr_init(r22561); mpfr_init(r22562); mpfr_init(r22563); mpfr_init(r22564); mpfr_init(r22565); mpfr_init(r22566); mpfr_init(r22567); mpfr_init(r22568); mpfr_init_set_str(r22569, "1.6580614592888172e+277", 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); } double f_dm(double x, double cos, double sin) { ; mpfr_set_d(r22558, x, MPFR_RNDN); mpfr_mul(r22559, r22557, r22558, MPFR_RNDN); mpfr_cos(r22560, r22559, MPFR_RNDN); mpfr_set_d(r22561, sin, MPFR_RNDN); mpfr_set_d(r22562, cos, MPFR_RNDN); mpfr_mul(r22563, r22561, r22562, MPFR_RNDN); mpfr_mul(r22564, r22563, r22558, MPFR_RNDN); mpfr_abs(r22565, r22564, MPFR_RNDN); mpfr_abs(r22566, r22565, MPFR_RNDN); mpfr_mul(r22567, r22566, r22566, MPFR_RNDN); mpfr_div(r22568, r22560, r22567, MPFR_RNDN); ; mpfr_set_si(r22570, mpfr_cmp(r22568, r22569) <= 0, MPFR_RNDN); mpfr_mul(r22571, r22558, r22557, MPFR_RNDN); mpfr_cos(r22572, r22571, MPFR_RNDN); mpfr_sqrt(r22573, r22572, MPFR_RNDN); mpfr_sqrt(r22574, r22560, MPFR_RNDN); mpfr_mul(r22575, r22558, r22561, MPFR_RNDN); mpfr_mul(r22576, r22562, r22575, MPFR_RNDN); mpfr_abs(r22577, r22576, MPFR_RNDN); mpfr_pow(r22578, r22577, r22557, MPFR_RNDN); mpfr_div(r22579, r22574, r22578, MPFR_RNDN); mpfr_mul(r22580, r22573, r22579, MPFR_RNDN); if (mpfr_get_si(r22570, MPFR_RNDN)) { mpfr_set(r22581, r22568, MPFR_RNDN); } else { mpfr_set(r22581, r22580, MPFR_RNDN); }; return mpfr_get_d(r22581, MPFR_RNDN); }