#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix A"; double f_if(float a, float k, float m) { float r27213 = a; float r27214 = k; float r27215 = m; float r27216 = pow(r27214, r27215); float r27217 = r27213 * r27216; float r27218 = 1; float r27219 = 10; float r27220 = r27219 * r27214; float r27221 = r27218 + r27220; float r27222 = r27214 * r27214; float r27223 = r27221 + r27222; float r27224 = r27217 / r27223; return r27224; } double f_id(double a, double k, double m) { double r27225 = a; double r27226 = k; double r27227 = m; double r27228 = pow(r27226, r27227); double r27229 = r27225 * r27228; double r27230 = 1; double r27231 = 10; double r27232 = r27231 * r27226; double r27233 = r27230 + r27232; double r27234 = r27226 * r27226; double r27235 = r27233 + r27234; double r27236 = r27229 / r27235; return r27236; } double f_of(float a, float k, float m) { float r27237 = k; float r27238 = m; float r27239 = pow(r27237, r27238); float r27240 = a; float r27241 = r27239 * r27240; float r27242 = 10; float r27243 = r27242 + r27237; float r27244 = 1; float r27245 = fma(r27237, r27243, r27244); float r27246 = sqrt(r27245); float r27247 = r27241 / r27246; float r27248 = r27242 * r27237; float r27249 = r27244 + r27248; float r27250 = r27237 * r27237; float r27251 = r27249 + r27250; float r27252 = sqrt(r27251); float r27253 = r27247 / r27252; return r27253; } double f_od(double a, double k, double m) { double r27254 = k; double r27255 = m; double r27256 = pow(r27254, r27255); double r27257 = a; double r27258 = r27256 * r27257; double r27259 = 10; double r27260 = r27259 + r27254; double r27261 = 1; double r27262 = fma(r27254, r27260, r27261); double r27263 = sqrt(r27262); double r27264 = r27258 / r27263; double r27265 = r27259 * r27254; double r27266 = r27261 + r27265; double r27267 = r27254 * r27254; double r27268 = r27266 + r27267; double r27269 = sqrt(r27268); double r27270 = r27264 / r27269; return r27270; } 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 r27271, r27272, r27273, r27274, r27275, r27276, r27277, r27278, r27279, r27280, r27281, r27282; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r27271); mpfr_init(r27272); mpfr_init(r27273); mpfr_init(r27274); mpfr_init(r27275); mpfr_init_set_str(r27276, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27277, "10", 10, MPFR_RNDN); mpfr_init(r27278); mpfr_init(r27279); mpfr_init(r27280); mpfr_init(r27281); mpfr_init(r27282); } double f_im(double a, double k, double m) { mpfr_set_d(r27271, a, MPFR_RNDN); mpfr_set_d(r27272, k, MPFR_RNDN); mpfr_set_d(r27273, m, MPFR_RNDN); mpfr_pow(r27274, r27272, r27273, MPFR_RNDN); mpfr_mul(r27275, r27271, r27274, MPFR_RNDN); ; ; mpfr_mul(r27278, r27277, r27272, MPFR_RNDN); mpfr_add(r27279, r27276, r27278, MPFR_RNDN); mpfr_mul(r27280, r27272, r27272, MPFR_RNDN); mpfr_add(r27281, r27279, r27280, MPFR_RNDN); mpfr_div(r27282, r27275, r27281, MPFR_RNDN); return mpfr_get_d(r27282, MPFR_RNDN); } static mpfr_t r27283, r27284, r27285, r27286, r27287, r27288, r27289, r27290, r27291, r27292, r27293, r27294, r27295, r27296, r27297, r27298, r27299; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r27283); mpfr_init(r27284); mpfr_init(r27285); mpfr_init(r27286); mpfr_init(r27287); mpfr_init_set_str(r27288, "10", 10, MPFR_RNDN); mpfr_init(r27289); mpfr_init_set_str(r27290, "1", 10, MPFR_RNDN); mpfr_init(r27291); mpfr_init(r27292); mpfr_init(r27293); mpfr_init(r27294); mpfr_init(r27295); mpfr_init(r27296); mpfr_init(r27297); mpfr_init(r27298); mpfr_init(r27299); } double f_fm(double a, double k, double m) { mpfr_set_d(r27283, k, MPFR_RNDN); mpfr_set_d(r27284, m, MPFR_RNDN); mpfr_pow(r27285, r27283, r27284, MPFR_RNDN); mpfr_set_d(r27286, a, MPFR_RNDN); mpfr_mul(r27287, r27285, r27286, MPFR_RNDN); ; mpfr_add(r27289, r27288, r27283, MPFR_RNDN); ; mpfr_fma(r27291, r27283, r27289, r27290, MPFR_RNDN); mpfr_sqrt(r27292, r27291, MPFR_RNDN); mpfr_div(r27293, r27287, r27292, MPFR_RNDN); mpfr_mul(r27294, r27288, r27283, MPFR_RNDN); mpfr_add(r27295, r27290, r27294, MPFR_RNDN); mpfr_mul(r27296, r27283, r27283, MPFR_RNDN); mpfr_add(r27297, r27295, r27296, MPFR_RNDN); mpfr_sqrt(r27298, r27297, MPFR_RNDN); mpfr_div(r27299, r27293, r27298, MPFR_RNDN); return mpfr_get_d(r27299, MPFR_RNDN); } static mpfr_t r27300, r27301, r27302, r27303, r27304, r27305, r27306, r27307, r27308, r27309, r27310, r27311, r27312, r27313, r27314, r27315, r27316; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r27300); mpfr_init(r27301); mpfr_init(r27302); mpfr_init(r27303); mpfr_init(r27304); mpfr_init_set_str(r27305, "10", 10, MPFR_RNDN); mpfr_init(r27306); mpfr_init_set_str(r27307, "1", 10, MPFR_RNDN); mpfr_init(r27308); mpfr_init(r27309); mpfr_init(r27310); mpfr_init(r27311); mpfr_init(r27312); mpfr_init(r27313); mpfr_init(r27314); mpfr_init(r27315); mpfr_init(r27316); } double f_dm(double a, double k, double m) { mpfr_set_d(r27300, k, MPFR_RNDN); mpfr_set_d(r27301, m, MPFR_RNDN); mpfr_pow(r27302, r27300, r27301, MPFR_RNDN); mpfr_set_d(r27303, a, MPFR_RNDN); mpfr_mul(r27304, r27302, r27303, MPFR_RNDN); ; mpfr_add(r27306, r27305, r27300, MPFR_RNDN); ; mpfr_fma(r27308, r27300, r27306, r27307, MPFR_RNDN); mpfr_sqrt(r27309, r27308, MPFR_RNDN); mpfr_div(r27310, r27304, r27309, MPFR_RNDN); mpfr_mul(r27311, r27305, r27300, MPFR_RNDN); mpfr_add(r27312, r27307, r27311, MPFR_RNDN); mpfr_mul(r27313, r27300, r27300, MPFR_RNDN); mpfr_add(r27314, r27312, r27313, MPFR_RNDN); mpfr_sqrt(r27315, r27314, MPFR_RNDN); mpfr_div(r27316, r27310, r27315, MPFR_RNDN); return mpfr_get_d(r27316, MPFR_RNDN); }