#include #include #include #include #include char *name = "Migdal et al, Equation (64)"; double f_if(float a1, float a2, float th) { float r25254 = th; float r25255 = cos(r25254); float r25256 = 2; float r25257 = sqrt(r25256); float r25258 = r25255 / r25257; float r25259 = a1; float r25260 = r25259 * r25259; float r25261 = r25258 * r25260; float r25262 = a2; float r25263 = r25262 * r25262; float r25264 = r25258 * r25263; float r25265 = r25261 + r25264; return r25265; } double f_id(double a1, double a2, double th) { double r25266 = th; double r25267 = cos(r25266); double r25268 = 2; double r25269 = sqrt(r25268); double r25270 = r25267 / r25269; double r25271 = a1; double r25272 = r25271 * r25271; double r25273 = r25270 * r25272; double r25274 = a2; double r25275 = r25274 * r25274; double r25276 = r25270 * r25275; double r25277 = r25273 + r25276; return r25277; } double f_of(float a1, float a2, float th) { float r25278 = th; float r25279 = cos(r25278); float r25280 = 3; float r25281 = pow(r25279, r25280); float r25282 = 2; float r25283 = sqrt(r25282); float r25284 = r25282 * r25283; float r25285 = r25281 / r25284; float r25286 = cbrt(r25285); float r25287 = a1; float r25288 = r25287 * r25287; float r25289 = r25286 * r25288; float r25290 = r25279 / r25283; float r25291 = a2; float r25292 = r25291 * r25291; float r25293 = r25290 * r25292; float r25294 = r25289 + r25293; return r25294; } double f_od(double a1, double a2, double th) { double r25295 = th; double r25296 = cos(r25295); double r25297 = 3; double r25298 = pow(r25296, r25297); double r25299 = 2; double r25300 = sqrt(r25299); double r25301 = r25299 * r25300; double r25302 = r25298 / r25301; double r25303 = cbrt(r25302); double r25304 = a1; double r25305 = r25304 * r25304; double r25306 = r25303 * r25305; double r25307 = r25296 / r25300; double r25308 = a2; double r25309 = r25308 * r25308; double r25310 = r25307 * r25309; double r25311 = r25306 + r25310; return r25311; } 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 r25312, r25313, r25314, r25315, r25316, r25317, r25318, r25319, r25320, r25321, r25322, r25323; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25312); mpfr_init(r25313); mpfr_init_set_str(r25314, "2", 10, MPFR_RNDN); mpfr_init(r25315); mpfr_init(r25316); mpfr_init(r25317); mpfr_init(r25318); mpfr_init(r25319); mpfr_init(r25320); mpfr_init(r25321); mpfr_init(r25322); mpfr_init(r25323); } double f_im(double a1, double a2, double th) { mpfr_set_d(r25312, th, MPFR_RNDN); mpfr_cos(r25313, r25312, MPFR_RNDN); ; mpfr_sqrt(r25315, r25314, MPFR_RNDN); mpfr_div(r25316, r25313, r25315, MPFR_RNDN); mpfr_set_d(r25317, a1, MPFR_RNDN); mpfr_mul(r25318, r25317, r25317, MPFR_RNDN); mpfr_mul(r25319, r25316, r25318, MPFR_RNDN); mpfr_set_d(r25320, a2, MPFR_RNDN); mpfr_mul(r25321, r25320, r25320, MPFR_RNDN); mpfr_mul(r25322, r25316, r25321, MPFR_RNDN); mpfr_add(r25323, r25319, r25322, MPFR_RNDN); return mpfr_get_d(r25323, MPFR_RNDN); } static mpfr_t r25324, r25325, r25326, r25327, r25328, r25329, r25330, r25331, r25332, r25333, r25334, r25335, r25336, r25337, r25338, r25339, r25340; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25324); mpfr_init(r25325); mpfr_init_set_str(r25326, "3", 10, MPFR_RNDN); mpfr_init(r25327); mpfr_init_set_str(r25328, "2", 10, MPFR_RNDN); mpfr_init(r25329); mpfr_init(r25330); mpfr_init(r25331); mpfr_init(r25332); mpfr_init(r25333); mpfr_init(r25334); mpfr_init(r25335); mpfr_init(r25336); mpfr_init(r25337); mpfr_init(r25338); mpfr_init(r25339); mpfr_init(r25340); } double f_fm(double a1, double a2, double th) { mpfr_set_d(r25324, th, MPFR_RNDN); mpfr_cos(r25325, r25324, MPFR_RNDN); ; mpfr_pow(r25327, r25325, r25326, MPFR_RNDN); ; mpfr_sqrt(r25329, r25328, MPFR_RNDN); mpfr_mul(r25330, r25328, r25329, MPFR_RNDN); mpfr_div(r25331, r25327, r25330, MPFR_RNDN); mpfr_cbrt(r25332, r25331, MPFR_RNDN); mpfr_set_d(r25333, a1, MPFR_RNDN); mpfr_mul(r25334, r25333, r25333, MPFR_RNDN); mpfr_mul(r25335, r25332, r25334, MPFR_RNDN); mpfr_div(r25336, r25325, r25329, MPFR_RNDN); mpfr_set_d(r25337, a2, MPFR_RNDN); mpfr_mul(r25338, r25337, r25337, MPFR_RNDN); mpfr_mul(r25339, r25336, r25338, MPFR_RNDN); mpfr_add(r25340, r25335, r25339, MPFR_RNDN); return mpfr_get_d(r25340, MPFR_RNDN); } static mpfr_t r25341, r25342, r25343, r25344, r25345, r25346, r25347, r25348, r25349, r25350, r25351, r25352, r25353, r25354, r25355, r25356, r25357; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25341); mpfr_init(r25342); mpfr_init_set_str(r25343, "3", 10, MPFR_RNDN); mpfr_init(r25344); mpfr_init_set_str(r25345, "2", 10, MPFR_RNDN); mpfr_init(r25346); mpfr_init(r25347); mpfr_init(r25348); mpfr_init(r25349); mpfr_init(r25350); mpfr_init(r25351); mpfr_init(r25352); mpfr_init(r25353); mpfr_init(r25354); mpfr_init(r25355); mpfr_init(r25356); mpfr_init(r25357); } double f_dm(double a1, double a2, double th) { mpfr_set_d(r25341, th, MPFR_RNDN); mpfr_cos(r25342, r25341, MPFR_RNDN); ; mpfr_pow(r25344, r25342, r25343, MPFR_RNDN); ; mpfr_sqrt(r25346, r25345, MPFR_RNDN); mpfr_mul(r25347, r25345, r25346, MPFR_RNDN); mpfr_div(r25348, r25344, r25347, MPFR_RNDN); mpfr_cbrt(r25349, r25348, MPFR_RNDN); mpfr_set_d(r25350, a1, MPFR_RNDN); mpfr_mul(r25351, r25350, r25350, MPFR_RNDN); mpfr_mul(r25352, r25349, r25351, MPFR_RNDN); mpfr_div(r25353, r25342, r25346, MPFR_RNDN); mpfr_set_d(r25354, a2, MPFR_RNDN); mpfr_mul(r25355, r25354, r25354, MPFR_RNDN); mpfr_mul(r25356, r25353, r25355, MPFR_RNDN); mpfr_add(r25357, r25352, r25356, MPFR_RNDN); return mpfr_get_d(r25357, MPFR_RNDN); }