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