#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r24703 = 2; float r24704 = sqrt(r24703); float r24705 = 4; float r24706 = r24704 / r24705; float r24707 = 1; float r24708 = 3; float r24709 = v; float r24710 = r24709 * r24709; float r24711 = r24708 * r24710; float r24712 = r24707 - r24711; float r24713 = sqrt(r24712); float r24714 = r24706 * r24713; float r24715 = r24707 - r24710; float r24716 = r24714 * r24715; return r24716; } double f_id(double v) { double r24717 = 2; double r24718 = sqrt(r24717); double r24719 = 4; double r24720 = r24718 / r24719; double r24721 = 1; double r24722 = 3; double r24723 = v; double r24724 = r24723 * r24723; double r24725 = r24722 * r24724; double r24726 = r24721 - r24725; double r24727 = sqrt(r24726); double r24728 = r24720 * r24727; double r24729 = r24721 - r24724; double r24730 = r24728 * r24729; return r24730; } double f_of(float v) { float r24731 = 2; float r24732 = sqrt(r24731); float r24733 = 4; float r24734 = r24732 / r24733; float r24735 = 1; float r24736 = 3; float r24737 = v; float r24738 = r24737 * r24737; float r24739 = r24736 * r24738; float r24740 = r24735 - r24739; float r24741 = sqrt(r24740); float r24742 = r24734 * r24741; float r24743 = log(r24742); float r24744 = r24735 - r24738; float r24745 = log(r24744); float r24746 = r24743 + r24745; float r24747 = exp(r24746); return r24747; } double f_od(double v) { double r24748 = 2; double r24749 = sqrt(r24748); double r24750 = 4; double r24751 = r24749 / r24750; double r24752 = 1; double r24753 = 3; double r24754 = v; double r24755 = r24754 * r24754; double r24756 = r24753 * r24755; double r24757 = r24752 - r24756; double r24758 = sqrt(r24757); double r24759 = r24751 * r24758; double r24760 = log(r24759); double r24761 = r24752 - r24755; double r24762 = log(r24761); double r24763 = r24760 + r24762; double r24764 = exp(r24763); return r24764; } 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 r24765, r24766, r24767, r24768, r24769, r24770, r24771, r24772, r24773, r24774, r24775, r24776, r24777, r24778; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r24765, "2", 10, MPFR_RNDN); mpfr_init(r24766); mpfr_init_set_str(r24767, "4", 10, MPFR_RNDN); mpfr_init(r24768); mpfr_init_set_str(r24769, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24770, "3", 10, MPFR_RNDN); mpfr_init(r24771); mpfr_init(r24772); mpfr_init(r24773); mpfr_init(r24774); mpfr_init(r24775); mpfr_init(r24776); mpfr_init(r24777); mpfr_init(r24778); } double f_im(double v) { ; mpfr_sqrt(r24766, r24765, MPFR_RNDN); ; mpfr_div(r24768, r24766, r24767, MPFR_RNDN); ; ; mpfr_set_d(r24771, v, MPFR_RNDN); mpfr_mul(r24772, r24771, r24771, MPFR_RNDN); mpfr_mul(r24773, r24770, r24772, MPFR_RNDN); mpfr_sub(r24774, r24769, r24773, MPFR_RNDN); mpfr_sqrt(r24775, r24774, MPFR_RNDN); mpfr_mul(r24776, r24768, r24775, MPFR_RNDN); mpfr_sub(r24777, r24769, r24772, MPFR_RNDN); mpfr_mul(r24778, r24776, r24777, MPFR_RNDN); return mpfr_get_d(r24778, MPFR_RNDN); } static mpfr_t r24779, r24780, r24781, r24782, r24783, r24784, r24785, r24786, r24787, r24788, r24789, r24790, r24791, r24792, r24793, r24794, r24795; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r24779, "2", 10, MPFR_RNDN); mpfr_init(r24780); mpfr_init_set_str(r24781, "4", 10, MPFR_RNDN); mpfr_init(r24782); mpfr_init_set_str(r24783, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24784, "3", 10, MPFR_RNDN); mpfr_init(r24785); mpfr_init(r24786); mpfr_init(r24787); mpfr_init(r24788); mpfr_init(r24789); mpfr_init(r24790); mpfr_init(r24791); mpfr_init(r24792); mpfr_init(r24793); mpfr_init(r24794); mpfr_init(r24795); } double f_fm(double v) { ; mpfr_sqrt(r24780, r24779, MPFR_RNDN); ; mpfr_div(r24782, r24780, r24781, MPFR_RNDN); ; ; mpfr_set_d(r24785, v, MPFR_RNDN); mpfr_mul(r24786, r24785, r24785, MPFR_RNDN); mpfr_mul(r24787, r24784, r24786, MPFR_RNDN); mpfr_sub(r24788, r24783, r24787, MPFR_RNDN); mpfr_sqrt(r24789, r24788, MPFR_RNDN); mpfr_mul(r24790, r24782, r24789, MPFR_RNDN); mpfr_log(r24791, r24790, MPFR_RNDN); mpfr_sub(r24792, r24783, r24786, MPFR_RNDN); mpfr_log(r24793, r24792, MPFR_RNDN); mpfr_add(r24794, r24791, r24793, MPFR_RNDN); mpfr_exp(r24795, r24794, MPFR_RNDN); return mpfr_get_d(r24795, MPFR_RNDN); } static mpfr_t r24796, r24797, r24798, r24799, r24800, r24801, r24802, r24803, r24804, r24805, r24806, r24807, r24808, r24809, r24810, r24811, r24812; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r24796, "2", 10, MPFR_RNDN); mpfr_init(r24797); mpfr_init_set_str(r24798, "4", 10, MPFR_RNDN); mpfr_init(r24799); mpfr_init_set_str(r24800, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24801, "3", 10, MPFR_RNDN); mpfr_init(r24802); mpfr_init(r24803); mpfr_init(r24804); mpfr_init(r24805); mpfr_init(r24806); mpfr_init(r24807); mpfr_init(r24808); mpfr_init(r24809); mpfr_init(r24810); mpfr_init(r24811); mpfr_init(r24812); } double f_dm(double v) { ; mpfr_sqrt(r24797, r24796, MPFR_RNDN); ; mpfr_div(r24799, r24797, r24798, MPFR_RNDN); ; ; mpfr_set_d(r24802, v, MPFR_RNDN); mpfr_mul(r24803, r24802, r24802, MPFR_RNDN); mpfr_mul(r24804, r24801, r24803, MPFR_RNDN); mpfr_sub(r24805, r24800, r24804, MPFR_RNDN); mpfr_sqrt(r24806, r24805, MPFR_RNDN); mpfr_mul(r24807, r24799, r24806, MPFR_RNDN); mpfr_log(r24808, r24807, MPFR_RNDN); mpfr_sub(r24809, r24800, r24803, MPFR_RNDN); mpfr_log(r24810, r24809, MPFR_RNDN); mpfr_add(r24811, r24808, r24810, MPFR_RNDN); mpfr_exp(r24812, r24811, MPFR_RNDN); return mpfr_get_d(r24812, MPFR_RNDN); }