#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r25453 = J; float r25454 = l; float r25455 = exp(r25454); float r25456 = -r25454; float r25457 = exp(r25456); float r25458 = r25455 - r25457; float r25459 = r25453 * r25458; float r25460 = K; float r25461 = 2; float r25462 = r25460 / r25461; float r25463 = cos(r25462); float r25464 = r25459 * r25463; float r25465 = U; float r25466 = r25464 + r25465; return r25466; } double f_id(double J, double l, double K, double U) { double r25467 = J; double r25468 = l; double r25469 = exp(r25468); double r25470 = -r25468; double r25471 = exp(r25470); double r25472 = r25469 - r25471; double r25473 = r25467 * r25472; double r25474 = K; double r25475 = 2; double r25476 = r25474 / r25475; double r25477 = cos(r25476); double r25478 = r25473 * r25477; double r25479 = U; double r25480 = r25478 + r25479; return r25480; } double f_of(float J, float l, float K, float U) { float r25481 = K; float r25482 = 2; float r25483 = r25481 / r25482; float r25484 = cos(r25483); float r25485 = J; float r25486 = r25484 * r25485; float r25487 = l; float r25488 = 1/3; float r25489 = r25487 * r25488; float r25490 = fma(r25489, r25487, r25482); float r25491 = 5; float r25492 = pow(r25487, r25491); float r25493 = 1/60; float r25494 = r25492 * r25493; float r25495 = fma(r25487, r25490, r25494); float r25496 = r25486 * r25495; float r25497 = U; float r25498 = r25496 + r25497; return r25498; } double f_od(double J, double l, double K, double U) { double r25499 = K; double r25500 = 2; double r25501 = r25499 / r25500; double r25502 = cos(r25501); double r25503 = J; double r25504 = r25502 * r25503; double r25505 = l; double r25506 = 1/3; double r25507 = r25505 * r25506; double r25508 = fma(r25507, r25505, r25500); double r25509 = 5; double r25510 = pow(r25505, r25509); double r25511 = 1/60; double r25512 = r25510 * r25511; double r25513 = fma(r25505, r25508, r25512); double r25514 = r25504 * r25513; double r25515 = U; double r25516 = r25514 + r25515; return r25516; } 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 r25517, r25518, r25519, r25520, r25521, r25522, r25523, r25524, r25525, r25526, r25527, r25528, r25529, r25530; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); 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_set_str(r25525, "2", 10, MPFR_RNDN); mpfr_init(r25526); mpfr_init(r25527); mpfr_init(r25528); mpfr_init(r25529); mpfr_init(r25530); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r25517, J, MPFR_RNDN); mpfr_set_d(r25518, l, MPFR_RNDN); mpfr_exp(r25519, r25518, MPFR_RNDN); mpfr_neg(r25520, r25518, MPFR_RNDN); mpfr_exp(r25521, r25520, MPFR_RNDN); mpfr_sub(r25522, r25519, r25521, MPFR_RNDN); mpfr_mul(r25523, r25517, r25522, MPFR_RNDN); mpfr_set_d(r25524, K, MPFR_RNDN); ; mpfr_div(r25526, r25524, r25525, MPFR_RNDN); mpfr_cos(r25527, r25526, MPFR_RNDN); mpfr_mul(r25528, r25523, r25527, MPFR_RNDN); mpfr_set_d(r25529, U, MPFR_RNDN); mpfr_add(r25530, r25528, r25529, MPFR_RNDN); return mpfr_get_d(r25530, MPFR_RNDN); } static mpfr_t r25531, r25532, r25533, r25534, r25535, r25536, r25537, r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545, r25546, r25547, r25548; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25531); mpfr_init_set_str(r25532, "2", 10, MPFR_RNDN); mpfr_init(r25533); mpfr_init(r25534); mpfr_init(r25535); mpfr_init(r25536); mpfr_init(r25537); mpfr_init_set_str(r25538, "1/3", 10, MPFR_RNDN); mpfr_init(r25539); mpfr_init(r25540); mpfr_init_set_str(r25541, "5", 10, MPFR_RNDN); mpfr_init(r25542); mpfr_init_set_str(r25543, "1/60", 10, MPFR_RNDN); mpfr_init(r25544); mpfr_init(r25545); mpfr_init(r25546); mpfr_init(r25547); mpfr_init(r25548); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r25531, K, MPFR_RNDN); ; mpfr_div(r25533, r25531, r25532, MPFR_RNDN); mpfr_cos(r25534, r25533, MPFR_RNDN); mpfr_set_d(r25535, J, MPFR_RNDN); mpfr_mul(r25536, r25534, r25535, MPFR_RNDN); mpfr_set_d(r25537, l, MPFR_RNDN); ; mpfr_mul(r25539, r25537, r25538, MPFR_RNDN); mpfr_fma(r25540, r25539, r25537, r25532, MPFR_RNDN); ; mpfr_pow(r25542, r25537, r25541, MPFR_RNDN); ; mpfr_mul(r25544, r25542, r25543, MPFR_RNDN); mpfr_fma(r25545, r25537, r25540, r25544, MPFR_RNDN); mpfr_mul(r25546, r25536, r25545, MPFR_RNDN); mpfr_set_d(r25547, U, MPFR_RNDN); mpfr_add(r25548, r25546, r25547, MPFR_RNDN); return mpfr_get_d(r25548, MPFR_RNDN); } static mpfr_t r25549, r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564, r25565, r25566; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25549); mpfr_init_set_str(r25550, "2", 10, MPFR_RNDN); mpfr_init(r25551); mpfr_init(r25552); mpfr_init(r25553); mpfr_init(r25554); mpfr_init(r25555); mpfr_init_set_str(r25556, "1/3", 10, MPFR_RNDN); mpfr_init(r25557); mpfr_init(r25558); mpfr_init_set_str(r25559, "5", 10, MPFR_RNDN); mpfr_init(r25560); mpfr_init_set_str(r25561, "1/60", 10, MPFR_RNDN); mpfr_init(r25562); mpfr_init(r25563); mpfr_init(r25564); mpfr_init(r25565); mpfr_init(r25566); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r25549, K, MPFR_RNDN); ; mpfr_div(r25551, r25549, r25550, MPFR_RNDN); mpfr_cos(r25552, r25551, MPFR_RNDN); mpfr_set_d(r25553, J, MPFR_RNDN); mpfr_mul(r25554, r25552, r25553, MPFR_RNDN); mpfr_set_d(r25555, l, MPFR_RNDN); ; mpfr_mul(r25557, r25555, r25556, MPFR_RNDN); mpfr_fma(r25558, r25557, r25555, r25550, MPFR_RNDN); ; mpfr_pow(r25560, r25555, r25559, MPFR_RNDN); ; mpfr_mul(r25562, r25560, r25561, MPFR_RNDN); mpfr_fma(r25563, r25555, r25558, r25562, MPFR_RNDN); mpfr_mul(r25564, r25554, r25563, MPFR_RNDN); mpfr_set_d(r25565, U, MPFR_RNDN); mpfr_add(r25566, r25564, r25565, MPFR_RNDN); return mpfr_get_d(r25566, MPFR_RNDN); }