#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10-)"; double f_if(float t, float l, float k) { float r25406 = 2; float r25407 = t; float r25408 = 3; float r25409 = pow(r25407, r25408); float r25410 = l; float r25411 = r25410 * r25410; float r25412 = r25409 / r25411; float r25413 = k; float r25414 = sin(r25413); float r25415 = r25412 * r25414; float r25416 = tan(r25413); float r25417 = r25415 * r25416; float r25418 = 1; float r25419 = r25413 / r25407; float r25420 = pow(r25419, r25406); float r25421 = r25418 + r25420; float r25422 = r25421 - r25418; float r25423 = r25417 * r25422; float r25424 = r25406 / r25423; return r25424; } double f_id(double t, double l, double k) { double r25425 = 2; double r25426 = t; double r25427 = 3; double r25428 = pow(r25426, r25427); double r25429 = l; double r25430 = r25429 * r25429; double r25431 = r25428 / r25430; double r25432 = k; double r25433 = sin(r25432); double r25434 = r25431 * r25433; double r25435 = tan(r25432); double r25436 = r25434 * r25435; double r25437 = 1; double r25438 = r25432 / r25426; double r25439 = pow(r25438, r25425); double r25440 = r25437 + r25439; double r25441 = r25440 - r25437; double r25442 = r25436 * r25441; double r25443 = r25425 / r25442; return r25443; } double f_of(float t, float l, float k) { float r25444 = 2; float r25445 = 1; float r25446 = r25444 / r25445; float r25447 = k; float r25448 = l; float r25449 = r25447 / r25448; float r25450 = sin(r25447); float r25451 = r25449 * r25450; float r25452 = t; float r25453 = cos(r25447); float r25454 = r25452 / r25453; float r25455 = r25451 * r25454; float r25456 = r25451 * r25455; float r25457 = r25446 / r25456; return r25457; } double f_od(double t, double l, double k) { double r25458 = 2; double r25459 = 1; double r25460 = r25458 / r25459; double r25461 = k; double r25462 = l; double r25463 = r25461 / r25462; double r25464 = sin(r25461); double r25465 = r25463 * r25464; double r25466 = t; double r25467 = cos(r25461); double r25468 = r25466 / r25467; double r25469 = r25465 * r25468; double r25470 = r25465 * r25469; double r25471 = r25460 / r25470; return r25471; } 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 r25472, r25473, r25474, r25475, r25476, r25477, r25478, r25479, r25480, r25481, r25482, r25483, r25484, r25485, r25486, r25487, r25488, r25489, r25490; void setup_mpfr_f_im() { mpfr_set_default_prec(4176); mpfr_init_set_str(r25472, "2", 10, MPFR_RNDN); mpfr_init(r25473); mpfr_init_set_str(r25474, "3", 10, MPFR_RNDN); mpfr_init(r25475); mpfr_init(r25476); mpfr_init(r25477); mpfr_init(r25478); mpfr_init(r25479); mpfr_init(r25480); mpfr_init(r25481); mpfr_init(r25482); mpfr_init(r25483); mpfr_init_set_str(r25484, "1", 10, MPFR_RNDN); mpfr_init(r25485); mpfr_init(r25486); mpfr_init(r25487); mpfr_init(r25488); mpfr_init(r25489); mpfr_init(r25490); } double f_im(double t, double l, double k) { ; mpfr_set_d(r25473, t, MPFR_RNDN); ; mpfr_pow(r25475, r25473, r25474, MPFR_RNDN); mpfr_set_d(r25476, l, MPFR_RNDN); mpfr_mul(r25477, r25476, r25476, MPFR_RNDN); mpfr_div(r25478, r25475, r25477, MPFR_RNDN); mpfr_set_d(r25479, k, MPFR_RNDN); mpfr_sin(r25480, r25479, MPFR_RNDN); mpfr_mul(r25481, r25478, r25480, MPFR_RNDN); mpfr_tan(r25482, r25479, MPFR_RNDN); mpfr_mul(r25483, r25481, r25482, MPFR_RNDN); ; mpfr_div(r25485, r25479, r25473, MPFR_RNDN); mpfr_pow(r25486, r25485, r25472, MPFR_RNDN); mpfr_add(r25487, r25484, r25486, MPFR_RNDN); mpfr_sub(r25488, r25487, r25484, MPFR_RNDN); mpfr_mul(r25489, r25483, r25488, MPFR_RNDN); mpfr_div(r25490, r25472, r25489, MPFR_RNDN); return mpfr_get_d(r25490, MPFR_RNDN); } static mpfr_t r25491, r25492, r25493, r25494, r25495, r25496, r25497, r25498, r25499, r25500, r25501, r25502, r25503, r25504; void setup_mpfr_f_fm() { mpfr_set_default_prec(4176); mpfr_init_set_str(r25491, "2", 10, MPFR_RNDN); mpfr_init_set_str(r25492, "1", 10, MPFR_RNDN); mpfr_init(r25493); mpfr_init(r25494); mpfr_init(r25495); mpfr_init(r25496); mpfr_init(r25497); mpfr_init(r25498); mpfr_init(r25499); mpfr_init(r25500); mpfr_init(r25501); mpfr_init(r25502); mpfr_init(r25503); mpfr_init(r25504); } double f_fm(double t, double l, double k) { ; ; mpfr_div(r25493, r25491, r25492, MPFR_RNDN); mpfr_set_d(r25494, k, MPFR_RNDN); mpfr_set_d(r25495, l, MPFR_RNDN); mpfr_div(r25496, r25494, r25495, MPFR_RNDN); mpfr_sin(r25497, r25494, MPFR_RNDN); mpfr_mul(r25498, r25496, r25497, MPFR_RNDN); mpfr_set_d(r25499, t, MPFR_RNDN); mpfr_cos(r25500, r25494, MPFR_RNDN); mpfr_div(r25501, r25499, r25500, MPFR_RNDN); mpfr_mul(r25502, r25498, r25501, MPFR_RNDN); mpfr_mul(r25503, r25498, r25502, MPFR_RNDN); mpfr_div(r25504, r25493, r25503, MPFR_RNDN); return mpfr_get_d(r25504, MPFR_RNDN); } static mpfr_t r25505, r25506, r25507, r25508, r25509, r25510, r25511, r25512, r25513, r25514, r25515, r25516, r25517, r25518; void setup_mpfr_f_dm() { mpfr_set_default_prec(4176); mpfr_init_set_str(r25505, "2", 10, MPFR_RNDN); mpfr_init_set_str(r25506, "1", 10, MPFR_RNDN); mpfr_init(r25507); mpfr_init(r25508); mpfr_init(r25509); mpfr_init(r25510); mpfr_init(r25511); mpfr_init(r25512); mpfr_init(r25513); mpfr_init(r25514); mpfr_init(r25515); mpfr_init(r25516); mpfr_init(r25517); mpfr_init(r25518); } double f_dm(double t, double l, double k) { ; ; mpfr_div(r25507, r25505, r25506, MPFR_RNDN); mpfr_set_d(r25508, k, MPFR_RNDN); mpfr_set_d(r25509, l, MPFR_RNDN); mpfr_div(r25510, r25508, r25509, MPFR_RNDN); mpfr_sin(r25511, r25508, MPFR_RNDN); mpfr_mul(r25512, r25510, r25511, MPFR_RNDN); mpfr_set_d(r25513, t, MPFR_RNDN); mpfr_cos(r25514, r25508, MPFR_RNDN); mpfr_div(r25515, r25513, r25514, MPFR_RNDN); mpfr_mul(r25516, r25512, r25515, MPFR_RNDN); mpfr_mul(r25517, r25512, r25516, MPFR_RNDN); mpfr_div(r25518, r25507, r25517, MPFR_RNDN); return mpfr_get_d(r25518, MPFR_RNDN); }