#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r22172 = 1; float r22173 = 2; float r22174 = r22172 / r22173; float r22175 = l; float r22176 = r22173 * r22175; float r22177 = Om; float r22178 = r22176 / r22177; float r22179 = pow(r22178, r22173); float r22180 = kx; float r22181 = sin(r22180); float r22182 = pow(r22181, r22173); float r22183 = ky; float r22184 = sin(r22183); float r22185 = pow(r22184, r22173); float r22186 = r22182 + r22185; float r22187 = r22179 * r22186; float r22188 = r22172 + r22187; float r22189 = sqrt(r22188); float r22190 = r22172 / r22189; float r22191 = r22172 + r22190; float r22192 = r22174 * r22191; float r22193 = sqrt(r22192); return r22193; } double f_id(double l, double Om, double kx, double ky) { double r22194 = 1; double r22195 = 2; double r22196 = r22194 / r22195; double r22197 = l; double r22198 = r22195 * r22197; double r22199 = Om; double r22200 = r22198 / r22199; double r22201 = pow(r22200, r22195); double r22202 = kx; double r22203 = sin(r22202); double r22204 = pow(r22203, r22195); double r22205 = ky; double r22206 = sin(r22205); double r22207 = pow(r22206, r22195); double r22208 = r22204 + r22207; double r22209 = r22201 * r22208; double r22210 = r22194 + r22209; double r22211 = sqrt(r22210); double r22212 = r22194 / r22211; double r22213 = r22194 + r22212; double r22214 = r22196 * r22213; double r22215 = sqrt(r22214); return r22215; } double f_of(float l, float Om, float kx, float ky) { float r22216 = 1; float r22217 = 2; float r22218 = r22216 / r22217; float r22219 = l; float r22220 = r22217 * r22219; float r22221 = Om; float r22222 = r22220 / r22221; float r22223 = pow(r22222, r22217); float r22224 = kx; float r22225 = sin(r22224); float r22226 = pow(r22225, r22217); float r22227 = ky; float r22228 = sin(r22227); float r22229 = pow(r22228, r22217); float r22230 = r22226 + r22229; float r22231 = r22223 * r22230; float r22232 = r22216 + r22231; float r22233 = sqrt(r22232); float r22234 = r22216 / r22233; float r22235 = r22216 + r22234; float r22236 = r22218 * r22235; float r22237 = sqrt(r22236); return r22237; } double f_od(double l, double Om, double kx, double ky) { double r22238 = 1; double r22239 = 2; double r22240 = r22238 / r22239; double r22241 = l; double r22242 = r22239 * r22241; double r22243 = Om; double r22244 = r22242 / r22243; double r22245 = pow(r22244, r22239); double r22246 = kx; double r22247 = sin(r22246); double r22248 = pow(r22247, r22239); double r22249 = ky; double r22250 = sin(r22249); double r22251 = pow(r22250, r22239); double r22252 = r22248 + r22251; double r22253 = r22245 * r22252; double r22254 = r22238 + r22253; double r22255 = sqrt(r22254); double r22256 = r22238 / r22255; double r22257 = r22238 + r22256; double r22258 = r22240 * r22257; double r22259 = sqrt(r22258); return r22259; } 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 r22260, r22261, r22262, r22263, r22264, r22265, r22266, r22267, r22268, r22269, r22270, r22271, r22272, r22273, r22274, r22275, r22276, r22277, r22278, r22279, r22280, r22281; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r22260, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22261, "2", 10, MPFR_RNDN); mpfr_init(r22262); mpfr_init(r22263); mpfr_init(r22264); mpfr_init(r22265); mpfr_init(r22266); mpfr_init(r22267); mpfr_init(r22268); mpfr_init(r22269); mpfr_init(r22270); mpfr_init(r22271); mpfr_init(r22272); mpfr_init(r22273); mpfr_init(r22274); mpfr_init(r22275); mpfr_init(r22276); mpfr_init(r22277); mpfr_init(r22278); mpfr_init(r22279); mpfr_init(r22280); mpfr_init(r22281); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22262, r22260, r22261, MPFR_RNDN); mpfr_set_d(r22263, l, MPFR_RNDN); mpfr_mul(r22264, r22261, r22263, MPFR_RNDN); mpfr_set_d(r22265, Om, MPFR_RNDN); mpfr_div(r22266, r22264, r22265, MPFR_RNDN); mpfr_pow(r22267, r22266, r22261, MPFR_RNDN); mpfr_set_d(r22268, kx, MPFR_RNDN); mpfr_sin(r22269, r22268, MPFR_RNDN); mpfr_pow(r22270, r22269, r22261, MPFR_RNDN); mpfr_set_d(r22271, ky, MPFR_RNDN); mpfr_sin(r22272, r22271, MPFR_RNDN); mpfr_pow(r22273, r22272, r22261, MPFR_RNDN); mpfr_add(r22274, r22270, r22273, MPFR_RNDN); mpfr_mul(r22275, r22267, r22274, MPFR_RNDN); mpfr_add(r22276, r22260, r22275, MPFR_RNDN); mpfr_sqrt(r22277, r22276, MPFR_RNDN); mpfr_div(r22278, r22260, r22277, MPFR_RNDN); mpfr_add(r22279, r22260, r22278, MPFR_RNDN); mpfr_mul(r22280, r22262, r22279, MPFR_RNDN); mpfr_sqrt(r22281, r22280, MPFR_RNDN); return mpfr_get_d(r22281, MPFR_RNDN); } static mpfr_t r22282, r22283, r22284, r22285, r22286, r22287, r22288, r22289, r22290, r22291, r22292, r22293, r22294, r22295, r22296, r22297, r22298, r22299, r22300, r22301, r22302, r22303; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r22282, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22283, "2", 10, MPFR_RNDN); mpfr_init(r22284); mpfr_init(r22285); mpfr_init(r22286); mpfr_init(r22287); mpfr_init(r22288); mpfr_init(r22289); mpfr_init(r22290); mpfr_init(r22291); mpfr_init(r22292); mpfr_init(r22293); mpfr_init(r22294); mpfr_init(r22295); mpfr_init(r22296); mpfr_init(r22297); mpfr_init(r22298); mpfr_init(r22299); mpfr_init(r22300); mpfr_init(r22301); mpfr_init(r22302); mpfr_init(r22303); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22284, r22282, r22283, MPFR_RNDN); mpfr_set_d(r22285, l, MPFR_RNDN); mpfr_mul(r22286, r22283, r22285, MPFR_RNDN); mpfr_set_d(r22287, Om, MPFR_RNDN); mpfr_div(r22288, r22286, r22287, MPFR_RNDN); mpfr_pow(r22289, r22288, r22283, MPFR_RNDN); mpfr_set_d(r22290, kx, MPFR_RNDN); mpfr_sin(r22291, r22290, MPFR_RNDN); mpfr_pow(r22292, r22291, r22283, MPFR_RNDN); mpfr_set_d(r22293, ky, MPFR_RNDN); mpfr_sin(r22294, r22293, MPFR_RNDN); mpfr_pow(r22295, r22294, r22283, MPFR_RNDN); mpfr_add(r22296, r22292, r22295, MPFR_RNDN); mpfr_mul(r22297, r22289, r22296, MPFR_RNDN); mpfr_add(r22298, r22282, r22297, MPFR_RNDN); mpfr_sqrt(r22299, r22298, MPFR_RNDN); mpfr_div(r22300, r22282, r22299, MPFR_RNDN); mpfr_add(r22301, r22282, r22300, MPFR_RNDN); mpfr_mul(r22302, r22284, r22301, MPFR_RNDN); mpfr_sqrt(r22303, r22302, MPFR_RNDN); return mpfr_get_d(r22303, MPFR_RNDN); } static mpfr_t r22304, r22305, r22306, r22307, r22308, r22309, r22310, r22311, r22312, r22313, r22314, r22315, r22316, r22317, r22318, r22319, r22320, r22321, r22322, r22323, r22324, r22325; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r22304, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22305, "2", 10, MPFR_RNDN); mpfr_init(r22306); mpfr_init(r22307); mpfr_init(r22308); mpfr_init(r22309); mpfr_init(r22310); mpfr_init(r22311); mpfr_init(r22312); mpfr_init(r22313); mpfr_init(r22314); mpfr_init(r22315); mpfr_init(r22316); mpfr_init(r22317); mpfr_init(r22318); mpfr_init(r22319); mpfr_init(r22320); mpfr_init(r22321); mpfr_init(r22322); mpfr_init(r22323); mpfr_init(r22324); mpfr_init(r22325); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22306, r22304, r22305, MPFR_RNDN); mpfr_set_d(r22307, l, MPFR_RNDN); mpfr_mul(r22308, r22305, r22307, MPFR_RNDN); mpfr_set_d(r22309, Om, MPFR_RNDN); mpfr_div(r22310, r22308, r22309, MPFR_RNDN); mpfr_pow(r22311, r22310, r22305, MPFR_RNDN); mpfr_set_d(r22312, kx, MPFR_RNDN); mpfr_sin(r22313, r22312, MPFR_RNDN); mpfr_pow(r22314, r22313, r22305, MPFR_RNDN); mpfr_set_d(r22315, ky, MPFR_RNDN); mpfr_sin(r22316, r22315, MPFR_RNDN); mpfr_pow(r22317, r22316, r22305, MPFR_RNDN); mpfr_add(r22318, r22314, r22317, MPFR_RNDN); mpfr_mul(r22319, r22311, r22318, MPFR_RNDN); mpfr_add(r22320, r22304, r22319, MPFR_RNDN); mpfr_sqrt(r22321, r22320, MPFR_RNDN); mpfr_div(r22322, r22304, r22321, MPFR_RNDN); mpfr_add(r22323, r22304, r22322, MPFR_RNDN); mpfr_mul(r22324, r22306, r22323, MPFR_RNDN); mpfr_sqrt(r22325, r22324, MPFR_RNDN); return mpfr_get_d(r22325, MPFR_RNDN); }