#include #include #include #include #include char *name = "NMSE problem 3.3.7"; double f_if(float x) { float r15945 = x; float r15946 = exp(r15945); float r15947 = 2.0f; float r15948 = r15946 - r15947; float r15949 = -r15945; float r15950 = exp(r15949); float r15951 = r15948 + r15950; return r15951; } double f_id(double x) { double r15952 = x; double r15953 = exp(r15952); double r15954 = 2.0; double r15955 = r15953 - r15954; double r15956 = -r15952; double r15957 = exp(r15956); double r15958 = r15955 + r15957; return r15958; } double f_of(float x) { float r15959 = x; float r15960 = -0.12694145739078522f; bool r15961 = r15959 <= r15960; float r15962 = exp(r15959); float r15963 = r15962 * r15962; float r15964 = 2.0f; float r15965 = r15964 * r15964; float r15966 = r15963 - r15965; float r15967 = r15966 * r15962; float r15968 = r15962 + r15964; float r15969 = 1.0f; float r15970 = r15968 * r15969; float r15971 = r15967 + r15970; float r15972 = r15968 * r15962; float r15973 = r15971 / r15972; float r15974 = 6.0f; float r15975 = pow(r15959, r15974); float r15976 = 0.0027777778450399637f; float r15977 = 4.0f; float r15978 = pow(r15959, r15977); float r15979 = 0.0833333358168602f; float r15980 = r15959 * r15959; float r15981 = fma(r15978, r15979, r15980); float r15982 = fma(r15975, r15976, r15981); float r15983 = r15961 ? r15973 : r15982; return r15983; } double f_od(double x) { double r15984 = x; double r15985 = -0.12694145739078522; bool r15986 = r15984 <= r15985; double r15987 = exp(r15984); double r15988 = r15987 * r15987; double r15989 = 2.0; double r15990 = r15989 * r15989; double r15991 = r15988 - r15990; double r15992 = r15991 * r15987; double r15993 = r15987 + r15989; double r15994 = 1.0; double r15995 = r15993 * r15994; double r15996 = r15992 + r15995; double r15997 = r15993 * r15987; double r15998 = r15996 / r15997; double r15999 = 6.0; double r16000 = pow(r15984, r15999); double r16001 = 0.0027777778450399637; double r16002 = 4.0; double r16003 = pow(r15984, r16002); double r16004 = 0.0833333358168602; double r16005 = r15984 * r15984; double r16006 = fma(r16003, r16004, r16005); double r16007 = fma(r16000, r16001, r16006); double r16008 = r15986 ? r15998 : r16007; return r16008; } 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 r16009, r16010, r16011, r16012, r16013, r16014, r16015; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16009); mpfr_init(r16010); mpfr_init_set_str(r16011, "2", 10, MPFR_RNDN); mpfr_init(r16012); mpfr_init(r16013); mpfr_init(r16014); mpfr_init(r16015); } double f_im(double x) { mpfr_set_d(r16009, x, MPFR_RNDN); mpfr_exp(r16010, r16009, MPFR_RNDN); ; mpfr_sub(r16012, r16010, r16011, MPFR_RNDN); mpfr_neg(r16013, r16009, MPFR_RNDN); mpfr_exp(r16014, r16013, MPFR_RNDN); mpfr_add(r16015, r16012, r16014, MPFR_RNDN); return mpfr_get_d(r16015, MPFR_RNDN); } static mpfr_t r16016, r16017, r16018, r16019, r16020, r16021, r16022, r16023, r16024, r16025, r16026, r16027, r16028, r16029, r16030, r16031, r16032, r16033, r16034, r16035, r16036, r16037, r16038, r16039, r16040; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16016); mpfr_init_set_str(r16017, "-0.12694146f0", 10, MPFR_RNDN); mpfr_init(r16018); mpfr_init(r16019); mpfr_init(r16020); mpfr_init_set_str(r16021, "2", 10, MPFR_RNDN); mpfr_init(r16022); mpfr_init(r16023); mpfr_init(r16024); mpfr_init(r16025); mpfr_init_set_str(r16026, "1", 10, MPFR_RNDN); mpfr_init(r16027); mpfr_init(r16028); mpfr_init(r16029); mpfr_init(r16030); mpfr_init_set_str(r16031, "6", 10, MPFR_RNDN); mpfr_init(r16032); mpfr_init_set_str(r16033, "1/360", 10, MPFR_RNDN); mpfr_init_set_str(r16034, "4", 10, MPFR_RNDN); mpfr_init(r16035); mpfr_init_set_str(r16036, "1/12", 10, MPFR_RNDN); mpfr_init(r16037); mpfr_init(r16038); mpfr_init(r16039); mpfr_init(r16040); } double f_fm(double x) { mpfr_set_d(r16016, x, MPFR_RNDN); ; mpfr_set_si(r16018, mpfr_cmp(r16016, r16017) <= 0, MPFR_RNDN); mpfr_exp(r16019, r16016, MPFR_RNDN); mpfr_sqr(r16020, r16019, MPFR_RNDN); ; mpfr_sqr(r16022, r16021, MPFR_RNDN); mpfr_sub(r16023, r16020, r16022, MPFR_RNDN); mpfr_mul(r16024, r16023, r16019, MPFR_RNDN); mpfr_add(r16025, r16019, r16021, MPFR_RNDN); ; mpfr_mul(r16027, r16025, r16026, MPFR_RNDN); mpfr_add(r16028, r16024, r16027, MPFR_RNDN); mpfr_mul(r16029, r16025, r16019, MPFR_RNDN); mpfr_div(r16030, r16028, r16029, MPFR_RNDN); ; mpfr_pow(r16032, r16016, r16031, MPFR_RNDN); ; ; mpfr_pow(r16035, r16016, r16034, MPFR_RNDN); ; mpfr_sqr(r16037, r16016, MPFR_RNDN); mpfr_fma(r16038, r16035, r16036, r16037, MPFR_RNDN); mpfr_fma(r16039, r16032, r16033, r16038, MPFR_RNDN); if (mpfr_get_si(r16018, MPFR_RNDN)) { mpfr_set(r16040, r16030, MPFR_RNDN); } else { mpfr_set(r16040, r16039, MPFR_RNDN); }; return mpfr_get_d(r16040, MPFR_RNDN); } static mpfr_t r16041, r16042, r16043, r16044, r16045, r16046, r16047, r16048, r16049, r16050, r16051, r16052, r16053, r16054, r16055, r16056, r16057, r16058, r16059, r16060, r16061, r16062, r16063, r16064, r16065; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16041); mpfr_init_set_str(r16042, "-0.12694146f0", 10, MPFR_RNDN); mpfr_init(r16043); mpfr_init(r16044); mpfr_init(r16045); mpfr_init_set_str(r16046, "2", 10, MPFR_RNDN); mpfr_init(r16047); mpfr_init(r16048); mpfr_init(r16049); mpfr_init(r16050); mpfr_init_set_str(r16051, "1", 10, MPFR_RNDN); mpfr_init(r16052); mpfr_init(r16053); mpfr_init(r16054); mpfr_init(r16055); mpfr_init_set_str(r16056, "6", 10, MPFR_RNDN); mpfr_init(r16057); mpfr_init_set_str(r16058, "1/360", 10, MPFR_RNDN); mpfr_init_set_str(r16059, "4", 10, MPFR_RNDN); mpfr_init(r16060); mpfr_init_set_str(r16061, "1/12", 10, MPFR_RNDN); mpfr_init(r16062); mpfr_init(r16063); mpfr_init(r16064); mpfr_init(r16065); } double f_dm(double x) { mpfr_set_d(r16041, x, MPFR_RNDN); ; mpfr_set_si(r16043, mpfr_cmp(r16041, r16042) <= 0, MPFR_RNDN); mpfr_exp(r16044, r16041, MPFR_RNDN); mpfr_sqr(r16045, r16044, MPFR_RNDN); ; mpfr_sqr(r16047, r16046, MPFR_RNDN); mpfr_sub(r16048, r16045, r16047, MPFR_RNDN); mpfr_mul(r16049, r16048, r16044, MPFR_RNDN); mpfr_add(r16050, r16044, r16046, MPFR_RNDN); ; mpfr_mul(r16052, r16050, r16051, MPFR_RNDN); mpfr_add(r16053, r16049, r16052, MPFR_RNDN); mpfr_mul(r16054, r16050, r16044, MPFR_RNDN); mpfr_div(r16055, r16053, r16054, MPFR_RNDN); ; mpfr_pow(r16057, r16041, r16056, MPFR_RNDN); ; ; mpfr_pow(r16060, r16041, r16059, MPFR_RNDN); ; mpfr_sqr(r16062, r16041, MPFR_RNDN); mpfr_fma(r16063, r16060, r16061, r16062, MPFR_RNDN); mpfr_fma(r16064, r16057, r16058, r16063, MPFR_RNDN); if (mpfr_get_si(r16043, MPFR_RNDN)) { mpfr_set(r16065, r16055, MPFR_RNDN); } else { mpfr_set(r16065, r16064, MPFR_RNDN); }; return mpfr_get_d(r16065, MPFR_RNDN); }