#include #include #include #include #include char *name = "NMSE Section 6.1 mentioned, A"; double f_if(float x, float eps) { float r21145 = 1; float r21146 = eps; float r21147 = r21145 / r21146; float r21148 = r21145 + r21147; float r21149 = r21145 - r21146; float r21150 = x; float r21151 = r21149 * r21150; float r21152 = -r21151; float r21153 = exp(r21152); float r21154 = r21148 * r21153; float r21155 = r21147 - r21145; float r21156 = r21145 + r21146; float r21157 = r21156 * r21150; float r21158 = -r21157; float r21159 = exp(r21158); float r21160 = r21155 * r21159; float r21161 = r21154 - r21160; float r21162 = 2; float r21163 = r21161 / r21162; return r21163; } double f_id(double x, double eps) { double r21164 = 1; double r21165 = eps; double r21166 = r21164 / r21165; double r21167 = r21164 + r21166; double r21168 = r21164 - r21165; double r21169 = x; double r21170 = r21168 * r21169; double r21171 = -r21170; double r21172 = exp(r21171); double r21173 = r21167 * r21172; double r21174 = r21166 - r21164; double r21175 = r21164 + r21165; double r21176 = r21175 * r21169; double r21177 = -r21176; double r21178 = exp(r21177); double r21179 = r21174 * r21178; double r21180 = r21173 - r21179; double r21181 = 2; double r21182 = r21180 / r21181; return r21182; } double f_of(float x, float eps) { float r21183 = 2/3; float r21184 = x; float r21185 = 3; float r21186 = pow(r21184, r21185); float r21187 = r21183 / r21186; float r21188 = 27/4; float r21189 = r21188 / r21184; float r21190 = 9/2; float r21191 = r21189 - r21190; float r21192 = r21187 * r21191; float r21193 = r21184 / r21183; float r21194 = r21193 * r21186; float r21195 = r21194 / r21187; float r21196 = r21192 / r21195; float r21197 = pow(r21187, r21185); float r21198 = r21196 + r21197; float r21199 = cbrt(r21198); float r21200 = 2; float r21201 = r21199 / r21200; float r21202 = -2.0703859980849916e+22; bool r21203 = r21201 <= r21202; float r21204 = r21183 * r21186; float r21205 = r21200 + r21204; float r21206 = pow(r21184, r21200); float r21207 = r21205 - r21206; float r21208 = r21207 / r21200; float r21209 = 5.842364092482766e-05; bool r21210 = r21201 <= r21209; float r21211 = 1; float r21212 = eps; float r21213 = r21211 / r21212; float r21214 = r21211 + r21213; float r21215 = r21211 - r21212; float r21216 = r21215 * r21184; float r21217 = -r21216; float r21218 = exp(r21217); float r21219 = sqrt(r21218); float r21220 = r21214 * r21219; float r21221 = r21220 * r21219; float r21222 = r21213 - r21211; float r21223 = r21211 + r21212; float r21224 = r21223 * r21184; float r21225 = -r21224; float r21226 = exp(r21225); float r21227 = r21222 * r21226; float r21228 = r21221 - r21227; float r21229 = r21228 / r21200; float r21230 = r21210 ? r21229 : r21208; float r21231 = r21203 ? r21208 : r21230; return r21231; } double f_od(double x, double eps) { double r21232 = 2/3; double r21233 = x; double r21234 = 3; double r21235 = pow(r21233, r21234); double r21236 = r21232 / r21235; double r21237 = 27/4; double r21238 = r21237 / r21233; double r21239 = 9/2; double r21240 = r21238 - r21239; double r21241 = r21236 * r21240; double r21242 = r21233 / r21232; double r21243 = r21242 * r21235; double r21244 = r21243 / r21236; double r21245 = r21241 / r21244; double r21246 = pow(r21236, r21234); double r21247 = r21245 + r21246; double r21248 = cbrt(r21247); double r21249 = 2; double r21250 = r21248 / r21249; double r21251 = -2.0703859980849916e+22; bool r21252 = r21250 <= r21251; double r21253 = r21232 * r21235; double r21254 = r21249 + r21253; double r21255 = pow(r21233, r21249); double r21256 = r21254 - r21255; double r21257 = r21256 / r21249; double r21258 = 5.842364092482766e-05; bool r21259 = r21250 <= r21258; double r21260 = 1; double r21261 = eps; double r21262 = r21260 / r21261; double r21263 = r21260 + r21262; double r21264 = r21260 - r21261; double r21265 = r21264 * r21233; double r21266 = -r21265; double r21267 = exp(r21266); double r21268 = sqrt(r21267); double r21269 = r21263 * r21268; double r21270 = r21269 * r21268; double r21271 = r21262 - r21260; double r21272 = r21260 + r21261; double r21273 = r21272 * r21233; double r21274 = -r21273; double r21275 = exp(r21274); double r21276 = r21271 * r21275; double r21277 = r21270 - r21276; double r21278 = r21277 / r21249; double r21279 = r21259 ? r21278 : r21257; double r21280 = r21252 ? r21257 : r21279; return r21280; } 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 r21281, r21282, r21283, r21284, r21285, r21286, r21287, r21288, r21289, r21290, r21291, r21292, r21293, r21294, r21295, r21296, r21297, r21298, r21299; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r21281, "1", 10, MPFR_RNDN); mpfr_init(r21282); mpfr_init(r21283); mpfr_init(r21284); mpfr_init(r21285); mpfr_init(r21286); mpfr_init(r21287); mpfr_init(r21288); mpfr_init(r21289); mpfr_init(r21290); mpfr_init(r21291); mpfr_init(r21292); mpfr_init(r21293); mpfr_init(r21294); mpfr_init(r21295); mpfr_init(r21296); mpfr_init(r21297); mpfr_init_set_str(r21298, "2", 10, MPFR_RNDN); mpfr_init(r21299); } double f_im(double x, double eps) { ; mpfr_set_d(r21282, eps, MPFR_RNDN); mpfr_div(r21283, r21281, r21282, MPFR_RNDN); mpfr_add(r21284, r21281, r21283, MPFR_RNDN); mpfr_sub(r21285, r21281, r21282, MPFR_RNDN); mpfr_set_d(r21286, x, MPFR_RNDN); mpfr_mul(r21287, r21285, r21286, MPFR_RNDN); mpfr_neg(r21288, r21287, MPFR_RNDN); mpfr_exp(r21289, r21288, MPFR_RNDN); mpfr_mul(r21290, r21284, r21289, MPFR_RNDN); mpfr_sub(r21291, r21283, r21281, MPFR_RNDN); mpfr_add(r21292, r21281, r21282, MPFR_RNDN); mpfr_mul(r21293, r21292, r21286, MPFR_RNDN); mpfr_neg(r21294, r21293, MPFR_RNDN); mpfr_exp(r21295, r21294, MPFR_RNDN); mpfr_mul(r21296, r21291, r21295, MPFR_RNDN); mpfr_sub(r21297, r21290, r21296, MPFR_RNDN); ; mpfr_div(r21299, r21297, r21298, MPFR_RNDN); return mpfr_get_d(r21299, MPFR_RNDN); } static mpfr_t r21300, r21301, r21302, r21303, r21304, r21305, r21306, r21307, r21308, r21309, r21310, r21311, r21312, r21313, r21314, r21315, r21316, r21317, r21318, r21319, r21320, r21321, r21322, r21323, r21324, r21325, r21326, r21327, r21328, r21329, r21330, r21331, r21332, r21333, r21334, r21335, r21336, r21337, r21338, r21339, r21340, r21341, r21342, r21343, r21344, r21345, r21346, r21347, r21348; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r21300, "2/3", 10, MPFR_RNDN); mpfr_init(r21301); mpfr_init_set_str(r21302, "3", 10, MPFR_RNDN); mpfr_init(r21303); mpfr_init(r21304); mpfr_init_set_str(r21305, "27/4", 10, MPFR_RNDN); mpfr_init(r21306); mpfr_init_set_str(r21307, "9/2", 10, MPFR_RNDN); mpfr_init(r21308); mpfr_init(r21309); mpfr_init(r21310); mpfr_init(r21311); mpfr_init(r21312); mpfr_init(r21313); mpfr_init(r21314); mpfr_init(r21315); mpfr_init(r21316); mpfr_init_set_str(r21317, "2", 10, MPFR_RNDN); mpfr_init(r21318); mpfr_init_set_str(r21319, "-2.0703859980849916e+22", 10, MPFR_RNDN); mpfr_init(r21320); mpfr_init(r21321); mpfr_init(r21322); mpfr_init(r21323); mpfr_init(r21324); mpfr_init(r21325); mpfr_init_set_str(r21326, "5.842364092482766e-05", 10, MPFR_RNDN); mpfr_init(r21327); mpfr_init_set_str(r21328, "1", 10, MPFR_RNDN); mpfr_init(r21329); mpfr_init(r21330); mpfr_init(r21331); mpfr_init(r21332); mpfr_init(r21333); mpfr_init(r21334); mpfr_init(r21335); mpfr_init(r21336); mpfr_init(r21337); mpfr_init(r21338); mpfr_init(r21339); mpfr_init(r21340); mpfr_init(r21341); mpfr_init(r21342); mpfr_init(r21343); mpfr_init(r21344); mpfr_init(r21345); mpfr_init(r21346); mpfr_init(r21347); mpfr_init(r21348); } double f_fm(double x, double eps) { ; mpfr_set_d(r21301, x, MPFR_RNDN); ; mpfr_pow(r21303, r21301, r21302, MPFR_RNDN); mpfr_div(r21304, r21300, r21303, MPFR_RNDN); ; mpfr_div(r21306, r21305, r21301, MPFR_RNDN); ; mpfr_sub(r21308, r21306, r21307, MPFR_RNDN); mpfr_mul(r21309, r21304, r21308, MPFR_RNDN); mpfr_div(r21310, r21301, r21300, MPFR_RNDN); mpfr_mul(r21311, r21310, r21303, MPFR_RNDN); mpfr_div(r21312, r21311, r21304, MPFR_RNDN); mpfr_div(r21313, r21309, r21312, MPFR_RNDN); mpfr_pow(r21314, r21304, r21302, MPFR_RNDN); mpfr_add(r21315, r21313, r21314, MPFR_RNDN); mpfr_cbrt(r21316, r21315, MPFR_RNDN); ; mpfr_div(r21318, r21316, r21317, MPFR_RNDN); ; mpfr_set_si(r21320, mpfr_cmp(r21318, r21319) <= 0, MPFR_RNDN); mpfr_mul(r21321, r21300, r21303, MPFR_RNDN); mpfr_add(r21322, r21317, r21321, MPFR_RNDN); mpfr_pow(r21323, r21301, r21317, MPFR_RNDN); mpfr_sub(r21324, r21322, r21323, MPFR_RNDN); mpfr_div(r21325, r21324, r21317, MPFR_RNDN); ; mpfr_set_si(r21327, mpfr_cmp(r21318, r21326) <= 0, MPFR_RNDN); ; mpfr_set_d(r21329, eps, MPFR_RNDN); mpfr_div(r21330, r21328, r21329, MPFR_RNDN); mpfr_add(r21331, r21328, r21330, MPFR_RNDN); mpfr_sub(r21332, r21328, r21329, MPFR_RNDN); mpfr_mul(r21333, r21332, r21301, MPFR_RNDN); mpfr_neg(r21334, r21333, MPFR_RNDN); mpfr_exp(r21335, r21334, MPFR_RNDN); mpfr_sqrt(r21336, r21335, MPFR_RNDN); mpfr_mul(r21337, r21331, r21336, MPFR_RNDN); mpfr_mul(r21338, r21337, r21336, MPFR_RNDN); mpfr_sub(r21339, r21330, r21328, MPFR_RNDN); mpfr_add(r21340, r21328, r21329, MPFR_RNDN); mpfr_mul(r21341, r21340, r21301, MPFR_RNDN); mpfr_neg(r21342, r21341, MPFR_RNDN); mpfr_exp(r21343, r21342, MPFR_RNDN); mpfr_mul(r21344, r21339, r21343, MPFR_RNDN); mpfr_sub(r21345, r21338, r21344, MPFR_RNDN); mpfr_div(r21346, r21345, r21317, MPFR_RNDN); if (mpfr_get_si(r21327, MPFR_RNDN)) { mpfr_set(r21347, r21346, MPFR_RNDN); } else { mpfr_set(r21347, r21325, MPFR_RNDN); }; if (mpfr_get_si(r21320, MPFR_RNDN)) { mpfr_set(r21348, r21325, MPFR_RNDN); } else { mpfr_set(r21348, r21347, MPFR_RNDN); }; return mpfr_get_d(r21348, MPFR_RNDN); } static mpfr_t r21349, r21350, r21351, r21352, r21353, r21354, r21355, r21356, r21357, r21358, r21359, r21360, r21361, r21362, r21363, r21364, r21365, r21366, r21367, r21368, r21369, r21370, r21371, r21372, r21373, r21374, r21375, r21376, r21377, r21378, r21379, r21380, r21381, r21382, r21383, r21384, r21385, r21386, r21387, r21388, r21389, r21390, r21391, r21392, r21393, r21394, r21395, r21396, r21397; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r21349, "2/3", 10, MPFR_RNDN); mpfr_init(r21350); mpfr_init_set_str(r21351, "3", 10, MPFR_RNDN); mpfr_init(r21352); mpfr_init(r21353); mpfr_init_set_str(r21354, "27/4", 10, MPFR_RNDN); mpfr_init(r21355); mpfr_init_set_str(r21356, "9/2", 10, MPFR_RNDN); mpfr_init(r21357); mpfr_init(r21358); mpfr_init(r21359); mpfr_init(r21360); mpfr_init(r21361); mpfr_init(r21362); mpfr_init(r21363); mpfr_init(r21364); mpfr_init(r21365); mpfr_init_set_str(r21366, "2", 10, MPFR_RNDN); mpfr_init(r21367); mpfr_init_set_str(r21368, "-2.0703859980849916e+22", 10, MPFR_RNDN); mpfr_init(r21369); mpfr_init(r21370); mpfr_init(r21371); mpfr_init(r21372); mpfr_init(r21373); mpfr_init(r21374); mpfr_init_set_str(r21375, "5.842364092482766e-05", 10, MPFR_RNDN); mpfr_init(r21376); mpfr_init_set_str(r21377, "1", 10, MPFR_RNDN); mpfr_init(r21378); mpfr_init(r21379); mpfr_init(r21380); mpfr_init(r21381); mpfr_init(r21382); mpfr_init(r21383); mpfr_init(r21384); mpfr_init(r21385); mpfr_init(r21386); mpfr_init(r21387); mpfr_init(r21388); mpfr_init(r21389); mpfr_init(r21390); mpfr_init(r21391); mpfr_init(r21392); mpfr_init(r21393); mpfr_init(r21394); mpfr_init(r21395); mpfr_init(r21396); mpfr_init(r21397); } double f_dm(double x, double eps) { ; mpfr_set_d(r21350, x, MPFR_RNDN); ; mpfr_pow(r21352, r21350, r21351, MPFR_RNDN); mpfr_div(r21353, r21349, r21352, MPFR_RNDN); ; mpfr_div(r21355, r21354, r21350, MPFR_RNDN); ; mpfr_sub(r21357, r21355, r21356, MPFR_RNDN); mpfr_mul(r21358, r21353, r21357, MPFR_RNDN); mpfr_div(r21359, r21350, r21349, MPFR_RNDN); mpfr_mul(r21360, r21359, r21352, MPFR_RNDN); mpfr_div(r21361, r21360, r21353, MPFR_RNDN); mpfr_div(r21362, r21358, r21361, MPFR_RNDN); mpfr_pow(r21363, r21353, r21351, MPFR_RNDN); mpfr_add(r21364, r21362, r21363, MPFR_RNDN); mpfr_cbrt(r21365, r21364, MPFR_RNDN); ; mpfr_div(r21367, r21365, r21366, MPFR_RNDN); ; mpfr_set_si(r21369, mpfr_cmp(r21367, r21368) <= 0, MPFR_RNDN); mpfr_mul(r21370, r21349, r21352, MPFR_RNDN); mpfr_add(r21371, r21366, r21370, MPFR_RNDN); mpfr_pow(r21372, r21350, r21366, MPFR_RNDN); mpfr_sub(r21373, r21371, r21372, MPFR_RNDN); mpfr_div(r21374, r21373, r21366, MPFR_RNDN); ; mpfr_set_si(r21376, mpfr_cmp(r21367, r21375) <= 0, MPFR_RNDN); ; mpfr_set_d(r21378, eps, MPFR_RNDN); mpfr_div(r21379, r21377, r21378, MPFR_RNDN); mpfr_add(r21380, r21377, r21379, MPFR_RNDN); mpfr_sub(r21381, r21377, r21378, MPFR_RNDN); mpfr_mul(r21382, r21381, r21350, MPFR_RNDN); mpfr_neg(r21383, r21382, MPFR_RNDN); mpfr_exp(r21384, r21383, MPFR_RNDN); mpfr_sqrt(r21385, r21384, MPFR_RNDN); mpfr_mul(r21386, r21380, r21385, MPFR_RNDN); mpfr_mul(r21387, r21386, r21385, MPFR_RNDN); mpfr_sub(r21388, r21379, r21377, MPFR_RNDN); mpfr_add(r21389, r21377, r21378, MPFR_RNDN); mpfr_mul(r21390, r21389, r21350, MPFR_RNDN); mpfr_neg(r21391, r21390, MPFR_RNDN); mpfr_exp(r21392, r21391, MPFR_RNDN); mpfr_mul(r21393, r21388, r21392, MPFR_RNDN); mpfr_sub(r21394, r21387, r21393, MPFR_RNDN); mpfr_div(r21395, r21394, r21366, MPFR_RNDN); if (mpfr_get_si(r21376, MPFR_RNDN)) { mpfr_set(r21396, r21395, MPFR_RNDN); } else { mpfr_set(r21396, r21374, MPFR_RNDN); }; if (mpfr_get_si(r21369, MPFR_RNDN)) { mpfr_set(r21397, r21374, MPFR_RNDN); } else { mpfr_set(r21397, r21396, MPFR_RNDN); }; return mpfr_get_d(r21397, MPFR_RNDN); }