#include #include #include #include #include char *name = "NMSE problem 3.4.5"; double f_if(float x) { float r14899 = x; float r14900 = sin(r14899); float r14901 = r14899 - r14900; float r14902 = tan(r14899); float r14903 = r14899 - r14902; float r14904 = r14901 / r14903; return r14904; } double f_id(double x) { double r14905 = x; double r14906 = sin(r14905); double r14907 = r14905 - r14906; double r14908 = tan(r14905); double r14909 = r14905 - r14908; double r14910 = r14907 / r14909; return r14910; } double f_of(float x) { float r14911 = x; float r14912 = -8.61646249982272e-12f; bool r14913 = r14911 <= r14912; float r14914 = sin(r14911); float r14915 = r14911 - r14914; float r14916 = tan(r14911); float r14917 = r14911 - r14916; float r14918 = r14915 / r14917; float r14919 = 8.584780819621077f; bool r14920 = r14911 <= r14919; float r14921 = 0.225f; float r14922 = r14911 * r14921; float r14923 = r14922 * r14911; float r14924 = 0.5f; float r14925 = r14923 - r14924; float r14926 = 0.009642857142857142f; float r14927 = 4.0f; float r14928 = pow(r14911, r14927); float r14929 = r14926 * r14928; float r14930 = r14925 - r14929; float r14931 = r14911 / r14917; float r14932 = r14914 / r14917; float r14933 = r14931 - r14932; float r14934 = cbrt(r14933); float r14935 = r14934 * (r14934 * r14934); float r14936 = cbrt(r14935); float r14937 = r14936 * (r14936 * r14936); float r14938 = r14920 ? r14930 : r14937; float r14939 = r14913 ? r14918 : r14938; return r14939; } double f_od(double x) { double r14940 = x; double r14941 = -8.61646249982272e-12; bool r14942 = r14940 <= r14941; double r14943 = sin(r14940); double r14944 = r14940 - r14943; double r14945 = tan(r14940); double r14946 = r14940 - r14945; double r14947 = r14944 / r14946; double r14948 = 8.584780819621077; bool r14949 = r14940 <= r14948; double r14950 = 0.225; double r14951 = r14940 * r14950; double r14952 = r14951 * r14940; double r14953 = 0.5; double r14954 = r14952 - r14953; double r14955 = 0.009642857142857142; double r14956 = 4.0; double r14957 = pow(r14940, r14956); double r14958 = r14955 * r14957; double r14959 = r14954 - r14958; double r14960 = r14940 / r14946; double r14961 = r14943 / r14946; double r14962 = r14960 - r14961; double r14963 = cbrt(r14962); double r14964 = r14963 * (r14963 * r14963); double r14965 = cbrt(r14964); double r14966 = r14965 * (r14965 * r14965); double r14967 = r14949 ? r14959 : r14966; double r14968 = r14942 ? r14947 : r14967; return r14968; } 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 r14969, r14970, r14971, r14972, r14973, r14974; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r14969); mpfr_init(r14970); mpfr_init(r14971); mpfr_init(r14972); mpfr_init(r14973); mpfr_init(r14974); } double f_im(double x) { mpfr_set_d(r14969, x, MPFR_RNDN); mpfr_sin(r14970, r14969, MPFR_RNDN); mpfr_sub(r14971, r14969, r14970, MPFR_RNDN); mpfr_tan(r14972, r14969, MPFR_RNDN); mpfr_sub(r14973, r14969, r14972, MPFR_RNDN); mpfr_div(r14974, r14971, r14973, MPFR_RNDN); return mpfr_get_d(r14974, MPFR_RNDN); } static mpfr_t r14975, r14976, r14977, r14978, r14979, r14980, r14981, r14982, r14983, r14984, r14985, r14986, r14987, r14988, r14989, r14990, r14991, r14992, r14993, r14994, r14995, r14996, r14997, r14998, r14999, r15000, r15001, r15002, r15003; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r14975); mpfr_init_set_str(r14976, "-8.61646249982272e-12", 10, MPFR_RNDN); mpfr_init(r14977); mpfr_init(r14978); mpfr_init(r14979); mpfr_init(r14980); mpfr_init(r14981); mpfr_init(r14982); mpfr_init_set_str(r14983, "8.584780819621077", 10, MPFR_RNDN); mpfr_init(r14984); mpfr_init_set_str(r14985, "9/40", 10, MPFR_RNDN); mpfr_init(r14986); mpfr_init(r14987); mpfr_init_set_str(r14988, "1/2", 10, MPFR_RNDN); mpfr_init(r14989); mpfr_init_set_str(r14990, "27/2800", 10, MPFR_RNDN); mpfr_init_set_str(r14991, "4", 10, MPFR_RNDN); mpfr_init(r14992); mpfr_init(r14993); mpfr_init(r14994); mpfr_init(r14995); mpfr_init(r14996); mpfr_init(r14997); mpfr_init(r14998); mpfr_init(r14999); mpfr_init(r15000); mpfr_init(r15001); mpfr_init(r15002); mpfr_init(r15003); } double f_fm(double x) { mpfr_set_d(r14975, x, MPFR_RNDN); ; mpfr_set_si(r14977, mpfr_cmp(r14975, r14976) <= 0, MPFR_RNDN); mpfr_sin(r14978, r14975, MPFR_RNDN); mpfr_sub(r14979, r14975, r14978, MPFR_RNDN); mpfr_tan(r14980, r14975, MPFR_RNDN); mpfr_sub(r14981, r14975, r14980, MPFR_RNDN); mpfr_div(r14982, r14979, r14981, MPFR_RNDN); ; mpfr_set_si(r14984, mpfr_cmp(r14975, r14983) <= 0, MPFR_RNDN); ; mpfr_mul(r14986, r14975, r14985, MPFR_RNDN); mpfr_mul(r14987, r14986, r14975, MPFR_RNDN); ; mpfr_sub(r14989, r14987, r14988, MPFR_RNDN); ; ; mpfr_pow(r14992, r14975, r14991, MPFR_RNDN); mpfr_mul(r14993, r14990, r14992, MPFR_RNDN); mpfr_sub(r14994, r14989, r14993, MPFR_RNDN); mpfr_div(r14995, r14975, r14981, MPFR_RNDN); mpfr_div(r14996, r14978, r14981, MPFR_RNDN); mpfr_sub(r14997, r14995, r14996, MPFR_RNDN); mpfr_cbrt(r14998, r14997, MPFR_RNDN); mpfr_mul(r14999, r14998, r14998, MPFR_RNDN); mpfr_mul(r14999, r14999, r14998, MPFR_RNDN); mpfr_cbrt(r15000, r14999, MPFR_RNDN); mpfr_mul(r15001, r15000, r15000, MPFR_RNDN); mpfr_mul(r15001, r15001, r15000, MPFR_RNDN); if (mpfr_get_si(r14984, MPFR_RNDN)) { mpfr_set(r15002, r14994, MPFR_RNDN); } else { mpfr_set(r15002, r15001, MPFR_RNDN); }; if (mpfr_get_si(r14977, MPFR_RNDN)) { mpfr_set(r15003, r14982, MPFR_RNDN); } else { mpfr_set(r15003, r15002, MPFR_RNDN); }; return mpfr_get_d(r15003, MPFR_RNDN); } static mpfr_t r15004, r15005, r15006, r15007, r15008, r15009, r15010, r15011, r15012, r15013, r15014, r15015, r15016, r15017, r15018, r15019, r15020, r15021, r15022, r15023, r15024, r15025, r15026, r15027, r15028, r15029, r15030, r15031, r15032; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15004); mpfr_init_set_str(r15005, "-8.61646249982272e-12", 10, MPFR_RNDN); mpfr_init(r15006); mpfr_init(r15007); mpfr_init(r15008); mpfr_init(r15009); mpfr_init(r15010); mpfr_init(r15011); mpfr_init_set_str(r15012, "8.584780819621077", 10, MPFR_RNDN); mpfr_init(r15013); mpfr_init_set_str(r15014, "9/40", 10, MPFR_RNDN); mpfr_init(r15015); mpfr_init(r15016); mpfr_init_set_str(r15017, "1/2", 10, MPFR_RNDN); mpfr_init(r15018); mpfr_init_set_str(r15019, "27/2800", 10, MPFR_RNDN); mpfr_init_set_str(r15020, "4", 10, MPFR_RNDN); mpfr_init(r15021); mpfr_init(r15022); mpfr_init(r15023); mpfr_init(r15024); mpfr_init(r15025); mpfr_init(r15026); mpfr_init(r15027); mpfr_init(r15028); mpfr_init(r15029); mpfr_init(r15030); mpfr_init(r15031); mpfr_init(r15032); } double f_dm(double x) { mpfr_set_d(r15004, x, MPFR_RNDN); ; mpfr_set_si(r15006, mpfr_cmp(r15004, r15005) <= 0, MPFR_RNDN); mpfr_sin(r15007, r15004, MPFR_RNDN); mpfr_sub(r15008, r15004, r15007, MPFR_RNDN); mpfr_tan(r15009, r15004, MPFR_RNDN); mpfr_sub(r15010, r15004, r15009, MPFR_RNDN); mpfr_div(r15011, r15008, r15010, MPFR_RNDN); ; mpfr_set_si(r15013, mpfr_cmp(r15004, r15012) <= 0, MPFR_RNDN); ; mpfr_mul(r15015, r15004, r15014, MPFR_RNDN); mpfr_mul(r15016, r15015, r15004, MPFR_RNDN); ; mpfr_sub(r15018, r15016, r15017, MPFR_RNDN); ; ; mpfr_pow(r15021, r15004, r15020, MPFR_RNDN); mpfr_mul(r15022, r15019, r15021, MPFR_RNDN); mpfr_sub(r15023, r15018, r15022, MPFR_RNDN); mpfr_div(r15024, r15004, r15010, MPFR_RNDN); mpfr_div(r15025, r15007, r15010, MPFR_RNDN); mpfr_sub(r15026, r15024, r15025, MPFR_RNDN); mpfr_cbrt(r15027, r15026, MPFR_RNDN); mpfr_mul(r15028, r15027, r15027, MPFR_RNDN); mpfr_mul(r15028, r15028, r15027, MPFR_RNDN); mpfr_cbrt(r15029, r15028, MPFR_RNDN); mpfr_mul(r15030, r15029, r15029, MPFR_RNDN); mpfr_mul(r15030, r15030, r15029, MPFR_RNDN); if (mpfr_get_si(r15013, MPFR_RNDN)) { mpfr_set(r15031, r15023, MPFR_RNDN); } else { mpfr_set(r15031, r15030, MPFR_RNDN); }; if (mpfr_get_si(r15006, MPFR_RNDN)) { mpfr_set(r15032, r15011, MPFR_RNDN); } else { mpfr_set(r15032, r15031, MPFR_RNDN); }; return mpfr_get_d(r15032, MPFR_RNDN); }