#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r21206 = 1; float r21207 = 2; float r21208 = r21206 / r21207; float r21209 = l; float r21210 = r21207 * r21209; float r21211 = Om; float r21212 = r21210 / r21211; float r21213 = pow(r21212, r21207); float r21214 = kx; float r21215 = sin(r21214); float r21216 = pow(r21215, r21207); float r21217 = ky; float r21218 = sin(r21217); float r21219 = pow(r21218, r21207); float r21220 = r21216 + r21219; float r21221 = r21213 * r21220; float r21222 = r21206 + r21221; float r21223 = sqrt(r21222); float r21224 = r21206 / r21223; float r21225 = r21206 + r21224; float r21226 = r21208 * r21225; float r21227 = sqrt(r21226); return r21227; } double f_id(double l, double Om, double kx, double ky) { double r21228 = 1; double r21229 = 2; double r21230 = r21228 / r21229; double r21231 = l; double r21232 = r21229 * r21231; double r21233 = Om; double r21234 = r21232 / r21233; double r21235 = pow(r21234, r21229); double r21236 = kx; double r21237 = sin(r21236); double r21238 = pow(r21237, r21229); double r21239 = ky; double r21240 = sin(r21239); double r21241 = pow(r21240, r21229); double r21242 = r21238 + r21241; double r21243 = r21235 * r21242; double r21244 = r21228 + r21243; double r21245 = sqrt(r21244); double r21246 = r21228 / r21245; double r21247 = r21228 + r21246; double r21248 = r21230 * r21247; double r21249 = sqrt(r21248); return r21249; } double f_of(float l, float Om, float kx, float ky) { float r21250 = 1; float r21251 = 2; float r21252 = r21250 / r21251; float r21253 = l; float r21254 = r21251 * r21253; float r21255 = Om; float r21256 = r21254 / r21255; float r21257 = kx; float r21258 = sin(r21257); float r21259 = pow(r21258, r21251); float r21260 = ky; float r21261 = sin(r21260); float r21262 = pow(r21261, r21251); float r21263 = r21259 + r21262; float r21264 = r21256 * r21263; float r21265 = r21256 * r21264; float r21266 = r21250 + r21265; float r21267 = sqrt(r21266); float r21268 = r21250 / r21267; float r21269 = r21250 + r21268; float r21270 = r21252 * r21269; float r21271 = sqrt(r21270); return r21271; } double f_od(double l, double Om, double kx, double ky) { double r21272 = 1; double r21273 = 2; double r21274 = r21272 / r21273; double r21275 = l; double r21276 = r21273 * r21275; double r21277 = Om; double r21278 = r21276 / r21277; double r21279 = kx; double r21280 = sin(r21279); double r21281 = pow(r21280, r21273); double r21282 = ky; double r21283 = sin(r21282); double r21284 = pow(r21283, r21273); double r21285 = r21281 + r21284; double r21286 = r21278 * r21285; double r21287 = r21278 * r21286; double r21288 = r21272 + r21287; double r21289 = sqrt(r21288); double r21290 = r21272 / r21289; double r21291 = r21272 + r21290; double r21292 = r21274 * r21291; double r21293 = sqrt(r21292); return r21293; } 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 r21294, r21295, r21296, r21297, r21298, r21299, r21300, r21301, r21302, r21303, r21304, r21305, r21306, r21307, r21308, r21309, r21310, r21311, r21312, r21313, r21314, r21315; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21294, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21295, "2", 10, MPFR_RNDN); mpfr_init(r21296); mpfr_init(r21297); mpfr_init(r21298); mpfr_init(r21299); mpfr_init(r21300); mpfr_init(r21301); mpfr_init(r21302); mpfr_init(r21303); mpfr_init(r21304); mpfr_init(r21305); mpfr_init(r21306); mpfr_init(r21307); mpfr_init(r21308); mpfr_init(r21309); mpfr_init(r21310); mpfr_init(r21311); mpfr_init(r21312); mpfr_init(r21313); mpfr_init(r21314); mpfr_init(r21315); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21296, r21294, r21295, MPFR_RNDN); mpfr_set_d(r21297, l, MPFR_RNDN); mpfr_mul(r21298, r21295, r21297, MPFR_RNDN); mpfr_set_d(r21299, Om, MPFR_RNDN); mpfr_div(r21300, r21298, r21299, MPFR_RNDN); mpfr_pow(r21301, r21300, r21295, MPFR_RNDN); mpfr_set_d(r21302, kx, MPFR_RNDN); mpfr_sin(r21303, r21302, MPFR_RNDN); mpfr_pow(r21304, r21303, r21295, MPFR_RNDN); mpfr_set_d(r21305, ky, MPFR_RNDN); mpfr_sin(r21306, r21305, MPFR_RNDN); mpfr_pow(r21307, r21306, r21295, MPFR_RNDN); mpfr_add(r21308, r21304, r21307, MPFR_RNDN); mpfr_mul(r21309, r21301, r21308, MPFR_RNDN); mpfr_add(r21310, r21294, r21309, MPFR_RNDN); mpfr_sqrt(r21311, r21310, MPFR_RNDN); mpfr_div(r21312, r21294, r21311, MPFR_RNDN); mpfr_add(r21313, r21294, r21312, MPFR_RNDN); mpfr_mul(r21314, r21296, r21313, MPFR_RNDN); mpfr_sqrt(r21315, r21314, MPFR_RNDN); return mpfr_get_d(r21315, MPFR_RNDN); } static mpfr_t r21316, r21317, r21318, r21319, r21320, r21321, r21322, r21323, r21324, r21325, r21326, r21327, r21328, r21329, r21330, r21331, r21332, r21333, r21334, r21335, r21336, r21337; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21316, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21317, "2", 10, MPFR_RNDN); mpfr_init(r21318); mpfr_init(r21319); mpfr_init(r21320); mpfr_init(r21321); mpfr_init(r21322); mpfr_init(r21323); mpfr_init(r21324); mpfr_init(r21325); mpfr_init(r21326); mpfr_init(r21327); mpfr_init(r21328); mpfr_init(r21329); mpfr_init(r21330); mpfr_init(r21331); mpfr_init(r21332); mpfr_init(r21333); mpfr_init(r21334); mpfr_init(r21335); mpfr_init(r21336); mpfr_init(r21337); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21318, r21316, r21317, MPFR_RNDN); mpfr_set_d(r21319, l, MPFR_RNDN); mpfr_mul(r21320, r21317, r21319, MPFR_RNDN); mpfr_set_d(r21321, Om, MPFR_RNDN); mpfr_div(r21322, r21320, r21321, MPFR_RNDN); mpfr_set_d(r21323, kx, MPFR_RNDN); mpfr_sin(r21324, r21323, MPFR_RNDN); mpfr_pow(r21325, r21324, r21317, MPFR_RNDN); mpfr_set_d(r21326, ky, MPFR_RNDN); mpfr_sin(r21327, r21326, MPFR_RNDN); mpfr_pow(r21328, r21327, r21317, MPFR_RNDN); mpfr_add(r21329, r21325, r21328, MPFR_RNDN); mpfr_mul(r21330, r21322, r21329, MPFR_RNDN); mpfr_mul(r21331, r21322, r21330, MPFR_RNDN); mpfr_add(r21332, r21316, r21331, MPFR_RNDN); mpfr_sqrt(r21333, r21332, MPFR_RNDN); mpfr_div(r21334, r21316, r21333, MPFR_RNDN); mpfr_add(r21335, r21316, r21334, MPFR_RNDN); mpfr_mul(r21336, r21318, r21335, MPFR_RNDN); mpfr_sqrt(r21337, r21336, MPFR_RNDN); return mpfr_get_d(r21337, MPFR_RNDN); } static mpfr_t r21338, r21339, r21340, r21341, r21342, r21343, r21344, r21345, r21346, r21347, r21348, r21349, r21350, r21351, r21352, r21353, r21354, r21355, r21356, r21357, r21358, r21359; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21338, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21339, "2", 10, MPFR_RNDN); mpfr_init(r21340); mpfr_init(r21341); mpfr_init(r21342); mpfr_init(r21343); mpfr_init(r21344); mpfr_init(r21345); mpfr_init(r21346); mpfr_init(r21347); mpfr_init(r21348); mpfr_init(r21349); mpfr_init(r21350); mpfr_init(r21351); mpfr_init(r21352); mpfr_init(r21353); mpfr_init(r21354); mpfr_init(r21355); mpfr_init(r21356); mpfr_init(r21357); mpfr_init(r21358); mpfr_init(r21359); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21340, r21338, r21339, MPFR_RNDN); mpfr_set_d(r21341, l, MPFR_RNDN); mpfr_mul(r21342, r21339, r21341, MPFR_RNDN); mpfr_set_d(r21343, Om, MPFR_RNDN); mpfr_div(r21344, r21342, r21343, MPFR_RNDN); mpfr_set_d(r21345, kx, MPFR_RNDN); mpfr_sin(r21346, r21345, MPFR_RNDN); mpfr_pow(r21347, r21346, r21339, MPFR_RNDN); mpfr_set_d(r21348, ky, MPFR_RNDN); mpfr_sin(r21349, r21348, MPFR_RNDN); mpfr_pow(r21350, r21349, r21339, MPFR_RNDN); mpfr_add(r21351, r21347, r21350, MPFR_RNDN); mpfr_mul(r21352, r21344, r21351, MPFR_RNDN); mpfr_mul(r21353, r21344, r21352, MPFR_RNDN); mpfr_add(r21354, r21338, r21353, MPFR_RNDN); mpfr_sqrt(r21355, r21354, MPFR_RNDN); mpfr_div(r21356, r21338, r21355, MPFR_RNDN); mpfr_add(r21357, r21338, r21356, MPFR_RNDN); mpfr_mul(r21358, r21340, r21357, MPFR_RNDN); mpfr_sqrt(r21359, r21358, MPFR_RNDN); return mpfr_get_d(r21359, MPFR_RNDN); }