#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r22060 = 1; float r22061 = 2; float r22062 = r22060 / r22061; float r22063 = l; float r22064 = r22061 * r22063; float r22065 = Om; float r22066 = r22064 / r22065; float r22067 = pow(r22066, r22061); float r22068 = kx; float r22069 = sin(r22068); float r22070 = pow(r22069, r22061); float r22071 = ky; float r22072 = sin(r22071); float r22073 = pow(r22072, r22061); float r22074 = r22070 + r22073; float r22075 = r22067 * r22074; float r22076 = r22060 + r22075; float r22077 = sqrt(r22076); float r22078 = r22060 / r22077; float r22079 = r22060 + r22078; float r22080 = r22062 * r22079; float r22081 = sqrt(r22080); return r22081; } double f_id(double l, double Om, double kx, double ky) { double r22082 = 1; double r22083 = 2; double r22084 = r22082 / r22083; double r22085 = l; double r22086 = r22083 * r22085; double r22087 = Om; double r22088 = r22086 / r22087; double r22089 = pow(r22088, r22083); double r22090 = kx; double r22091 = sin(r22090); double r22092 = pow(r22091, r22083); double r22093 = ky; double r22094 = sin(r22093); double r22095 = pow(r22094, r22083); double r22096 = r22092 + r22095; double r22097 = r22089 * r22096; double r22098 = r22082 + r22097; double r22099 = sqrt(r22098); double r22100 = r22082 / r22099; double r22101 = r22082 + r22100; double r22102 = r22084 * r22101; double r22103 = sqrt(r22102); return r22103; } double f_of(float l, float Om, float kx, float ky) { float r22104 = 1; float r22105 = 2; float r22106 = r22104 / r22105; float r22107 = l; float r22108 = Om; float r22109 = r22108 / r22105; float r22110 = r22107 / r22109; float r22111 = kx; float r22112 = sin(r22111); float r22113 = ky; float r22114 = sin(r22113); float r22115 = hypot(r22112, r22114); float r22116 = r22110 * r22115; float r22117 = hypot(r22104, r22116); float r22118 = expm1(r22117); float r22119 = log1p(r22118); float r22120 = r22104 / r22119; float r22121 = r22104 + r22120; float r22122 = r22106 * r22121; float r22123 = sqrt(r22122); return r22123; } double f_od(double l, double Om, double kx, double ky) { double r22124 = 1; double r22125 = 2; double r22126 = r22124 / r22125; double r22127 = l; double r22128 = Om; double r22129 = r22128 / r22125; double r22130 = r22127 / r22129; double r22131 = kx; double r22132 = sin(r22131); double r22133 = ky; double r22134 = sin(r22133); double r22135 = hypot(r22132, r22134); double r22136 = r22130 * r22135; double r22137 = hypot(r22124, r22136); double r22138 = expm1(r22137); double r22139 = log1p(r22138); double r22140 = r22124 / r22139; double r22141 = r22124 + r22140; double r22142 = r22126 * r22141; double r22143 = sqrt(r22142); return r22143; } 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 r22144, r22145, r22146, r22147, r22148, r22149, r22150, r22151, r22152, r22153, r22154, r22155, r22156, r22157, r22158, r22159, r22160, r22161, r22162, r22163, r22164, r22165; void setup_mpfr_f_im() { mpfr_set_default_prec(400); 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); mpfr_init(r22164); mpfr_init(r22165); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22146, r22144, r22145, MPFR_RNDN); mpfr_set_d(r22147, l, MPFR_RNDN); mpfr_mul(r22148, r22145, r22147, MPFR_RNDN); mpfr_set_d(r22149, Om, MPFR_RNDN); mpfr_div(r22150, r22148, r22149, MPFR_RNDN); mpfr_pow(r22151, r22150, r22145, MPFR_RNDN); mpfr_set_d(r22152, kx, MPFR_RNDN); mpfr_sin(r22153, r22152, MPFR_RNDN); mpfr_pow(r22154, r22153, r22145, MPFR_RNDN); mpfr_set_d(r22155, ky, MPFR_RNDN); mpfr_sin(r22156, r22155, MPFR_RNDN); mpfr_pow(r22157, r22156, r22145, MPFR_RNDN); mpfr_add(r22158, r22154, r22157, MPFR_RNDN); mpfr_mul(r22159, r22151, r22158, MPFR_RNDN); mpfr_add(r22160, r22144, r22159, MPFR_RNDN); mpfr_sqrt(r22161, r22160, MPFR_RNDN); mpfr_div(r22162, r22144, r22161, MPFR_RNDN); mpfr_add(r22163, r22144, r22162, MPFR_RNDN); mpfr_mul(r22164, r22146, r22163, MPFR_RNDN); mpfr_sqrt(r22165, r22164, MPFR_RNDN); return mpfr_get_d(r22165, MPFR_RNDN); } static mpfr_t r22166, r22167, r22168, r22169, r22170, r22171, r22172, r22173, r22174, r22175, r22176, r22177, r22178, r22179, r22180, r22181, r22182, r22183, r22184, r22185; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22166, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22167, "2", 10, MPFR_RNDN); mpfr_init(r22168); mpfr_init(r22169); mpfr_init(r22170); mpfr_init(r22171); mpfr_init(r22172); mpfr_init(r22173); mpfr_init(r22174); mpfr_init(r22175); mpfr_init(r22176); mpfr_init(r22177); mpfr_init(r22178); mpfr_init(r22179); mpfr_init(r22180); mpfr_init(r22181); mpfr_init(r22182); mpfr_init(r22183); mpfr_init(r22184); mpfr_init(r22185); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22168, r22166, r22167, MPFR_RNDN); mpfr_set_d(r22169, l, MPFR_RNDN); mpfr_set_d(r22170, Om, MPFR_RNDN); mpfr_div(r22171, r22170, r22167, MPFR_RNDN); mpfr_div(r22172, r22169, r22171, MPFR_RNDN); mpfr_set_d(r22173, kx, MPFR_RNDN); mpfr_sin(r22174, r22173, MPFR_RNDN); mpfr_set_d(r22175, ky, MPFR_RNDN); mpfr_sin(r22176, r22175, MPFR_RNDN); mpfr_hypot(r22177, r22174, r22176, MPFR_RNDN); mpfr_mul(r22178, r22172, r22177, MPFR_RNDN); mpfr_hypot(r22179, r22166, r22178, MPFR_RNDN); mpfr_expm1(r22180, r22179, MPFR_RNDN); mpfr_log1p(r22181, r22180, MPFR_RNDN); mpfr_div(r22182, r22166, r22181, MPFR_RNDN); mpfr_add(r22183, r22166, r22182, MPFR_RNDN); mpfr_mul(r22184, r22168, r22183, MPFR_RNDN); mpfr_sqrt(r22185, r22184, MPFR_RNDN); return mpfr_get_d(r22185, MPFR_RNDN); } static mpfr_t r22186, r22187, r22188, r22189, r22190, r22191, r22192, r22193, r22194, r22195, r22196, r22197, r22198, r22199, r22200, r22201, r22202, r22203, r22204, r22205; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22186, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22187, "2", 10, MPFR_RNDN); mpfr_init(r22188); mpfr_init(r22189); mpfr_init(r22190); mpfr_init(r22191); mpfr_init(r22192); mpfr_init(r22193); mpfr_init(r22194); mpfr_init(r22195); mpfr_init(r22196); mpfr_init(r22197); mpfr_init(r22198); mpfr_init(r22199); mpfr_init(r22200); mpfr_init(r22201); mpfr_init(r22202); mpfr_init(r22203); mpfr_init(r22204); mpfr_init(r22205); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22188, r22186, r22187, MPFR_RNDN); mpfr_set_d(r22189, l, MPFR_RNDN); mpfr_set_d(r22190, Om, MPFR_RNDN); mpfr_div(r22191, r22190, r22187, MPFR_RNDN); mpfr_div(r22192, r22189, r22191, MPFR_RNDN); mpfr_set_d(r22193, kx, MPFR_RNDN); mpfr_sin(r22194, r22193, MPFR_RNDN); mpfr_set_d(r22195, ky, MPFR_RNDN); mpfr_sin(r22196, r22195, MPFR_RNDN); mpfr_hypot(r22197, r22194, r22196, MPFR_RNDN); mpfr_mul(r22198, r22192, r22197, MPFR_RNDN); mpfr_hypot(r22199, r22186, r22198, MPFR_RNDN); mpfr_expm1(r22200, r22199, MPFR_RNDN); mpfr_log1p(r22201, r22200, MPFR_RNDN); mpfr_div(r22202, r22186, r22201, MPFR_RNDN); mpfr_add(r22203, r22186, r22202, MPFR_RNDN); mpfr_mul(r22204, r22188, r22203, MPFR_RNDN); mpfr_sqrt(r22205, r22204, MPFR_RNDN); return mpfr_get_d(r22205, MPFR_RNDN); }