#include #include #include #include #include char *name = "NMSE example 3.3"; double f_if(float x, float eps) { float r17324 = x; float r17325 = eps; float r17326 = r17324 + r17325; float r17327 = sin(r17326); float r17328 = sin(r17324); float r17329 = r17327 - r17328; return r17329; } double f_id(double x, double eps) { double r17330 = x; double r17331 = eps; double r17332 = r17330 + r17331; double r17333 = sin(r17332); double r17334 = sin(r17330); double r17335 = r17333 - r17334; return r17335; } double f_of(float x, float eps) { float r17336 = eps; float r17337 = -1.5876557880878384e-13f; bool r17338 = r17336 <= r17337; float r17339 = x; float r17340 = sin(r17339); float r17341 = cos(r17336); float r17342 = r17340 * r17341; float r17343 = r17342 * r17342; float r17344 = cos(r17339); float r17345 = sin(r17336); float r17346 = r17344 * r17345; float r17347 = r17346 - r17340; float r17348 = r17347 * r17347; float r17349 = r17343 - r17348; float r17350 = r17342 - r17347; float r17351 = r17349 / r17350; float r17352 = 1.9313832044564974e-20f; bool r17353 = r17336 <= r17352; float r17354 = 0.5f; float r17355 = r17336 + r17339; float r17356 = r17339 * r17336; float r17357 = r17355 * r17356; float r17358 = r17354 * r17357; float r17359 = r17336 - r17358; float r17360 = r17353 ? r17359 : r17351; float r17361 = r17338 ? r17351 : r17360; return r17361; } double f_od(double x, double eps) { double r17362 = eps; double r17363 = -1.5876557880878384e-13; bool r17364 = r17362 <= r17363; double r17365 = x; double r17366 = sin(r17365); double r17367 = cos(r17362); double r17368 = r17366 * r17367; double r17369 = r17368 * r17368; double r17370 = cos(r17365); double r17371 = sin(r17362); double r17372 = r17370 * r17371; double r17373 = r17372 - r17366; double r17374 = r17373 * r17373; double r17375 = r17369 - r17374; double r17376 = r17368 - r17373; double r17377 = r17375 / r17376; double r17378 = 1.9313832044564974e-20; bool r17379 = r17362 <= r17378; double r17380 = 0.5; double r17381 = r17362 + r17365; double r17382 = r17365 * r17362; double r17383 = r17381 * r17382; double r17384 = r17380 * r17383; double r17385 = r17362 - r17384; double r17386 = r17379 ? r17385 : r17377; double r17387 = r17364 ? r17377 : r17386; return r17387; } 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 r17388, r17389, r17390, r17391, r17392, r17393; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17388); mpfr_init(r17389); mpfr_init(r17390); mpfr_init(r17391); mpfr_init(r17392); mpfr_init(r17393); } double f_im(double x, double eps) { mpfr_set_d(r17388, x, MPFR_RNDN); mpfr_set_d(r17389, eps, MPFR_RNDN); mpfr_add(r17390, r17388, r17389, MPFR_RNDN); mpfr_sin(r17391, r17390, MPFR_RNDN); mpfr_sin(r17392, r17388, MPFR_RNDN); mpfr_sub(r17393, r17391, r17392, MPFR_RNDN); return mpfr_get_d(r17393, MPFR_RNDN); } static mpfr_t r17394, r17395, r17396, r17397, r17398, r17399, r17400, r17401, r17402, r17403, r17404, r17405, r17406, r17407, r17408, r17409, r17410, r17411, r17412, r17413, r17414, r17415, r17416, r17417, r17418, r17419; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17394); mpfr_init_set_str(r17395, "-1.5876558f-13", 10, MPFR_RNDN); mpfr_init(r17396); mpfr_init(r17397); mpfr_init(r17398); mpfr_init(r17399); mpfr_init(r17400); mpfr_init(r17401); mpfr_init(r17402); mpfr_init(r17403); mpfr_init(r17404); mpfr_init(r17405); mpfr_init(r17406); mpfr_init(r17407); mpfr_init(r17408); mpfr_init(r17409); mpfr_init_set_str(r17410, "1.9313832f-20", 10, MPFR_RNDN); mpfr_init(r17411); mpfr_init_set_str(r17412, "1/2", 10, MPFR_RNDN); mpfr_init(r17413); mpfr_init(r17414); mpfr_init(r17415); mpfr_init(r17416); mpfr_init(r17417); mpfr_init(r17418); mpfr_init(r17419); } double f_fm(double x, double eps) { mpfr_set_d(r17394, eps, MPFR_RNDN); ; mpfr_set_si(r17396, mpfr_cmp(r17394, r17395) <= 0, MPFR_RNDN); mpfr_set_d(r17397, x, MPFR_RNDN); mpfr_sin(r17398, r17397, MPFR_RNDN); mpfr_cos(r17399, r17394, MPFR_RNDN); mpfr_mul(r17400, r17398, r17399, MPFR_RNDN); mpfr_sqr(r17401, r17400, MPFR_RNDN); mpfr_cos(r17402, r17397, MPFR_RNDN); mpfr_sin(r17403, r17394, MPFR_RNDN); mpfr_mul(r17404, r17402, r17403, MPFR_RNDN); mpfr_sub(r17405, r17404, r17398, MPFR_RNDN); mpfr_sqr(r17406, r17405, MPFR_RNDN); mpfr_sub(r17407, r17401, r17406, MPFR_RNDN); mpfr_sub(r17408, r17400, r17405, MPFR_RNDN); mpfr_div(r17409, r17407, r17408, MPFR_RNDN); ; mpfr_set_si(r17411, mpfr_cmp(r17394, r17410) <= 0, MPFR_RNDN); ; mpfr_add(r17413, r17394, r17397, MPFR_RNDN); mpfr_mul(r17414, r17397, r17394, MPFR_RNDN); mpfr_mul(r17415, r17413, r17414, MPFR_RNDN); mpfr_mul(r17416, r17412, r17415, MPFR_RNDN); mpfr_sub(r17417, r17394, r17416, MPFR_RNDN); if (mpfr_get_si(r17411, MPFR_RNDN)) { mpfr_set(r17418, r17417, MPFR_RNDN); } else { mpfr_set(r17418, r17409, MPFR_RNDN); }; if (mpfr_get_si(r17396, MPFR_RNDN)) { mpfr_set(r17419, r17409, MPFR_RNDN); } else { mpfr_set(r17419, r17418, MPFR_RNDN); }; return mpfr_get_d(r17419, MPFR_RNDN); } static mpfr_t r17420, r17421, r17422, r17423, r17424, r17425, r17426, r17427, r17428, r17429, r17430, r17431, r17432, r17433, r17434, r17435, r17436, r17437, r17438, r17439, r17440, r17441, r17442, r17443, r17444, r17445; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17420); mpfr_init_set_str(r17421, "-1.5876558f-13", 10, MPFR_RNDN); mpfr_init(r17422); mpfr_init(r17423); mpfr_init(r17424); mpfr_init(r17425); mpfr_init(r17426); mpfr_init(r17427); mpfr_init(r17428); mpfr_init(r17429); mpfr_init(r17430); mpfr_init(r17431); mpfr_init(r17432); mpfr_init(r17433); mpfr_init(r17434); mpfr_init(r17435); mpfr_init_set_str(r17436, "1.9313832f-20", 10, MPFR_RNDN); mpfr_init(r17437); mpfr_init_set_str(r17438, "1/2", 10, MPFR_RNDN); mpfr_init(r17439); mpfr_init(r17440); mpfr_init(r17441); mpfr_init(r17442); mpfr_init(r17443); mpfr_init(r17444); mpfr_init(r17445); } double f_dm(double x, double eps) { mpfr_set_d(r17420, eps, MPFR_RNDN); ; mpfr_set_si(r17422, mpfr_cmp(r17420, r17421) <= 0, MPFR_RNDN); mpfr_set_d(r17423, x, MPFR_RNDN); mpfr_sin(r17424, r17423, MPFR_RNDN); mpfr_cos(r17425, r17420, MPFR_RNDN); mpfr_mul(r17426, r17424, r17425, MPFR_RNDN); mpfr_sqr(r17427, r17426, MPFR_RNDN); mpfr_cos(r17428, r17423, MPFR_RNDN); mpfr_sin(r17429, r17420, MPFR_RNDN); mpfr_mul(r17430, r17428, r17429, MPFR_RNDN); mpfr_sub(r17431, r17430, r17424, MPFR_RNDN); mpfr_sqr(r17432, r17431, MPFR_RNDN); mpfr_sub(r17433, r17427, r17432, MPFR_RNDN); mpfr_sub(r17434, r17426, r17431, MPFR_RNDN); mpfr_div(r17435, r17433, r17434, MPFR_RNDN); ; mpfr_set_si(r17437, mpfr_cmp(r17420, r17436) <= 0, MPFR_RNDN); ; mpfr_add(r17439, r17420, r17423, MPFR_RNDN); mpfr_mul(r17440, r17423, r17420, MPFR_RNDN); mpfr_mul(r17441, r17439, r17440, MPFR_RNDN); mpfr_mul(r17442, r17438, r17441, MPFR_RNDN); mpfr_sub(r17443, r17420, r17442, MPFR_RNDN); if (mpfr_get_si(r17437, MPFR_RNDN)) { mpfr_set(r17444, r17443, MPFR_RNDN); } else { mpfr_set(r17444, r17435, MPFR_RNDN); }; if (mpfr_get_si(r17422, MPFR_RNDN)) { mpfr_set(r17445, r17435, MPFR_RNDN); } else { mpfr_set(r17445, r17444, MPFR_RNDN); }; return mpfr_get_d(r17445, MPFR_RNDN); }