#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r23167 = J; float r23168 = l; float r23169 = exp(r23168); float r23170 = -r23168; float r23171 = exp(r23170); float r23172 = r23169 - r23171; float r23173 = r23167 * r23172; float r23174 = K; float r23175 = 2; float r23176 = r23174 / r23175; float r23177 = cos(r23176); float r23178 = r23173 * r23177; float r23179 = U; float r23180 = r23178 + r23179; return r23180; } double f_id(double J, double l, double K, double U) { double r23181 = J; double r23182 = l; double r23183 = exp(r23182); double r23184 = -r23182; double r23185 = exp(r23184); double r23186 = r23183 - r23185; double r23187 = r23181 * r23186; double r23188 = K; double r23189 = 2; double r23190 = r23188 / r23189; double r23191 = cos(r23190); double r23192 = r23187 * r23191; double r23193 = U; double r23194 = r23192 + r23193; return r23194; } double f_of(float J, float l, float K, float U) { float r23195 = K; float r23196 = 2; float r23197 = r23195 / r23196; float r23198 = cos(r23197); float r23199 = J; float r23200 = r23198 * r23199; float r23201 = l; float r23202 = 1/3; float r23203 = r23201 * r23202; float r23204 = fma(r23203, r23201, r23196); float r23205 = 5; float r23206 = pow(r23201, r23205); float r23207 = 1/60; float r23208 = r23206 * r23207; float r23209 = fma(r23201, r23204, r23208); float r23210 = U; float r23211 = fma(r23200, r23209, r23210); return r23211; } double f_od(double J, double l, double K, double U) { double r23212 = K; double r23213 = 2; double r23214 = r23212 / r23213; double r23215 = cos(r23214); double r23216 = J; double r23217 = r23215 * r23216; double r23218 = l; double r23219 = 1/3; double r23220 = r23218 * r23219; double r23221 = fma(r23220, r23218, r23213); double r23222 = 5; double r23223 = pow(r23218, r23222); double r23224 = 1/60; double r23225 = r23223 * r23224; double r23226 = fma(r23218, r23221, r23225); double r23227 = U; double r23228 = fma(r23217, r23226, r23227); return r23228; } 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 r23229, r23230, r23231, r23232, r23233, r23234, r23235, r23236, r23237, r23238, r23239, r23240, r23241, r23242; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r23229); mpfr_init(r23230); mpfr_init(r23231); mpfr_init(r23232); mpfr_init(r23233); mpfr_init(r23234); mpfr_init(r23235); mpfr_init(r23236); mpfr_init_set_str(r23237, "2", 10, MPFR_RNDN); mpfr_init(r23238); mpfr_init(r23239); mpfr_init(r23240); mpfr_init(r23241); mpfr_init(r23242); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r23229, J, MPFR_RNDN); mpfr_set_d(r23230, l, MPFR_RNDN); mpfr_exp(r23231, r23230, MPFR_RNDN); mpfr_neg(r23232, r23230, MPFR_RNDN); mpfr_exp(r23233, r23232, MPFR_RNDN); mpfr_sub(r23234, r23231, r23233, MPFR_RNDN); mpfr_mul(r23235, r23229, r23234, MPFR_RNDN); mpfr_set_d(r23236, K, MPFR_RNDN); ; mpfr_div(r23238, r23236, r23237, MPFR_RNDN); mpfr_cos(r23239, r23238, MPFR_RNDN); mpfr_mul(r23240, r23235, r23239, MPFR_RNDN); mpfr_set_d(r23241, U, MPFR_RNDN); mpfr_add(r23242, r23240, r23241, MPFR_RNDN); return mpfr_get_d(r23242, MPFR_RNDN); } static mpfr_t r23243, r23244, r23245, r23246, r23247, r23248, r23249, r23250, r23251, r23252, r23253, r23254, r23255, r23256, r23257, r23258, r23259; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r23243); mpfr_init_set_str(r23244, "2", 10, MPFR_RNDN); mpfr_init(r23245); mpfr_init(r23246); mpfr_init(r23247); mpfr_init(r23248); mpfr_init(r23249); mpfr_init_set_str(r23250, "1/3", 10, MPFR_RNDN); mpfr_init(r23251); mpfr_init(r23252); mpfr_init_set_str(r23253, "5", 10, MPFR_RNDN); mpfr_init(r23254); mpfr_init_set_str(r23255, "1/60", 10, MPFR_RNDN); mpfr_init(r23256); mpfr_init(r23257); mpfr_init(r23258); mpfr_init(r23259); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r23243, K, MPFR_RNDN); ; mpfr_div(r23245, r23243, r23244, MPFR_RNDN); mpfr_cos(r23246, r23245, MPFR_RNDN); mpfr_set_d(r23247, J, MPFR_RNDN); mpfr_mul(r23248, r23246, r23247, MPFR_RNDN); mpfr_set_d(r23249, l, MPFR_RNDN); ; mpfr_mul(r23251, r23249, r23250, MPFR_RNDN); mpfr_fma(r23252, r23251, r23249, r23244, MPFR_RNDN); ; mpfr_pow(r23254, r23249, r23253, MPFR_RNDN); ; mpfr_mul(r23256, r23254, r23255, MPFR_RNDN); mpfr_fma(r23257, r23249, r23252, r23256, MPFR_RNDN); mpfr_set_d(r23258, U, MPFR_RNDN); mpfr_fma(r23259, r23248, r23257, r23258, MPFR_RNDN); return mpfr_get_d(r23259, MPFR_RNDN); } static mpfr_t r23260, r23261, r23262, r23263, r23264, r23265, r23266, r23267, r23268, r23269, r23270, r23271, r23272, r23273, r23274, r23275, r23276; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r23260); mpfr_init_set_str(r23261, "2", 10, MPFR_RNDN); mpfr_init(r23262); mpfr_init(r23263); mpfr_init(r23264); mpfr_init(r23265); mpfr_init(r23266); mpfr_init_set_str(r23267, "1/3", 10, MPFR_RNDN); mpfr_init(r23268); mpfr_init(r23269); mpfr_init_set_str(r23270, "5", 10, MPFR_RNDN); mpfr_init(r23271); mpfr_init_set_str(r23272, "1/60", 10, MPFR_RNDN); mpfr_init(r23273); mpfr_init(r23274); mpfr_init(r23275); mpfr_init(r23276); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r23260, K, MPFR_RNDN); ; mpfr_div(r23262, r23260, r23261, MPFR_RNDN); mpfr_cos(r23263, r23262, MPFR_RNDN); mpfr_set_d(r23264, J, MPFR_RNDN); mpfr_mul(r23265, r23263, r23264, MPFR_RNDN); mpfr_set_d(r23266, l, MPFR_RNDN); ; mpfr_mul(r23268, r23266, r23267, MPFR_RNDN); mpfr_fma(r23269, r23268, r23266, r23261, MPFR_RNDN); ; mpfr_pow(r23271, r23266, r23270, MPFR_RNDN); ; mpfr_mul(r23273, r23271, r23272, MPFR_RNDN); mpfr_fma(r23274, r23266, r23269, r23273, MPFR_RNDN); mpfr_set_d(r23275, U, MPFR_RNDN); mpfr_fma(r23276, r23265, r23274, r23275, MPFR_RNDN); return mpfr_get_d(r23276, MPFR_RNDN); }