#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r27074 = 4; float r27075 = 3; float r27076 = atan2(1.0, 0.0); float r27077 = r27075 * r27076; float r27078 = 1; float r27079 = v; float r27080 = r27079 * r27079; float r27081 = r27078 - r27080; float r27082 = r27077 * r27081; float r27083 = 2; float r27084 = 6; float r27085 = r27084 * r27080; float r27086 = r27083 - r27085; float r27087 = sqrt(r27086); float r27088 = r27082 * r27087; float r27089 = r27074 / r27088; return r27089; } double f_id(double v) { double r27090 = 4; double r27091 = 3; double r27092 = atan2(1.0, 0.0); double r27093 = r27091 * r27092; double r27094 = 1; double r27095 = v; double r27096 = r27095 * r27095; double r27097 = r27094 - r27096; double r27098 = r27093 * r27097; double r27099 = 2; double r27100 = 6; double r27101 = r27100 * r27096; double r27102 = r27099 - r27101; double r27103 = sqrt(r27102); double r27104 = r27098 * r27103; double r27105 = r27090 / r27104; return r27105; } double f_of(float v) { float r27106 = 4; float r27107 = atan2(1.0, 0.0); float r27108 = 3; float r27109 = r27107 * r27108; float r27110 = r27106 / r27109; float r27111 = 1; float r27112 = v; float r27113 = r27112 * r27112; float r27114 = r27111 - r27113; float r27115 = r27110 / r27114; float r27116 = 6; float r27117 = -r27116; float r27118 = 2; float r27119 = fma(r27117, r27113, r27118); float r27120 = sqrt(r27119); float r27121 = r27115 / r27120; float r27122 = pow(r27121, r27108); float r27123 = cbrt(r27122); return r27123; } double f_od(double v) { double r27124 = 4; double r27125 = atan2(1.0, 0.0); double r27126 = 3; double r27127 = r27125 * r27126; double r27128 = r27124 / r27127; double r27129 = 1; double r27130 = v; double r27131 = r27130 * r27130; double r27132 = r27129 - r27131; double r27133 = r27128 / r27132; double r27134 = 6; double r27135 = -r27134; double r27136 = 2; double r27137 = fma(r27135, r27131, r27136); double r27138 = sqrt(r27137); double r27139 = r27133 / r27138; double r27140 = pow(r27139, r27126); double r27141 = cbrt(r27140); return r27141; } 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 r27142, r27143, r27144, r27145, r27146, r27147, r27148, r27149, r27150, r27151, r27152, r27153, r27154, r27155, r27156, r27157; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r27142, "4", 10, MPFR_RNDN); mpfr_init_set_str(r27143, "3", 10, MPFR_RNDN); mpfr_init(r27144); mpfr_init(r27145); mpfr_init_set_str(r27146, "1", 10, MPFR_RNDN); mpfr_init(r27147); mpfr_init(r27148); mpfr_init(r27149); mpfr_init(r27150); mpfr_init_set_str(r27151, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27152, "6", 10, MPFR_RNDN); mpfr_init(r27153); mpfr_init(r27154); mpfr_init(r27155); mpfr_init(r27156); mpfr_init(r27157); } double f_im(double v) { ; ; mpfr_const_pi(r27144, MPFR_RNDN); mpfr_mul(r27145, r27143, r27144, MPFR_RNDN); ; mpfr_set_d(r27147, v, MPFR_RNDN); mpfr_mul(r27148, r27147, r27147, MPFR_RNDN); mpfr_sub(r27149, r27146, r27148, MPFR_RNDN); mpfr_mul(r27150, r27145, r27149, MPFR_RNDN); ; ; mpfr_mul(r27153, r27152, r27148, MPFR_RNDN); mpfr_sub(r27154, r27151, r27153, MPFR_RNDN); mpfr_sqrt(r27155, r27154, MPFR_RNDN); mpfr_mul(r27156, r27150, r27155, MPFR_RNDN); mpfr_div(r27157, r27142, r27156, MPFR_RNDN); return mpfr_get_d(r27157, MPFR_RNDN); } static mpfr_t r27158, r27159, r27160, r27161, r27162, r27163, r27164, r27165, r27166, r27167, r27168, r27169, r27170, r27171, r27172, r27173, r27174, r27175; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27158, "4", 10, MPFR_RNDN); mpfr_init(r27159); mpfr_init_set_str(r27160, "3", 10, MPFR_RNDN); mpfr_init(r27161); mpfr_init(r27162); mpfr_init_set_str(r27163, "1", 10, MPFR_RNDN); mpfr_init(r27164); mpfr_init(r27165); mpfr_init(r27166); mpfr_init(r27167); mpfr_init_set_str(r27168, "6", 10, MPFR_RNDN); mpfr_init(r27169); mpfr_init_set_str(r27170, "2", 10, MPFR_RNDN); mpfr_init(r27171); mpfr_init(r27172); mpfr_init(r27173); mpfr_init(r27174); mpfr_init(r27175); } double f_fm(double v) { ; mpfr_const_pi(r27159, MPFR_RNDN); ; mpfr_mul(r27161, r27159, r27160, MPFR_RNDN); mpfr_div(r27162, r27158, r27161, MPFR_RNDN); ; mpfr_set_d(r27164, v, MPFR_RNDN); mpfr_mul(r27165, r27164, r27164, MPFR_RNDN); mpfr_sub(r27166, r27163, r27165, MPFR_RNDN); mpfr_div(r27167, r27162, r27166, MPFR_RNDN); ; mpfr_neg(r27169, r27168, MPFR_RNDN); ; mpfr_fma(r27171, r27169, r27165, r27170, MPFR_RNDN); mpfr_sqrt(r27172, r27171, MPFR_RNDN); mpfr_div(r27173, r27167, r27172, MPFR_RNDN); mpfr_pow(r27174, r27173, r27160, MPFR_RNDN); mpfr_cbrt(r27175, r27174, MPFR_RNDN); return mpfr_get_d(r27175, MPFR_RNDN); } static mpfr_t r27176, r27177, r27178, r27179, r27180, r27181, r27182, r27183, r27184, r27185, r27186, r27187, r27188, r27189, r27190, r27191, r27192, r27193; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27176, "4", 10, MPFR_RNDN); mpfr_init(r27177); mpfr_init_set_str(r27178, "3", 10, MPFR_RNDN); mpfr_init(r27179); mpfr_init(r27180); mpfr_init_set_str(r27181, "1", 10, MPFR_RNDN); mpfr_init(r27182); mpfr_init(r27183); mpfr_init(r27184); mpfr_init(r27185); mpfr_init_set_str(r27186, "6", 10, MPFR_RNDN); mpfr_init(r27187); mpfr_init_set_str(r27188, "2", 10, MPFR_RNDN); mpfr_init(r27189); mpfr_init(r27190); mpfr_init(r27191); mpfr_init(r27192); mpfr_init(r27193); } double f_dm(double v) { ; mpfr_const_pi(r27177, MPFR_RNDN); ; mpfr_mul(r27179, r27177, r27178, MPFR_RNDN); mpfr_div(r27180, r27176, r27179, MPFR_RNDN); ; mpfr_set_d(r27182, v, MPFR_RNDN); mpfr_mul(r27183, r27182, r27182, MPFR_RNDN); mpfr_sub(r27184, r27181, r27183, MPFR_RNDN); mpfr_div(r27185, r27180, r27184, MPFR_RNDN); ; mpfr_neg(r27187, r27186, MPFR_RNDN); ; mpfr_fma(r27189, r27187, r27183, r27188, MPFR_RNDN); mpfr_sqrt(r27190, r27189, MPFR_RNDN); mpfr_div(r27191, r27185, r27190, MPFR_RNDN); mpfr_pow(r27192, r27191, r27178, MPFR_RNDN); mpfr_cbrt(r27193, r27192, MPFR_RNDN); return mpfr_get_d(r27193, MPFR_RNDN); }