#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r22995 = J; float r22996 = l; float r22997 = exp(r22996); float r22998 = -r22996; float r22999 = exp(r22998); float r23000 = r22997 - r22999; float r23001 = r22995 * r23000; float r23002 = K; float r23003 = 2; float r23004 = r23002 / r23003; float r23005 = cos(r23004); float r23006 = r23001 * r23005; float r23007 = U; float r23008 = r23006 + r23007; return r23008; } double f_id(double J, double l, double K, double U) { double r23009 = J; double r23010 = l; double r23011 = exp(r23010); double r23012 = -r23010; double r23013 = exp(r23012); double r23014 = r23011 - r23013; double r23015 = r23009 * r23014; double r23016 = K; double r23017 = 2; double r23018 = r23016 / r23017; double r23019 = cos(r23018); double r23020 = r23015 * r23019; double r23021 = U; double r23022 = r23020 + r23021; return r23022; } double f_of(float J, float l, float K, float U) { float r23023 = K; float r23024 = 2; float r23025 = r23023 / r23024; float r23026 = cos(r23025); float r23027 = J; float r23028 = r23026 * r23027; float r23029 = l; float r23030 = 1/3; float r23031 = r23029 * r23030; float r23032 = fma(r23031, r23029, r23024); float r23033 = 5; float r23034 = pow(r23029, r23033); float r23035 = 1/60; float r23036 = r23034 * r23035; float r23037 = fma(r23029, r23032, r23036); float r23038 = U; float r23039 = fma(r23028, r23037, r23038); return r23039; } double f_od(double J, double l, double K, double U) { double r23040 = K; double r23041 = 2; double r23042 = r23040 / r23041; double r23043 = cos(r23042); double r23044 = J; double r23045 = r23043 * r23044; double r23046 = l; double r23047 = 1/3; double r23048 = r23046 * r23047; double r23049 = fma(r23048, r23046, r23041); double r23050 = 5; double r23051 = pow(r23046, r23050); double r23052 = 1/60; double r23053 = r23051 * r23052; double r23054 = fma(r23046, r23049, r23053); double r23055 = U; double r23056 = fma(r23045, r23054, r23055); return r23056; } 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 r23057, r23058, r23059, r23060, r23061, r23062, r23063, r23064, r23065, r23066, r23067, r23068, r23069, r23070; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r23057); mpfr_init(r23058); mpfr_init(r23059); mpfr_init(r23060); mpfr_init(r23061); mpfr_init(r23062); mpfr_init(r23063); mpfr_init(r23064); mpfr_init_set_str(r23065, "2", 10, MPFR_RNDN); mpfr_init(r23066); mpfr_init(r23067); mpfr_init(r23068); mpfr_init(r23069); mpfr_init(r23070); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r23057, J, MPFR_RNDN); mpfr_set_d(r23058, l, MPFR_RNDN); mpfr_exp(r23059, r23058, MPFR_RNDN); mpfr_neg(r23060, r23058, MPFR_RNDN); mpfr_exp(r23061, r23060, MPFR_RNDN); mpfr_sub(r23062, r23059, r23061, MPFR_RNDN); mpfr_mul(r23063, r23057, r23062, MPFR_RNDN); mpfr_set_d(r23064, K, MPFR_RNDN); ; mpfr_div(r23066, r23064, r23065, MPFR_RNDN); mpfr_cos(r23067, r23066, MPFR_RNDN); mpfr_mul(r23068, r23063, r23067, MPFR_RNDN); mpfr_set_d(r23069, U, MPFR_RNDN); mpfr_add(r23070, r23068, r23069, MPFR_RNDN); return mpfr_get_d(r23070, MPFR_RNDN); } static mpfr_t r23071, r23072, r23073, r23074, r23075, r23076, r23077, r23078, r23079, r23080, r23081, r23082, r23083, r23084, r23085, r23086, r23087; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r23071); mpfr_init_set_str(r23072, "2", 10, MPFR_RNDN); mpfr_init(r23073); mpfr_init(r23074); mpfr_init(r23075); mpfr_init(r23076); mpfr_init(r23077); mpfr_init_set_str(r23078, "1/3", 10, MPFR_RNDN); mpfr_init(r23079); mpfr_init(r23080); mpfr_init_set_str(r23081, "5", 10, MPFR_RNDN); mpfr_init(r23082); mpfr_init_set_str(r23083, "1/60", 10, MPFR_RNDN); mpfr_init(r23084); mpfr_init(r23085); mpfr_init(r23086); mpfr_init(r23087); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r23071, K, MPFR_RNDN); ; mpfr_div(r23073, r23071, r23072, MPFR_RNDN); mpfr_cos(r23074, r23073, MPFR_RNDN); mpfr_set_d(r23075, J, MPFR_RNDN); mpfr_mul(r23076, r23074, r23075, MPFR_RNDN); mpfr_set_d(r23077, l, MPFR_RNDN); ; mpfr_mul(r23079, r23077, r23078, MPFR_RNDN); mpfr_fma(r23080, r23079, r23077, r23072, MPFR_RNDN); ; mpfr_pow(r23082, r23077, r23081, MPFR_RNDN); ; mpfr_mul(r23084, r23082, r23083, MPFR_RNDN); mpfr_fma(r23085, r23077, r23080, r23084, MPFR_RNDN); mpfr_set_d(r23086, U, MPFR_RNDN); mpfr_fma(r23087, r23076, r23085, r23086, MPFR_RNDN); return mpfr_get_d(r23087, MPFR_RNDN); } static mpfr_t r23088, r23089, r23090, r23091, r23092, r23093, r23094, r23095, r23096, r23097, r23098, r23099, r23100, r23101, r23102, r23103, r23104; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r23088); mpfr_init_set_str(r23089, "2", 10, MPFR_RNDN); mpfr_init(r23090); mpfr_init(r23091); mpfr_init(r23092); mpfr_init(r23093); mpfr_init(r23094); mpfr_init_set_str(r23095, "1/3", 10, MPFR_RNDN); mpfr_init(r23096); mpfr_init(r23097); mpfr_init_set_str(r23098, "5", 10, MPFR_RNDN); mpfr_init(r23099); mpfr_init_set_str(r23100, "1/60", 10, MPFR_RNDN); mpfr_init(r23101); mpfr_init(r23102); mpfr_init(r23103); mpfr_init(r23104); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r23088, K, MPFR_RNDN); ; mpfr_div(r23090, r23088, r23089, MPFR_RNDN); mpfr_cos(r23091, r23090, MPFR_RNDN); mpfr_set_d(r23092, J, MPFR_RNDN); mpfr_mul(r23093, r23091, r23092, MPFR_RNDN); mpfr_set_d(r23094, l, MPFR_RNDN); ; mpfr_mul(r23096, r23094, r23095, MPFR_RNDN); mpfr_fma(r23097, r23096, r23094, r23089, MPFR_RNDN); ; mpfr_pow(r23099, r23094, r23098, MPFR_RNDN); ; mpfr_mul(r23101, r23099, r23100, MPFR_RNDN); mpfr_fma(r23102, r23094, r23097, r23101, MPFR_RNDN); mpfr_set_d(r23103, U, MPFR_RNDN); mpfr_fma(r23104, r23093, r23102, r23103, MPFR_RNDN); return mpfr_get_d(r23104, MPFR_RNDN); }