#include #include #include #include #include char *name = "2nthrt (problem 3.4.6)"; double f_if(float x, float n) { float r25426 = x; float r25427 = 1; float r25428 = r25426 + r25427; float r25429 = n; float r25430 = r25427 / r25429; float r25431 = pow(r25428, r25430); float r25432 = pow(r25426, r25430); float r25433 = r25431 - r25432; return r25433; } double f_id(double x, double n) { double r25434 = x; double r25435 = 1; double r25436 = r25434 + r25435; double r25437 = n; double r25438 = r25435 / r25437; double r25439 = pow(r25436, r25438); double r25440 = pow(r25434, r25438); double r25441 = r25439 - r25440; return r25441; } double f_of(float x, float n) { float r25442 = x; float r25443 = 2.2675200116336946e-05; bool r25444 = r25442 <= r25443; float r25445 = 1; float r25446 = r25445 + r25442; float r25447 = n; float r25448 = r25445 / r25447; float r25449 = pow(r25446, r25448); float r25450 = r25449 - r25445; float r25451 = log(r25442); float r25452 = 1/2; float r25453 = r25452 / r25447; float r25454 = r25451 * r25453; float r25455 = r25447 / r25451; float r25456 = r25454 / r25455; float r25457 = r25450 - r25456; float r25458 = r25451 / r25447; float r25459 = r25457 - r25458; float r25460 = r25442 + r25445; float r25461 = pow(r25460, r25448); float r25462 = pow(r25442, r25448); float r25463 = r25461 - r25462; float r25464 = cbrt(r25463); float r25465 = r25464 * r25464; float r25466 = r25465 * r25464; float r25467 = r25444 ? r25459 : r25466; return r25467; } double f_od(double x, double n) { double r25468 = x; double r25469 = 2.2675200116336946e-05; bool r25470 = r25468 <= r25469; double r25471 = 1; double r25472 = r25471 + r25468; double r25473 = n; double r25474 = r25471 / r25473; double r25475 = pow(r25472, r25474); double r25476 = r25475 - r25471; double r25477 = log(r25468); double r25478 = 1/2; double r25479 = r25478 / r25473; double r25480 = r25477 * r25479; double r25481 = r25473 / r25477; double r25482 = r25480 / r25481; double r25483 = r25476 - r25482; double r25484 = r25477 / r25473; double r25485 = r25483 - r25484; double r25486 = r25468 + r25471; double r25487 = pow(r25486, r25474); double r25488 = pow(r25468, r25474); double r25489 = r25487 - r25488; double r25490 = cbrt(r25489); double r25491 = r25490 * r25490; double r25492 = r25491 * r25490; double r25493 = r25470 ? r25485 : r25492; return r25493; } 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 r25494, r25495, r25496, r25497, r25498, r25499, r25500, r25501; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25494); mpfr_init_set_str(r25495, "1", 10, MPFR_RNDN); mpfr_init(r25496); mpfr_init(r25497); mpfr_init(r25498); mpfr_init(r25499); mpfr_init(r25500); mpfr_init(r25501); } double f_im(double x, double n) { mpfr_set_d(r25494, x, MPFR_RNDN); ; mpfr_add(r25496, r25494, r25495, MPFR_RNDN); mpfr_set_d(r25497, n, MPFR_RNDN); mpfr_div(r25498, r25495, r25497, MPFR_RNDN); mpfr_pow(r25499, r25496, r25498, MPFR_RNDN); mpfr_pow(r25500, r25494, r25498, MPFR_RNDN); mpfr_sub(r25501, r25499, r25500, MPFR_RNDN); return mpfr_get_d(r25501, MPFR_RNDN); } static mpfr_t r25502, r25503, r25504, r25505, r25506, r25507, r25508, r25509, r25510, r25511, r25512, r25513, r25514, r25515, r25516, r25517, r25518, r25519, r25520, r25521, r25522, r25523, r25524, r25525, r25526, r25527; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25502); mpfr_init_set_str(r25503, "2.2675200116336946e-05", 10, MPFR_RNDN); mpfr_init(r25504); mpfr_init_set_str(r25505, "1", 10, MPFR_RNDN); mpfr_init(r25506); mpfr_init(r25507); mpfr_init(r25508); mpfr_init(r25509); mpfr_init(r25510); mpfr_init(r25511); mpfr_init_set_str(r25512, "1/2", 10, MPFR_RNDN); mpfr_init(r25513); mpfr_init(r25514); mpfr_init(r25515); mpfr_init(r25516); mpfr_init(r25517); mpfr_init(r25518); mpfr_init(r25519); mpfr_init(r25520); mpfr_init(r25521); mpfr_init(r25522); mpfr_init(r25523); mpfr_init(r25524); mpfr_init(r25525); mpfr_init(r25526); mpfr_init(r25527); } double f_fm(double x, double n) { mpfr_set_d(r25502, x, MPFR_RNDN); ; mpfr_set_si(r25504, mpfr_cmp(r25502, r25503) <= 0, MPFR_RNDN); ; mpfr_add(r25506, r25505, r25502, MPFR_RNDN); mpfr_set_d(r25507, n, MPFR_RNDN); mpfr_div(r25508, r25505, r25507, MPFR_RNDN); mpfr_pow(r25509, r25506, r25508, MPFR_RNDN); mpfr_sub(r25510, r25509, r25505, MPFR_RNDN); mpfr_log(r25511, r25502, MPFR_RNDN); ; mpfr_div(r25513, r25512, r25507, MPFR_RNDN); mpfr_mul(r25514, r25511, r25513, MPFR_RNDN); mpfr_div(r25515, r25507, r25511, MPFR_RNDN); mpfr_div(r25516, r25514, r25515, MPFR_RNDN); mpfr_sub(r25517, r25510, r25516, MPFR_RNDN); mpfr_div(r25518, r25511, r25507, MPFR_RNDN); mpfr_sub(r25519, r25517, r25518, MPFR_RNDN); mpfr_add(r25520, r25502, r25505, MPFR_RNDN); mpfr_pow(r25521, r25520, r25508, MPFR_RNDN); mpfr_pow(r25522, r25502, r25508, MPFR_RNDN); mpfr_sub(r25523, r25521, r25522, MPFR_RNDN); mpfr_cbrt(r25524, r25523, MPFR_RNDN); mpfr_mul(r25525, r25524, r25524, MPFR_RNDN); mpfr_mul(r25526, r25525, r25524, MPFR_RNDN); if (mpfr_get_si(r25504, MPFR_RNDN)) { mpfr_set(r25527, r25519, MPFR_RNDN); } else { mpfr_set(r25527, r25526, MPFR_RNDN); }; return mpfr_get_d(r25527, MPFR_RNDN); } static mpfr_t r25528, r25529, r25530, r25531, r25532, r25533, r25534, r25535, r25536, r25537, r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545, r25546, r25547, r25548, r25549, r25550, r25551, r25552, r25553; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25528); mpfr_init_set_str(r25529, "2.2675200116336946e-05", 10, MPFR_RNDN); mpfr_init(r25530); mpfr_init_set_str(r25531, "1", 10, MPFR_RNDN); mpfr_init(r25532); mpfr_init(r25533); mpfr_init(r25534); mpfr_init(r25535); mpfr_init(r25536); mpfr_init(r25537); mpfr_init_set_str(r25538, "1/2", 10, MPFR_RNDN); mpfr_init(r25539); mpfr_init(r25540); mpfr_init(r25541); mpfr_init(r25542); mpfr_init(r25543); mpfr_init(r25544); mpfr_init(r25545); mpfr_init(r25546); mpfr_init(r25547); mpfr_init(r25548); mpfr_init(r25549); mpfr_init(r25550); mpfr_init(r25551); mpfr_init(r25552); mpfr_init(r25553); } double f_dm(double x, double n) { mpfr_set_d(r25528, x, MPFR_RNDN); ; mpfr_set_si(r25530, mpfr_cmp(r25528, r25529) <= 0, MPFR_RNDN); ; mpfr_add(r25532, r25531, r25528, MPFR_RNDN); mpfr_set_d(r25533, n, MPFR_RNDN); mpfr_div(r25534, r25531, r25533, MPFR_RNDN); mpfr_pow(r25535, r25532, r25534, MPFR_RNDN); mpfr_sub(r25536, r25535, r25531, MPFR_RNDN); mpfr_log(r25537, r25528, MPFR_RNDN); ; mpfr_div(r25539, r25538, r25533, MPFR_RNDN); mpfr_mul(r25540, r25537, r25539, MPFR_RNDN); mpfr_div(r25541, r25533, r25537, MPFR_RNDN); mpfr_div(r25542, r25540, r25541, MPFR_RNDN); mpfr_sub(r25543, r25536, r25542, MPFR_RNDN); mpfr_div(r25544, r25537, r25533, MPFR_RNDN); mpfr_sub(r25545, r25543, r25544, MPFR_RNDN); mpfr_add(r25546, r25528, r25531, MPFR_RNDN); mpfr_pow(r25547, r25546, r25534, MPFR_RNDN); mpfr_pow(r25548, r25528, r25534, MPFR_RNDN); mpfr_sub(r25549, r25547, r25548, MPFR_RNDN); mpfr_cbrt(r25550, r25549, MPFR_RNDN); mpfr_mul(r25551, r25550, r25550, MPFR_RNDN); mpfr_mul(r25552, r25551, r25550, MPFR_RNDN); if (mpfr_get_si(r25530, MPFR_RNDN)) { mpfr_set(r25553, r25545, MPFR_RNDN); } else { mpfr_set(r25553, r25552, MPFR_RNDN); }; return mpfr_get_d(r25553, MPFR_RNDN); }