#include #include #include #include #include char *name = "NMSE problem 3.4.5"; double f_if(float x) { float r14839 = x; float r14840 = sin(r14839); float r14841 = r14839 - r14840; float r14842 = tan(r14839); float r14843 = r14839 - r14842; float r14844 = r14841 / r14843; return r14844; } double f_id(double x) { double r14845 = x; double r14846 = sin(r14845); double r14847 = r14845 - r14846; double r14848 = tan(r14845); double r14849 = r14845 - r14848; double r14850 = r14847 / r14849; return r14850; } double f_of(float x) { float r14851 = x; float r14852 = -1.7643470284436324e-06f; bool r14853 = r14851 <= r14852; float r14854 = tan(r14851); float r14855 = r14851 - r14854; float r14856 = r14851 / r14855; float r14857 = sin(r14851); float r14858 = r14857 / r14855; float r14859 = r14856 - r14858; float r14860 = 5.107891217442523e+18f; bool r14861 = r14851 <= r14860; float r14862 = 0.225f; float r14863 = r14851 * r14851; float r14864 = r14862 * r14863; float r14865 = exp(r14864); float r14866 = 0.009642857142857142f; float r14867 = 4.0f; float r14868 = pow(r14851, r14867); float r14869 = r14866 * r14868; float r14870 = 0.5f; float r14871 = r14869 + r14870; float r14872 = exp(r14871); float r14873 = r14865 / r14872; float r14874 = log(r14873); float r14875 = r14861 ? r14874 : r14859; float r14876 = r14853 ? r14859 : r14875; return r14876; } double f_od(double x) { double r14877 = x; double r14878 = -1.7643470284436324e-06; bool r14879 = r14877 <= r14878; double r14880 = tan(r14877); double r14881 = r14877 - r14880; double r14882 = r14877 / r14881; double r14883 = sin(r14877); double r14884 = r14883 / r14881; double r14885 = r14882 - r14884; double r14886 = 5.107891217442523e+18; bool r14887 = r14877 <= r14886; double r14888 = 0.225; double r14889 = r14877 * r14877; double r14890 = r14888 * r14889; double r14891 = exp(r14890); double r14892 = 0.009642857142857142; double r14893 = 4.0; double r14894 = pow(r14877, r14893); double r14895 = r14892 * r14894; double r14896 = 0.5; double r14897 = r14895 + r14896; double r14898 = exp(r14897); double r14899 = r14891 / r14898; double r14900 = log(r14899); double r14901 = r14887 ? r14900 : r14885; double r14902 = r14879 ? r14885 : r14901; return r14902; } 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 r14903, r14904, r14905, r14906, r14907, r14908; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r14903); mpfr_init(r14904); mpfr_init(r14905); mpfr_init(r14906); mpfr_init(r14907); mpfr_init(r14908); } double f_im(double x) { mpfr_set_d(r14903, x, MPFR_RNDN); mpfr_sin(r14904, r14903, MPFR_RNDN); mpfr_sub(r14905, r14903, r14904, MPFR_RNDN); mpfr_tan(r14906, r14903, MPFR_RNDN); mpfr_sub(r14907, r14903, r14906, MPFR_RNDN); mpfr_div(r14908, r14905, r14907, MPFR_RNDN); return mpfr_get_d(r14908, MPFR_RNDN); } static mpfr_t r14909, r14910, r14911, r14912, r14913, r14914, r14915, r14916, r14917, r14918, r14919, r14920, r14921, r14922, r14923, r14924, r14925, r14926, r14927, r14928, r14929, r14930, r14931, r14932, r14933, r14934; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r14909); mpfr_init_set_str(r14910, "-1.7643470284436324e-06", 10, MPFR_RNDN); mpfr_init(r14911); mpfr_init(r14912); mpfr_init(r14913); mpfr_init(r14914); mpfr_init(r14915); mpfr_init(r14916); mpfr_init(r14917); mpfr_init_set_str(r14918, "5.107891217442523e+18", 10, MPFR_RNDN); mpfr_init(r14919); mpfr_init_set_str(r14920, "9/40", 10, MPFR_RNDN); mpfr_init(r14921); mpfr_init(r14922); mpfr_init(r14923); mpfr_init_set_str(r14924, "27/2800", 10, MPFR_RNDN); mpfr_init_set_str(r14925, "4", 10, MPFR_RNDN); mpfr_init(r14926); mpfr_init(r14927); mpfr_init_set_str(r14928, "1/2", 10, MPFR_RNDN); mpfr_init(r14929); mpfr_init(r14930); mpfr_init(r14931); mpfr_init(r14932); mpfr_init(r14933); mpfr_init(r14934); } double f_fm(double x) { mpfr_set_d(r14909, x, MPFR_RNDN); ; mpfr_set_si(r14911, mpfr_cmp(r14909, r14910) <= 0, MPFR_RNDN); mpfr_tan(r14912, r14909, MPFR_RNDN); mpfr_sub(r14913, r14909, r14912, MPFR_RNDN); mpfr_div(r14914, r14909, r14913, MPFR_RNDN); mpfr_sin(r14915, r14909, MPFR_RNDN); mpfr_div(r14916, r14915, r14913, MPFR_RNDN); mpfr_sub(r14917, r14914, r14916, MPFR_RNDN); ; mpfr_set_si(r14919, mpfr_cmp(r14909, r14918) <= 0, MPFR_RNDN); ; mpfr_sqr(r14921, r14909, MPFR_RNDN); mpfr_mul(r14922, r14920, r14921, MPFR_RNDN); mpfr_exp(r14923, r14922, MPFR_RNDN); ; ; mpfr_pow(r14926, r14909, r14925, MPFR_RNDN); mpfr_mul(r14927, r14924, r14926, MPFR_RNDN); ; mpfr_add(r14929, r14927, r14928, MPFR_RNDN); mpfr_exp(r14930, r14929, MPFR_RNDN); mpfr_div(r14931, r14923, r14930, MPFR_RNDN); mpfr_log(r14932, r14931, MPFR_RNDN); if (mpfr_get_si(r14919, MPFR_RNDN)) { mpfr_set(r14933, r14932, MPFR_RNDN); } else { mpfr_set(r14933, r14917, MPFR_RNDN); }; if (mpfr_get_si(r14911, MPFR_RNDN)) { mpfr_set(r14934, r14917, MPFR_RNDN); } else { mpfr_set(r14934, r14933, MPFR_RNDN); }; return mpfr_get_d(r14934, MPFR_RNDN); } static mpfr_t r14935, r14936, r14937, r14938, r14939, r14940, r14941, r14942, r14943, r14944, r14945, r14946, r14947, r14948, r14949, r14950, r14951, r14952, r14953, r14954, r14955, r14956, r14957, r14958, r14959, r14960; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r14935); mpfr_init_set_str(r14936, "-1.7643470284436324e-06", 10, MPFR_RNDN); mpfr_init(r14937); mpfr_init(r14938); mpfr_init(r14939); mpfr_init(r14940); mpfr_init(r14941); mpfr_init(r14942); mpfr_init(r14943); mpfr_init_set_str(r14944, "5.107891217442523e+18", 10, MPFR_RNDN); mpfr_init(r14945); mpfr_init_set_str(r14946, "9/40", 10, MPFR_RNDN); mpfr_init(r14947); mpfr_init(r14948); mpfr_init(r14949); mpfr_init_set_str(r14950, "27/2800", 10, MPFR_RNDN); mpfr_init_set_str(r14951, "4", 10, MPFR_RNDN); mpfr_init(r14952); mpfr_init(r14953); mpfr_init_set_str(r14954, "1/2", 10, MPFR_RNDN); mpfr_init(r14955); mpfr_init(r14956); mpfr_init(r14957); mpfr_init(r14958); mpfr_init(r14959); mpfr_init(r14960); } double f_dm(double x) { mpfr_set_d(r14935, x, MPFR_RNDN); ; mpfr_set_si(r14937, mpfr_cmp(r14935, r14936) <= 0, MPFR_RNDN); mpfr_tan(r14938, r14935, MPFR_RNDN); mpfr_sub(r14939, r14935, r14938, MPFR_RNDN); mpfr_div(r14940, r14935, r14939, MPFR_RNDN); mpfr_sin(r14941, r14935, MPFR_RNDN); mpfr_div(r14942, r14941, r14939, MPFR_RNDN); mpfr_sub(r14943, r14940, r14942, MPFR_RNDN); ; mpfr_set_si(r14945, mpfr_cmp(r14935, r14944) <= 0, MPFR_RNDN); ; mpfr_sqr(r14947, r14935, MPFR_RNDN); mpfr_mul(r14948, r14946, r14947, MPFR_RNDN); mpfr_exp(r14949, r14948, MPFR_RNDN); ; ; mpfr_pow(r14952, r14935, r14951, MPFR_RNDN); mpfr_mul(r14953, r14950, r14952, MPFR_RNDN); ; mpfr_add(r14955, r14953, r14954, MPFR_RNDN); mpfr_exp(r14956, r14955, MPFR_RNDN); mpfr_div(r14957, r14949, r14956, MPFR_RNDN); mpfr_log(r14958, r14957, MPFR_RNDN); if (mpfr_get_si(r14945, MPFR_RNDN)) { mpfr_set(r14959, r14958, MPFR_RNDN); } else { mpfr_set(r14959, r14943, MPFR_RNDN); }; if (mpfr_get_si(r14937, MPFR_RNDN)) { mpfr_set(r14960, r14943, MPFR_RNDN); } else { mpfr_set(r14960, r14959, MPFR_RNDN); }; return mpfr_get_d(r14960, MPFR_RNDN); }