#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r23804 = 2; float r23805 = sqrt(r23804); float r23806 = 4; float r23807 = r23805 / r23806; float r23808 = 1; float r23809 = 3; float r23810 = v; float r23811 = r23810 * r23810; float r23812 = r23809 * r23811; float r23813 = r23808 - r23812; float r23814 = sqrt(r23813); float r23815 = r23807 * r23814; float r23816 = r23808 - r23811; float r23817 = r23815 * r23816; return r23817; } double f_id(double v) { double r23818 = 2; double r23819 = sqrt(r23818); double r23820 = 4; double r23821 = r23819 / r23820; double r23822 = 1; double r23823 = 3; double r23824 = v; double r23825 = r23824 * r23824; double r23826 = r23823 * r23825; double r23827 = r23822 - r23826; double r23828 = sqrt(r23827); double r23829 = r23821 * r23828; double r23830 = r23822 - r23825; double r23831 = r23829 * r23830; return r23831; } double f_of(float v) { float r23832 = 2; float r23833 = sqrt(r23832); float r23834 = 4; float r23835 = r23833 / r23834; float r23836 = 1; float r23837 = 3; float r23838 = v; float r23839 = r23838 * r23838; float r23840 = r23837 * r23839; float r23841 = r23836 - r23840; float r23842 = sqrt(r23841); float r23843 = r23835 * r23842; float r23844 = r23836 - r23839; float r23845 = r23843 * r23844; float r23846 = expm1(r23845); float r23847 = log1p(r23846); return r23847; } double f_od(double v) { double r23848 = 2; double r23849 = sqrt(r23848); double r23850 = 4; double r23851 = r23849 / r23850; double r23852 = 1; double r23853 = 3; double r23854 = v; double r23855 = r23854 * r23854; double r23856 = r23853 * r23855; double r23857 = r23852 - r23856; double r23858 = sqrt(r23857); double r23859 = r23851 * r23858; double r23860 = r23852 - r23855; double r23861 = r23859 * r23860; double r23862 = expm1(r23861); double r23863 = log1p(r23862); return r23863; } 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 r23864, r23865, r23866, r23867, r23868, r23869, r23870, r23871, r23872, r23873, r23874, r23875, r23876, r23877; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r23864, "2", 10, MPFR_RNDN); mpfr_init(r23865); mpfr_init_set_str(r23866, "4", 10, MPFR_RNDN); mpfr_init(r23867); mpfr_init_set_str(r23868, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23869, "3", 10, MPFR_RNDN); mpfr_init(r23870); mpfr_init(r23871); mpfr_init(r23872); mpfr_init(r23873); mpfr_init(r23874); mpfr_init(r23875); mpfr_init(r23876); mpfr_init(r23877); } double f_im(double v) { ; mpfr_sqrt(r23865, r23864, MPFR_RNDN); ; mpfr_div(r23867, r23865, r23866, MPFR_RNDN); ; ; mpfr_set_d(r23870, v, MPFR_RNDN); mpfr_mul(r23871, r23870, r23870, MPFR_RNDN); mpfr_mul(r23872, r23869, r23871, MPFR_RNDN); mpfr_sub(r23873, r23868, r23872, MPFR_RNDN); mpfr_sqrt(r23874, r23873, MPFR_RNDN); mpfr_mul(r23875, r23867, r23874, MPFR_RNDN); mpfr_sub(r23876, r23868, r23871, MPFR_RNDN); mpfr_mul(r23877, r23875, r23876, MPFR_RNDN); return mpfr_get_d(r23877, MPFR_RNDN); } static mpfr_t r23878, r23879, r23880, r23881, r23882, r23883, r23884, r23885, r23886, r23887, r23888, r23889, r23890, r23891, r23892, r23893; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23878, "2", 10, MPFR_RNDN); mpfr_init(r23879); mpfr_init_set_str(r23880, "4", 10, MPFR_RNDN); mpfr_init(r23881); mpfr_init_set_str(r23882, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23883, "3", 10, MPFR_RNDN); mpfr_init(r23884); mpfr_init(r23885); mpfr_init(r23886); mpfr_init(r23887); mpfr_init(r23888); mpfr_init(r23889); mpfr_init(r23890); mpfr_init(r23891); mpfr_init(r23892); mpfr_init(r23893); } double f_fm(double v) { ; mpfr_sqrt(r23879, r23878, MPFR_RNDN); ; mpfr_div(r23881, r23879, r23880, MPFR_RNDN); ; ; mpfr_set_d(r23884, v, MPFR_RNDN); mpfr_mul(r23885, r23884, r23884, MPFR_RNDN); mpfr_mul(r23886, r23883, r23885, MPFR_RNDN); mpfr_sub(r23887, r23882, r23886, MPFR_RNDN); mpfr_sqrt(r23888, r23887, MPFR_RNDN); mpfr_mul(r23889, r23881, r23888, MPFR_RNDN); mpfr_sub(r23890, r23882, r23885, MPFR_RNDN); mpfr_mul(r23891, r23889, r23890, MPFR_RNDN); mpfr_expm1(r23892, r23891, MPFR_RNDN); mpfr_log1p(r23893, r23892, MPFR_RNDN); return mpfr_get_d(r23893, MPFR_RNDN); } static mpfr_t r23894, r23895, r23896, r23897, r23898, r23899, r23900, r23901, r23902, r23903, r23904, r23905, r23906, r23907, r23908, r23909; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23894, "2", 10, MPFR_RNDN); mpfr_init(r23895); mpfr_init_set_str(r23896, "4", 10, MPFR_RNDN); mpfr_init(r23897); mpfr_init_set_str(r23898, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23899, "3", 10, MPFR_RNDN); mpfr_init(r23900); mpfr_init(r23901); mpfr_init(r23902); mpfr_init(r23903); mpfr_init(r23904); mpfr_init(r23905); mpfr_init(r23906); mpfr_init(r23907); mpfr_init(r23908); mpfr_init(r23909); } double f_dm(double v) { ; mpfr_sqrt(r23895, r23894, MPFR_RNDN); ; mpfr_div(r23897, r23895, r23896, MPFR_RNDN); ; ; mpfr_set_d(r23900, v, MPFR_RNDN); mpfr_mul(r23901, r23900, r23900, MPFR_RNDN); mpfr_mul(r23902, r23899, r23901, MPFR_RNDN); mpfr_sub(r23903, r23898, r23902, MPFR_RNDN); mpfr_sqrt(r23904, r23903, MPFR_RNDN); mpfr_mul(r23905, r23897, r23904, MPFR_RNDN); mpfr_sub(r23906, r23898, r23901, MPFR_RNDN); mpfr_mul(r23907, r23905, r23906, MPFR_RNDN); mpfr_expm1(r23908, r23907, MPFR_RNDN); mpfr_log1p(r23909, r23908, MPFR_RNDN); return mpfr_get_d(r23909, MPFR_RNDN); }