#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r22089 = 1; float r22090 = 2; float r22091 = r22089 / r22090; float r22092 = l; float r22093 = r22090 * r22092; float r22094 = Om; float r22095 = r22093 / r22094; float r22096 = pow(r22095, r22090); float r22097 = kx; float r22098 = sin(r22097); float r22099 = pow(r22098, r22090); float r22100 = ky; float r22101 = sin(r22100); float r22102 = pow(r22101, r22090); float r22103 = r22099 + r22102; float r22104 = r22096 * r22103; float r22105 = r22089 + r22104; float r22106 = sqrt(r22105); float r22107 = r22089 / r22106; float r22108 = r22089 + r22107; float r22109 = r22091 * r22108; float r22110 = sqrt(r22109); return r22110; } double f_id(double l, double Om, double kx, double ky) { double r22111 = 1; double r22112 = 2; double r22113 = r22111 / r22112; double r22114 = l; double r22115 = r22112 * r22114; double r22116 = Om; double r22117 = r22115 / r22116; double r22118 = pow(r22117, r22112); double r22119 = kx; double r22120 = sin(r22119); double r22121 = pow(r22120, r22112); double r22122 = ky; double r22123 = sin(r22122); double r22124 = pow(r22123, r22112); double r22125 = r22121 + r22124; double r22126 = r22118 * r22125; double r22127 = r22111 + r22126; double r22128 = sqrt(r22127); double r22129 = r22111 / r22128; double r22130 = r22111 + r22129; double r22131 = r22113 * r22130; double r22132 = sqrt(r22131); return r22132; } double f_of(float l, float Om, float kx, float ky) { float r22133 = 1; float r22134 = 2; float r22135 = r22133 / r22134; float r22136 = l; float r22137 = Om; float r22138 = r22137 / r22134; float r22139 = r22136 / r22138; float r22140 = kx; float r22141 = sin(r22140); float r22142 = ky; float r22143 = sin(r22142); float r22144 = hypot(r22141, r22143); float r22145 = r22139 * r22144; float r22146 = hypot(r22133, r22145); float r22147 = expm1(r22146); float r22148 = log1p(r22147); float r22149 = r22133 / r22148; float r22150 = r22133 + r22149; float r22151 = r22135 * r22150; float r22152 = sqrt(r22151); return r22152; } double f_od(double l, double Om, double kx, double ky) { double r22153 = 1; double r22154 = 2; double r22155 = r22153 / r22154; double r22156 = l; double r22157 = Om; double r22158 = r22157 / r22154; double r22159 = r22156 / r22158; double r22160 = kx; double r22161 = sin(r22160); double r22162 = ky; double r22163 = sin(r22162); double r22164 = hypot(r22161, r22163); double r22165 = r22159 * r22164; double r22166 = hypot(r22153, r22165); double r22167 = expm1(r22166); double r22168 = log1p(r22167); double r22169 = r22153 / r22168; double r22170 = r22153 + r22169; double r22171 = r22155 * r22170; double r22172 = sqrt(r22171); return r22172; } 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 r22173, r22174, r22175, r22176, r22177, r22178, r22179, r22180, r22181, r22182, r22183, r22184, r22185, r22186, r22187, r22188, r22189, r22190, r22191, r22192, r22193, r22194; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22173, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22174, "2", 10, MPFR_RNDN); 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); mpfr_init(r22186); mpfr_init(r22187); mpfr_init(r22188); mpfr_init(r22189); mpfr_init(r22190); mpfr_init(r22191); mpfr_init(r22192); mpfr_init(r22193); mpfr_init(r22194); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22175, r22173, r22174, MPFR_RNDN); mpfr_set_d(r22176, l, MPFR_RNDN); mpfr_mul(r22177, r22174, r22176, MPFR_RNDN); mpfr_set_d(r22178, Om, MPFR_RNDN); mpfr_div(r22179, r22177, r22178, MPFR_RNDN); mpfr_pow(r22180, r22179, r22174, MPFR_RNDN); mpfr_set_d(r22181, kx, MPFR_RNDN); mpfr_sin(r22182, r22181, MPFR_RNDN); mpfr_pow(r22183, r22182, r22174, MPFR_RNDN); mpfr_set_d(r22184, ky, MPFR_RNDN); mpfr_sin(r22185, r22184, MPFR_RNDN); mpfr_pow(r22186, r22185, r22174, MPFR_RNDN); mpfr_add(r22187, r22183, r22186, MPFR_RNDN); mpfr_mul(r22188, r22180, r22187, MPFR_RNDN); mpfr_add(r22189, r22173, r22188, MPFR_RNDN); mpfr_sqrt(r22190, r22189, MPFR_RNDN); mpfr_div(r22191, r22173, r22190, MPFR_RNDN); mpfr_add(r22192, r22173, r22191, MPFR_RNDN); mpfr_mul(r22193, r22175, r22192, MPFR_RNDN); mpfr_sqrt(r22194, r22193, MPFR_RNDN); return mpfr_get_d(r22194, MPFR_RNDN); } static mpfr_t r22195, r22196, r22197, r22198, r22199, r22200, r22201, r22202, r22203, r22204, r22205, r22206, r22207, r22208, r22209, r22210, r22211, r22212, r22213, r22214; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22195, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22196, "2", 10, MPFR_RNDN); 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); mpfr_init(r22206); mpfr_init(r22207); mpfr_init(r22208); mpfr_init(r22209); mpfr_init(r22210); mpfr_init(r22211); mpfr_init(r22212); mpfr_init(r22213); mpfr_init(r22214); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22197, r22195, r22196, MPFR_RNDN); mpfr_set_d(r22198, l, MPFR_RNDN); mpfr_set_d(r22199, Om, MPFR_RNDN); mpfr_div(r22200, r22199, r22196, MPFR_RNDN); mpfr_div(r22201, r22198, r22200, MPFR_RNDN); mpfr_set_d(r22202, kx, MPFR_RNDN); mpfr_sin(r22203, r22202, MPFR_RNDN); mpfr_set_d(r22204, ky, MPFR_RNDN); mpfr_sin(r22205, r22204, MPFR_RNDN); mpfr_hypot(r22206, r22203, r22205, MPFR_RNDN); mpfr_mul(r22207, r22201, r22206, MPFR_RNDN); mpfr_hypot(r22208, r22195, r22207, MPFR_RNDN); mpfr_expm1(r22209, r22208, MPFR_RNDN); mpfr_log1p(r22210, r22209, MPFR_RNDN); mpfr_div(r22211, r22195, r22210, MPFR_RNDN); mpfr_add(r22212, r22195, r22211, MPFR_RNDN); mpfr_mul(r22213, r22197, r22212, MPFR_RNDN); mpfr_sqrt(r22214, r22213, MPFR_RNDN); return mpfr_get_d(r22214, MPFR_RNDN); } static mpfr_t r22215, r22216, r22217, r22218, r22219, r22220, r22221, r22222, r22223, r22224, r22225, r22226, r22227, r22228, r22229, r22230, r22231, r22232, r22233, r22234; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22215, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22216, "2", 10, MPFR_RNDN); mpfr_init(r22217); mpfr_init(r22218); mpfr_init(r22219); mpfr_init(r22220); mpfr_init(r22221); mpfr_init(r22222); mpfr_init(r22223); mpfr_init(r22224); mpfr_init(r22225); mpfr_init(r22226); mpfr_init(r22227); mpfr_init(r22228); mpfr_init(r22229); mpfr_init(r22230); mpfr_init(r22231); mpfr_init(r22232); mpfr_init(r22233); mpfr_init(r22234); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22217, r22215, r22216, MPFR_RNDN); mpfr_set_d(r22218, l, MPFR_RNDN); mpfr_set_d(r22219, Om, MPFR_RNDN); mpfr_div(r22220, r22219, r22216, MPFR_RNDN); mpfr_div(r22221, r22218, r22220, MPFR_RNDN); mpfr_set_d(r22222, kx, MPFR_RNDN); mpfr_sin(r22223, r22222, MPFR_RNDN); mpfr_set_d(r22224, ky, MPFR_RNDN); mpfr_sin(r22225, r22224, MPFR_RNDN); mpfr_hypot(r22226, r22223, r22225, MPFR_RNDN); mpfr_mul(r22227, r22221, r22226, MPFR_RNDN); mpfr_hypot(r22228, r22215, r22227, MPFR_RNDN); mpfr_expm1(r22229, r22228, MPFR_RNDN); mpfr_log1p(r22230, r22229, MPFR_RNDN); mpfr_div(r22231, r22215, r22230, MPFR_RNDN); mpfr_add(r22232, r22215, r22231, MPFR_RNDN); mpfr_mul(r22233, r22217, r22232, MPFR_RNDN); mpfr_sqrt(r22234, r22233, MPFR_RNDN); return mpfr_get_d(r22234, MPFR_RNDN); }