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