#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r22018 = 1; float r22019 = 2; float r22020 = r22018 / r22019; float r22021 = l; float r22022 = r22019 * r22021; float r22023 = Om; float r22024 = r22022 / r22023; float r22025 = pow(r22024, r22019); float r22026 = kx; float r22027 = sin(r22026); float r22028 = pow(r22027, r22019); float r22029 = ky; float r22030 = sin(r22029); float r22031 = pow(r22030, r22019); float r22032 = r22028 + r22031; float r22033 = r22025 * r22032; float r22034 = r22018 + r22033; float r22035 = sqrt(r22034); float r22036 = r22018 / r22035; float r22037 = r22018 + r22036; float r22038 = r22020 * r22037; float r22039 = sqrt(r22038); return r22039; } double f_id(double l, double Om, double kx, double ky) { double r22040 = 1; double r22041 = 2; double r22042 = r22040 / r22041; double r22043 = l; double r22044 = r22041 * r22043; double r22045 = Om; double r22046 = r22044 / r22045; double r22047 = pow(r22046, r22041); double r22048 = kx; double r22049 = sin(r22048); double r22050 = pow(r22049, r22041); double r22051 = ky; double r22052 = sin(r22051); double r22053 = pow(r22052, r22041); double r22054 = r22050 + r22053; double r22055 = r22047 * r22054; double r22056 = r22040 + r22055; double r22057 = sqrt(r22056); double r22058 = r22040 / r22057; double r22059 = r22040 + r22058; double r22060 = r22042 * r22059; double r22061 = sqrt(r22060); return r22061; } double f_of(float l, float Om, float kx, float ky) { float r22062 = 1; float r22063 = 2; float r22064 = r22062 / r22063; float r22065 = l; float r22066 = Om; float r22067 = r22066 / r22063; float r22068 = r22065 / r22067; float r22069 = kx; float r22070 = sin(r22069); float r22071 = ky; float r22072 = sin(r22071); float r22073 = hypot(r22070, r22072); float r22074 = r22068 * r22073; float r22075 = hypot(r22062, r22074); float r22076 = expm1(r22075); float r22077 = log1p(r22076); float r22078 = r22062 / r22077; float r22079 = r22062 + r22078; float r22080 = r22064 * r22079; float r22081 = sqrt(r22080); return r22081; } double f_od(double l, double Om, double kx, double ky) { double r22082 = 1; double r22083 = 2; double r22084 = r22082 / r22083; double r22085 = l; double r22086 = Om; double r22087 = r22086 / r22083; double r22088 = r22085 / r22087; double r22089 = kx; double r22090 = sin(r22089); double r22091 = ky; double r22092 = sin(r22091); double r22093 = hypot(r22090, r22092); double r22094 = r22088 * r22093; double r22095 = hypot(r22082, r22094); double r22096 = expm1(r22095); double r22097 = log1p(r22096); double r22098 = r22082 / r22097; double r22099 = r22082 + r22098; double r22100 = r22084 * r22099; double r22101 = sqrt(r22100); return r22101; } 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 r22102, r22103, r22104, r22105, r22106, r22107, r22108, r22109, r22110, r22111, r22112, r22113, r22114, r22115, r22116, r22117, r22118, r22119, r22120, r22121, r22122, r22123; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r22102, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22103, "2", 10, MPFR_RNDN); mpfr_init(r22104); mpfr_init(r22105); mpfr_init(r22106); mpfr_init(r22107); mpfr_init(r22108); mpfr_init(r22109); mpfr_init(r22110); mpfr_init(r22111); mpfr_init(r22112); mpfr_init(r22113); mpfr_init(r22114); mpfr_init(r22115); mpfr_init(r22116); mpfr_init(r22117); mpfr_init(r22118); mpfr_init(r22119); mpfr_init(r22120); mpfr_init(r22121); mpfr_init(r22122); mpfr_init(r22123); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22104, r22102, r22103, MPFR_RNDN); mpfr_set_d(r22105, l, MPFR_RNDN); mpfr_mul(r22106, r22103, r22105, MPFR_RNDN); mpfr_set_d(r22107, Om, MPFR_RNDN); mpfr_div(r22108, r22106, r22107, MPFR_RNDN); mpfr_pow(r22109, r22108, r22103, MPFR_RNDN); mpfr_set_d(r22110, kx, MPFR_RNDN); mpfr_sin(r22111, r22110, MPFR_RNDN); mpfr_pow(r22112, r22111, r22103, MPFR_RNDN); mpfr_set_d(r22113, ky, MPFR_RNDN); mpfr_sin(r22114, r22113, MPFR_RNDN); mpfr_pow(r22115, r22114, r22103, MPFR_RNDN); mpfr_add(r22116, r22112, r22115, MPFR_RNDN); mpfr_mul(r22117, r22109, r22116, MPFR_RNDN); mpfr_add(r22118, r22102, r22117, MPFR_RNDN); mpfr_sqrt(r22119, r22118, MPFR_RNDN); mpfr_div(r22120, r22102, r22119, MPFR_RNDN); mpfr_add(r22121, r22102, r22120, MPFR_RNDN); mpfr_mul(r22122, r22104, r22121, MPFR_RNDN); mpfr_sqrt(r22123, r22122, MPFR_RNDN); return mpfr_get_d(r22123, MPFR_RNDN); } static mpfr_t r22124, r22125, r22126, r22127, r22128, r22129, r22130, r22131, r22132, r22133, r22134, r22135, r22136, r22137, r22138, r22139, r22140, r22141, r22142, r22143; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r22124, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22125, "2", 10, MPFR_RNDN); mpfr_init(r22126); mpfr_init(r22127); mpfr_init(r22128); mpfr_init(r22129); mpfr_init(r22130); mpfr_init(r22131); mpfr_init(r22132); mpfr_init(r22133); mpfr_init(r22134); mpfr_init(r22135); mpfr_init(r22136); mpfr_init(r22137); mpfr_init(r22138); mpfr_init(r22139); mpfr_init(r22140); mpfr_init(r22141); mpfr_init(r22142); mpfr_init(r22143); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22126, r22124, r22125, MPFR_RNDN); mpfr_set_d(r22127, l, MPFR_RNDN); mpfr_set_d(r22128, Om, MPFR_RNDN); mpfr_div(r22129, r22128, r22125, MPFR_RNDN); mpfr_div(r22130, r22127, r22129, MPFR_RNDN); mpfr_set_d(r22131, kx, MPFR_RNDN); mpfr_sin(r22132, r22131, MPFR_RNDN); mpfr_set_d(r22133, ky, MPFR_RNDN); mpfr_sin(r22134, r22133, MPFR_RNDN); mpfr_hypot(r22135, r22132, r22134, MPFR_RNDN); mpfr_mul(r22136, r22130, r22135, MPFR_RNDN); mpfr_hypot(r22137, r22124, r22136, MPFR_RNDN); mpfr_expm1(r22138, r22137, MPFR_RNDN); mpfr_log1p(r22139, r22138, MPFR_RNDN); mpfr_div(r22140, r22124, r22139, MPFR_RNDN); mpfr_add(r22141, r22124, r22140, MPFR_RNDN); mpfr_mul(r22142, r22126, r22141, MPFR_RNDN); mpfr_sqrt(r22143, r22142, MPFR_RNDN); return mpfr_get_d(r22143, MPFR_RNDN); } static mpfr_t r22144, r22145, r22146, r22147, r22148, r22149, r22150, r22151, r22152, r22153, r22154, r22155, r22156, r22157, r22158, r22159, r22160, r22161, r22162, r22163; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r22144, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22145, "2", 10, MPFR_RNDN); mpfr_init(r22146); mpfr_init(r22147); mpfr_init(r22148); mpfr_init(r22149); mpfr_init(r22150); mpfr_init(r22151); mpfr_init(r22152); mpfr_init(r22153); mpfr_init(r22154); mpfr_init(r22155); mpfr_init(r22156); mpfr_init(r22157); mpfr_init(r22158); mpfr_init(r22159); mpfr_init(r22160); mpfr_init(r22161); mpfr_init(r22162); mpfr_init(r22163); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22146, r22144, r22145, MPFR_RNDN); mpfr_set_d(r22147, l, MPFR_RNDN); mpfr_set_d(r22148, Om, MPFR_RNDN); mpfr_div(r22149, r22148, r22145, MPFR_RNDN); mpfr_div(r22150, r22147, r22149, MPFR_RNDN); mpfr_set_d(r22151, kx, MPFR_RNDN); mpfr_sin(r22152, r22151, MPFR_RNDN); mpfr_set_d(r22153, ky, MPFR_RNDN); mpfr_sin(r22154, r22153, MPFR_RNDN); mpfr_hypot(r22155, r22152, r22154, MPFR_RNDN); mpfr_mul(r22156, r22150, r22155, MPFR_RNDN); mpfr_hypot(r22157, r22144, r22156, MPFR_RNDN); mpfr_expm1(r22158, r22157, MPFR_RNDN); mpfr_log1p(r22159, r22158, MPFR_RNDN); mpfr_div(r22160, r22144, r22159, MPFR_RNDN); mpfr_add(r22161, r22144, r22160, MPFR_RNDN); mpfr_mul(r22162, r22146, r22161, MPFR_RNDN); mpfr_sqrt(r22163, r22162, MPFR_RNDN); return mpfr_get_d(r22163, MPFR_RNDN); }