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