#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r21216 = 4.0f; float r21217 = 3.0f; float r21218 = atan2(1.0, 0.0); float r21219 = r21217 * r21218; float r21220 = 1.0f; float r21221 = v; float r21222 = r21221 * r21221; float r21223 = r21220 - r21222; float r21224 = r21219 * r21223; float r21225 = 2.0f; float r21226 = 6.0f; float r21227 = r21226 * r21222; float r21228 = r21225 - r21227; float r21229 = sqrt(r21228); float r21230 = r21224 * r21229; float r21231 = r21216 / r21230; return r21231; } double f_id(double v) { double r21232 = 4.0; double r21233 = 3.0; double r21234 = atan2(1.0, 0.0); double r21235 = r21233 * r21234; double r21236 = 1.0; double r21237 = v; double r21238 = r21237 * r21237; double r21239 = r21236 - r21238; double r21240 = r21235 * r21239; double r21241 = 2.0; double r21242 = 6.0; double r21243 = r21242 * r21238; double r21244 = r21241 - r21243; double r21245 = sqrt(r21244); double r21246 = r21240 * r21245; double r21247 = r21232 / r21246; return r21247; } double f_of(float v) { float r21248 = 1.3333333333333333f; float r21249 = atan2(1.0, 0.0); float r21250 = r21248 / r21249; float r21251 = 2.0f; float r21252 = sqrt(r21251); float r21253 = r21250 / r21252; float r21254 = v; float r21255 = r21254 / r21252; float r21256 = 4.0f; float r21257 = r21256 / r21251; float r21258 = r21257 + r21248; float r21259 = r21255 * r21258; float r21260 = r21249 / r21254; float r21261 = r21259 / r21260; float r21262 = r21253 + r21261; float r21263 = pow(r21254, r21256); float r21264 = 18.0f; float r21265 = r21264 / r21249; float r21266 = r21263 * r21265; float r21267 = 5.0f; float r21268 = pow(r21252, r21267); float r21269 = r21266 / r21268; float r21270 = r21248 + r21257; float r21271 = r21263 * r21270; float r21272 = r21252 * r21249; float r21273 = r21271 / r21272; float r21274 = r21269 + r21273; float r21275 = r21262 + r21274; return r21275; } double f_od(double v) { double r21276 = 1.3333333333333333; double r21277 = atan2(1.0, 0.0); double r21278 = r21276 / r21277; double r21279 = 2.0; double r21280 = sqrt(r21279); double r21281 = r21278 / r21280; double r21282 = v; double r21283 = r21282 / r21280; double r21284 = 4.0; double r21285 = r21284 / r21279; double r21286 = r21285 + r21276; double r21287 = r21283 * r21286; double r21288 = r21277 / r21282; double r21289 = r21287 / r21288; double r21290 = r21281 + r21289; double r21291 = pow(r21282, r21284); double r21292 = 18.0; double r21293 = r21292 / r21277; double r21294 = r21291 * r21293; double r21295 = 5.0; double r21296 = pow(r21280, r21295); double r21297 = r21294 / r21296; double r21298 = r21276 + r21285; double r21299 = r21291 * r21298; double r21300 = r21280 * r21277; double r21301 = r21299 / r21300; double r21302 = r21297 + r21301; double r21303 = r21290 + r21302; return r21303; } 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 r21304, r21305, r21306, r21307, r21308, r21309, r21310, r21311, r21312, r21313, r21314, r21315, r21316, r21317, r21318, r21319; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21304, "4", 10, MPFR_RNDN); mpfr_init_set_str(r21305, "3", 10, MPFR_RNDN); mpfr_init(r21306); mpfr_init(r21307); mpfr_init_set_str(r21308, "1", 10, MPFR_RNDN); mpfr_init(r21309); mpfr_init(r21310); mpfr_init(r21311); mpfr_init(r21312); mpfr_init_set_str(r21313, "2", 10, MPFR_RNDN); mpfr_init_set_str(r21314, "6", 10, MPFR_RNDN); mpfr_init(r21315); mpfr_init(r21316); mpfr_init(r21317); mpfr_init(r21318); mpfr_init(r21319); } double f_im(double v) { ; ; mpfr_const_pi(r21306, MPFR_RNDN); mpfr_mul(r21307, r21305, r21306, MPFR_RNDN); ; mpfr_set_d(r21309, v, MPFR_RNDN); mpfr_sqr(r21310, r21309, MPFR_RNDN); mpfr_sub(r21311, r21308, r21310, MPFR_RNDN); mpfr_mul(r21312, r21307, r21311, MPFR_RNDN); ; ; mpfr_mul(r21315, r21314, r21310, MPFR_RNDN); mpfr_sub(r21316, r21313, r21315, MPFR_RNDN); mpfr_sqrt(r21317, r21316, MPFR_RNDN); mpfr_mul(r21318, r21312, r21317, MPFR_RNDN); mpfr_div(r21319, r21304, r21318, MPFR_RNDN); return mpfr_get_d(r21319, MPFR_RNDN); } static mpfr_t r21320, r21321, r21322, r21323, r21324, r21325, r21326, r21327, r21328, r21329, r21330, r21331, r21332, r21333, r21334, r21335, r21336, r21337, r21338, r21339, r21340, r21341, r21342, r21343, r21344, r21345, r21346, r21347; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21320, "4/3", 10, MPFR_RNDN); mpfr_init(r21321); mpfr_init(r21322); mpfr_init_set_str(r21323, "2", 10, MPFR_RNDN); mpfr_init(r21324); mpfr_init(r21325); mpfr_init(r21326); mpfr_init(r21327); mpfr_init_set_str(r21328, "4", 10, MPFR_RNDN); mpfr_init(r21329); mpfr_init(r21330); mpfr_init(r21331); mpfr_init(r21332); mpfr_init(r21333); mpfr_init(r21334); mpfr_init(r21335); mpfr_init_set_str(r21336, "18", 10, MPFR_RNDN); mpfr_init(r21337); mpfr_init(r21338); mpfr_init_set_str(r21339, "5", 10, MPFR_RNDN); mpfr_init(r21340); mpfr_init(r21341); mpfr_init(r21342); mpfr_init(r21343); mpfr_init(r21344); mpfr_init(r21345); mpfr_init(r21346); mpfr_init(r21347); } double f_fm(double v) { ; mpfr_const_pi(r21321, MPFR_RNDN); mpfr_div(r21322, r21320, r21321, MPFR_RNDN); ; mpfr_sqrt(r21324, r21323, MPFR_RNDN); mpfr_div(r21325, r21322, r21324, MPFR_RNDN); mpfr_set_d(r21326, v, MPFR_RNDN); mpfr_div(r21327, r21326, r21324, MPFR_RNDN); ; mpfr_div(r21329, r21328, r21323, MPFR_RNDN); mpfr_add(r21330, r21329, r21320, MPFR_RNDN); mpfr_mul(r21331, r21327, r21330, MPFR_RNDN); mpfr_div(r21332, r21321, r21326, MPFR_RNDN); mpfr_div(r21333, r21331, r21332, MPFR_RNDN); mpfr_add(r21334, r21325, r21333, MPFR_RNDN); mpfr_pow(r21335, r21326, r21328, MPFR_RNDN); ; mpfr_div(r21337, r21336, r21321, MPFR_RNDN); mpfr_mul(r21338, r21335, r21337, MPFR_RNDN); ; mpfr_pow(r21340, r21324, r21339, MPFR_RNDN); mpfr_div(r21341, r21338, r21340, MPFR_RNDN); mpfr_add(r21342, r21320, r21329, MPFR_RNDN); mpfr_mul(r21343, r21335, r21342, MPFR_RNDN); mpfr_mul(r21344, r21324, r21321, MPFR_RNDN); mpfr_div(r21345, r21343, r21344, MPFR_RNDN); mpfr_add(r21346, r21341, r21345, MPFR_RNDN); mpfr_add(r21347, r21334, r21346, MPFR_RNDN); return mpfr_get_d(r21347, MPFR_RNDN); } static mpfr_t r21348, r21349, r21350, r21351, r21352, r21353, r21354, r21355, r21356, r21357, r21358, r21359, r21360, r21361, r21362, r21363, r21364, r21365, r21366, r21367, r21368, r21369, r21370, r21371, r21372, r21373, r21374, r21375; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21348, "4/3", 10, MPFR_RNDN); mpfr_init(r21349); mpfr_init(r21350); mpfr_init_set_str(r21351, "2", 10, MPFR_RNDN); mpfr_init(r21352); mpfr_init(r21353); mpfr_init(r21354); mpfr_init(r21355); mpfr_init_set_str(r21356, "4", 10, MPFR_RNDN); mpfr_init(r21357); mpfr_init(r21358); mpfr_init(r21359); mpfr_init(r21360); mpfr_init(r21361); mpfr_init(r21362); mpfr_init(r21363); mpfr_init_set_str(r21364, "18", 10, MPFR_RNDN); mpfr_init(r21365); mpfr_init(r21366); mpfr_init_set_str(r21367, "5", 10, MPFR_RNDN); mpfr_init(r21368); mpfr_init(r21369); mpfr_init(r21370); mpfr_init(r21371); mpfr_init(r21372); mpfr_init(r21373); mpfr_init(r21374); mpfr_init(r21375); } double f_dm(double v) { ; mpfr_const_pi(r21349, MPFR_RNDN); mpfr_div(r21350, r21348, r21349, MPFR_RNDN); ; mpfr_sqrt(r21352, r21351, MPFR_RNDN); mpfr_div(r21353, r21350, r21352, MPFR_RNDN); mpfr_set_d(r21354, v, MPFR_RNDN); mpfr_div(r21355, r21354, r21352, MPFR_RNDN); ; mpfr_div(r21357, r21356, r21351, MPFR_RNDN); mpfr_add(r21358, r21357, r21348, MPFR_RNDN); mpfr_mul(r21359, r21355, r21358, MPFR_RNDN); mpfr_div(r21360, r21349, r21354, MPFR_RNDN); mpfr_div(r21361, r21359, r21360, MPFR_RNDN); mpfr_add(r21362, r21353, r21361, MPFR_RNDN); mpfr_pow(r21363, r21354, r21356, MPFR_RNDN); ; mpfr_div(r21365, r21364, r21349, MPFR_RNDN); mpfr_mul(r21366, r21363, r21365, MPFR_RNDN); ; mpfr_pow(r21368, r21352, r21367, MPFR_RNDN); mpfr_div(r21369, r21366, r21368, MPFR_RNDN); mpfr_add(r21370, r21348, r21357, MPFR_RNDN); mpfr_mul(r21371, r21363, r21370, MPFR_RNDN); mpfr_mul(r21372, r21352, r21349, MPFR_RNDN); mpfr_div(r21373, r21371, r21372, MPFR_RNDN); mpfr_add(r21374, r21369, r21373, MPFR_RNDN); mpfr_add(r21375, r21362, r21374, MPFR_RNDN); return mpfr_get_d(r21375, MPFR_RNDN); }