#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r23051 = J; float r23052 = l; float r23053 = exp(r23052); float r23054 = -r23052; float r23055 = exp(r23054); float r23056 = r23053 - r23055; float r23057 = r23051 * r23056; float r23058 = K; float r23059 = 2; float r23060 = r23058 / r23059; float r23061 = cos(r23060); float r23062 = r23057 * r23061; float r23063 = U; float r23064 = r23062 + r23063; return r23064; } double f_id(double J, double l, double K, double U) { double r23065 = J; double r23066 = l; double r23067 = exp(r23066); double r23068 = -r23066; double r23069 = exp(r23068); double r23070 = r23067 - r23069; double r23071 = r23065 * r23070; double r23072 = K; double r23073 = 2; double r23074 = r23072 / r23073; double r23075 = cos(r23074); double r23076 = r23071 * r23075; double r23077 = U; double r23078 = r23076 + r23077; return r23078; } double f_of(float J, float l, float K, float U) { float r23079 = K; float r23080 = 2; float r23081 = r23079 / r23080; float r23082 = cos(r23081); float r23083 = J; float r23084 = r23082 * r23083; float r23085 = l; float r23086 = 1/3; float r23087 = r23085 * r23086; float r23088 = fma(r23087, r23085, r23080); float r23089 = 5; float r23090 = pow(r23085, r23089); float r23091 = 1/60; float r23092 = r23090 * r23091; float r23093 = fma(r23085, r23088, r23092); float r23094 = U; float r23095 = fma(r23084, r23093, r23094); return r23095; } double f_od(double J, double l, double K, double U) { double r23096 = K; double r23097 = 2; double r23098 = r23096 / r23097; double r23099 = cos(r23098); double r23100 = J; double r23101 = r23099 * r23100; double r23102 = l; double r23103 = 1/3; double r23104 = r23102 * r23103; double r23105 = fma(r23104, r23102, r23097); double r23106 = 5; double r23107 = pow(r23102, r23106); double r23108 = 1/60; double r23109 = r23107 * r23108; double r23110 = fma(r23102, r23105, r23109); double r23111 = U; double r23112 = fma(r23101, r23110, r23111); return r23112; } 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 r23113, r23114, r23115, r23116, r23117, r23118, r23119, r23120, r23121, r23122, r23123, r23124, r23125, r23126; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r23113); mpfr_init(r23114); mpfr_init(r23115); mpfr_init(r23116); mpfr_init(r23117); mpfr_init(r23118); mpfr_init(r23119); mpfr_init(r23120); mpfr_init_set_str(r23121, "2", 10, MPFR_RNDN); mpfr_init(r23122); mpfr_init(r23123); mpfr_init(r23124); mpfr_init(r23125); mpfr_init(r23126); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r23113, J, MPFR_RNDN); mpfr_set_d(r23114, l, MPFR_RNDN); mpfr_exp(r23115, r23114, MPFR_RNDN); mpfr_neg(r23116, r23114, MPFR_RNDN); mpfr_exp(r23117, r23116, MPFR_RNDN); mpfr_sub(r23118, r23115, r23117, MPFR_RNDN); mpfr_mul(r23119, r23113, r23118, MPFR_RNDN); mpfr_set_d(r23120, K, MPFR_RNDN); ; mpfr_div(r23122, r23120, r23121, MPFR_RNDN); mpfr_cos(r23123, r23122, MPFR_RNDN); mpfr_mul(r23124, r23119, r23123, MPFR_RNDN); mpfr_set_d(r23125, U, MPFR_RNDN); mpfr_add(r23126, r23124, r23125, MPFR_RNDN); return mpfr_get_d(r23126, MPFR_RNDN); } static mpfr_t r23127, r23128, r23129, r23130, r23131, r23132, r23133, r23134, r23135, r23136, r23137, r23138, r23139, r23140, r23141, r23142, r23143; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r23127); mpfr_init_set_str(r23128, "2", 10, MPFR_RNDN); mpfr_init(r23129); mpfr_init(r23130); mpfr_init(r23131); mpfr_init(r23132); mpfr_init(r23133); mpfr_init_set_str(r23134, "1/3", 10, MPFR_RNDN); mpfr_init(r23135); mpfr_init(r23136); mpfr_init_set_str(r23137, "5", 10, MPFR_RNDN); mpfr_init(r23138); mpfr_init_set_str(r23139, "1/60", 10, MPFR_RNDN); mpfr_init(r23140); mpfr_init(r23141); mpfr_init(r23142); mpfr_init(r23143); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r23127, K, MPFR_RNDN); ; mpfr_div(r23129, r23127, r23128, MPFR_RNDN); mpfr_cos(r23130, r23129, MPFR_RNDN); mpfr_set_d(r23131, J, MPFR_RNDN); mpfr_mul(r23132, r23130, r23131, MPFR_RNDN); mpfr_set_d(r23133, l, MPFR_RNDN); ; mpfr_mul(r23135, r23133, r23134, MPFR_RNDN); mpfr_fma(r23136, r23135, r23133, r23128, MPFR_RNDN); ; mpfr_pow(r23138, r23133, r23137, MPFR_RNDN); ; mpfr_mul(r23140, r23138, r23139, MPFR_RNDN); mpfr_fma(r23141, r23133, r23136, r23140, MPFR_RNDN); mpfr_set_d(r23142, U, MPFR_RNDN); mpfr_fma(r23143, r23132, r23141, r23142, MPFR_RNDN); return mpfr_get_d(r23143, MPFR_RNDN); } static mpfr_t r23144, r23145, r23146, r23147, r23148, r23149, r23150, r23151, r23152, r23153, r23154, r23155, r23156, r23157, r23158, r23159, r23160; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r23144); mpfr_init_set_str(r23145, "2", 10, MPFR_RNDN); mpfr_init(r23146); mpfr_init(r23147); mpfr_init(r23148); mpfr_init(r23149); mpfr_init(r23150); mpfr_init_set_str(r23151, "1/3", 10, MPFR_RNDN); mpfr_init(r23152); mpfr_init(r23153); mpfr_init_set_str(r23154, "5", 10, MPFR_RNDN); mpfr_init(r23155); mpfr_init_set_str(r23156, "1/60", 10, MPFR_RNDN); mpfr_init(r23157); mpfr_init(r23158); mpfr_init(r23159); mpfr_init(r23160); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r23144, K, MPFR_RNDN); ; mpfr_div(r23146, r23144, r23145, MPFR_RNDN); mpfr_cos(r23147, r23146, MPFR_RNDN); mpfr_set_d(r23148, J, MPFR_RNDN); mpfr_mul(r23149, r23147, r23148, MPFR_RNDN); mpfr_set_d(r23150, l, MPFR_RNDN); ; mpfr_mul(r23152, r23150, r23151, MPFR_RNDN); mpfr_fma(r23153, r23152, r23150, r23145, MPFR_RNDN); ; mpfr_pow(r23155, r23150, r23154, MPFR_RNDN); ; mpfr_mul(r23157, r23155, r23156, MPFR_RNDN); mpfr_fma(r23158, r23150, r23153, r23157, MPFR_RNDN); mpfr_set_d(r23159, U, MPFR_RNDN); mpfr_fma(r23160, r23149, r23158, r23159, MPFR_RNDN); return mpfr_get_d(r23160, MPFR_RNDN); }