#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r23219 = J; float r23220 = l; float r23221 = exp(r23220); float r23222 = -r23220; float r23223 = exp(r23222); float r23224 = r23221 - r23223; float r23225 = r23219 * r23224; float r23226 = K; float r23227 = 2; float r23228 = r23226 / r23227; float r23229 = cos(r23228); float r23230 = r23225 * r23229; float r23231 = U; float r23232 = r23230 + r23231; return r23232; } double f_id(double J, double l, double K, double U) { double r23233 = J; double r23234 = l; double r23235 = exp(r23234); double r23236 = -r23234; double r23237 = exp(r23236); double r23238 = r23235 - r23237; double r23239 = r23233 * r23238; double r23240 = K; double r23241 = 2; double r23242 = r23240 / r23241; double r23243 = cos(r23242); double r23244 = r23239 * r23243; double r23245 = U; double r23246 = r23244 + r23245; return r23246; } double f_of(float J, float l, float K, float U) { float r23247 = K; float r23248 = 2; float r23249 = r23247 / r23248; float r23250 = cos(r23249); float r23251 = J; float r23252 = r23250 * r23251; float r23253 = l; float r23254 = 1/3; float r23255 = r23253 * r23254; float r23256 = fma(r23255, r23253, r23248); float r23257 = 5; float r23258 = pow(r23253, r23257); float r23259 = 1/60; float r23260 = r23258 * r23259; float r23261 = fma(r23253, r23256, r23260); float r23262 = U; float r23263 = fma(r23252, r23261, r23262); return r23263; } double f_od(double J, double l, double K, double U) { double r23264 = K; double r23265 = 2; double r23266 = r23264 / r23265; double r23267 = cos(r23266); double r23268 = J; double r23269 = r23267 * r23268; double r23270 = l; double r23271 = 1/3; double r23272 = r23270 * r23271; double r23273 = fma(r23272, r23270, r23265); double r23274 = 5; double r23275 = pow(r23270, r23274); double r23276 = 1/60; double r23277 = r23275 * r23276; double r23278 = fma(r23270, r23273, r23277); double r23279 = U; double r23280 = fma(r23269, r23278, r23279); return r23280; } 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 r23281, r23282, r23283, r23284, r23285, r23286, r23287, r23288, r23289, r23290, r23291, r23292, r23293, r23294; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r23281); mpfr_init(r23282); mpfr_init(r23283); mpfr_init(r23284); mpfr_init(r23285); mpfr_init(r23286); mpfr_init(r23287); mpfr_init(r23288); mpfr_init_set_str(r23289, "2", 10, MPFR_RNDN); mpfr_init(r23290); mpfr_init(r23291); mpfr_init(r23292); mpfr_init(r23293); mpfr_init(r23294); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r23281, J, MPFR_RNDN); mpfr_set_d(r23282, l, MPFR_RNDN); mpfr_exp(r23283, r23282, MPFR_RNDN); mpfr_neg(r23284, r23282, MPFR_RNDN); mpfr_exp(r23285, r23284, MPFR_RNDN); mpfr_sub(r23286, r23283, r23285, MPFR_RNDN); mpfr_mul(r23287, r23281, r23286, MPFR_RNDN); mpfr_set_d(r23288, K, MPFR_RNDN); ; mpfr_div(r23290, r23288, r23289, MPFR_RNDN); mpfr_cos(r23291, r23290, MPFR_RNDN); mpfr_mul(r23292, r23287, r23291, MPFR_RNDN); mpfr_set_d(r23293, U, MPFR_RNDN); mpfr_add(r23294, r23292, r23293, MPFR_RNDN); return mpfr_get_d(r23294, MPFR_RNDN); } static mpfr_t r23295, r23296, r23297, r23298, r23299, r23300, r23301, r23302, r23303, r23304, r23305, r23306, r23307, r23308, r23309, r23310, r23311; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r23295); mpfr_init_set_str(r23296, "2", 10, MPFR_RNDN); mpfr_init(r23297); mpfr_init(r23298); mpfr_init(r23299); mpfr_init(r23300); mpfr_init(r23301); mpfr_init_set_str(r23302, "1/3", 10, MPFR_RNDN); mpfr_init(r23303); mpfr_init(r23304); mpfr_init_set_str(r23305, "5", 10, MPFR_RNDN); mpfr_init(r23306); mpfr_init_set_str(r23307, "1/60", 10, MPFR_RNDN); mpfr_init(r23308); mpfr_init(r23309); mpfr_init(r23310); mpfr_init(r23311); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r23295, K, MPFR_RNDN); ; mpfr_div(r23297, r23295, r23296, MPFR_RNDN); mpfr_cos(r23298, r23297, MPFR_RNDN); mpfr_set_d(r23299, J, MPFR_RNDN); mpfr_mul(r23300, r23298, r23299, MPFR_RNDN); mpfr_set_d(r23301, l, MPFR_RNDN); ; mpfr_mul(r23303, r23301, r23302, MPFR_RNDN); mpfr_fma(r23304, r23303, r23301, r23296, MPFR_RNDN); ; mpfr_pow(r23306, r23301, r23305, MPFR_RNDN); ; mpfr_mul(r23308, r23306, r23307, MPFR_RNDN); mpfr_fma(r23309, r23301, r23304, r23308, MPFR_RNDN); mpfr_set_d(r23310, U, MPFR_RNDN); mpfr_fma(r23311, r23300, r23309, r23310, MPFR_RNDN); return mpfr_get_d(r23311, MPFR_RNDN); } static mpfr_t r23312, r23313, r23314, r23315, r23316, r23317, r23318, r23319, r23320, r23321, r23322, r23323, r23324, r23325, r23326, r23327, r23328; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r23312); mpfr_init_set_str(r23313, "2", 10, MPFR_RNDN); mpfr_init(r23314); mpfr_init(r23315); mpfr_init(r23316); mpfr_init(r23317); mpfr_init(r23318); mpfr_init_set_str(r23319, "1/3", 10, MPFR_RNDN); mpfr_init(r23320); mpfr_init(r23321); mpfr_init_set_str(r23322, "5", 10, MPFR_RNDN); mpfr_init(r23323); mpfr_init_set_str(r23324, "1/60", 10, MPFR_RNDN); mpfr_init(r23325); mpfr_init(r23326); mpfr_init(r23327); mpfr_init(r23328); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r23312, K, MPFR_RNDN); ; mpfr_div(r23314, r23312, r23313, MPFR_RNDN); mpfr_cos(r23315, r23314, MPFR_RNDN); mpfr_set_d(r23316, J, MPFR_RNDN); mpfr_mul(r23317, r23315, r23316, MPFR_RNDN); mpfr_set_d(r23318, l, MPFR_RNDN); ; mpfr_mul(r23320, r23318, r23319, MPFR_RNDN); mpfr_fma(r23321, r23320, r23318, r23313, MPFR_RNDN); ; mpfr_pow(r23323, r23318, r23322, MPFR_RNDN); ; mpfr_mul(r23325, r23323, r23324, MPFR_RNDN); mpfr_fma(r23326, r23318, r23321, r23325, MPFR_RNDN); mpfr_set_d(r23327, U, MPFR_RNDN); mpfr_fma(r23328, r23317, r23326, r23327, MPFR_RNDN); return mpfr_get_d(r23328, MPFR_RNDN); }