#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (32)"; double f_if(float K, float m, float n, float M, float l) { float r23760 = K; float r23761 = m; float r23762 = n; float r23763 = r23761 + r23762; float r23764 = r23760 * r23763; float r23765 = 2; float r23766 = r23764 / r23765; float r23767 = M; float r23768 = r23766 - r23767; float r23769 = cos(r23768); float r23770 = r23763 / r23765; float r23771 = r23770 - r23767; float r23772 = pow(r23771, r23765); float r23773 = -r23772; float r23774 = l; float r23775 = r23761 - r23762; float r23776 = fabs(r23775); float r23777 = r23774 - r23776; float r23778 = r23773 - r23777; float r23779 = exp(r23778); float r23780 = r23769 * r23779; return r23780; } double f_id(double K, double m, double n, double M, double l) { double r23781 = K; double r23782 = m; double r23783 = n; double r23784 = r23782 + r23783; double r23785 = r23781 * r23784; double r23786 = 2; double r23787 = r23785 / r23786; double r23788 = M; double r23789 = r23787 - r23788; double r23790 = cos(r23789); double r23791 = r23784 / r23786; double r23792 = r23791 - r23788; double r23793 = pow(r23792, r23786); double r23794 = -r23793; double r23795 = l; double r23796 = r23782 - r23783; double r23797 = fabs(r23796); double r23798 = r23795 - r23797; double r23799 = r23794 - r23798; double r23800 = exp(r23799); double r23801 = r23790 * r23800; return r23801; } double f_of(float __attribute__((unused)) K, float m, float n, float M, float l) { float r23802 = 1; float r23803 = m; float r23804 = n; float r23805 = r23803 + r23804; float r23806 = 2; float r23807 = r23805 / r23806; float r23808 = M; float r23809 = r23807 - r23808; float r23810 = pow(r23809, r23806); float r23811 = -r23810; float r23812 = l; float r23813 = r23803 - r23804; float r23814 = fabs(r23813); float r23815 = r23812 - r23814; float r23816 = r23811 - r23815; float r23817 = exp(r23816); float r23818 = exp(r23817); float r23819 = log(r23818); float r23820 = cbrt(r23819); float r23821 = r23820 * r23820; float r23822 = r23821 * r23820; float r23823 = r23802 * r23822; return r23823; } double f_od(double __attribute__((unused)) K, double m, double n, double M, double l) { double r23824 = 1; double r23825 = m; double r23826 = n; double r23827 = r23825 + r23826; double r23828 = 2; double r23829 = r23827 / r23828; double r23830 = M; double r23831 = r23829 - r23830; double r23832 = pow(r23831, r23828); double r23833 = -r23832; double r23834 = l; double r23835 = r23825 - r23826; double r23836 = fabs(r23835); double r23837 = r23834 - r23836; double r23838 = r23833 - r23837; double r23839 = exp(r23838); double r23840 = exp(r23839); double r23841 = log(r23840); double r23842 = cbrt(r23841); double r23843 = r23842 * r23842; double r23844 = r23843 * r23842; double r23845 = r23824 * r23844; return r23845; } 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 r23846, r23847, r23848, r23849, r23850, r23851, r23852, r23853, r23854, r23855, r23856, r23857, r23858, r23859, r23860, r23861, r23862, r23863, r23864, r23865, r23866; void setup_mpfr_f_im() { mpfr_set_default_prec(1616); mpfr_init(r23846); mpfr_init(r23847); mpfr_init(r23848); mpfr_init(r23849); mpfr_init(r23850); mpfr_init_set_str(r23851, "2", 10, MPFR_RNDN); mpfr_init(r23852); mpfr_init(r23853); mpfr_init(r23854); mpfr_init(r23855); mpfr_init(r23856); mpfr_init(r23857); mpfr_init(r23858); mpfr_init(r23859); mpfr_init(r23860); mpfr_init(r23861); mpfr_init(r23862); mpfr_init(r23863); mpfr_init(r23864); mpfr_init(r23865); mpfr_init(r23866); } double f_im(double K, double m, double n, double M, double l) { mpfr_set_d(r23846, K, MPFR_RNDN); mpfr_set_d(r23847, m, MPFR_RNDN); mpfr_set_d(r23848, n, MPFR_RNDN); mpfr_add(r23849, r23847, r23848, MPFR_RNDN); mpfr_mul(r23850, r23846, r23849, MPFR_RNDN); ; mpfr_div(r23852, r23850, r23851, MPFR_RNDN); mpfr_set_d(r23853, M, MPFR_RNDN); mpfr_sub(r23854, r23852, r23853, MPFR_RNDN); mpfr_cos(r23855, r23854, MPFR_RNDN); mpfr_div(r23856, r23849, r23851, MPFR_RNDN); mpfr_sub(r23857, r23856, r23853, MPFR_RNDN); mpfr_pow(r23858, r23857, r23851, MPFR_RNDN); mpfr_neg(r23859, r23858, MPFR_RNDN); mpfr_set_d(r23860, l, MPFR_RNDN); mpfr_sub(r23861, r23847, r23848, MPFR_RNDN); mpfr_abs(r23862, r23861, MPFR_RNDN); mpfr_sub(r23863, r23860, r23862, MPFR_RNDN); mpfr_sub(r23864, r23859, r23863, MPFR_RNDN); mpfr_exp(r23865, r23864, MPFR_RNDN); mpfr_mul(r23866, r23855, r23865, MPFR_RNDN); return mpfr_get_d(r23866, MPFR_RNDN); } static mpfr_t r23867, r23868, r23869, r23870, r23871, r23872, r23873, r23874, r23875, r23876, r23877, r23878, r23879, r23880, r23881, r23882, r23883, r23884, r23885, r23886, r23887, r23888; void setup_mpfr_f_fm() { mpfr_set_default_prec(1616); mpfr_init_set_str(r23867, "1", 10, MPFR_RNDN); mpfr_init(r23868); mpfr_init(r23869); mpfr_init(r23870); mpfr_init_set_str(r23871, "2", 10, MPFR_RNDN); mpfr_init(r23872); mpfr_init(r23873); mpfr_init(r23874); mpfr_init(r23875); mpfr_init(r23876); mpfr_init(r23877); mpfr_init(r23878); mpfr_init(r23879); mpfr_init(r23880); mpfr_init(r23881); mpfr_init(r23882); mpfr_init(r23883); mpfr_init(r23884); mpfr_init(r23885); mpfr_init(r23886); mpfr_init(r23887); mpfr_init(r23888); } double f_fm(double __attribute__((unused)) K, double m, double n, double M, double l) { ; mpfr_set_d(r23868, m, MPFR_RNDN); mpfr_set_d(r23869, n, MPFR_RNDN); mpfr_add(r23870, r23868, r23869, MPFR_RNDN); ; mpfr_div(r23872, r23870, r23871, MPFR_RNDN); mpfr_set_d(r23873, M, MPFR_RNDN); mpfr_sub(r23874, r23872, r23873, MPFR_RNDN); mpfr_pow(r23875, r23874, r23871, MPFR_RNDN); mpfr_neg(r23876, r23875, MPFR_RNDN); mpfr_set_d(r23877, l, MPFR_RNDN); mpfr_sub(r23878, r23868, r23869, MPFR_RNDN); mpfr_abs(r23879, r23878, MPFR_RNDN); mpfr_sub(r23880, r23877, r23879, MPFR_RNDN); mpfr_sub(r23881, r23876, r23880, MPFR_RNDN); mpfr_exp(r23882, r23881, MPFR_RNDN); mpfr_exp(r23883, r23882, MPFR_RNDN); mpfr_log(r23884, r23883, MPFR_RNDN); mpfr_cbrt(r23885, r23884, MPFR_RNDN); mpfr_mul(r23886, r23885, r23885, MPFR_RNDN); mpfr_mul(r23887, r23886, r23885, MPFR_RNDN); mpfr_mul(r23888, r23867, r23887, MPFR_RNDN); return mpfr_get_d(r23888, MPFR_RNDN); } static mpfr_t r23889, r23890, r23891, r23892, r23893, r23894, r23895, r23896, r23897, r23898, r23899, r23900, r23901, r23902, r23903, r23904, r23905, r23906, r23907, r23908, r23909, r23910; void setup_mpfr_f_dm() { mpfr_set_default_prec(1616); mpfr_init_set_str(r23889, "1", 10, MPFR_RNDN); mpfr_init(r23890); mpfr_init(r23891); mpfr_init(r23892); mpfr_init_set_str(r23893, "2", 10, MPFR_RNDN); mpfr_init(r23894); mpfr_init(r23895); mpfr_init(r23896); mpfr_init(r23897); mpfr_init(r23898); mpfr_init(r23899); mpfr_init(r23900); mpfr_init(r23901); mpfr_init(r23902); mpfr_init(r23903); mpfr_init(r23904); mpfr_init(r23905); mpfr_init(r23906); mpfr_init(r23907); mpfr_init(r23908); mpfr_init(r23909); mpfr_init(r23910); } double f_dm(double __attribute__((unused)) K, double m, double n, double M, double l) { ; mpfr_set_d(r23890, m, MPFR_RNDN); mpfr_set_d(r23891, n, MPFR_RNDN); mpfr_add(r23892, r23890, r23891, MPFR_RNDN); ; mpfr_div(r23894, r23892, r23893, MPFR_RNDN); mpfr_set_d(r23895, M, MPFR_RNDN); mpfr_sub(r23896, r23894, r23895, MPFR_RNDN); mpfr_pow(r23897, r23896, r23893, MPFR_RNDN); mpfr_neg(r23898, r23897, MPFR_RNDN); mpfr_set_d(r23899, l, MPFR_RNDN); mpfr_sub(r23900, r23890, r23891, MPFR_RNDN); mpfr_abs(r23901, r23900, MPFR_RNDN); mpfr_sub(r23902, r23899, r23901, MPFR_RNDN); mpfr_sub(r23903, r23898, r23902, MPFR_RNDN); mpfr_exp(r23904, r23903, MPFR_RNDN); mpfr_exp(r23905, r23904, MPFR_RNDN); mpfr_log(r23906, r23905, MPFR_RNDN); mpfr_cbrt(r23907, r23906, MPFR_RNDN); mpfr_mul(r23908, r23907, r23907, MPFR_RNDN); mpfr_mul(r23909, r23908, r23907, MPFR_RNDN); mpfr_mul(r23910, r23889, r23909, MPFR_RNDN); return mpfr_get_d(r23910, MPFR_RNDN); }