#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r24043 = 4; float r24044 = 3; float r24045 = atan2(1.0, 0.0); float r24046 = r24044 * r24045; float r24047 = 1; float r24048 = v; float r24049 = r24048 * r24048; float r24050 = r24047 - r24049; float r24051 = r24046 * r24050; float r24052 = 2; float r24053 = 6; float r24054 = r24053 * r24049; float r24055 = r24052 - r24054; float r24056 = sqrt(r24055); float r24057 = r24051 * r24056; float r24058 = r24043 / r24057; return r24058; } double f_id(double v) { double r24059 = 4; double r24060 = 3; double r24061 = atan2(1.0, 0.0); double r24062 = r24060 * r24061; double r24063 = 1; double r24064 = v; double r24065 = r24064 * r24064; double r24066 = r24063 - r24065; double r24067 = r24062 * r24066; double r24068 = 2; double r24069 = 6; double r24070 = r24069 * r24065; double r24071 = r24068 - r24070; double r24072 = sqrt(r24071); double r24073 = r24067 * r24072; double r24074 = r24059 / r24073; return r24074; } double f_of(float v) { float r24075 = 4; float r24076 = atan2(1.0, 0.0); float r24077 = 3; float r24078 = r24076 * r24077; float r24079 = r24075 / r24078; float r24080 = 1; float r24081 = v; float r24082 = r24081 * r24081; float r24083 = r24080 - r24082; float r24084 = r24079 / r24083; float r24085 = 6; float r24086 = -r24085; float r24087 = 2; float r24088 = fma(r24086, r24082, r24087); float r24089 = sqrt(r24088); float r24090 = r24084 / r24089; float r24091 = pow(r24090, r24077); float r24092 = cbrt(r24091); return r24092; } double f_od(double v) { double r24093 = 4; double r24094 = atan2(1.0, 0.0); double r24095 = 3; double r24096 = r24094 * r24095; double r24097 = r24093 / r24096; double r24098 = 1; double r24099 = v; double r24100 = r24099 * r24099; double r24101 = r24098 - r24100; double r24102 = r24097 / r24101; double r24103 = 6; double r24104 = -r24103; double r24105 = 2; double r24106 = fma(r24104, r24100, r24105); double r24107 = sqrt(r24106); double r24108 = r24102 / r24107; double r24109 = pow(r24108, r24095); double r24110 = cbrt(r24109); return r24110; } 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 r24111, r24112, r24113, r24114, r24115, r24116, r24117, r24118, r24119, r24120, r24121, r24122, r24123, r24124, r24125, r24126; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r24111, "4", 10, MPFR_RNDN); mpfr_init_set_str(r24112, "3", 10, MPFR_RNDN); mpfr_init(r24113); mpfr_init(r24114); mpfr_init_set_str(r24115, "1", 10, MPFR_RNDN); mpfr_init(r24116); mpfr_init(r24117); mpfr_init(r24118); mpfr_init(r24119); mpfr_init_set_str(r24120, "2", 10, MPFR_RNDN); mpfr_init_set_str(r24121, "6", 10, MPFR_RNDN); mpfr_init(r24122); mpfr_init(r24123); mpfr_init(r24124); mpfr_init(r24125); mpfr_init(r24126); } double f_im(double v) { ; ; mpfr_const_pi(r24113, MPFR_RNDN); mpfr_mul(r24114, r24112, r24113, MPFR_RNDN); ; mpfr_set_d(r24116, v, MPFR_RNDN); mpfr_mul(r24117, r24116, r24116, MPFR_RNDN); mpfr_sub(r24118, r24115, r24117, MPFR_RNDN); mpfr_mul(r24119, r24114, r24118, MPFR_RNDN); ; ; mpfr_mul(r24122, r24121, r24117, MPFR_RNDN); mpfr_sub(r24123, r24120, r24122, MPFR_RNDN); mpfr_sqrt(r24124, r24123, MPFR_RNDN); mpfr_mul(r24125, r24119, r24124, MPFR_RNDN); mpfr_div(r24126, r24111, r24125, MPFR_RNDN); return mpfr_get_d(r24126, MPFR_RNDN); } static mpfr_t r24127, r24128, r24129, r24130, r24131, r24132, r24133, r24134, r24135, r24136, r24137, r24138, r24139, r24140, r24141, r24142, r24143, r24144; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r24127, "4", 10, MPFR_RNDN); mpfr_init(r24128); mpfr_init_set_str(r24129, "3", 10, MPFR_RNDN); mpfr_init(r24130); mpfr_init(r24131); mpfr_init_set_str(r24132, "1", 10, MPFR_RNDN); mpfr_init(r24133); mpfr_init(r24134); mpfr_init(r24135); mpfr_init(r24136); mpfr_init_set_str(r24137, "6", 10, MPFR_RNDN); mpfr_init(r24138); mpfr_init_set_str(r24139, "2", 10, MPFR_RNDN); mpfr_init(r24140); mpfr_init(r24141); mpfr_init(r24142); mpfr_init(r24143); mpfr_init(r24144); } double f_fm(double v) { ; mpfr_const_pi(r24128, MPFR_RNDN); ; mpfr_mul(r24130, r24128, r24129, MPFR_RNDN); mpfr_div(r24131, r24127, r24130, MPFR_RNDN); ; mpfr_set_d(r24133, v, MPFR_RNDN); mpfr_mul(r24134, r24133, r24133, MPFR_RNDN); mpfr_sub(r24135, r24132, r24134, MPFR_RNDN); mpfr_div(r24136, r24131, r24135, MPFR_RNDN); ; mpfr_neg(r24138, r24137, MPFR_RNDN); ; mpfr_fma(r24140, r24138, r24134, r24139, MPFR_RNDN); mpfr_sqrt(r24141, r24140, MPFR_RNDN); mpfr_div(r24142, r24136, r24141, MPFR_RNDN); mpfr_pow(r24143, r24142, r24129, MPFR_RNDN); mpfr_cbrt(r24144, r24143, MPFR_RNDN); return mpfr_get_d(r24144, MPFR_RNDN); } static mpfr_t r24145, r24146, r24147, r24148, r24149, r24150, r24151, r24152, r24153, r24154, r24155, r24156, r24157, r24158, r24159, r24160, r24161, r24162; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r24145, "4", 10, MPFR_RNDN); mpfr_init(r24146); mpfr_init_set_str(r24147, "3", 10, MPFR_RNDN); mpfr_init(r24148); mpfr_init(r24149); mpfr_init_set_str(r24150, "1", 10, MPFR_RNDN); mpfr_init(r24151); mpfr_init(r24152); mpfr_init(r24153); mpfr_init(r24154); mpfr_init_set_str(r24155, "6", 10, MPFR_RNDN); mpfr_init(r24156); mpfr_init_set_str(r24157, "2", 10, MPFR_RNDN); mpfr_init(r24158); mpfr_init(r24159); mpfr_init(r24160); mpfr_init(r24161); mpfr_init(r24162); } double f_dm(double v) { ; mpfr_const_pi(r24146, MPFR_RNDN); ; mpfr_mul(r24148, r24146, r24147, MPFR_RNDN); mpfr_div(r24149, r24145, r24148, MPFR_RNDN); ; mpfr_set_d(r24151, v, MPFR_RNDN); mpfr_mul(r24152, r24151, r24151, MPFR_RNDN); mpfr_sub(r24153, r24150, r24152, MPFR_RNDN); mpfr_div(r24154, r24149, r24153, MPFR_RNDN); ; mpfr_neg(r24156, r24155, MPFR_RNDN); ; mpfr_fma(r24158, r24156, r24152, r24157, MPFR_RNDN); mpfr_sqrt(r24159, r24158, MPFR_RNDN); mpfr_div(r24160, r24154, r24159, MPFR_RNDN); mpfr_pow(r24161, r24160, r24147, MPFR_RNDN); mpfr_cbrt(r24162, r24161, MPFR_RNDN); return mpfr_get_d(r24162, MPFR_RNDN); }