#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r23239 = J; float r23240 = l; float r23241 = exp(r23240); float r23242 = -r23240; float r23243 = exp(r23242); float r23244 = r23241 - r23243; float r23245 = r23239 * r23244; float r23246 = K; float r23247 = 2; float r23248 = r23246 / r23247; float r23249 = cos(r23248); float r23250 = r23245 * r23249; float r23251 = U; float r23252 = r23250 + r23251; return r23252; } double f_id(double J, double l, double K, double U) { double r23253 = J; double r23254 = l; double r23255 = exp(r23254); double r23256 = -r23254; double r23257 = exp(r23256); double r23258 = r23255 - r23257; double r23259 = r23253 * r23258; double r23260 = K; double r23261 = 2; double r23262 = r23260 / r23261; double r23263 = cos(r23262); double r23264 = r23259 * r23263; double r23265 = U; double r23266 = r23264 + r23265; return r23266; } double f_of(float J, float l, float K, float U) { float r23267 = K; float r23268 = 1/2; float r23269 = r23267 * r23268; float r23270 = cos(r23269); float r23271 = l; float r23272 = r23271 * r23271; float r23273 = 1/3; float r23274 = r23272 * r23273; float r23275 = J; float r23276 = r23275 * r23271; float r23277 = 2; float r23278 = r23277 * r23275; float r23279 = r23278 * r23271; float r23280 = fma(r23274, r23276, r23279); float r23281 = U; float r23282 = fma(r23270, r23280, r23281); return r23282; } double f_od(double J, double l, double K, double U) { double r23283 = K; double r23284 = 1/2; double r23285 = r23283 * r23284; double r23286 = cos(r23285); double r23287 = l; double r23288 = r23287 * r23287; double r23289 = 1/3; double r23290 = r23288 * r23289; double r23291 = J; double r23292 = r23291 * r23287; double r23293 = 2; double r23294 = r23293 * r23291; double r23295 = r23294 * r23287; double r23296 = fma(r23290, r23292, r23295); double r23297 = U; double r23298 = fma(r23286, r23296, r23297); return r23298; } 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 r23299, r23300, r23301, r23302, r23303, r23304, r23305, r23306, r23307, r23308, r23309, r23310, r23311, r23312; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r23299); mpfr_init(r23300); mpfr_init(r23301); mpfr_init(r23302); mpfr_init(r23303); mpfr_init(r23304); mpfr_init(r23305); mpfr_init(r23306); mpfr_init_set_str(r23307, "2", 10, MPFR_RNDN); mpfr_init(r23308); mpfr_init(r23309); mpfr_init(r23310); mpfr_init(r23311); mpfr_init(r23312); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r23299, J, MPFR_RNDN); mpfr_set_d(r23300, l, MPFR_RNDN); mpfr_exp(r23301, r23300, MPFR_RNDN); mpfr_neg(r23302, r23300, MPFR_RNDN); mpfr_exp(r23303, r23302, MPFR_RNDN); mpfr_sub(r23304, r23301, r23303, MPFR_RNDN); mpfr_mul(r23305, r23299, r23304, MPFR_RNDN); mpfr_set_d(r23306, K, MPFR_RNDN); ; mpfr_div(r23308, r23306, r23307, MPFR_RNDN); mpfr_cos(r23309, r23308, MPFR_RNDN); mpfr_mul(r23310, r23305, r23309, MPFR_RNDN); mpfr_set_d(r23311, U, MPFR_RNDN); mpfr_add(r23312, r23310, r23311, MPFR_RNDN); return mpfr_get_d(r23312, MPFR_RNDN); } static mpfr_t r23313, r23314, r23315, r23316, r23317, r23318, r23319, r23320, r23321, r23322, r23323, r23324, r23325, r23326, r23327, r23328; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r23313); mpfr_init_set_str(r23314, "1/2", 10, MPFR_RNDN); 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(r23322); mpfr_init_set_str(r23323, "2", 10, MPFR_RNDN); mpfr_init(r23324); mpfr_init(r23325); mpfr_init(r23326); mpfr_init(r23327); mpfr_init(r23328); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r23313, K, MPFR_RNDN); ; mpfr_mul(r23315, r23313, r23314, MPFR_RNDN); mpfr_cos(r23316, r23315, MPFR_RNDN); mpfr_set_d(r23317, l, MPFR_RNDN); mpfr_mul(r23318, r23317, r23317, MPFR_RNDN); ; mpfr_mul(r23320, r23318, r23319, MPFR_RNDN); mpfr_set_d(r23321, J, MPFR_RNDN); mpfr_mul(r23322, r23321, r23317, MPFR_RNDN); ; mpfr_mul(r23324, r23323, r23321, MPFR_RNDN); mpfr_mul(r23325, r23324, r23317, MPFR_RNDN); mpfr_fma(r23326, r23320, r23322, r23325, MPFR_RNDN); mpfr_set_d(r23327, U, MPFR_RNDN); mpfr_fma(r23328, r23316, r23326, r23327, MPFR_RNDN); return mpfr_get_d(r23328, MPFR_RNDN); } static mpfr_t r23329, r23330, r23331, r23332, r23333, r23334, r23335, r23336, r23337, r23338, r23339, r23340, r23341, r23342, r23343, r23344; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r23329); mpfr_init_set_str(r23330, "1/2", 10, MPFR_RNDN); mpfr_init(r23331); mpfr_init(r23332); mpfr_init(r23333); mpfr_init(r23334); mpfr_init_set_str(r23335, "1/3", 10, MPFR_RNDN); mpfr_init(r23336); mpfr_init(r23337); mpfr_init(r23338); mpfr_init_set_str(r23339, "2", 10, MPFR_RNDN); mpfr_init(r23340); mpfr_init(r23341); mpfr_init(r23342); mpfr_init(r23343); mpfr_init(r23344); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r23329, K, MPFR_RNDN); ; mpfr_mul(r23331, r23329, r23330, MPFR_RNDN); mpfr_cos(r23332, r23331, MPFR_RNDN); mpfr_set_d(r23333, l, MPFR_RNDN); mpfr_mul(r23334, r23333, r23333, MPFR_RNDN); ; mpfr_mul(r23336, r23334, r23335, MPFR_RNDN); mpfr_set_d(r23337, J, MPFR_RNDN); mpfr_mul(r23338, r23337, r23333, MPFR_RNDN); ; mpfr_mul(r23340, r23339, r23337, MPFR_RNDN); mpfr_mul(r23341, r23340, r23333, MPFR_RNDN); mpfr_fma(r23342, r23336, r23338, r23341, MPFR_RNDN); mpfr_set_d(r23343, U, MPFR_RNDN); mpfr_fma(r23344, r23332, r23342, r23343, MPFR_RNDN); return mpfr_get_d(r23344, MPFR_RNDN); }