#include #include #include #include #include char *name = "Migdal et al, Equation (64)"; double f_if(float a1, float a2, float th) { float r25462 = th; float r25463 = cos(r25462); float r25464 = 2; float r25465 = sqrt(r25464); float r25466 = r25463 / r25465; float r25467 = a1; float r25468 = r25467 * r25467; float r25469 = r25466 * r25468; float r25470 = a2; float r25471 = r25470 * r25470; float r25472 = r25466 * r25471; float r25473 = r25469 + r25472; return r25473; } double f_id(double a1, double a2, double th) { double r25474 = th; double r25475 = cos(r25474); double r25476 = 2; double r25477 = sqrt(r25476); double r25478 = r25475 / r25477; double r25479 = a1; double r25480 = r25479 * r25479; double r25481 = r25478 * r25480; double r25482 = a2; double r25483 = r25482 * r25482; double r25484 = r25478 * r25483; double r25485 = r25481 + r25484; return r25485; } double f_of(float a1, float a2, float th) { float r25486 = th; float r25487 = cos(r25486); float r25488 = a1; float r25489 = 2; float r25490 = sqrt(r25489); float r25491 = cbrt(r25490); float r25492 = r25491 * r25491; float r25493 = r25488 / r25492; float r25494 = r25488 / r25491; float r25495 = r25493 * r25494; float r25496 = r25487 * r25495; float r25497 = r25487 / r25490; float r25498 = a2; float r25499 = r25498 * r25498; float r25500 = r25497 * r25499; float r25501 = r25496 + r25500; return r25501; } double f_od(double a1, double a2, double th) { double r25502 = th; double r25503 = cos(r25502); double r25504 = a1; double r25505 = 2; double r25506 = sqrt(r25505); double r25507 = cbrt(r25506); double r25508 = r25507 * r25507; double r25509 = r25504 / r25508; double r25510 = r25504 / r25507; double r25511 = r25509 * r25510; double r25512 = r25503 * r25511; double r25513 = r25503 / r25506; double r25514 = a2; double r25515 = r25514 * r25514; double r25516 = r25513 * r25515; double r25517 = r25512 + r25516; return r25517; } 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 r25518, r25519, r25520, r25521, r25522, r25523, r25524, r25525, r25526, r25527, r25528, r25529; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25518); mpfr_init(r25519); mpfr_init_set_str(r25520, "2", 10, MPFR_RNDN); mpfr_init(r25521); mpfr_init(r25522); mpfr_init(r25523); mpfr_init(r25524); mpfr_init(r25525); mpfr_init(r25526); mpfr_init(r25527); mpfr_init(r25528); mpfr_init(r25529); } double f_im(double a1, double a2, double th) { mpfr_set_d(r25518, th, MPFR_RNDN); mpfr_cos(r25519, r25518, MPFR_RNDN); ; mpfr_sqrt(r25521, r25520, MPFR_RNDN); mpfr_div(r25522, r25519, r25521, MPFR_RNDN); mpfr_set_d(r25523, a1, MPFR_RNDN); mpfr_mul(r25524, r25523, r25523, MPFR_RNDN); mpfr_mul(r25525, r25522, r25524, MPFR_RNDN); mpfr_set_d(r25526, a2, MPFR_RNDN); mpfr_mul(r25527, r25526, r25526, MPFR_RNDN); mpfr_mul(r25528, r25522, r25527, MPFR_RNDN); mpfr_add(r25529, r25525, r25528, MPFR_RNDN); return mpfr_get_d(r25529, MPFR_RNDN); } static mpfr_t r25530, r25531, r25532, r25533, r25534, r25535, r25536, r25537, r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25530); mpfr_init(r25531); mpfr_init(r25532); mpfr_init_set_str(r25533, "2", 10, MPFR_RNDN); mpfr_init(r25534); mpfr_init(r25535); mpfr_init(r25536); mpfr_init(r25537); mpfr_init(r25538); mpfr_init(r25539); mpfr_init(r25540); mpfr_init(r25541); mpfr_init(r25542); mpfr_init(r25543); mpfr_init(r25544); mpfr_init(r25545); } double f_fm(double a1, double a2, double th) { mpfr_set_d(r25530, th, MPFR_RNDN); mpfr_cos(r25531, r25530, MPFR_RNDN); mpfr_set_d(r25532, a1, MPFR_RNDN); ; mpfr_sqrt(r25534, r25533, MPFR_RNDN); mpfr_cbrt(r25535, r25534, MPFR_RNDN); mpfr_mul(r25536, r25535, r25535, MPFR_RNDN); mpfr_div(r25537, r25532, r25536, MPFR_RNDN); mpfr_div(r25538, r25532, r25535, MPFR_RNDN); mpfr_mul(r25539, r25537, r25538, MPFR_RNDN); mpfr_mul(r25540, r25531, r25539, MPFR_RNDN); mpfr_div(r25541, r25531, r25534, MPFR_RNDN); mpfr_set_d(r25542, a2, MPFR_RNDN); mpfr_mul(r25543, r25542, r25542, MPFR_RNDN); mpfr_mul(r25544, r25541, r25543, MPFR_RNDN); mpfr_add(r25545, r25540, r25544, MPFR_RNDN); return mpfr_get_d(r25545, MPFR_RNDN); } static mpfr_t r25546, r25547, r25548, r25549, r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25546); mpfr_init(r25547); mpfr_init(r25548); mpfr_init_set_str(r25549, "2", 10, MPFR_RNDN); mpfr_init(r25550); mpfr_init(r25551); mpfr_init(r25552); mpfr_init(r25553); mpfr_init(r25554); mpfr_init(r25555); mpfr_init(r25556); mpfr_init(r25557); mpfr_init(r25558); mpfr_init(r25559); mpfr_init(r25560); mpfr_init(r25561); } double f_dm(double a1, double a2, double th) { mpfr_set_d(r25546, th, MPFR_RNDN); mpfr_cos(r25547, r25546, MPFR_RNDN); mpfr_set_d(r25548, a1, MPFR_RNDN); ; mpfr_sqrt(r25550, r25549, MPFR_RNDN); mpfr_cbrt(r25551, r25550, MPFR_RNDN); mpfr_mul(r25552, r25551, r25551, MPFR_RNDN); mpfr_div(r25553, r25548, r25552, MPFR_RNDN); mpfr_div(r25554, r25548, r25551, MPFR_RNDN); mpfr_mul(r25555, r25553, r25554, MPFR_RNDN); mpfr_mul(r25556, r25547, r25555, MPFR_RNDN); mpfr_div(r25557, r25547, r25550, MPFR_RNDN); mpfr_set_d(r25558, a2, MPFR_RNDN); mpfr_mul(r25559, r25558, r25558, MPFR_RNDN); mpfr_mul(r25560, r25557, r25559, MPFR_RNDN); mpfr_add(r25561, r25556, r25560, MPFR_RNDN); return mpfr_get_d(r25561, MPFR_RNDN); }