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