#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 r21549 = 2; float r21550 = x; float r21551 = r21549 * r21550; float r21552 = cos(r21551); float r21553 = cos; float r21554 = pow(r21553, r21549); float r21555 = sin; float r21556 = pow(r21555, r21549); float r21557 = r21550 * r21556; float r21558 = r21557 * r21550; float r21559 = r21554 * r21558; float r21560 = r21552 / r21559; return r21560; } double f_id(double x, double cos, double sin) { double r21561 = 2; double r21562 = x; double r21563 = r21561 * r21562; double r21564 = cos(r21563); double r21565 = cos; double r21566 = pow(r21565, r21561); double r21567 = sin; double r21568 = pow(r21567, r21561); double r21569 = r21562 * r21568; double r21570 = r21569 * r21562; double r21571 = r21566 * r21570; double r21572 = r21564 / r21571; return r21572; } double f_of(float x, float cos, float sin) { float r21573 = 1; float r21574 = cos; float r21575 = sin; float r21576 = x; float r21577 = r21575 * r21576; float r21578 = fabs(r21577); float r21579 = r21574 * r21578; float r21580 = r21573 / r21579; float r21581 = 2; float r21582 = r21581 * r21576; float r21583 = cos(r21582); float r21584 = sqrt(r21578); float r21585 = r21574 * r21584; float r21586 = r21585 * r21584; float r21587 = r21583 / r21586; float r21588 = r21580 * r21587; return r21588; } double f_od(double x, double cos, double sin) { double r21589 = 1; double r21590 = cos; double r21591 = sin; double r21592 = x; double r21593 = r21591 * r21592; double r21594 = fabs(r21593); double r21595 = r21590 * r21594; double r21596 = r21589 / r21595; double r21597 = 2; double r21598 = r21597 * r21592; double r21599 = cos(r21598); double r21600 = sqrt(r21594); double r21601 = r21590 * r21600; double r21602 = r21601 * r21600; double r21603 = r21599 / r21602; double r21604 = r21596 * r21603; return r21604; } 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 r21605, r21606, r21607, r21608, r21609, r21610, r21611, r21612, r21613, r21614, r21615, r21616; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21605, "2", 10, MPFR_RNDN); mpfr_init(r21606); mpfr_init(r21607); mpfr_init(r21608); mpfr_init(r21609); mpfr_init(r21610); mpfr_init(r21611); mpfr_init(r21612); mpfr_init(r21613); mpfr_init(r21614); mpfr_init(r21615); mpfr_init(r21616); } double f_im(double x, double cos, double sin) { ; mpfr_set_d(r21606, x, MPFR_RNDN); mpfr_mul(r21607, r21605, r21606, MPFR_RNDN); mpfr_cos(r21608, r21607, MPFR_RNDN); mpfr_set_d(r21609, cos, MPFR_RNDN); mpfr_pow(r21610, r21609, r21605, MPFR_RNDN); mpfr_set_d(r21611, sin, MPFR_RNDN); mpfr_pow(r21612, r21611, r21605, MPFR_RNDN); mpfr_mul(r21613, r21606, r21612, MPFR_RNDN); mpfr_mul(r21614, r21613, r21606, MPFR_RNDN); mpfr_mul(r21615, r21610, r21614, MPFR_RNDN); mpfr_div(r21616, r21608, r21615, MPFR_RNDN); return mpfr_get_d(r21616, MPFR_RNDN); } static mpfr_t r21617, r21618, r21619, r21620, r21621, r21622, r21623, r21624, r21625, r21626, r21627, r21628, r21629, r21630, r21631, r21632; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21617, "1", 10, MPFR_RNDN); mpfr_init(r21618); mpfr_init(r21619); mpfr_init(r21620); mpfr_init(r21621); mpfr_init(r21622); mpfr_init(r21623); mpfr_init(r21624); mpfr_init_set_str(r21625, "2", 10, MPFR_RNDN); mpfr_init(r21626); mpfr_init(r21627); mpfr_init(r21628); mpfr_init(r21629); mpfr_init(r21630); mpfr_init(r21631); mpfr_init(r21632); } double f_fm(double x, double cos, double sin) { ; mpfr_set_d(r21618, cos, MPFR_RNDN); mpfr_set_d(r21619, sin, MPFR_RNDN); mpfr_set_d(r21620, x, MPFR_RNDN); mpfr_mul(r21621, r21619, r21620, MPFR_RNDN); mpfr_abs(r21622, r21621, MPFR_RNDN); mpfr_mul(r21623, r21618, r21622, MPFR_RNDN); mpfr_div(r21624, r21617, r21623, MPFR_RNDN); ; mpfr_mul(r21626, r21625, r21620, MPFR_RNDN); mpfr_cos(r21627, r21626, MPFR_RNDN); mpfr_sqrt(r21628, r21622, MPFR_RNDN); mpfr_mul(r21629, r21618, r21628, MPFR_RNDN); mpfr_mul(r21630, r21629, r21628, MPFR_RNDN); mpfr_div(r21631, r21627, r21630, MPFR_RNDN); mpfr_mul(r21632, r21624, r21631, MPFR_RNDN); return mpfr_get_d(r21632, MPFR_RNDN); } static mpfr_t r21633, r21634, r21635, r21636, r21637, r21638, r21639, r21640, r21641, r21642, r21643, r21644, r21645, r21646, r21647, r21648; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21633, "1", 10, MPFR_RNDN); mpfr_init(r21634); mpfr_init(r21635); mpfr_init(r21636); mpfr_init(r21637); mpfr_init(r21638); mpfr_init(r21639); mpfr_init(r21640); mpfr_init_set_str(r21641, "2", 10, MPFR_RNDN); mpfr_init(r21642); mpfr_init(r21643); mpfr_init(r21644); mpfr_init(r21645); mpfr_init(r21646); mpfr_init(r21647); mpfr_init(r21648); } double f_dm(double x, double cos, double sin) { ; mpfr_set_d(r21634, cos, MPFR_RNDN); mpfr_set_d(r21635, sin, MPFR_RNDN); mpfr_set_d(r21636, x, MPFR_RNDN); mpfr_mul(r21637, r21635, r21636, MPFR_RNDN); mpfr_abs(r21638, r21637, MPFR_RNDN); mpfr_mul(r21639, r21634, r21638, MPFR_RNDN); mpfr_div(r21640, r21633, r21639, MPFR_RNDN); ; mpfr_mul(r21642, r21641, r21636, MPFR_RNDN); mpfr_cos(r21643, r21642, MPFR_RNDN); mpfr_sqrt(r21644, r21638, MPFR_RNDN); mpfr_mul(r21645, r21634, r21644, MPFR_RNDN); mpfr_mul(r21646, r21645, r21644, MPFR_RNDN); mpfr_div(r21647, r21643, r21646, MPFR_RNDN); mpfr_mul(r21648, r21640, r21647, MPFR_RNDN); return mpfr_get_d(r21648, MPFR_RNDN); }