#include #include #include #include #include char *name = "NMSE problem 3.3.2"; double f_if(float x, float eps) { float r4619060 = x; float r4619061 = eps; float r4619062 = r4619060 + r4619061; float r4619063 = tan(r4619062); float r4619064 = tan(r4619060); float r4619065 = r4619063 - r4619064; return r4619065; } double f_id(double x, double eps) { double r4619066 = x; double r4619067 = eps; double r4619068 = r4619066 + r4619067; double r4619069 = tan(r4619068); double r4619070 = tan(r4619066); double r4619071 = r4619069 - r4619070; return r4619071; } double f_of(float x, float eps) { float r4619072 = eps; float r4619073 = -1.1324884610140404e-37f; bool r4619074 = r4619072 <= r4619073; float r4619075 = x; float r4619076 = cos(r4619075); float r4619077 = r4619075 + r4619072; float r4619078 = 1.0 / tan(r4619077); float r4619079 = cbrt(r4619078); float r4619080 = sin(r4619075); float r4619081 = cbrt(r4619080); float r4619082 = r4619079 * r4619081; float r4619083 = r4619082 * (r4619082 * r4619082); float r4619084 = r4619076 - r4619083; float r4619085 = r4619078 * r4619076; float r4619086 = r4619084 / r4619085; float r4619087 = -8.343087553984244e-186f; bool r4619088 = r4619072 <= r4619087; float r4619089 = r4619080 * r4619072; float r4619090 = 3.0f; float r4619091 = pow(r4619072, r4619090); float r4619092 = 4.0f; float r4619093 = pow(r4619076, r4619092); float r4619094 = r4619091 * r4619093; float r4619095 = pow(r4619080, r4619090); float r4619096 = r4619094 / r4619095; float r4619097 = r4619076 * r4619076; float r4619098 = r4619072 * r4619097; float r4619099 = r4619098 / r4619080; float r4619100 = 0.3333333333333333f; float r4619101 = r4619080 * r4619091; float r4619102 = r4619100 * r4619101; float r4619103 = 1.3333333333333333f; float r4619104 = r4619091 * r4619097; float r4619105 = r4619104 / r4619080; float r4619106 = r4619103 * r4619105; float r4619107 = r4619102 + r4619106; float r4619108 = r4619099 + r4619107; float r4619109 = r4619096 + r4619108; float r4619110 = r4619089 + r4619109; float r4619111 = r4619072 * r4619072; float r4619112 = pow(r4619076, r4619090); float r4619113 = r4619111 * r4619112; float r4619114 = r4619080 * r4619080; float r4619115 = r4619113 / r4619114; float r4619116 = r4619111 * r4619076; float r4619117 = r4619115 + r4619116; float r4619118 = r4619110 - r4619117; float r4619119 = r4619118 / r4619085; float r4619120 = -1.4036887096279737e-262f; bool r4619121 = r4619072 <= r4619120; float r4619122 = pow(r4619072, r4619092); float r4619123 = r4619075 * (r4619075 * r4619075); float r4619124 = r4619122 * r4619123; float r4619125 = r4619072 + r4619124; float r4619126 = r4619075 * r4619075; float r4619127 = r4619072 * (r4619072 * r4619072); float r4619128 = r4619126 * r4619127; float r4619129 = r4619125 + r4619128; float r4619130 = 4.0170806835946265e-42f; bool r4619131 = r4619072 <= r4619130; float r4619132 = r4619131 ? r4619119 : r4619086; float r4619133 = r4619121 ? r4619129 : r4619132; float r4619134 = r4619088 ? r4619119 : r4619133; float r4619135 = r4619074 ? r4619086 : r4619134; return r4619135; } double f_od(double x, double eps) { double r4619136 = eps; double r4619137 = -1.1324884610140404e-37; bool r4619138 = r4619136 <= r4619137; double r4619139 = x; double r4619140 = cos(r4619139); double r4619141 = r4619139 + r4619136; double r4619142 = 1.0 / tan(r4619141); double r4619143 = cbrt(r4619142); double r4619144 = sin(r4619139); double r4619145 = cbrt(r4619144); double r4619146 = r4619143 * r4619145; double r4619147 = r4619146 * (r4619146 * r4619146); double r4619148 = r4619140 - r4619147; double r4619149 = r4619142 * r4619140; double r4619150 = r4619148 / r4619149; double r4619151 = -8.343087553984244e-186; bool r4619152 = r4619136 <= r4619151; double r4619153 = r4619144 * r4619136; double r4619154 = 3.0; double r4619155 = pow(r4619136, r4619154); double r4619156 = 4.0; double r4619157 = pow(r4619140, r4619156); double r4619158 = r4619155 * r4619157; double r4619159 = pow(r4619144, r4619154); double r4619160 = r4619158 / r4619159; double r4619161 = r4619140 * r4619140; double r4619162 = r4619136 * r4619161; double r4619163 = r4619162 / r4619144; double r4619164 = 0.3333333333333333; double r4619165 = r4619144 * r4619155; double r4619166 = r4619164 * r4619165; double r4619167 = 1.3333333333333333; double r4619168 = r4619155 * r4619161; double r4619169 = r4619168 / r4619144; double r4619170 = r4619167 * r4619169; double r4619171 = r4619166 + r4619170; double r4619172 = r4619163 + r4619171; double r4619173 = r4619160 + r4619172; double r4619174 = r4619153 + r4619173; double r4619175 = r4619136 * r4619136; double r4619176 = pow(r4619140, r4619154); double r4619177 = r4619175 * r4619176; double r4619178 = r4619144 * r4619144; double r4619179 = r4619177 / r4619178; double r4619180 = r4619175 * r4619140; double r4619181 = r4619179 + r4619180; double r4619182 = r4619174 - r4619181; double r4619183 = r4619182 / r4619149; double r4619184 = -1.4036887096279737e-262; bool r4619185 = r4619136 <= r4619184; double r4619186 = pow(r4619136, r4619156); double r4619187 = r4619139 * (r4619139 * r4619139); double r4619188 = r4619186 * r4619187; double r4619189 = r4619136 + r4619188; double r4619190 = r4619139 * r4619139; double r4619191 = r4619136 * (r4619136 * r4619136); double r4619192 = r4619190 * r4619191; double r4619193 = r4619189 + r4619192; double r4619194 = 4.0170806835946265e-42; bool r4619195 = r4619136 <= r4619194; double r4619196 = r4619195 ? r4619183 : r4619150; double r4619197 = r4619185 ? r4619193 : r4619196; double r4619198 = r4619152 ? r4619183 : r4619197; double r4619199 = r4619138 ? r4619150 : r4619198; return r4619199; } 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 r4619200, r4619201, r4619202, r4619203, r4619204, r4619205; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r4619200); mpfr_init(r4619201); mpfr_init(r4619202); mpfr_init(r4619203); mpfr_init(r4619204); mpfr_init(r4619205); } double f_im(double x, double eps) { mpfr_set_d(r4619200, x, MPFR_RNDN); mpfr_set_d(r4619201, eps, MPFR_RNDN); mpfr_add(r4619202, r4619200, r4619201, MPFR_RNDN); mpfr_tan(r4619203, r4619202, MPFR_RNDN); mpfr_tan(r4619204, r4619200, MPFR_RNDN); mpfr_sub(r4619205, r4619203, r4619204, MPFR_RNDN); return mpfr_get_d(r4619205, MPFR_RNDN); } static mpfr_t r4619206, r4619207, r4619208, r4619209, r4619210, r4619211, r4619212, r4619213, r4619214, r4619215, r4619216, r4619217, r4619218, r4619219, r4619220, r4619221, r4619222, r4619223, r4619224, r4619225, r4619226, r4619227, r4619228, r4619229, r4619230, r4619231, r4619232, r4619233, r4619234, r4619235, r4619236, r4619237, r4619238, r4619239, r4619240, r4619241, r4619242, r4619243, r4619244, r4619245, r4619246, r4619247, r4619248, r4619249, r4619250, r4619251, r4619252, r4619253, r4619254, r4619255, r4619256, r4619257, r4619258, r4619259, r4619260, r4619261, r4619262, r4619263, r4619264, r4619265, r4619266, r4619267, r4619268, r4619269; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r4619206); mpfr_init_set_str(r4619207, "-1.1324884610140404e-37", 10, MPFR_RNDN); mpfr_init(r4619208); mpfr_init(r4619209); mpfr_init(r4619210); mpfr_init(r4619211); mpfr_init(r4619212); mpfr_init(r4619213); mpfr_init(r4619214); mpfr_init(r4619215); mpfr_init(r4619216); mpfr_init(r4619217); mpfr_init(r4619218); mpfr_init(r4619219); mpfr_init(r4619220); mpfr_init_set_str(r4619221, "-8.343087553984244e-186", 10, MPFR_RNDN); mpfr_init(r4619222); mpfr_init(r4619223); mpfr_init_set_str(r4619224, "3", 10, MPFR_RNDN); mpfr_init(r4619225); mpfr_init_set_str(r4619226, "4", 10, MPFR_RNDN); mpfr_init(r4619227); mpfr_init(r4619228); mpfr_init(r4619229); mpfr_init(r4619230); mpfr_init(r4619231); mpfr_init(r4619232); mpfr_init(r4619233); mpfr_init_set_str(r4619234, "1/3", 10, MPFR_RNDN); mpfr_init(r4619235); mpfr_init(r4619236); mpfr_init_set_str(r4619237, "4/3", 10, MPFR_RNDN); mpfr_init(r4619238); mpfr_init(r4619239); mpfr_init(r4619240); mpfr_init(r4619241); mpfr_init(r4619242); mpfr_init(r4619243); mpfr_init(r4619244); mpfr_init(r4619245); mpfr_init(r4619246); mpfr_init(r4619247); mpfr_init(r4619248); mpfr_init(r4619249); mpfr_init(r4619250); mpfr_init(r4619251); mpfr_init(r4619252); mpfr_init(r4619253); mpfr_init_set_str(r4619254, "-1.4036887096279737e-262", 10, MPFR_RNDN); mpfr_init(r4619255); mpfr_init(r4619256); mpfr_init(r4619257); mpfr_init(r4619258); mpfr_init(r4619259); mpfr_init(r4619260); mpfr_init(r4619261); mpfr_init(r4619262); mpfr_init(r4619263); mpfr_init_set_str(r4619264, "4.0170806835946265e-42", 10, MPFR_RNDN); mpfr_init(r4619265); mpfr_init(r4619266); mpfr_init(r4619267); mpfr_init(r4619268); mpfr_init(r4619269); } double f_fm(double x, double eps) { mpfr_set_d(r4619206, eps, MPFR_RNDN); ; mpfr_set_si(r4619208, mpfr_cmp(r4619206, r4619207) <= 0, MPFR_RNDN); mpfr_set_d(r4619209, x, MPFR_RNDN); mpfr_cos(r4619210, r4619209, MPFR_RNDN); mpfr_add(r4619211, r4619209, r4619206, MPFR_RNDN); mpfr_cot(r4619212, r4619211, MPFR_RNDN); mpfr_cbrt(r4619213, r4619212, MPFR_RNDN); mpfr_sin(r4619214, r4619209, MPFR_RNDN); mpfr_cbrt(r4619215, r4619214, MPFR_RNDN); mpfr_mul(r4619216, r4619213, r4619215, MPFR_RNDN); mpfr_mul(r4619217, r4619216, r4619216, MPFR_RNDN); mpfr_mul(r4619217, r4619217, r4619216, MPFR_RNDN); mpfr_sub(r4619218, r4619210, r4619217, MPFR_RNDN); mpfr_mul(r4619219, r4619212, r4619210, MPFR_RNDN); mpfr_div(r4619220, r4619218, r4619219, MPFR_RNDN); ; mpfr_set_si(r4619222, mpfr_cmp(r4619206, r4619221) <= 0, MPFR_RNDN); mpfr_mul(r4619223, r4619214, r4619206, MPFR_RNDN); ; mpfr_pow(r4619225, r4619206, r4619224, MPFR_RNDN); ; mpfr_pow(r4619227, r4619210, r4619226, MPFR_RNDN); mpfr_mul(r4619228, r4619225, r4619227, MPFR_RNDN); mpfr_pow(r4619229, r4619214, r4619224, MPFR_RNDN); mpfr_div(r4619230, r4619228, r4619229, MPFR_RNDN); mpfr_sqr(r4619231, r4619210, MPFR_RNDN); mpfr_mul(r4619232, r4619206, r4619231, MPFR_RNDN); mpfr_div(r4619233, r4619232, r4619214, MPFR_RNDN); ; mpfr_mul(r4619235, r4619214, r4619225, MPFR_RNDN); mpfr_mul(r4619236, r4619234, r4619235, MPFR_RNDN); ; mpfr_mul(r4619238, r4619225, r4619231, MPFR_RNDN); mpfr_div(r4619239, r4619238, r4619214, MPFR_RNDN); mpfr_mul(r4619240, r4619237, r4619239, MPFR_RNDN); mpfr_add(r4619241, r4619236, r4619240, MPFR_RNDN); mpfr_add(r4619242, r4619233, r4619241, MPFR_RNDN); mpfr_add(r4619243, r4619230, r4619242, MPFR_RNDN); mpfr_add(r4619244, r4619223, r4619243, MPFR_RNDN); mpfr_sqr(r4619245, r4619206, MPFR_RNDN); mpfr_pow(r4619246, r4619210, r4619224, MPFR_RNDN); mpfr_mul(r4619247, r4619245, r4619246, MPFR_RNDN); mpfr_sqr(r4619248, r4619214, MPFR_RNDN); mpfr_div(r4619249, r4619247, r4619248, MPFR_RNDN); mpfr_mul(r4619250, r4619245, r4619210, MPFR_RNDN); mpfr_add(r4619251, r4619249, r4619250, MPFR_RNDN); mpfr_sub(r4619252, r4619244, r4619251, MPFR_RNDN); mpfr_div(r4619253, r4619252, r4619219, MPFR_RNDN); ; mpfr_set_si(r4619255, mpfr_cmp(r4619206, r4619254) <= 0, MPFR_RNDN); mpfr_pow(r4619256, r4619206, r4619226, MPFR_RNDN); mpfr_mul(r4619257, r4619209, r4619209, MPFR_RNDN); mpfr_mul(r4619257, r4619257, r4619209, MPFR_RNDN); mpfr_mul(r4619258, r4619256, r4619257, MPFR_RNDN); mpfr_add(r4619259, r4619206, r4619258, MPFR_RNDN); mpfr_mul(r4619260, r4619209, r4619209, MPFR_RNDN); mpfr_mul(r4619261, r4619206, r4619206, MPFR_RNDN); mpfr_mul(r4619261, r4619261, r4619206, MPFR_RNDN); mpfr_mul(r4619262, r4619260, r4619261, MPFR_RNDN); mpfr_add(r4619263, r4619259, r4619262, MPFR_RNDN); ; mpfr_set_si(r4619265, mpfr_cmp(r4619206, r4619264) <= 0, MPFR_RNDN); if (mpfr_get_si(r4619265, MPFR_RNDN)) { mpfr_set(r4619266, r4619253, MPFR_RNDN); } else { mpfr_set(r4619266, r4619220, MPFR_RNDN); }; if (mpfr_get_si(r4619255, MPFR_RNDN)) { mpfr_set(r4619267, r4619263, MPFR_RNDN); } else { mpfr_set(r4619267, r4619266, MPFR_RNDN); }; if (mpfr_get_si(r4619222, MPFR_RNDN)) { mpfr_set(r4619268, r4619253, MPFR_RNDN); } else { mpfr_set(r4619268, r4619267, MPFR_RNDN); }; if (mpfr_get_si(r4619208, MPFR_RNDN)) { mpfr_set(r4619269, r4619220, MPFR_RNDN); } else { mpfr_set(r4619269, r4619268, MPFR_RNDN); }; return mpfr_get_d(r4619269, MPFR_RNDN); } static mpfr_t r4619270, r4619271, r4619272, r4619273, r4619274, r4619275, r4619276, r4619277, r4619278, r4619279, r4619280, r4619281, r4619282, r4619283, r4619284, r4619285, r4619286, r4619287, r4619288, r4619289, r4619290, r4619291, r4619292, r4619293, r4619294, r4619295, r4619296, r4619297, r4619298, r4619299, r4619300, r4619301, r4619302, r4619303, r4619304, r4619305, r4619306, r4619307, r4619308, r4619309, r4619310, r4619311, r4619312, r4619313, r4619314, r4619315, r4619316, r4619317, r4619318, r4619319, r4619320, r4619321, r4619322, r4619323, r4619324, r4619325, r4619326, r4619327, r4619328, r4619329, r4619330, r4619331, r4619332, r4619333; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r4619270); mpfr_init_set_str(r4619271, "-1.1324884610140404e-37", 10, MPFR_RNDN); mpfr_init(r4619272); mpfr_init(r4619273); mpfr_init(r4619274); mpfr_init(r4619275); mpfr_init(r4619276); mpfr_init(r4619277); mpfr_init(r4619278); mpfr_init(r4619279); mpfr_init(r4619280); mpfr_init(r4619281); mpfr_init(r4619282); mpfr_init(r4619283); mpfr_init(r4619284); mpfr_init_set_str(r4619285, "-8.343087553984244e-186", 10, MPFR_RNDN); mpfr_init(r4619286); mpfr_init(r4619287); mpfr_init_set_str(r4619288, "3", 10, MPFR_RNDN); mpfr_init(r4619289); mpfr_init_set_str(r4619290, "4", 10, MPFR_RNDN); mpfr_init(r4619291); mpfr_init(r4619292); mpfr_init(r4619293); mpfr_init(r4619294); mpfr_init(r4619295); mpfr_init(r4619296); mpfr_init(r4619297); mpfr_init_set_str(r4619298, "1/3", 10, MPFR_RNDN); mpfr_init(r4619299); mpfr_init(r4619300); mpfr_init_set_str(r4619301, "4/3", 10, MPFR_RNDN); mpfr_init(r4619302); mpfr_init(r4619303); mpfr_init(r4619304); mpfr_init(r4619305); mpfr_init(r4619306); mpfr_init(r4619307); mpfr_init(r4619308); mpfr_init(r4619309); mpfr_init(r4619310); mpfr_init(r4619311); mpfr_init(r4619312); mpfr_init(r4619313); mpfr_init(r4619314); mpfr_init(r4619315); mpfr_init(r4619316); mpfr_init(r4619317); mpfr_init_set_str(r4619318, "-1.4036887096279737e-262", 10, MPFR_RNDN); mpfr_init(r4619319); mpfr_init(r4619320); mpfr_init(r4619321); mpfr_init(r4619322); mpfr_init(r4619323); mpfr_init(r4619324); mpfr_init(r4619325); mpfr_init(r4619326); mpfr_init(r4619327); mpfr_init_set_str(r4619328, "4.0170806835946265e-42", 10, MPFR_RNDN); mpfr_init(r4619329); mpfr_init(r4619330); mpfr_init(r4619331); mpfr_init(r4619332); mpfr_init(r4619333); } double f_dm(double x, double eps) { mpfr_set_d(r4619270, eps, MPFR_RNDN); ; mpfr_set_si(r4619272, mpfr_cmp(r4619270, r4619271) <= 0, MPFR_RNDN); mpfr_set_d(r4619273, x, MPFR_RNDN); mpfr_cos(r4619274, r4619273, MPFR_RNDN); mpfr_add(r4619275, r4619273, r4619270, MPFR_RNDN); mpfr_cot(r4619276, r4619275, MPFR_RNDN); mpfr_cbrt(r4619277, r4619276, MPFR_RNDN); mpfr_sin(r4619278, r4619273, MPFR_RNDN); mpfr_cbrt(r4619279, r4619278, MPFR_RNDN); mpfr_mul(r4619280, r4619277, r4619279, MPFR_RNDN); mpfr_mul(r4619281, r4619280, r4619280, MPFR_RNDN); mpfr_mul(r4619281, r4619281, r4619280, MPFR_RNDN); mpfr_sub(r4619282, r4619274, r4619281, MPFR_RNDN); mpfr_mul(r4619283, r4619276, r4619274, MPFR_RNDN); mpfr_div(r4619284, r4619282, r4619283, MPFR_RNDN); ; mpfr_set_si(r4619286, mpfr_cmp(r4619270, r4619285) <= 0, MPFR_RNDN); mpfr_mul(r4619287, r4619278, r4619270, MPFR_RNDN); ; mpfr_pow(r4619289, r4619270, r4619288, MPFR_RNDN); ; mpfr_pow(r4619291, r4619274, r4619290, MPFR_RNDN); mpfr_mul(r4619292, r4619289, r4619291, MPFR_RNDN); mpfr_pow(r4619293, r4619278, r4619288, MPFR_RNDN); mpfr_div(r4619294, r4619292, r4619293, MPFR_RNDN); mpfr_sqr(r4619295, r4619274, MPFR_RNDN); mpfr_mul(r4619296, r4619270, r4619295, MPFR_RNDN); mpfr_div(r4619297, r4619296, r4619278, MPFR_RNDN); ; mpfr_mul(r4619299, r4619278, r4619289, MPFR_RNDN); mpfr_mul(r4619300, r4619298, r4619299, MPFR_RNDN); ; mpfr_mul(r4619302, r4619289, r4619295, MPFR_RNDN); mpfr_div(r4619303, r4619302, r4619278, MPFR_RNDN); mpfr_mul(r4619304, r4619301, r4619303, MPFR_RNDN); mpfr_add(r4619305, r4619300, r4619304, MPFR_RNDN); mpfr_add(r4619306, r4619297, r4619305, MPFR_RNDN); mpfr_add(r4619307, r4619294, r4619306, MPFR_RNDN); mpfr_add(r4619308, r4619287, r4619307, MPFR_RNDN); mpfr_sqr(r4619309, r4619270, MPFR_RNDN); mpfr_pow(r4619310, r4619274, r4619288, MPFR_RNDN); mpfr_mul(r4619311, r4619309, r4619310, MPFR_RNDN); mpfr_sqr(r4619312, r4619278, MPFR_RNDN); mpfr_div(r4619313, r4619311, r4619312, MPFR_RNDN); mpfr_mul(r4619314, r4619309, r4619274, MPFR_RNDN); mpfr_add(r4619315, r4619313, r4619314, MPFR_RNDN); mpfr_sub(r4619316, r4619308, r4619315, MPFR_RNDN); mpfr_div(r4619317, r4619316, r4619283, MPFR_RNDN); ; mpfr_set_si(r4619319, mpfr_cmp(r4619270, r4619318) <= 0, MPFR_RNDN); mpfr_pow(r4619320, r4619270, r4619290, MPFR_RNDN); mpfr_mul(r4619321, r4619273, r4619273, MPFR_RNDN); mpfr_mul(r4619321, r4619321, r4619273, MPFR_RNDN); mpfr_mul(r4619322, r4619320, r4619321, MPFR_RNDN); mpfr_add(r4619323, r4619270, r4619322, MPFR_RNDN); mpfr_mul(r4619324, r4619273, r4619273, MPFR_RNDN); mpfr_mul(r4619325, r4619270, r4619270, MPFR_RNDN); mpfr_mul(r4619325, r4619325, r4619270, MPFR_RNDN); mpfr_mul(r4619326, r4619324, r4619325, MPFR_RNDN); mpfr_add(r4619327, r4619323, r4619326, MPFR_RNDN); ; mpfr_set_si(r4619329, mpfr_cmp(r4619270, r4619328) <= 0, MPFR_RNDN); if (mpfr_get_si(r4619329, MPFR_RNDN)) { mpfr_set(r4619330, r4619317, MPFR_RNDN); } else { mpfr_set(r4619330, r4619284, MPFR_RNDN); }; if (mpfr_get_si(r4619319, MPFR_RNDN)) { mpfr_set(r4619331, r4619327, MPFR_RNDN); } else { mpfr_set(r4619331, r4619330, MPFR_RNDN); }; if (mpfr_get_si(r4619286, MPFR_RNDN)) { mpfr_set(r4619332, r4619317, MPFR_RNDN); } else { mpfr_set(r4619332, r4619331, MPFR_RNDN); }; if (mpfr_get_si(r4619272, MPFR_RNDN)) { mpfr_set(r4619333, r4619284, MPFR_RNDN); } else { mpfr_set(r4619333, r4619332, MPFR_RNDN); }; return mpfr_get_d(r4619333, MPFR_RNDN); }