#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r21184 = 4.0f; float r21185 = 3.0f; float r21186 = atan2(1.0, 0.0); float r21187 = r21185 * r21186; float r21188 = 1.0f; float r21189 = v; float r21190 = r21189 * r21189; float r21191 = r21188 - r21190; float r21192 = r21187 * r21191; float r21193 = 2.0f; float r21194 = 6.0f; float r21195 = r21194 * r21190; float r21196 = r21193 - r21195; float r21197 = sqrt(r21196); float r21198 = r21192 * r21197; float r21199 = r21184 / r21198; return r21199; } double f_id(double v) { double r21200 = 4.0; double r21201 = 3.0; double r21202 = atan2(1.0, 0.0); double r21203 = r21201 * r21202; double r21204 = 1.0; double r21205 = v; double r21206 = r21205 * r21205; double r21207 = r21204 - r21206; double r21208 = r21203 * r21207; double r21209 = 2.0; double r21210 = 6.0; double r21211 = r21210 * r21206; double r21212 = r21209 - r21211; double r21213 = sqrt(r21212); double r21214 = r21208 * r21213; double r21215 = r21200 / r21214; return r21215; } double f_of(float v) { float r21216 = 1.3333333333333333f; float r21217 = atan2(1.0, 0.0); float r21218 = r21216 / r21217; float r21219 = 2.0f; float r21220 = sqrt(r21219); float r21221 = r21218 / r21220; float r21222 = v; float r21223 = r21222 / r21220; float r21224 = 4.0f; float r21225 = r21224 / r21219; float r21226 = r21225 + r21216; float r21227 = r21223 * r21226; float r21228 = r21217 / r21222; float r21229 = r21227 / r21228; float r21230 = r21221 + r21229; float r21231 = pow(r21222, r21224); float r21232 = 18.0f; float r21233 = r21232 / r21217; float r21234 = r21231 * r21233; float r21235 = 5.0f; float r21236 = pow(r21220, r21235); float r21237 = r21234 / r21236; float r21238 = r21216 + r21225; float r21239 = r21231 * r21238; float r21240 = r21220 * r21217; float r21241 = r21239 / r21240; float r21242 = r21237 + r21241; float r21243 = r21230 + r21242; return r21243; } double f_od(double v) { double r21244 = 1.3333333333333333; double r21245 = atan2(1.0, 0.0); double r21246 = r21244 / r21245; double r21247 = 2.0; double r21248 = sqrt(r21247); double r21249 = r21246 / r21248; double r21250 = v; double r21251 = r21250 / r21248; double r21252 = 4.0; double r21253 = r21252 / r21247; double r21254 = r21253 + r21244; double r21255 = r21251 * r21254; double r21256 = r21245 / r21250; double r21257 = r21255 / r21256; double r21258 = r21249 + r21257; double r21259 = pow(r21250, r21252); double r21260 = 18.0; double r21261 = r21260 / r21245; double r21262 = r21259 * r21261; double r21263 = 5.0; double r21264 = pow(r21248, r21263); double r21265 = r21262 / r21264; double r21266 = r21244 + r21253; double r21267 = r21259 * r21266; double r21268 = r21248 * r21245; double r21269 = r21267 / r21268; double r21270 = r21265 + r21269; double r21271 = r21258 + r21270; return r21271; } 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 r21272, r21273, r21274, r21275, r21276, r21277, r21278, r21279, r21280, r21281, r21282, r21283, r21284, r21285, r21286, r21287; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21272, "4", 10, MPFR_RNDN); mpfr_init_set_str(r21273, "3", 10, MPFR_RNDN); mpfr_init(r21274); mpfr_init(r21275); mpfr_init_set_str(r21276, "1", 10, MPFR_RNDN); mpfr_init(r21277); mpfr_init(r21278); mpfr_init(r21279); mpfr_init(r21280); mpfr_init_set_str(r21281, "2", 10, MPFR_RNDN); mpfr_init_set_str(r21282, "6", 10, MPFR_RNDN); mpfr_init(r21283); mpfr_init(r21284); mpfr_init(r21285); mpfr_init(r21286); mpfr_init(r21287); } double f_im(double v) { ; ; mpfr_const_pi(r21274, MPFR_RNDN); mpfr_mul(r21275, r21273, r21274, MPFR_RNDN); ; mpfr_set_d(r21277, v, MPFR_RNDN); mpfr_sqr(r21278, r21277, MPFR_RNDN); mpfr_sub(r21279, r21276, r21278, MPFR_RNDN); mpfr_mul(r21280, r21275, r21279, MPFR_RNDN); ; ; mpfr_mul(r21283, r21282, r21278, MPFR_RNDN); mpfr_sub(r21284, r21281, r21283, MPFR_RNDN); mpfr_sqrt(r21285, r21284, MPFR_RNDN); mpfr_mul(r21286, r21280, r21285, MPFR_RNDN); mpfr_div(r21287, r21272, r21286, MPFR_RNDN); return mpfr_get_d(r21287, MPFR_RNDN); } static mpfr_t r21288, r21289, r21290, r21291, r21292, r21293, r21294, r21295, r21296, r21297, r21298, r21299, r21300, r21301, r21302, r21303, r21304, r21305, r21306, r21307, r21308, r21309, r21310, r21311, r21312, r21313, r21314, r21315; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21288, "4/3", 10, MPFR_RNDN); mpfr_init(r21289); mpfr_init(r21290); mpfr_init_set_str(r21291, "2", 10, MPFR_RNDN); mpfr_init(r21292); mpfr_init(r21293); mpfr_init(r21294); mpfr_init(r21295); mpfr_init_set_str(r21296, "4", 10, MPFR_RNDN); mpfr_init(r21297); mpfr_init(r21298); mpfr_init(r21299); mpfr_init(r21300); mpfr_init(r21301); mpfr_init(r21302); mpfr_init(r21303); mpfr_init_set_str(r21304, "18", 10, MPFR_RNDN); mpfr_init(r21305); mpfr_init(r21306); mpfr_init_set_str(r21307, "5", 10, MPFR_RNDN); mpfr_init(r21308); mpfr_init(r21309); mpfr_init(r21310); mpfr_init(r21311); mpfr_init(r21312); mpfr_init(r21313); mpfr_init(r21314); mpfr_init(r21315); } double f_fm(double v) { ; mpfr_const_pi(r21289, MPFR_RNDN); mpfr_div(r21290, r21288, r21289, MPFR_RNDN); ; mpfr_sqrt(r21292, r21291, MPFR_RNDN); mpfr_div(r21293, r21290, r21292, MPFR_RNDN); mpfr_set_d(r21294, v, MPFR_RNDN); mpfr_div(r21295, r21294, r21292, MPFR_RNDN); ; mpfr_div(r21297, r21296, r21291, MPFR_RNDN); mpfr_add(r21298, r21297, r21288, MPFR_RNDN); mpfr_mul(r21299, r21295, r21298, MPFR_RNDN); mpfr_div(r21300, r21289, r21294, MPFR_RNDN); mpfr_div(r21301, r21299, r21300, MPFR_RNDN); mpfr_add(r21302, r21293, r21301, MPFR_RNDN); mpfr_pow(r21303, r21294, r21296, MPFR_RNDN); ; mpfr_div(r21305, r21304, r21289, MPFR_RNDN); mpfr_mul(r21306, r21303, r21305, MPFR_RNDN); ; mpfr_pow(r21308, r21292, r21307, MPFR_RNDN); mpfr_div(r21309, r21306, r21308, MPFR_RNDN); mpfr_add(r21310, r21288, r21297, MPFR_RNDN); mpfr_mul(r21311, r21303, r21310, MPFR_RNDN); mpfr_mul(r21312, r21292, r21289, MPFR_RNDN); mpfr_div(r21313, r21311, r21312, MPFR_RNDN); mpfr_add(r21314, r21309, r21313, MPFR_RNDN); mpfr_add(r21315, r21302, r21314, MPFR_RNDN); return mpfr_get_d(r21315, MPFR_RNDN); } static mpfr_t r21316, r21317, r21318, r21319, r21320, r21321, r21322, r21323, r21324, r21325, r21326, r21327, r21328, r21329, r21330, r21331, r21332, r21333, r21334, r21335, r21336, r21337, r21338, r21339, r21340, r21341, r21342, r21343; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21316, "4/3", 10, MPFR_RNDN); mpfr_init(r21317); mpfr_init(r21318); mpfr_init_set_str(r21319, "2", 10, MPFR_RNDN); mpfr_init(r21320); mpfr_init(r21321); mpfr_init(r21322); mpfr_init(r21323); mpfr_init_set_str(r21324, "4", 10, MPFR_RNDN); mpfr_init(r21325); mpfr_init(r21326); mpfr_init(r21327); mpfr_init(r21328); mpfr_init(r21329); mpfr_init(r21330); mpfr_init(r21331); mpfr_init_set_str(r21332, "18", 10, MPFR_RNDN); mpfr_init(r21333); mpfr_init(r21334); mpfr_init_set_str(r21335, "5", 10, MPFR_RNDN); mpfr_init(r21336); mpfr_init(r21337); mpfr_init(r21338); mpfr_init(r21339); mpfr_init(r21340); mpfr_init(r21341); mpfr_init(r21342); mpfr_init(r21343); } double f_dm(double v) { ; mpfr_const_pi(r21317, MPFR_RNDN); mpfr_div(r21318, r21316, r21317, MPFR_RNDN); ; mpfr_sqrt(r21320, r21319, MPFR_RNDN); mpfr_div(r21321, r21318, r21320, MPFR_RNDN); mpfr_set_d(r21322, v, MPFR_RNDN); mpfr_div(r21323, r21322, r21320, MPFR_RNDN); ; mpfr_div(r21325, r21324, r21319, MPFR_RNDN); mpfr_add(r21326, r21325, r21316, MPFR_RNDN); mpfr_mul(r21327, r21323, r21326, MPFR_RNDN); mpfr_div(r21328, r21317, r21322, MPFR_RNDN); mpfr_div(r21329, r21327, r21328, MPFR_RNDN); mpfr_add(r21330, r21321, r21329, MPFR_RNDN); mpfr_pow(r21331, r21322, r21324, MPFR_RNDN); ; mpfr_div(r21333, r21332, r21317, MPFR_RNDN); mpfr_mul(r21334, r21331, r21333, MPFR_RNDN); ; mpfr_pow(r21336, r21320, r21335, MPFR_RNDN); mpfr_div(r21337, r21334, r21336, MPFR_RNDN); mpfr_add(r21338, r21316, r21325, MPFR_RNDN); mpfr_mul(r21339, r21331, r21338, MPFR_RNDN); mpfr_mul(r21340, r21320, r21317, MPFR_RNDN); mpfr_div(r21341, r21339, r21340, MPFR_RNDN); mpfr_add(r21342, r21337, r21341, MPFR_RNDN); mpfr_add(r21343, r21330, r21342, MPFR_RNDN); return mpfr_get_d(r21343, MPFR_RNDN); }