#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 r22345 = 2; float r22346 = x; float r22347 = r22345 * r22346; float r22348 = cos(r22347); float r22349 = cos; float r22350 = pow(r22349, r22345); float r22351 = sin; float r22352 = pow(r22351, r22345); float r22353 = r22346 * r22352; float r22354 = r22353 * r22346; float r22355 = r22350 * r22354; float r22356 = r22348 / r22355; return r22356; } double f_id(double x, double cos, double sin) { double r22357 = 2; double r22358 = x; double r22359 = r22357 * r22358; double r22360 = cos(r22359); double r22361 = cos; double r22362 = pow(r22361, r22357); double r22363 = sin; double r22364 = pow(r22363, r22357); double r22365 = r22358 * r22364; double r22366 = r22365 * r22358; double r22367 = r22362 * r22366; double r22368 = r22360 / r22367; return r22368; } double f_of(float x, float cos, float sin) { float r22369 = 2; float r22370 = x; float r22371 = r22369 * r22370; float r22372 = cos(r22371); float r22373 = cos; float r22374 = pow(r22373, r22369); float r22375 = sin; float r22376 = pow(r22375, r22369); float r22377 = r22370 * r22376; float r22378 = r22374 * r22377; float r22379 = r22378 * r22370; float r22380 = r22372 / r22379; float r22381 = +inf.0; bool r22382 = r22380 <= r22381; float r22383 = r22370 * r22375; float r22384 = r22373 * r22383; float r22385 = fabs(r22384); float r22386 = pow(r22385, r22369); float r22387 = r22372 / r22386; float r22388 = r22370 * r22373; float r22389 = r22388 * r22375; float r22390 = fabs(r22389); float r22391 = r22373 * r22375; float r22392 = r22370 * r22391; float r22393 = fabs(r22392); float r22394 = r22390 * r22393; float r22395 = r22372 / r22394; float r22396 = r22382 ? r22387 : r22395; return r22396; } double f_od(double x, double cos, double sin) { double r22397 = 2; double r22398 = x; double r22399 = r22397 * r22398; double r22400 = cos(r22399); double r22401 = cos; double r22402 = pow(r22401, r22397); double r22403 = sin; double r22404 = pow(r22403, r22397); double r22405 = r22398 * r22404; double r22406 = r22402 * r22405; double r22407 = r22406 * r22398; double r22408 = r22400 / r22407; double r22409 = +inf.0; bool r22410 = r22408 <= r22409; double r22411 = r22398 * r22403; double r22412 = r22401 * r22411; double r22413 = fabs(r22412); double r22414 = pow(r22413, r22397); double r22415 = r22400 / r22414; double r22416 = r22398 * r22401; double r22417 = r22416 * r22403; double r22418 = fabs(r22417); double r22419 = r22401 * r22403; double r22420 = r22398 * r22419; double r22421 = fabs(r22420); double r22422 = r22418 * r22421; double r22423 = r22400 / r22422; double r22424 = r22410 ? r22415 : r22423; return r22424; } 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 r22425, r22426, r22427, r22428, r22429, r22430, r22431, r22432, r22433, r22434, r22435, r22436; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r22425, "2", 10, MPFR_RNDN); mpfr_init(r22426); mpfr_init(r22427); mpfr_init(r22428); mpfr_init(r22429); mpfr_init(r22430); mpfr_init(r22431); mpfr_init(r22432); mpfr_init(r22433); mpfr_init(r22434); mpfr_init(r22435); mpfr_init(r22436); } double f_im(double x, double cos, double sin) { ; mpfr_set_d(r22426, x, MPFR_RNDN); mpfr_mul(r22427, r22425, r22426, MPFR_RNDN); mpfr_cos(r22428, r22427, MPFR_RNDN); mpfr_set_d(r22429, cos, MPFR_RNDN); mpfr_pow(r22430, r22429, r22425, MPFR_RNDN); mpfr_set_d(r22431, sin, MPFR_RNDN); mpfr_pow(r22432, r22431, r22425, MPFR_RNDN); mpfr_mul(r22433, r22426, r22432, MPFR_RNDN); mpfr_mul(r22434, r22433, r22426, MPFR_RNDN); mpfr_mul(r22435, r22430, r22434, MPFR_RNDN); mpfr_div(r22436, r22428, r22435, MPFR_RNDN); return mpfr_get_d(r22436, MPFR_RNDN); } static mpfr_t r22437, r22438, r22439, r22440, r22441, r22442, r22443, r22444, r22445, r22446, r22447, r22448, r22449, r22450, r22451, r22452, r22453, r22454, r22455, r22456, r22457, r22458, r22459, r22460, r22461, r22462, r22463, r22464; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22437, "2", 10, MPFR_RNDN); mpfr_init(r22438); mpfr_init(r22439); mpfr_init(r22440); mpfr_init(r22441); mpfr_init(r22442); mpfr_init(r22443); mpfr_init(r22444); mpfr_init(r22445); mpfr_init(r22446); mpfr_init(r22447); mpfr_init(r22448); mpfr_init_set_str(r22449, "+inf.0", 10, MPFR_RNDN); mpfr_init(r22450); mpfr_init(r22451); mpfr_init(r22452); mpfr_init(r22453); mpfr_init(r22454); mpfr_init(r22455); mpfr_init(r22456); mpfr_init(r22457); mpfr_init(r22458); mpfr_init(r22459); mpfr_init(r22460); mpfr_init(r22461); mpfr_init(r22462); mpfr_init(r22463); mpfr_init(r22464); } double f_fm(double x, double cos, double sin) { ; mpfr_set_d(r22438, x, MPFR_RNDN); mpfr_mul(r22439, r22437, r22438, MPFR_RNDN); mpfr_cos(r22440, r22439, MPFR_RNDN); mpfr_set_d(r22441, cos, MPFR_RNDN); mpfr_pow(r22442, r22441, r22437, MPFR_RNDN); mpfr_set_d(r22443, sin, MPFR_RNDN); mpfr_pow(r22444, r22443, r22437, MPFR_RNDN); mpfr_mul(r22445, r22438, r22444, MPFR_RNDN); mpfr_mul(r22446, r22442, r22445, MPFR_RNDN); mpfr_mul(r22447, r22446, r22438, MPFR_RNDN); mpfr_div(r22448, r22440, r22447, MPFR_RNDN); ; mpfr_set_si(r22450, mpfr_cmp(r22448, r22449) <= 0, MPFR_RNDN); mpfr_mul(r22451, r22438, r22443, MPFR_RNDN); mpfr_mul(r22452, r22441, r22451, MPFR_RNDN); mpfr_abs(r22453, r22452, MPFR_RNDN); mpfr_pow(r22454, r22453, r22437, MPFR_RNDN); mpfr_div(r22455, r22440, r22454, MPFR_RNDN); mpfr_mul(r22456, r22438, r22441, MPFR_RNDN); mpfr_mul(r22457, r22456, r22443, MPFR_RNDN); mpfr_abs(r22458, r22457, MPFR_RNDN); mpfr_mul(r22459, r22441, r22443, MPFR_RNDN); mpfr_mul(r22460, r22438, r22459, MPFR_RNDN); mpfr_abs(r22461, r22460, MPFR_RNDN); mpfr_mul(r22462, r22458, r22461, MPFR_RNDN); mpfr_div(r22463, r22440, r22462, MPFR_RNDN); if (mpfr_get_si(r22450, MPFR_RNDN)) { mpfr_set(r22464, r22455, MPFR_RNDN); } else { mpfr_set(r22464, r22463, MPFR_RNDN); }; return mpfr_get_d(r22464, MPFR_RNDN); } static mpfr_t r22465, r22466, r22467, r22468, r22469, r22470, r22471, r22472, r22473, r22474, r22475, r22476, r22477, r22478, r22479, r22480, r22481, r22482, r22483, r22484, r22485, r22486, r22487, r22488, r22489, r22490, r22491, r22492; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22465, "2", 10, MPFR_RNDN); mpfr_init(r22466); mpfr_init(r22467); mpfr_init(r22468); mpfr_init(r22469); mpfr_init(r22470); mpfr_init(r22471); mpfr_init(r22472); mpfr_init(r22473); mpfr_init(r22474); mpfr_init(r22475); mpfr_init(r22476); mpfr_init_set_str(r22477, "+inf.0", 10, MPFR_RNDN); mpfr_init(r22478); mpfr_init(r22479); mpfr_init(r22480); mpfr_init(r22481); mpfr_init(r22482); mpfr_init(r22483); mpfr_init(r22484); mpfr_init(r22485); mpfr_init(r22486); mpfr_init(r22487); mpfr_init(r22488); mpfr_init(r22489); mpfr_init(r22490); mpfr_init(r22491); mpfr_init(r22492); } double f_dm(double x, double cos, double sin) { ; mpfr_set_d(r22466, x, MPFR_RNDN); mpfr_mul(r22467, r22465, r22466, MPFR_RNDN); mpfr_cos(r22468, r22467, MPFR_RNDN); mpfr_set_d(r22469, cos, MPFR_RNDN); mpfr_pow(r22470, r22469, r22465, MPFR_RNDN); mpfr_set_d(r22471, sin, MPFR_RNDN); mpfr_pow(r22472, r22471, r22465, MPFR_RNDN); mpfr_mul(r22473, r22466, r22472, MPFR_RNDN); mpfr_mul(r22474, r22470, r22473, MPFR_RNDN); mpfr_mul(r22475, r22474, r22466, MPFR_RNDN); mpfr_div(r22476, r22468, r22475, MPFR_RNDN); ; mpfr_set_si(r22478, mpfr_cmp(r22476, r22477) <= 0, MPFR_RNDN); mpfr_mul(r22479, r22466, r22471, MPFR_RNDN); mpfr_mul(r22480, r22469, r22479, MPFR_RNDN); mpfr_abs(r22481, r22480, MPFR_RNDN); mpfr_pow(r22482, r22481, r22465, MPFR_RNDN); mpfr_div(r22483, r22468, r22482, MPFR_RNDN); mpfr_mul(r22484, r22466, r22469, MPFR_RNDN); mpfr_mul(r22485, r22484, r22471, MPFR_RNDN); mpfr_abs(r22486, r22485, MPFR_RNDN); mpfr_mul(r22487, r22469, r22471, MPFR_RNDN); mpfr_mul(r22488, r22466, r22487, MPFR_RNDN); mpfr_abs(r22489, r22488, MPFR_RNDN); mpfr_mul(r22490, r22486, r22489, MPFR_RNDN); mpfr_div(r22491, r22468, r22490, MPFR_RNDN); if (mpfr_get_si(r22478, MPFR_RNDN)) { mpfr_set(r22492, r22483, MPFR_RNDN); } else { mpfr_set(r22492, r22491, MPFR_RNDN); }; return mpfr_get_d(r22492, MPFR_RNDN); }