#include #include #include #include #include char *name = "NMSE example 3.4"; double f_if(float x) { float r17562 = 1.0f; float r17563 = x; float r17564 = cos(r17563); float r17565 = r17562 - r17564; float r17566 = sin(r17563); float r17567 = r17565 / r17566; return r17567; } double f_id(double x) { double r17568 = 1.0; double r17569 = x; double r17570 = cos(r17569); double r17571 = r17568 - r17570; double r17572 = sin(r17569); double r17573 = r17571 / r17572; return r17573; } double f_of(float x) { float r17574 = x; float r17575 = -8.61646249982272e-12f; bool r17576 = r17574 <= r17575; float r17577 = sin(r17574); float r17578 = r17577 * r17577; float r17579 = 1.0f; float r17580 = cos(r17574); float r17581 = r17579 + r17580; float r17582 = r17578 / r17581; float r17583 = r17582 / r17577; float r17584 = 8.584780819621077f; bool r17585 = r17574 <= r17584; float r17586 = 0.5f; float r17587 = r17574 * r17586; float r17588 = 0.004166666666666667f; float r17589 = 5.0f; float r17590 = pow(r17574, r17589); float r17591 = r17588 * r17590; float r17592 = r17574 * (r17574 * r17574); float r17593 = 0.041666666666666664f; float r17594 = r17592 * r17593; float r17595 = r17591 + r17594; float r17596 = r17587 + r17595; float r17597 = r17579 - r17580; float r17598 = r17579 / r17577; float r17599 = r17597 * r17598; float r17600 = r17585 ? r17596 : r17599; float r17601 = r17576 ? r17583 : r17600; return r17601; } double f_od(double x) { double r17602 = x; double r17603 = -8.61646249982272e-12; bool r17604 = r17602 <= r17603; double r17605 = sin(r17602); double r17606 = r17605 * r17605; double r17607 = 1.0; double r17608 = cos(r17602); double r17609 = r17607 + r17608; double r17610 = r17606 / r17609; double r17611 = r17610 / r17605; double r17612 = 8.584780819621077; bool r17613 = r17602 <= r17612; double r17614 = 0.5; double r17615 = r17602 * r17614; double r17616 = 0.004166666666666667; double r17617 = 5.0; double r17618 = pow(r17602, r17617); double r17619 = r17616 * r17618; double r17620 = r17602 * (r17602 * r17602); double r17621 = 0.041666666666666664; double r17622 = r17620 * r17621; double r17623 = r17619 + r17622; double r17624 = r17615 + r17623; double r17625 = r17607 - r17608; double r17626 = r17607 / r17605; double r17627 = r17625 * r17626; double r17628 = r17613 ? r17624 : r17627; double r17629 = r17604 ? r17611 : r17628; return r17629; } 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 r17630, r17631, r17632, r17633, r17634, r17635; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r17630, "1", 10, MPFR_RNDN); mpfr_init(r17631); mpfr_init(r17632); mpfr_init(r17633); mpfr_init(r17634); mpfr_init(r17635); } double f_im(double x) { ; mpfr_set_d(r17631, x, MPFR_RNDN); mpfr_cos(r17632, r17631, MPFR_RNDN); mpfr_sub(r17633, r17630, r17632, MPFR_RNDN); mpfr_sin(r17634, r17631, MPFR_RNDN); mpfr_div(r17635, r17633, r17634, MPFR_RNDN); return mpfr_get_d(r17635, MPFR_RNDN); } static mpfr_t r17636, r17637, r17638, r17639, r17640, r17641, r17642, r17643, r17644, r17645, r17646, r17647, r17648, r17649, r17650, r17651, r17652, r17653, r17654, r17655, r17656, r17657, r17658, r17659, r17660, r17661, r17662, r17663; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17636); mpfr_init_set_str(r17637, "-8.61646249982272e-12", 10, MPFR_RNDN); mpfr_init(r17638); mpfr_init(r17639); mpfr_init(r17640); mpfr_init_set_str(r17641, "1", 10, MPFR_RNDN); mpfr_init(r17642); mpfr_init(r17643); mpfr_init(r17644); mpfr_init(r17645); mpfr_init_set_str(r17646, "8.584780819621077", 10, MPFR_RNDN); mpfr_init(r17647); mpfr_init_set_str(r17648, "1/2", 10, MPFR_RNDN); mpfr_init(r17649); mpfr_init_set_str(r17650, "1/240", 10, MPFR_RNDN); mpfr_init_set_str(r17651, "5", 10, MPFR_RNDN); mpfr_init(r17652); mpfr_init(r17653); mpfr_init(r17654); mpfr_init_set_str(r17655, "1/24", 10, MPFR_RNDN); mpfr_init(r17656); mpfr_init(r17657); mpfr_init(r17658); mpfr_init(r17659); mpfr_init(r17660); mpfr_init(r17661); mpfr_init(r17662); mpfr_init(r17663); } double f_fm(double x) { mpfr_set_d(r17636, x, MPFR_RNDN); ; mpfr_set_si(r17638, mpfr_cmp(r17636, r17637) <= 0, MPFR_RNDN); mpfr_sin(r17639, r17636, MPFR_RNDN); mpfr_sqr(r17640, r17639, MPFR_RNDN); ; mpfr_cos(r17642, r17636, MPFR_RNDN); mpfr_add(r17643, r17641, r17642, MPFR_RNDN); mpfr_div(r17644, r17640, r17643, MPFR_RNDN); mpfr_div(r17645, r17644, r17639, MPFR_RNDN); ; mpfr_set_si(r17647, mpfr_cmp(r17636, r17646) <= 0, MPFR_RNDN); ; mpfr_mul(r17649, r17636, r17648, MPFR_RNDN); ; ; mpfr_pow(r17652, r17636, r17651, MPFR_RNDN); mpfr_mul(r17653, r17650, r17652, MPFR_RNDN); mpfr_mul(r17654, r17636, r17636, MPFR_RNDN); mpfr_mul(r17654, r17654, r17636, MPFR_RNDN); ; mpfr_mul(r17656, r17654, r17655, MPFR_RNDN); mpfr_add(r17657, r17653, r17656, MPFR_RNDN); mpfr_add(r17658, r17649, r17657, MPFR_RNDN); mpfr_sub(r17659, r17641, r17642, MPFR_RNDN); mpfr_div(r17660, r17641, r17639, MPFR_RNDN); mpfr_mul(r17661, r17659, r17660, MPFR_RNDN); if (mpfr_get_si(r17647, MPFR_RNDN)) { mpfr_set(r17662, r17658, MPFR_RNDN); } else { mpfr_set(r17662, r17661, MPFR_RNDN); }; if (mpfr_get_si(r17638, MPFR_RNDN)) { mpfr_set(r17663, r17645, MPFR_RNDN); } else { mpfr_set(r17663, r17662, MPFR_RNDN); }; return mpfr_get_d(r17663, MPFR_RNDN); } static mpfr_t r17664, r17665, r17666, r17667, r17668, r17669, r17670, r17671, r17672, r17673, r17674, r17675, r17676, r17677, r17678, r17679, r17680, r17681, r17682, r17683, r17684, r17685, r17686, r17687, r17688, r17689, r17690, r17691; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17664); mpfr_init_set_str(r17665, "-8.61646249982272e-12", 10, MPFR_RNDN); mpfr_init(r17666); mpfr_init(r17667); mpfr_init(r17668); mpfr_init_set_str(r17669, "1", 10, MPFR_RNDN); mpfr_init(r17670); mpfr_init(r17671); mpfr_init(r17672); mpfr_init(r17673); mpfr_init_set_str(r17674, "8.584780819621077", 10, MPFR_RNDN); mpfr_init(r17675); mpfr_init_set_str(r17676, "1/2", 10, MPFR_RNDN); mpfr_init(r17677); mpfr_init_set_str(r17678, "1/240", 10, MPFR_RNDN); mpfr_init_set_str(r17679, "5", 10, MPFR_RNDN); mpfr_init(r17680); mpfr_init(r17681); mpfr_init(r17682); mpfr_init_set_str(r17683, "1/24", 10, MPFR_RNDN); mpfr_init(r17684); mpfr_init(r17685); mpfr_init(r17686); mpfr_init(r17687); mpfr_init(r17688); mpfr_init(r17689); mpfr_init(r17690); mpfr_init(r17691); } double f_dm(double x) { mpfr_set_d(r17664, x, MPFR_RNDN); ; mpfr_set_si(r17666, mpfr_cmp(r17664, r17665) <= 0, MPFR_RNDN); mpfr_sin(r17667, r17664, MPFR_RNDN); mpfr_sqr(r17668, r17667, MPFR_RNDN); ; mpfr_cos(r17670, r17664, MPFR_RNDN); mpfr_add(r17671, r17669, r17670, MPFR_RNDN); mpfr_div(r17672, r17668, r17671, MPFR_RNDN); mpfr_div(r17673, r17672, r17667, MPFR_RNDN); ; mpfr_set_si(r17675, mpfr_cmp(r17664, r17674) <= 0, MPFR_RNDN); ; mpfr_mul(r17677, r17664, r17676, MPFR_RNDN); ; ; mpfr_pow(r17680, r17664, r17679, MPFR_RNDN); mpfr_mul(r17681, r17678, r17680, MPFR_RNDN); mpfr_mul(r17682, r17664, r17664, MPFR_RNDN); mpfr_mul(r17682, r17682, r17664, MPFR_RNDN); ; mpfr_mul(r17684, r17682, r17683, MPFR_RNDN); mpfr_add(r17685, r17681, r17684, MPFR_RNDN); mpfr_add(r17686, r17677, r17685, MPFR_RNDN); mpfr_sub(r17687, r17669, r17670, MPFR_RNDN); mpfr_div(r17688, r17669, r17667, MPFR_RNDN); mpfr_mul(r17689, r17687, r17688, MPFR_RNDN); if (mpfr_get_si(r17675, MPFR_RNDN)) { mpfr_set(r17690, r17686, MPFR_RNDN); } else { mpfr_set(r17690, r17689, MPFR_RNDN); }; if (mpfr_get_si(r17666, MPFR_RNDN)) { mpfr_set(r17691, r17673, MPFR_RNDN); } else { mpfr_set(r17691, r17690, MPFR_RNDN); }; return mpfr_get_d(r17691, MPFR_RNDN); }