#include #include #include #include #include char *name = "Harley's example"; double f_if(float c_p, float c_n, float t, float s) { float r27652 = 1; float r27653 = s; float r27654 = -r27653; float r27655 = exp(r27654); float r27656 = r27652 + r27655; float r27657 = r27652 / r27656; float r27658 = c_p; float r27659 = pow(r27657, r27658); float r27660 = r27652 - r27657; float r27661 = c_n; float r27662 = pow(r27660, r27661); float r27663 = r27659 * r27662; float r27664 = t; float r27665 = -r27664; float r27666 = exp(r27665); float r27667 = r27652 + r27666; float r27668 = r27652 / r27667; float r27669 = pow(r27668, r27658); float r27670 = r27652 - r27668; float r27671 = pow(r27670, r27661); float r27672 = r27669 * r27671; float r27673 = r27663 / r27672; return r27673; } double f_id(double c_p, double c_n, double t, double s) { double r27674 = 1; double r27675 = s; double r27676 = -r27675; double r27677 = exp(r27676); double r27678 = r27674 + r27677; double r27679 = r27674 / r27678; double r27680 = c_p; double r27681 = pow(r27679, r27680); double r27682 = r27674 - r27679; double r27683 = c_n; double r27684 = pow(r27682, r27683); double r27685 = r27681 * r27684; double r27686 = t; double r27687 = -r27686; double r27688 = exp(r27687); double r27689 = r27674 + r27688; double r27690 = r27674 / r27689; double r27691 = pow(r27690, r27680); double r27692 = r27674 - r27690; double r27693 = pow(r27692, r27683); double r27694 = r27691 * r27693; double r27695 = r27685 / r27694; return r27695; } double f_of(float c_p, float c_n, float t, float s) { float r27696 = 1; float r27697 = s; float r27698 = -r27697; float r27699 = exp(r27698); float r27700 = r27696 + r27699; float r27701 = r27696 / r27700; float r27702 = r27696 - r27701; float r27703 = c_n; float r27704 = pow(r27702, r27703); float r27705 = t; float r27706 = -r27705; float r27707 = exp(r27706); float r27708 = r27696 + r27707; float r27709 = r27696 / r27708; float r27710 = r27696 - r27709; float r27711 = pow(r27710, r27703); float r27712 = c_p; float r27713 = 1/2; float r27714 = log(r27713); float r27715 = fma(r27713, r27705, r27714); float r27716 = fma(r27712, r27715, r27696); float r27717 = pow(r27701, r27712); float r27718 = r27716 / r27717; float r27719 = r27711 * r27718; float r27720 = r27704 / r27719; return r27720; } double f_od(double c_p, double c_n, double t, double s) { double r27721 = 1; double r27722 = s; double r27723 = -r27722; double r27724 = exp(r27723); double r27725 = r27721 + r27724; double r27726 = r27721 / r27725; double r27727 = r27721 - r27726; double r27728 = c_n; double r27729 = pow(r27727, r27728); double r27730 = t; double r27731 = -r27730; double r27732 = exp(r27731); double r27733 = r27721 + r27732; double r27734 = r27721 / r27733; double r27735 = r27721 - r27734; double r27736 = pow(r27735, r27728); double r27737 = c_p; double r27738 = 1/2; double r27739 = log(r27738); double r27740 = fma(r27738, r27730, r27739); double r27741 = fma(r27737, r27740, r27721); double r27742 = pow(r27726, r27737); double r27743 = r27741 / r27742; double r27744 = r27736 * r27743; double r27745 = r27729 / r27744; return r27745; } 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 r27746, r27747, r27748, r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756, r27757, r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767; void setup_mpfr_f_im() { mpfr_set_default_prec(2640); mpfr_init_set_str(r27746, "1", 10, MPFR_RNDN); 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); mpfr_init(r27761); mpfr_init(r27762); mpfr_init(r27763); mpfr_init(r27764); mpfr_init(r27765); mpfr_init(r27766); mpfr_init(r27767); } double f_im(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r27747, s, MPFR_RNDN); mpfr_neg(r27748, r27747, MPFR_RNDN); mpfr_exp(r27749, r27748, MPFR_RNDN); mpfr_add(r27750, r27746, r27749, MPFR_RNDN); mpfr_div(r27751, r27746, r27750, MPFR_RNDN); mpfr_set_d(r27752, c_p, MPFR_RNDN); mpfr_pow(r27753, r27751, r27752, MPFR_RNDN); mpfr_sub(r27754, r27746, r27751, MPFR_RNDN); mpfr_set_d(r27755, c_n, MPFR_RNDN); mpfr_pow(r27756, r27754, r27755, MPFR_RNDN); mpfr_mul(r27757, r27753, r27756, MPFR_RNDN); mpfr_set_d(r27758, t, MPFR_RNDN); mpfr_neg(r27759, r27758, MPFR_RNDN); mpfr_exp(r27760, r27759, MPFR_RNDN); mpfr_add(r27761, r27746, r27760, MPFR_RNDN); mpfr_div(r27762, r27746, r27761, MPFR_RNDN); mpfr_pow(r27763, r27762, r27752, MPFR_RNDN); mpfr_sub(r27764, r27746, r27762, MPFR_RNDN); mpfr_pow(r27765, r27764, r27755, MPFR_RNDN); mpfr_mul(r27766, r27763, r27765, MPFR_RNDN); mpfr_div(r27767, r27757, r27766, MPFR_RNDN); return mpfr_get_d(r27767, MPFR_RNDN); } static mpfr_t r27768, r27769, r27770, r27771, r27772, r27773, r27774, r27775, r27776, r27777, r27778, r27779, r27780, r27781, r27782, r27783, r27784, r27785, r27786, r27787, r27788, r27789, r27790, r27791, r27792; void setup_mpfr_f_fm() { mpfr_set_default_prec(2640); mpfr_init_set_str(r27768, "1", 10, MPFR_RNDN); 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); mpfr_init_set_str(r27785, "1/2", 10, MPFR_RNDN); mpfr_init(r27786); mpfr_init(r27787); mpfr_init(r27788); mpfr_init(r27789); mpfr_init(r27790); mpfr_init(r27791); mpfr_init(r27792); } double f_fm(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r27769, s, MPFR_RNDN); mpfr_neg(r27770, r27769, MPFR_RNDN); mpfr_exp(r27771, r27770, MPFR_RNDN); mpfr_add(r27772, r27768, r27771, MPFR_RNDN); mpfr_div(r27773, r27768, r27772, MPFR_RNDN); mpfr_sub(r27774, r27768, r27773, MPFR_RNDN); mpfr_set_d(r27775, c_n, MPFR_RNDN); mpfr_pow(r27776, r27774, r27775, MPFR_RNDN); mpfr_set_d(r27777, t, MPFR_RNDN); mpfr_neg(r27778, r27777, MPFR_RNDN); mpfr_exp(r27779, r27778, MPFR_RNDN); mpfr_add(r27780, r27768, r27779, MPFR_RNDN); mpfr_div(r27781, r27768, r27780, MPFR_RNDN); mpfr_sub(r27782, r27768, r27781, MPFR_RNDN); mpfr_pow(r27783, r27782, r27775, MPFR_RNDN); mpfr_set_d(r27784, c_p, MPFR_RNDN); ; mpfr_log(r27786, r27785, MPFR_RNDN); mpfr_fma(r27787, r27785, r27777, r27786, MPFR_RNDN); mpfr_fma(r27788, r27784, r27787, r27768, MPFR_RNDN); mpfr_pow(r27789, r27773, r27784, MPFR_RNDN); mpfr_div(r27790, r27788, r27789, MPFR_RNDN); mpfr_mul(r27791, r27783, r27790, MPFR_RNDN); mpfr_div(r27792, r27776, r27791, MPFR_RNDN); return mpfr_get_d(r27792, MPFR_RNDN); } static mpfr_t r27793, r27794, r27795, r27796, r27797, r27798, r27799, r27800, r27801, r27802, r27803, r27804, r27805, r27806, r27807, r27808, r27809, r27810, r27811, r27812, r27813, r27814, r27815, r27816, r27817; void setup_mpfr_f_dm() { mpfr_set_default_prec(2640); mpfr_init_set_str(r27793, "1", 10, MPFR_RNDN); mpfr_init(r27794); mpfr_init(r27795); mpfr_init(r27796); mpfr_init(r27797); mpfr_init(r27798); mpfr_init(r27799); mpfr_init(r27800); mpfr_init(r27801); mpfr_init(r27802); mpfr_init(r27803); mpfr_init(r27804); mpfr_init(r27805); mpfr_init(r27806); mpfr_init(r27807); mpfr_init(r27808); mpfr_init(r27809); mpfr_init_set_str(r27810, "1/2", 10, MPFR_RNDN); mpfr_init(r27811); mpfr_init(r27812); mpfr_init(r27813); mpfr_init(r27814); mpfr_init(r27815); mpfr_init(r27816); mpfr_init(r27817); } double f_dm(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r27794, s, MPFR_RNDN); mpfr_neg(r27795, r27794, MPFR_RNDN); mpfr_exp(r27796, r27795, MPFR_RNDN); mpfr_add(r27797, r27793, r27796, MPFR_RNDN); mpfr_div(r27798, r27793, r27797, MPFR_RNDN); mpfr_sub(r27799, r27793, r27798, MPFR_RNDN); mpfr_set_d(r27800, c_n, MPFR_RNDN); mpfr_pow(r27801, r27799, r27800, MPFR_RNDN); mpfr_set_d(r27802, t, MPFR_RNDN); mpfr_neg(r27803, r27802, MPFR_RNDN); mpfr_exp(r27804, r27803, MPFR_RNDN); mpfr_add(r27805, r27793, r27804, MPFR_RNDN); mpfr_div(r27806, r27793, r27805, MPFR_RNDN); mpfr_sub(r27807, r27793, r27806, MPFR_RNDN); mpfr_pow(r27808, r27807, r27800, MPFR_RNDN); mpfr_set_d(r27809, c_p, MPFR_RNDN); ; mpfr_log(r27811, r27810, MPFR_RNDN); mpfr_fma(r27812, r27810, r27802, r27811, MPFR_RNDN); mpfr_fma(r27813, r27809, r27812, r27793, MPFR_RNDN); mpfr_pow(r27814, r27798, r27809, MPFR_RNDN); mpfr_div(r27815, r27813, r27814, MPFR_RNDN); mpfr_mul(r27816, r27808, r27815, MPFR_RNDN); mpfr_div(r27817, r27801, r27816, MPFR_RNDN); return mpfr_get_d(r27817, MPFR_RNDN); }