#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix A"; double f_if(float a, float k, float m) { float r27177 = a; float r27178 = k; float r27179 = m; float r27180 = pow(r27178, r27179); float r27181 = r27177 * r27180; float r27182 = 1; float r27183 = 10; float r27184 = r27183 * r27178; float r27185 = r27182 + r27184; float r27186 = r27178 * r27178; float r27187 = r27185 + r27186; float r27188 = r27181 / r27187; return r27188; } double f_id(double a, double k, double m) { double r27189 = a; double r27190 = k; double r27191 = m; double r27192 = pow(r27190, r27191); double r27193 = r27189 * r27192; double r27194 = 1; double r27195 = 10; double r27196 = r27195 * r27190; double r27197 = r27194 + r27196; double r27198 = r27190 * r27190; double r27199 = r27197 + r27198; double r27200 = r27193 / r27199; return r27200; } double f_of(float a, float k, float m) { float r27201 = a; float r27202 = k; float r27203 = 10; float r27204 = r27202 + r27203; float r27205 = r27202 * r27204; float r27206 = 1; float r27207 = r27205 + r27206; float r27208 = sqrt(r27207); float r27209 = r27201 / r27208; float r27210 = m; float r27211 = pow(r27202, r27210); float r27212 = r27203 + r27202; float r27213 = r27202 * r27212; float r27214 = r27206 + r27213; float r27215 = sqrt(r27214); float r27216 = r27211 / r27215; float r27217 = r27209 * r27216; return r27217; } double f_od(double a, double k, double m) { double r27218 = a; double r27219 = k; double r27220 = 10; double r27221 = r27219 + r27220; double r27222 = r27219 * r27221; double r27223 = 1; double r27224 = r27222 + r27223; double r27225 = sqrt(r27224); double r27226 = r27218 / r27225; double r27227 = m; double r27228 = pow(r27219, r27227); double r27229 = r27220 + r27219; double r27230 = r27219 * r27229; double r27231 = r27223 + r27230; double r27232 = sqrt(r27231); double r27233 = r27228 / r27232; double r27234 = r27226 * r27233; return r27234; } 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 r27235, r27236, r27237, r27238, r27239, r27240, r27241, r27242, r27243, r27244, r27245, r27246; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r27235); mpfr_init(r27236); mpfr_init(r27237); mpfr_init(r27238); mpfr_init(r27239); mpfr_init_set_str(r27240, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27241, "10", 10, MPFR_RNDN); mpfr_init(r27242); mpfr_init(r27243); mpfr_init(r27244); mpfr_init(r27245); mpfr_init(r27246); } double f_im(double a, double k, double m) { mpfr_set_d(r27235, a, MPFR_RNDN); mpfr_set_d(r27236, k, MPFR_RNDN); mpfr_set_d(r27237, m, MPFR_RNDN); mpfr_pow(r27238, r27236, r27237, MPFR_RNDN); mpfr_mul(r27239, r27235, r27238, MPFR_RNDN); ; ; mpfr_mul(r27242, r27241, r27236, MPFR_RNDN); mpfr_add(r27243, r27240, r27242, MPFR_RNDN); mpfr_mul(r27244, r27236, r27236, MPFR_RNDN); mpfr_add(r27245, r27243, r27244, MPFR_RNDN); mpfr_div(r27246, r27239, r27245, MPFR_RNDN); return mpfr_get_d(r27246, MPFR_RNDN); } static mpfr_t r27247, r27248, r27249, r27250, r27251, r27252, r27253, r27254, r27255, r27256, r27257, r27258, r27259, r27260, r27261, r27262, r27263; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r27247); mpfr_init(r27248); mpfr_init_set_str(r27249, "10", 10, MPFR_RNDN); mpfr_init(r27250); mpfr_init(r27251); mpfr_init_set_str(r27252, "1", 10, MPFR_RNDN); mpfr_init(r27253); mpfr_init(r27254); mpfr_init(r27255); mpfr_init(r27256); mpfr_init(r27257); mpfr_init(r27258); mpfr_init(r27259); mpfr_init(r27260); mpfr_init(r27261); mpfr_init(r27262); mpfr_init(r27263); } double f_fm(double a, double k, double m) { mpfr_set_d(r27247, a, MPFR_RNDN); mpfr_set_d(r27248, k, MPFR_RNDN); ; mpfr_add(r27250, r27248, r27249, MPFR_RNDN); mpfr_mul(r27251, r27248, r27250, MPFR_RNDN); ; mpfr_add(r27253, r27251, r27252, MPFR_RNDN); mpfr_sqrt(r27254, r27253, MPFR_RNDN); mpfr_div(r27255, r27247, r27254, MPFR_RNDN); mpfr_set_d(r27256, m, MPFR_RNDN); mpfr_pow(r27257, r27248, r27256, MPFR_RNDN); mpfr_add(r27258, r27249, r27248, MPFR_RNDN); mpfr_mul(r27259, r27248, r27258, MPFR_RNDN); mpfr_add(r27260, r27252, r27259, MPFR_RNDN); mpfr_sqrt(r27261, r27260, MPFR_RNDN); mpfr_div(r27262, r27257, r27261, MPFR_RNDN); mpfr_mul(r27263, r27255, r27262, MPFR_RNDN); return mpfr_get_d(r27263, MPFR_RNDN); } static mpfr_t r27264, r27265, r27266, r27267, r27268, r27269, r27270, r27271, r27272, r27273, r27274, r27275, r27276, r27277, r27278, r27279, r27280; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r27264); mpfr_init(r27265); mpfr_init_set_str(r27266, "10", 10, MPFR_RNDN); mpfr_init(r27267); mpfr_init(r27268); mpfr_init_set_str(r27269, "1", 10, MPFR_RNDN); mpfr_init(r27270); mpfr_init(r27271); mpfr_init(r27272); mpfr_init(r27273); mpfr_init(r27274); mpfr_init(r27275); mpfr_init(r27276); mpfr_init(r27277); mpfr_init(r27278); mpfr_init(r27279); mpfr_init(r27280); } double f_dm(double a, double k, double m) { mpfr_set_d(r27264, a, MPFR_RNDN); mpfr_set_d(r27265, k, MPFR_RNDN); ; mpfr_add(r27267, r27265, r27266, MPFR_RNDN); mpfr_mul(r27268, r27265, r27267, MPFR_RNDN); ; mpfr_add(r27270, r27268, r27269, MPFR_RNDN); mpfr_sqrt(r27271, r27270, MPFR_RNDN); mpfr_div(r27272, r27264, r27271, MPFR_RNDN); mpfr_set_d(r27273, m, MPFR_RNDN); mpfr_pow(r27274, r27265, r27273, MPFR_RNDN); mpfr_add(r27275, r27266, r27265, MPFR_RNDN); mpfr_mul(r27276, r27265, r27275, MPFR_RNDN); mpfr_add(r27277, r27269, r27276, MPFR_RNDN); mpfr_sqrt(r27278, r27277, MPFR_RNDN); mpfr_div(r27279, r27274, r27278, MPFR_RNDN); mpfr_mul(r27280, r27272, r27279, MPFR_RNDN); return mpfr_get_d(r27280, MPFR_RNDN); }