#include #include #include #include #include char *name = "expax (section 3.5)"; double f_if(float a, float x) { float r10642 = a; float r10643 = x; float r10644 = r10642 * r10643; float r10645 = exp(r10644); float r10646 = 1; float r10647 = r10645 - r10646; return r10647; } double f_id(double a, double x) { double r10648 = a; double r10649 = x; double r10650 = r10648 * r10649; double r10651 = exp(r10650); double r10652 = 1; double r10653 = r10651 - r10652; return r10653; } double f_of(float a, float x) { float r10654 = a; float r10655 = x; float r10656 = r10654 * r10655; float r10657 = exp(r10656); float r10658 = 1; float r10659 = r10657 - r10658; float r10660 = -5.2697850804707055e-06; bool r10661 = r10659 <= r10660; float r10662 = cbrt(r10659); float r10663 = r10662 * r10662; float r10664 = cbrt(r10663); float r10665 = sqrt(r10657); float r10666 = r10665 + r10658; float r10667 = r10665 - r10658; float r10668 = r10666 * r10667; float r10669 = cbrt(r10668); float r10670 = cbrt(r10669); float r10671 = r10664 * r10670; float r10672 = r10663 * r10671; float r10673 = r10656 * r10656; float r10674 = 1/6; float r10675 = r10655 * r10674; float r10676 = r10675 * r10654; float r10677 = 1/2; float r10678 = r10676 + r10677; float r10679 = r10673 * r10678; float r10680 = r10679 + r10656; float r10681 = r10661 ? r10672 : r10680; return r10681; } double f_od(double a, double x) { double r10682 = a; double r10683 = x; double r10684 = r10682 * r10683; double r10685 = exp(r10684); double r10686 = 1; double r10687 = r10685 - r10686; double r10688 = -5.2697850804707055e-06; bool r10689 = r10687 <= r10688; double r10690 = cbrt(r10687); double r10691 = r10690 * r10690; double r10692 = cbrt(r10691); double r10693 = sqrt(r10685); double r10694 = r10693 + r10686; double r10695 = r10693 - r10686; double r10696 = r10694 * r10695; double r10697 = cbrt(r10696); double r10698 = cbrt(r10697); double r10699 = r10692 * r10698; double r10700 = r10691 * r10699; double r10701 = r10684 * r10684; double r10702 = 1/6; double r10703 = r10683 * r10702; double r10704 = r10703 * r10682; double r10705 = 1/2; double r10706 = r10704 + r10705; double r10707 = r10701 * r10706; double r10708 = r10707 + r10684; double r10709 = r10689 ? r10700 : r10708; return r10709; } 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 r10710, r10711, r10712, r10713, r10714, r10715; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r10710); mpfr_init(r10711); mpfr_init(r10712); mpfr_init(r10713); mpfr_init_set_str(r10714, "1", 10, MPFR_RNDN); mpfr_init(r10715); } double f_im(double a, double x) { mpfr_set_d(r10710, a, MPFR_RNDN); mpfr_set_d(r10711, x, MPFR_RNDN); mpfr_mul(r10712, r10710, r10711, MPFR_RNDN); mpfr_exp(r10713, r10712, MPFR_RNDN); ; mpfr_sub(r10715, r10713, r10714, MPFR_RNDN); return mpfr_get_d(r10715, MPFR_RNDN); } static mpfr_t r10716, r10717, r10718, r10719, r10720, r10721, r10722, r10723, r10724, r10725, r10726, r10727, r10728, r10729, r10730, r10731, r10732, r10733, r10734, r10735, r10736, r10737, r10738, r10739, r10740, r10741, r10742, r10743; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r10716); mpfr_init(r10717); mpfr_init(r10718); mpfr_init(r10719); mpfr_init_set_str(r10720, "1", 10, MPFR_RNDN); mpfr_init(r10721); mpfr_init_set_str(r10722, "-5.2697850804707055e-06", 10, MPFR_RNDN); mpfr_init(r10723); mpfr_init(r10724); mpfr_init(r10725); mpfr_init(r10726); mpfr_init(r10727); mpfr_init(r10728); mpfr_init(r10729); mpfr_init(r10730); mpfr_init(r10731); mpfr_init(r10732); mpfr_init(r10733); mpfr_init(r10734); mpfr_init(r10735); mpfr_init_set_str(r10736, "1/6", 10, MPFR_RNDN); mpfr_init(r10737); mpfr_init(r10738); mpfr_init_set_str(r10739, "1/2", 10, MPFR_RNDN); mpfr_init(r10740); mpfr_init(r10741); mpfr_init(r10742); mpfr_init(r10743); } double f_fm(double a, double x) { mpfr_set_d(r10716, a, MPFR_RNDN); mpfr_set_d(r10717, x, MPFR_RNDN); mpfr_mul(r10718, r10716, r10717, MPFR_RNDN); mpfr_exp(r10719, r10718, MPFR_RNDN); ; mpfr_sub(r10721, r10719, r10720, MPFR_RNDN); ; mpfr_set_si(r10723, mpfr_cmp(r10721, r10722) <= 0, MPFR_RNDN); mpfr_cbrt(r10724, r10721, MPFR_RNDN); mpfr_mul(r10725, r10724, r10724, MPFR_RNDN); mpfr_cbrt(r10726, r10725, MPFR_RNDN); mpfr_sqrt(r10727, r10719, MPFR_RNDN); mpfr_add(r10728, r10727, r10720, MPFR_RNDN); mpfr_sub(r10729, r10727, r10720, MPFR_RNDN); mpfr_mul(r10730, r10728, r10729, MPFR_RNDN); mpfr_cbrt(r10731, r10730, MPFR_RNDN); mpfr_cbrt(r10732, r10731, MPFR_RNDN); mpfr_mul(r10733, r10726, r10732, MPFR_RNDN); mpfr_mul(r10734, r10725, r10733, MPFR_RNDN); mpfr_mul(r10735, r10718, r10718, MPFR_RNDN); ; mpfr_mul(r10737, r10717, r10736, MPFR_RNDN); mpfr_mul(r10738, r10737, r10716, MPFR_RNDN); ; mpfr_add(r10740, r10738, r10739, MPFR_RNDN); mpfr_mul(r10741, r10735, r10740, MPFR_RNDN); mpfr_add(r10742, r10741, r10718, MPFR_RNDN); if (mpfr_get_si(r10723, MPFR_RNDN)) { mpfr_set(r10743, r10734, MPFR_RNDN); } else { mpfr_set(r10743, r10742, MPFR_RNDN); }; return mpfr_get_d(r10743, MPFR_RNDN); } static mpfr_t r10744, r10745, r10746, r10747, r10748, r10749, r10750, r10751, r10752, r10753, r10754, r10755, r10756, r10757, r10758, r10759, r10760, r10761, r10762, r10763, r10764, r10765, r10766, r10767, r10768, r10769, r10770, r10771; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r10744); mpfr_init(r10745); mpfr_init(r10746); mpfr_init(r10747); mpfr_init_set_str(r10748, "1", 10, MPFR_RNDN); mpfr_init(r10749); mpfr_init_set_str(r10750, "-5.2697850804707055e-06", 10, MPFR_RNDN); mpfr_init(r10751); mpfr_init(r10752); mpfr_init(r10753); mpfr_init(r10754); mpfr_init(r10755); mpfr_init(r10756); mpfr_init(r10757); mpfr_init(r10758); mpfr_init(r10759); mpfr_init(r10760); mpfr_init(r10761); mpfr_init(r10762); mpfr_init(r10763); mpfr_init_set_str(r10764, "1/6", 10, MPFR_RNDN); mpfr_init(r10765); mpfr_init(r10766); mpfr_init_set_str(r10767, "1/2", 10, MPFR_RNDN); mpfr_init(r10768); mpfr_init(r10769); mpfr_init(r10770); mpfr_init(r10771); } double f_dm(double a, double x) { mpfr_set_d(r10744, a, MPFR_RNDN); mpfr_set_d(r10745, x, MPFR_RNDN); mpfr_mul(r10746, r10744, r10745, MPFR_RNDN); mpfr_exp(r10747, r10746, MPFR_RNDN); ; mpfr_sub(r10749, r10747, r10748, MPFR_RNDN); ; mpfr_set_si(r10751, mpfr_cmp(r10749, r10750) <= 0, MPFR_RNDN); mpfr_cbrt(r10752, r10749, MPFR_RNDN); mpfr_mul(r10753, r10752, r10752, MPFR_RNDN); mpfr_cbrt(r10754, r10753, MPFR_RNDN); mpfr_sqrt(r10755, r10747, MPFR_RNDN); mpfr_add(r10756, r10755, r10748, MPFR_RNDN); mpfr_sub(r10757, r10755, r10748, MPFR_RNDN); mpfr_mul(r10758, r10756, r10757, MPFR_RNDN); mpfr_cbrt(r10759, r10758, MPFR_RNDN); mpfr_cbrt(r10760, r10759, MPFR_RNDN); mpfr_mul(r10761, r10754, r10760, MPFR_RNDN); mpfr_mul(r10762, r10753, r10761, MPFR_RNDN); mpfr_mul(r10763, r10746, r10746, MPFR_RNDN); ; mpfr_mul(r10765, r10745, r10764, MPFR_RNDN); mpfr_mul(r10766, r10765, r10744, MPFR_RNDN); ; mpfr_add(r10768, r10766, r10767, MPFR_RNDN); mpfr_mul(r10769, r10763, r10768, MPFR_RNDN); mpfr_add(r10770, r10769, r10746, MPFR_RNDN); if (mpfr_get_si(r10751, MPFR_RNDN)) { mpfr_set(r10771, r10762, MPFR_RNDN); } else { mpfr_set(r10771, r10770, MPFR_RNDN); }; return mpfr_get_d(r10771, MPFR_RNDN); }