#include #include #include #include #include char *name = "Harley's example"; double f_if(float c_p, float c_n, float t, float s) { float r24727 = 1; float r24728 = s; float r24729 = -r24728; float r24730 = exp(r24729); float r24731 = r24727 + r24730; float r24732 = r24727 / r24731; float r24733 = c_p; float r24734 = pow(r24732, r24733); float r24735 = r24727 - r24732; float r24736 = c_n; float r24737 = pow(r24735, r24736); float r24738 = r24734 * r24737; float r24739 = t; float r24740 = -r24739; float r24741 = exp(r24740); float r24742 = r24727 + r24741; float r24743 = r24727 / r24742; float r24744 = pow(r24743, r24733); float r24745 = r24727 - r24743; float r24746 = pow(r24745, r24736); float r24747 = r24744 * r24746; float r24748 = r24738 / r24747; return r24748; } double f_id(double c_p, double c_n, double t, double s) { double r24749 = 1; double r24750 = s; double r24751 = -r24750; double r24752 = exp(r24751); double r24753 = r24749 + r24752; double r24754 = r24749 / r24753; double r24755 = c_p; double r24756 = pow(r24754, r24755); double r24757 = r24749 - r24754; double r24758 = c_n; double r24759 = pow(r24757, r24758); double r24760 = r24756 * r24759; double r24761 = t; double r24762 = -r24761; double r24763 = exp(r24762); double r24764 = r24749 + r24763; double r24765 = r24749 / r24764; double r24766 = pow(r24765, r24755); double r24767 = r24749 - r24765; double r24768 = pow(r24767, r24758); double r24769 = r24766 * r24768; double r24770 = r24760 / r24769; return r24770; } double f_of(float c_p, float c_n, float t, float s) { float r24771 = -1; float r24772 = 1; float r24773 = s; float r24774 = -r24773; float r24775 = exp(r24774); float r24776 = r24772 + r24775; float r24777 = r24771 / r24776; float r24778 = log1p(r24777); float r24779 = c_n; float r24780 = r24778 * r24779; float r24781 = t; float r24782 = -r24781; float r24783 = exp(r24782); float r24784 = r24783 + r24772; float r24785 = r24771 / r24784; float r24786 = log1p(r24785); float r24787 = log1p(r24775); float r24788 = log1p(r24783); float r24789 = r24787 - r24788; float r24790 = c_p; float r24791 = r24789 * r24790; float r24792 = fma(r24786, r24779, r24791); float r24793 = r24780 - r24792; float r24794 = exp(r24793); return r24794; } double f_od(double c_p, double c_n, double t, double s) { double r24795 = -1; double r24796 = 1; double r24797 = s; double r24798 = -r24797; double r24799 = exp(r24798); double r24800 = r24796 + r24799; double r24801 = r24795 / r24800; double r24802 = log1p(r24801); double r24803 = c_n; double r24804 = r24802 * r24803; double r24805 = t; double r24806 = -r24805; double r24807 = exp(r24806); double r24808 = r24807 + r24796; double r24809 = r24795 / r24808; double r24810 = log1p(r24809); double r24811 = log1p(r24799); double r24812 = log1p(r24807); double r24813 = r24811 - r24812; double r24814 = c_p; double r24815 = r24813 * r24814; double r24816 = fma(r24810, r24803, r24815); double r24817 = r24804 - r24816; double r24818 = exp(r24817); return r24818; } 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 r24819, r24820, r24821, r24822, r24823, r24824, r24825, r24826, r24827, r24828, r24829, r24830, r24831, r24832, r24833, r24834, r24835, r24836, r24837, r24838, r24839, r24840; void setup_mpfr_f_im() { mpfr_set_default_prec(1872); mpfr_init_set_str(r24819, "1", 10, MPFR_RNDN); mpfr_init(r24820); mpfr_init(r24821); mpfr_init(r24822); mpfr_init(r24823); mpfr_init(r24824); mpfr_init(r24825); mpfr_init(r24826); mpfr_init(r24827); mpfr_init(r24828); mpfr_init(r24829); mpfr_init(r24830); mpfr_init(r24831); mpfr_init(r24832); mpfr_init(r24833); mpfr_init(r24834); mpfr_init(r24835); mpfr_init(r24836); mpfr_init(r24837); mpfr_init(r24838); mpfr_init(r24839); mpfr_init(r24840); } double f_im(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r24820, s, MPFR_RNDN); mpfr_neg(r24821, r24820, MPFR_RNDN); mpfr_exp(r24822, r24821, MPFR_RNDN); mpfr_add(r24823, r24819, r24822, MPFR_RNDN); mpfr_div(r24824, r24819, r24823, MPFR_RNDN); mpfr_set_d(r24825, c_p, MPFR_RNDN); mpfr_pow(r24826, r24824, r24825, MPFR_RNDN); mpfr_sub(r24827, r24819, r24824, MPFR_RNDN); mpfr_set_d(r24828, c_n, MPFR_RNDN); mpfr_pow(r24829, r24827, r24828, MPFR_RNDN); mpfr_mul(r24830, r24826, r24829, MPFR_RNDN); mpfr_set_d(r24831, t, MPFR_RNDN); mpfr_neg(r24832, r24831, MPFR_RNDN); mpfr_exp(r24833, r24832, MPFR_RNDN); mpfr_add(r24834, r24819, r24833, MPFR_RNDN); mpfr_div(r24835, r24819, r24834, MPFR_RNDN); mpfr_pow(r24836, r24835, r24825, MPFR_RNDN); mpfr_sub(r24837, r24819, r24835, MPFR_RNDN); mpfr_pow(r24838, r24837, r24828, MPFR_RNDN); mpfr_mul(r24839, r24836, r24838, MPFR_RNDN); mpfr_div(r24840, r24830, r24839, MPFR_RNDN); return mpfr_get_d(r24840, MPFR_RNDN); } static mpfr_t r24841, r24842, r24843, r24844, r24845, r24846, r24847, r24848, r24849, r24850, r24851, r24852, r24853, r24854, r24855, r24856, r24857, r24858, r24859, r24860, r24861, r24862, r24863, r24864; void setup_mpfr_f_fm() { mpfr_set_default_prec(1872); mpfr_init_set_str(r24841, "-1", 10, MPFR_RNDN); mpfr_init_set_str(r24842, "1", 10, MPFR_RNDN); mpfr_init(r24843); mpfr_init(r24844); mpfr_init(r24845); mpfr_init(r24846); mpfr_init(r24847); mpfr_init(r24848); mpfr_init(r24849); mpfr_init(r24850); mpfr_init(r24851); mpfr_init(r24852); mpfr_init(r24853); mpfr_init(r24854); mpfr_init(r24855); mpfr_init(r24856); mpfr_init(r24857); mpfr_init(r24858); mpfr_init(r24859); mpfr_init(r24860); mpfr_init(r24861); mpfr_init(r24862); mpfr_init(r24863); mpfr_init(r24864); } double f_fm(double c_p, double c_n, double t, double s) { ; ; mpfr_set_d(r24843, s, MPFR_RNDN); mpfr_neg(r24844, r24843, MPFR_RNDN); mpfr_exp(r24845, r24844, MPFR_RNDN); mpfr_add(r24846, r24842, r24845, MPFR_RNDN); mpfr_div(r24847, r24841, r24846, MPFR_RNDN); mpfr_log1p(r24848, r24847, MPFR_RNDN); mpfr_set_d(r24849, c_n, MPFR_RNDN); mpfr_mul(r24850, r24848, r24849, MPFR_RNDN); mpfr_set_d(r24851, t, MPFR_RNDN); mpfr_neg(r24852, r24851, MPFR_RNDN); mpfr_exp(r24853, r24852, MPFR_RNDN); mpfr_add(r24854, r24853, r24842, MPFR_RNDN); mpfr_div(r24855, r24841, r24854, MPFR_RNDN); mpfr_log1p(r24856, r24855, MPFR_RNDN); mpfr_log1p(r24857, r24845, MPFR_RNDN); mpfr_log1p(r24858, r24853, MPFR_RNDN); mpfr_sub(r24859, r24857, r24858, MPFR_RNDN); mpfr_set_d(r24860, c_p, MPFR_RNDN); mpfr_mul(r24861, r24859, r24860, MPFR_RNDN); mpfr_fma(r24862, r24856, r24849, r24861, MPFR_RNDN); mpfr_sub(r24863, r24850, r24862, MPFR_RNDN); mpfr_exp(r24864, r24863, MPFR_RNDN); return mpfr_get_d(r24864, MPFR_RNDN); } static mpfr_t r24865, r24866, r24867, r24868, r24869, r24870, r24871, r24872, r24873, r24874, r24875, r24876, r24877, r24878, r24879, r24880, r24881, r24882, r24883, r24884, r24885, r24886, r24887, r24888; void setup_mpfr_f_dm() { mpfr_set_default_prec(1872); mpfr_init_set_str(r24865, "-1", 10, MPFR_RNDN); mpfr_init_set_str(r24866, "1", 10, MPFR_RNDN); mpfr_init(r24867); mpfr_init(r24868); mpfr_init(r24869); mpfr_init(r24870); mpfr_init(r24871); mpfr_init(r24872); mpfr_init(r24873); mpfr_init(r24874); mpfr_init(r24875); mpfr_init(r24876); mpfr_init(r24877); mpfr_init(r24878); mpfr_init(r24879); mpfr_init(r24880); mpfr_init(r24881); mpfr_init(r24882); mpfr_init(r24883); mpfr_init(r24884); mpfr_init(r24885); mpfr_init(r24886); mpfr_init(r24887); mpfr_init(r24888); } double f_dm(double c_p, double c_n, double t, double s) { ; ; mpfr_set_d(r24867, s, MPFR_RNDN); mpfr_neg(r24868, r24867, MPFR_RNDN); mpfr_exp(r24869, r24868, MPFR_RNDN); mpfr_add(r24870, r24866, r24869, MPFR_RNDN); mpfr_div(r24871, r24865, r24870, MPFR_RNDN); mpfr_log1p(r24872, r24871, MPFR_RNDN); mpfr_set_d(r24873, c_n, MPFR_RNDN); mpfr_mul(r24874, r24872, r24873, MPFR_RNDN); mpfr_set_d(r24875, t, MPFR_RNDN); mpfr_neg(r24876, r24875, MPFR_RNDN); mpfr_exp(r24877, r24876, MPFR_RNDN); mpfr_add(r24878, r24877, r24866, MPFR_RNDN); mpfr_div(r24879, r24865, r24878, MPFR_RNDN); mpfr_log1p(r24880, r24879, MPFR_RNDN); mpfr_log1p(r24881, r24869, MPFR_RNDN); mpfr_log1p(r24882, r24877, MPFR_RNDN); mpfr_sub(r24883, r24881, r24882, MPFR_RNDN); mpfr_set_d(r24884, c_p, MPFR_RNDN); mpfr_mul(r24885, r24883, r24884, MPFR_RNDN); mpfr_fma(r24886, r24880, r24873, r24885, MPFR_RNDN); mpfr_sub(r24887, r24874, r24886, MPFR_RNDN); mpfr_exp(r24888, r24887, MPFR_RNDN); return mpfr_get_d(r24888, MPFR_RNDN); }