#include #include #include #include #include char *name = "NMSE example 3.3"; double f_if(float x, float eps) { float r17460 = x; float r17461 = eps; float r17462 = r17460 + r17461; float r17463 = sin(r17462); float r17464 = sin(r17460); float r17465 = r17463 - r17464; return r17465; } double f_id(double x, double eps) { double r17466 = x; double r17467 = eps; double r17468 = r17466 + r17467; double r17469 = sin(r17468); double r17470 = sin(r17466); double r17471 = r17469 - r17470; return r17471; } double f_of(float x, float eps) { float r17472 = eps; float r17473 = -0.0036671459674835205f; bool r17474 = r17472 <= r17473; float r17475 = x; float r17476 = sin(r17475); float r17477 = cos(r17472); float r17478 = r17476 * r17477; float r17479 = cos(r17475); float r17480 = sin(r17472); float r17481 = r17479 * r17480; float r17482 = 1.0f; float r17483 = pow(r17481, r17482); float r17484 = r17483 - r17476; float r17485 = r17478 + r17484; float r17486 = 3.4589021424835664e-08f; bool r17487 = r17472 <= r17486; float r17488 = 0.1666666716337204f; float r17489 = r17472 * (r17472 * r17472); float r17490 = r17488 * r17489; float r17491 = r17472 - r17490; float r17492 = r17491 * r17479; float r17493 = 0.5f; float r17494 = r17476 * r17493; float r17495 = r17472 * r17472; float r17496 = r17494 * r17495; float r17497 = r17492 - r17496; float r17498 = r17487 ? r17497 : r17485; float r17499 = r17474 ? r17485 : r17498; return r17499; } double f_od(double x, double eps) { double r17500 = eps; double r17501 = -0.0036671459674835205; bool r17502 = r17500 <= r17501; double r17503 = x; double r17504 = sin(r17503); double r17505 = cos(r17500); double r17506 = r17504 * r17505; double r17507 = cos(r17503); double r17508 = sin(r17500); double r17509 = r17507 * r17508; double r17510 = 1.0; double r17511 = pow(r17509, r17510); double r17512 = r17511 - r17504; double r17513 = r17506 + r17512; double r17514 = 3.4589021424835664e-08; bool r17515 = r17500 <= r17514; double r17516 = 0.1666666716337204; double r17517 = r17500 * (r17500 * r17500); double r17518 = r17516 * r17517; double r17519 = r17500 - r17518; double r17520 = r17519 * r17507; double r17521 = 0.5; double r17522 = r17504 * r17521; double r17523 = r17500 * r17500; double r17524 = r17522 * r17523; double r17525 = r17520 - r17524; double r17526 = r17515 ? r17525 : r17513; double r17527 = r17502 ? r17513 : r17526; return r17527; } 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 r17528, r17529, r17530, r17531, r17532, r17533; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17528); mpfr_init(r17529); mpfr_init(r17530); mpfr_init(r17531); mpfr_init(r17532); mpfr_init(r17533); } double f_im(double x, double eps) { mpfr_set_d(r17528, x, MPFR_RNDN); mpfr_set_d(r17529, eps, MPFR_RNDN); mpfr_add(r17530, r17528, r17529, MPFR_RNDN); mpfr_sin(r17531, r17530, MPFR_RNDN); mpfr_sin(r17532, r17528, MPFR_RNDN); mpfr_sub(r17533, r17531, r17532, MPFR_RNDN); return mpfr_get_d(r17533, MPFR_RNDN); } static mpfr_t r17534, r17535, r17536, r17537, r17538, r17539, r17540, r17541, r17542, r17543, r17544, r17545, r17546, r17547, r17548, r17549, r17550, r17551, r17552, r17553, r17554, r17555, r17556, r17557, r17558, r17559, r17560, r17561; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17534); mpfr_init_set_str(r17535, "-0.003667146f0", 10, MPFR_RNDN); mpfr_init(r17536); mpfr_init(r17537); mpfr_init(r17538); mpfr_init(r17539); mpfr_init(r17540); mpfr_init(r17541); mpfr_init(r17542); mpfr_init(r17543); mpfr_init_set_str(r17544, "1", 10, MPFR_RNDN); mpfr_init(r17545); mpfr_init(r17546); mpfr_init(r17547); mpfr_init_set_str(r17548, "3.458902f-08", 10, MPFR_RNDN); mpfr_init(r17549); mpfr_init_set_str(r17550, "1/6", 10, MPFR_RNDN); mpfr_init(r17551); mpfr_init(r17552); mpfr_init(r17553); mpfr_init(r17554); mpfr_init_set_str(r17555, "1/2", 10, MPFR_RNDN); mpfr_init(r17556); mpfr_init(r17557); mpfr_init(r17558); mpfr_init(r17559); mpfr_init(r17560); mpfr_init(r17561); } double f_fm(double x, double eps) { mpfr_set_d(r17534, eps, MPFR_RNDN); ; mpfr_set_si(r17536, mpfr_cmp(r17534, r17535) <= 0, MPFR_RNDN); mpfr_set_d(r17537, x, MPFR_RNDN); mpfr_sin(r17538, r17537, MPFR_RNDN); mpfr_cos(r17539, r17534, MPFR_RNDN); mpfr_mul(r17540, r17538, r17539, MPFR_RNDN); mpfr_cos(r17541, r17537, MPFR_RNDN); mpfr_sin(r17542, r17534, MPFR_RNDN); mpfr_mul(r17543, r17541, r17542, MPFR_RNDN); ; mpfr_pow(r17545, r17543, r17544, MPFR_RNDN); mpfr_sub(r17546, r17545, r17538, MPFR_RNDN); mpfr_add(r17547, r17540, r17546, MPFR_RNDN); ; mpfr_set_si(r17549, mpfr_cmp(r17534, r17548) <= 0, MPFR_RNDN); ; mpfr_mul(r17551, r17534, r17534, MPFR_RNDN); mpfr_mul(r17551, r17551, r17534, MPFR_RNDN); mpfr_mul(r17552, r17550, r17551, MPFR_RNDN); mpfr_sub(r17553, r17534, r17552, MPFR_RNDN); mpfr_mul(r17554, r17553, r17541, MPFR_RNDN); ; mpfr_mul(r17556, r17538, r17555, MPFR_RNDN); mpfr_sqr(r17557, r17534, MPFR_RNDN); mpfr_mul(r17558, r17556, r17557, MPFR_RNDN); mpfr_sub(r17559, r17554, r17558, MPFR_RNDN); if (mpfr_get_si(r17549, MPFR_RNDN)) { mpfr_set(r17560, r17559, MPFR_RNDN); } else { mpfr_set(r17560, r17547, MPFR_RNDN); }; if (mpfr_get_si(r17536, MPFR_RNDN)) { mpfr_set(r17561, r17547, MPFR_RNDN); } else { mpfr_set(r17561, r17560, MPFR_RNDN); }; return mpfr_get_d(r17561, MPFR_RNDN); } static mpfr_t r17562, r17563, r17564, r17565, r17566, r17567, r17568, r17569, r17570, r17571, r17572, r17573, r17574, r17575, r17576, r17577, r17578, r17579, r17580, r17581, r17582, r17583, r17584, r17585, r17586, r17587, r17588, r17589; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17562); mpfr_init_set_str(r17563, "-0.003667146f0", 10, MPFR_RNDN); mpfr_init(r17564); mpfr_init(r17565); mpfr_init(r17566); mpfr_init(r17567); mpfr_init(r17568); mpfr_init(r17569); mpfr_init(r17570); mpfr_init(r17571); mpfr_init_set_str(r17572, "1", 10, MPFR_RNDN); mpfr_init(r17573); mpfr_init(r17574); mpfr_init(r17575); mpfr_init_set_str(r17576, "3.458902f-08", 10, MPFR_RNDN); mpfr_init(r17577); mpfr_init_set_str(r17578, "1/6", 10, MPFR_RNDN); mpfr_init(r17579); mpfr_init(r17580); mpfr_init(r17581); mpfr_init(r17582); mpfr_init_set_str(r17583, "1/2", 10, MPFR_RNDN); mpfr_init(r17584); mpfr_init(r17585); mpfr_init(r17586); mpfr_init(r17587); mpfr_init(r17588); mpfr_init(r17589); } double f_dm(double x, double eps) { mpfr_set_d(r17562, eps, MPFR_RNDN); ; mpfr_set_si(r17564, mpfr_cmp(r17562, r17563) <= 0, MPFR_RNDN); mpfr_set_d(r17565, x, MPFR_RNDN); mpfr_sin(r17566, r17565, MPFR_RNDN); mpfr_cos(r17567, r17562, MPFR_RNDN); mpfr_mul(r17568, r17566, r17567, MPFR_RNDN); mpfr_cos(r17569, r17565, MPFR_RNDN); mpfr_sin(r17570, r17562, MPFR_RNDN); mpfr_mul(r17571, r17569, r17570, MPFR_RNDN); ; mpfr_pow(r17573, r17571, r17572, MPFR_RNDN); mpfr_sub(r17574, r17573, r17566, MPFR_RNDN); mpfr_add(r17575, r17568, r17574, MPFR_RNDN); ; mpfr_set_si(r17577, mpfr_cmp(r17562, r17576) <= 0, MPFR_RNDN); ; mpfr_mul(r17579, r17562, r17562, MPFR_RNDN); mpfr_mul(r17579, r17579, r17562, MPFR_RNDN); mpfr_mul(r17580, r17578, r17579, MPFR_RNDN); mpfr_sub(r17581, r17562, r17580, MPFR_RNDN); mpfr_mul(r17582, r17581, r17569, MPFR_RNDN); ; mpfr_mul(r17584, r17566, r17583, MPFR_RNDN); mpfr_sqr(r17585, r17562, MPFR_RNDN); mpfr_mul(r17586, r17584, r17585, MPFR_RNDN); mpfr_sub(r17587, r17582, r17586, MPFR_RNDN); if (mpfr_get_si(r17577, MPFR_RNDN)) { mpfr_set(r17588, r17587, MPFR_RNDN); } else { mpfr_set(r17588, r17575, MPFR_RNDN); }; if (mpfr_get_si(r17564, MPFR_RNDN)) { mpfr_set(r17589, r17575, MPFR_RNDN); } else { mpfr_set(r17589, r17588, MPFR_RNDN); }; return mpfr_get_d(r17589, MPFR_RNDN); }