#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r24756 = 2; float r24757 = sqrt(r24756); float r24758 = 4; float r24759 = r24757 / r24758; float r24760 = 1; float r24761 = 3; float r24762 = v; float r24763 = r24762 * r24762; float r24764 = r24761 * r24763; float r24765 = r24760 - r24764; float r24766 = sqrt(r24765); float r24767 = r24759 * r24766; float r24768 = r24760 - r24763; float r24769 = r24767 * r24768; return r24769; } double f_id(double v) { double r24770 = 2; double r24771 = sqrt(r24770); double r24772 = 4; double r24773 = r24771 / r24772; double r24774 = 1; double r24775 = 3; double r24776 = v; double r24777 = r24776 * r24776; double r24778 = r24775 * r24777; double r24779 = r24774 - r24778; double r24780 = sqrt(r24779); double r24781 = r24773 * r24780; double r24782 = r24774 - r24777; double r24783 = r24781 * r24782; return r24783; } double f_of(float v) { float r24784 = 2; float r24785 = sqrt(r24784); float r24786 = 4; float r24787 = r24785 / r24786; float r24788 = log(r24787); float r24789 = 1; float r24790 = 3; float r24791 = v; float r24792 = r24791 * r24791; float r24793 = r24790 * r24792; float r24794 = r24789 - r24793; float r24795 = sqrt(r24794); float r24796 = log(r24795); float r24797 = r24788 + r24796; float r24798 = exp(r24797); float r24799 = r24789 - r24792; float r24800 = r24798 * r24799; return r24800; } double f_od(double v) { double r24801 = 2; double r24802 = sqrt(r24801); double r24803 = 4; double r24804 = r24802 / r24803; double r24805 = log(r24804); double r24806 = 1; double r24807 = 3; double r24808 = v; double r24809 = r24808 * r24808; double r24810 = r24807 * r24809; double r24811 = r24806 - r24810; double r24812 = sqrt(r24811); double r24813 = log(r24812); double r24814 = r24805 + r24813; double r24815 = exp(r24814); double r24816 = r24806 - r24809; double r24817 = r24815 * r24816; return r24817; } 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 r24818, r24819, r24820, r24821, r24822, r24823, r24824, r24825, r24826, r24827, r24828, r24829, r24830, r24831; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r24818, "2", 10, MPFR_RNDN); mpfr_init(r24819); mpfr_init_set_str(r24820, "4", 10, MPFR_RNDN); mpfr_init(r24821); mpfr_init_set_str(r24822, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24823, "3", 10, MPFR_RNDN); mpfr_init(r24824); mpfr_init(r24825); mpfr_init(r24826); mpfr_init(r24827); mpfr_init(r24828); mpfr_init(r24829); mpfr_init(r24830); mpfr_init(r24831); } double f_im(double v) { ; mpfr_sqrt(r24819, r24818, MPFR_RNDN); ; mpfr_div(r24821, r24819, r24820, MPFR_RNDN); ; ; mpfr_set_d(r24824, v, MPFR_RNDN); mpfr_mul(r24825, r24824, r24824, MPFR_RNDN); mpfr_mul(r24826, r24823, r24825, MPFR_RNDN); mpfr_sub(r24827, r24822, r24826, MPFR_RNDN); mpfr_sqrt(r24828, r24827, MPFR_RNDN); mpfr_mul(r24829, r24821, r24828, MPFR_RNDN); mpfr_sub(r24830, r24822, r24825, MPFR_RNDN); mpfr_mul(r24831, r24829, r24830, MPFR_RNDN); return mpfr_get_d(r24831, MPFR_RNDN); } static mpfr_t r24832, r24833, r24834, r24835, r24836, r24837, r24838, r24839, r24840, r24841, r24842, r24843, r24844, r24845, r24846, r24847, r24848; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r24832, "2", 10, MPFR_RNDN); mpfr_init(r24833); mpfr_init_set_str(r24834, "4", 10, MPFR_RNDN); mpfr_init(r24835); mpfr_init(r24836); mpfr_init_set_str(r24837, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24838, "3", 10, MPFR_RNDN); mpfr_init(r24839); mpfr_init(r24840); mpfr_init(r24841); mpfr_init(r24842); mpfr_init(r24843); mpfr_init(r24844); mpfr_init(r24845); mpfr_init(r24846); mpfr_init(r24847); mpfr_init(r24848); } double f_fm(double v) { ; mpfr_sqrt(r24833, r24832, MPFR_RNDN); ; mpfr_div(r24835, r24833, r24834, MPFR_RNDN); mpfr_log(r24836, r24835, MPFR_RNDN); ; ; mpfr_set_d(r24839, v, MPFR_RNDN); mpfr_mul(r24840, r24839, r24839, MPFR_RNDN); mpfr_mul(r24841, r24838, r24840, MPFR_RNDN); mpfr_sub(r24842, r24837, r24841, MPFR_RNDN); mpfr_sqrt(r24843, r24842, MPFR_RNDN); mpfr_log(r24844, r24843, MPFR_RNDN); mpfr_add(r24845, r24836, r24844, MPFR_RNDN); mpfr_exp(r24846, r24845, MPFR_RNDN); mpfr_sub(r24847, r24837, r24840, MPFR_RNDN); mpfr_mul(r24848, r24846, r24847, MPFR_RNDN); return mpfr_get_d(r24848, MPFR_RNDN); } static mpfr_t r24849, r24850, r24851, r24852, r24853, r24854, r24855, r24856, r24857, r24858, r24859, r24860, r24861, r24862, r24863, r24864, r24865; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r24849, "2", 10, MPFR_RNDN); mpfr_init(r24850); mpfr_init_set_str(r24851, "4", 10, MPFR_RNDN); mpfr_init(r24852); mpfr_init(r24853); mpfr_init_set_str(r24854, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24855, "3", 10, MPFR_RNDN); 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); mpfr_init(r24865); } double f_dm(double v) { ; mpfr_sqrt(r24850, r24849, MPFR_RNDN); ; mpfr_div(r24852, r24850, r24851, MPFR_RNDN); mpfr_log(r24853, r24852, MPFR_RNDN); ; ; mpfr_set_d(r24856, v, MPFR_RNDN); mpfr_mul(r24857, r24856, r24856, MPFR_RNDN); mpfr_mul(r24858, r24855, r24857, MPFR_RNDN); mpfr_sub(r24859, r24854, r24858, MPFR_RNDN); mpfr_sqrt(r24860, r24859, MPFR_RNDN); mpfr_log(r24861, r24860, MPFR_RNDN); mpfr_add(r24862, r24853, r24861, MPFR_RNDN); mpfr_exp(r24863, r24862, MPFR_RNDN); mpfr_sub(r24864, r24854, r24857, MPFR_RNDN); mpfr_mul(r24865, r24863, r24864, MPFR_RNDN); return mpfr_get_d(r24865, MPFR_RNDN); }