#include #include #include #include #include char *name = "Harley's example"; double f_if(float c_p, float c_n, float t, float s) { float r27614 = 1; float r27615 = s; float r27616 = -r27615; float r27617 = exp(r27616); float r27618 = r27614 + r27617; float r27619 = r27614 / r27618; float r27620 = c_p; float r27621 = pow(r27619, r27620); float r27622 = r27614 - r27619; float r27623 = c_n; float r27624 = pow(r27622, r27623); float r27625 = r27621 * r27624; float r27626 = t; float r27627 = -r27626; float r27628 = exp(r27627); float r27629 = r27614 + r27628; float r27630 = r27614 / r27629; float r27631 = pow(r27630, r27620); float r27632 = r27614 - r27630; float r27633 = pow(r27632, r27623); float r27634 = r27631 * r27633; float r27635 = r27625 / r27634; return r27635; } double f_id(double c_p, double c_n, double t, double s) { double r27636 = 1; double r27637 = s; double r27638 = -r27637; double r27639 = exp(r27638); double r27640 = r27636 + r27639; double r27641 = r27636 / r27640; double r27642 = c_p; double r27643 = pow(r27641, r27642); double r27644 = r27636 - r27641; double r27645 = c_n; double r27646 = pow(r27644, r27645); double r27647 = r27643 * r27646; double r27648 = t; double r27649 = -r27648; double r27650 = exp(r27649); double r27651 = r27636 + r27650; double r27652 = r27636 / r27651; double r27653 = pow(r27652, r27642); double r27654 = r27636 - r27652; double r27655 = pow(r27654, r27645); double r27656 = r27653 * r27655; double r27657 = r27647 / r27656; return r27657; } double f_of(float c_p, float c_n, float t, float s) { float r27658 = 1; float r27659 = -r27658; float r27660 = s; float r27661 = -r27660; float r27662 = exp(r27661); float r27663 = r27658 + r27662; float r27664 = r27659 / r27663; float r27665 = log1p(r27664); float r27666 = t; float r27667 = -r27666; float r27668 = exp(r27667); float r27669 = r27668 + r27658; float r27670 = r27659 / r27669; float r27671 = log1p(r27670); float r27672 = r27665 - r27671; float r27673 = c_n; float r27674 = c_p; float r27675 = log1p(r27668); float r27676 = log1p(r27662); float r27677 = r27675 - r27676; float r27678 = r27674 * r27677; float r27679 = fma(r27672, r27673, r27678); float r27680 = exp(r27679); return r27680; } double f_od(double c_p, double c_n, double t, double s) { double r27681 = 1; double r27682 = -r27681; double r27683 = s; double r27684 = -r27683; double r27685 = exp(r27684); double r27686 = r27681 + r27685; double r27687 = r27682 / r27686; double r27688 = log1p(r27687); double r27689 = t; double r27690 = -r27689; double r27691 = exp(r27690); double r27692 = r27691 + r27681; double r27693 = r27682 / r27692; double r27694 = log1p(r27693); double r27695 = r27688 - r27694; double r27696 = c_n; double r27697 = c_p; double r27698 = log1p(r27691); double r27699 = log1p(r27685); double r27700 = r27698 - r27699; double r27701 = r27697 * r27700; double r27702 = fma(r27695, r27696, r27701); double r27703 = exp(r27702); return r27703; } 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 r27704, r27705, r27706, r27707, r27708, r27709, r27710, r27711, r27712, r27713, r27714, r27715, r27716, r27717, r27718, r27719, r27720, r27721, r27722, r27723, r27724, r27725; void setup_mpfr_f_im() { mpfr_set_default_prec(2128); mpfr_init_set_str(r27704, "1", 10, MPFR_RNDN); mpfr_init(r27705); mpfr_init(r27706); mpfr_init(r27707); mpfr_init(r27708); mpfr_init(r27709); mpfr_init(r27710); mpfr_init(r27711); mpfr_init(r27712); mpfr_init(r27713); mpfr_init(r27714); mpfr_init(r27715); mpfr_init(r27716); mpfr_init(r27717); mpfr_init(r27718); mpfr_init(r27719); mpfr_init(r27720); mpfr_init(r27721); mpfr_init(r27722); mpfr_init(r27723); mpfr_init(r27724); mpfr_init(r27725); } double f_im(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r27705, s, MPFR_RNDN); mpfr_neg(r27706, r27705, MPFR_RNDN); mpfr_exp(r27707, r27706, MPFR_RNDN); mpfr_add(r27708, r27704, r27707, MPFR_RNDN); mpfr_div(r27709, r27704, r27708, MPFR_RNDN); mpfr_set_d(r27710, c_p, MPFR_RNDN); mpfr_pow(r27711, r27709, r27710, MPFR_RNDN); mpfr_sub(r27712, r27704, r27709, MPFR_RNDN); mpfr_set_d(r27713, c_n, MPFR_RNDN); mpfr_pow(r27714, r27712, r27713, MPFR_RNDN); mpfr_mul(r27715, r27711, r27714, MPFR_RNDN); mpfr_set_d(r27716, t, MPFR_RNDN); mpfr_neg(r27717, r27716, MPFR_RNDN); mpfr_exp(r27718, r27717, MPFR_RNDN); mpfr_add(r27719, r27704, r27718, MPFR_RNDN); mpfr_div(r27720, r27704, r27719, MPFR_RNDN); mpfr_pow(r27721, r27720, r27710, MPFR_RNDN); mpfr_sub(r27722, r27704, r27720, MPFR_RNDN); mpfr_pow(r27723, r27722, r27713, MPFR_RNDN); mpfr_mul(r27724, r27721, r27723, MPFR_RNDN); mpfr_div(r27725, r27715, r27724, MPFR_RNDN); return mpfr_get_d(r27725, MPFR_RNDN); } static mpfr_t r27726, r27727, r27728, r27729, r27730, r27731, r27732, r27733, r27734, r27735, r27736, r27737, r27738, r27739, r27740, r27741, r27742, r27743, r27744, r27745, r27746, r27747, r27748; void setup_mpfr_f_fm() { mpfr_set_default_prec(2128); mpfr_init_set_str(r27726, "1", 10, MPFR_RNDN); mpfr_init(r27727); mpfr_init(r27728); mpfr_init(r27729); mpfr_init(r27730); mpfr_init(r27731); mpfr_init(r27732); mpfr_init(r27733); mpfr_init(r27734); mpfr_init(r27735); mpfr_init(r27736); mpfr_init(r27737); mpfr_init(r27738); mpfr_init(r27739); mpfr_init(r27740); mpfr_init(r27741); mpfr_init(r27742); mpfr_init(r27743); mpfr_init(r27744); mpfr_init(r27745); mpfr_init(r27746); mpfr_init(r27747); mpfr_init(r27748); } double f_fm(double c_p, double c_n, double t, double s) { ; mpfr_neg(r27727, r27726, MPFR_RNDN); mpfr_set_d(r27728, s, MPFR_RNDN); mpfr_neg(r27729, r27728, MPFR_RNDN); mpfr_exp(r27730, r27729, MPFR_RNDN); mpfr_add(r27731, r27726, r27730, MPFR_RNDN); mpfr_div(r27732, r27727, r27731, MPFR_RNDN); mpfr_log1p(r27733, r27732, MPFR_RNDN); mpfr_set_d(r27734, t, MPFR_RNDN); mpfr_neg(r27735, r27734, MPFR_RNDN); mpfr_exp(r27736, r27735, MPFR_RNDN); mpfr_add(r27737, r27736, r27726, MPFR_RNDN); mpfr_div(r27738, r27727, r27737, MPFR_RNDN); mpfr_log1p(r27739, r27738, MPFR_RNDN); mpfr_sub(r27740, r27733, r27739, MPFR_RNDN); mpfr_set_d(r27741, c_n, MPFR_RNDN); mpfr_set_d(r27742, c_p, MPFR_RNDN); mpfr_log1p(r27743, r27736, MPFR_RNDN); mpfr_log1p(r27744, r27730, MPFR_RNDN); mpfr_sub(r27745, r27743, r27744, MPFR_RNDN); mpfr_mul(r27746, r27742, r27745, MPFR_RNDN); mpfr_fma(r27747, r27740, r27741, r27746, MPFR_RNDN); mpfr_exp(r27748, r27747, MPFR_RNDN); return mpfr_get_d(r27748, MPFR_RNDN); } static mpfr_t r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756, r27757, r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767, r27768, r27769, r27770, r27771; void setup_mpfr_f_dm() { mpfr_set_default_prec(2128); mpfr_init_set_str(r27749, "1", 10, MPFR_RNDN); mpfr_init(r27750); mpfr_init(r27751); mpfr_init(r27752); mpfr_init(r27753); mpfr_init(r27754); mpfr_init(r27755); mpfr_init(r27756); mpfr_init(r27757); mpfr_init(r27758); mpfr_init(r27759); mpfr_init(r27760); mpfr_init(r27761); mpfr_init(r27762); mpfr_init(r27763); mpfr_init(r27764); mpfr_init(r27765); mpfr_init(r27766); mpfr_init(r27767); mpfr_init(r27768); mpfr_init(r27769); mpfr_init(r27770); mpfr_init(r27771); } double f_dm(double c_p, double c_n, double t, double s) { ; mpfr_neg(r27750, r27749, MPFR_RNDN); mpfr_set_d(r27751, s, MPFR_RNDN); mpfr_neg(r27752, r27751, MPFR_RNDN); mpfr_exp(r27753, r27752, MPFR_RNDN); mpfr_add(r27754, r27749, r27753, MPFR_RNDN); mpfr_div(r27755, r27750, r27754, MPFR_RNDN); mpfr_log1p(r27756, r27755, MPFR_RNDN); mpfr_set_d(r27757, t, MPFR_RNDN); mpfr_neg(r27758, r27757, MPFR_RNDN); mpfr_exp(r27759, r27758, MPFR_RNDN); mpfr_add(r27760, r27759, r27749, MPFR_RNDN); mpfr_div(r27761, r27750, r27760, MPFR_RNDN); mpfr_log1p(r27762, r27761, MPFR_RNDN); mpfr_sub(r27763, r27756, r27762, MPFR_RNDN); mpfr_set_d(r27764, c_n, MPFR_RNDN); mpfr_set_d(r27765, c_p, MPFR_RNDN); mpfr_log1p(r27766, r27759, MPFR_RNDN); mpfr_log1p(r27767, r27753, MPFR_RNDN); mpfr_sub(r27768, r27766, r27767, MPFR_RNDN); mpfr_mul(r27769, r27765, r27768, MPFR_RNDN); mpfr_fma(r27770, r27763, r27764, r27769, MPFR_RNDN); mpfr_exp(r27771, r27770, MPFR_RNDN); return mpfr_get_d(r27771, MPFR_RNDN); }