#include #include #include #include #include char *name = "2isqrt (example 3.6)"; double f_if(float x) { float r27232 = 1; float r27233 = x; float r27234 = sqrt(r27233); float r27235 = r27232 / r27234; float r27236 = r27233 + r27232; float r27237 = sqrt(r27236); float r27238 = r27232 / r27237; float r27239 = r27235 - r27238; return r27239; } double f_id(double x) { double r27240 = 1; double r27241 = x; double r27242 = sqrt(r27241); double r27243 = r27240 / r27242; double r27244 = r27241 + r27240; double r27245 = sqrt(r27244); double r27246 = r27240 / r27245; double r27247 = r27243 - r27246; return r27247; } double f_of(float x) { float r27248 = 1; float r27249 = x; float r27250 = sqrt(r27249); float r27251 = r27248 / r27250; float r27252 = r27249 + r27248; float r27253 = sqrt(r27252); float r27254 = r27248 / r27253; float r27255 = r27251 - r27254; float r27256 = 0.0048865480386920505; bool r27257 = r27255 <= r27256; float r27258 = r27252 / r27250; float r27259 = r27258 / r27249; float r27260 = r27249 / r27253; float r27261 = fma(r27249, r27259, r27260); float r27262 = r27248 / r27261; float r27263 = r27251 * r27251; float r27264 = r27254 * r27254; float r27265 = r27251 * r27254; float r27266 = r27264 - r27265; float r27267 = r27263 + r27266; float r27268 = r27262 * r27267; float r27269 = 1/2; float r27270 = -r27269; float r27271 = pow(r27249, r27270); float r27272 = r27271 - r27254; float r27273 = r27257 ? r27268 : r27272; return r27273; } double f_od(double x) { double r27274 = 1; double r27275 = x; double r27276 = sqrt(r27275); double r27277 = r27274 / r27276; double r27278 = r27275 + r27274; double r27279 = sqrt(r27278); double r27280 = r27274 / r27279; double r27281 = r27277 - r27280; double r27282 = 0.0048865480386920505; bool r27283 = r27281 <= r27282; double r27284 = r27278 / r27276; double r27285 = r27284 / r27275; double r27286 = r27275 / r27279; double r27287 = fma(r27275, r27285, r27286); double r27288 = r27274 / r27287; double r27289 = r27277 * r27277; double r27290 = r27280 * r27280; double r27291 = r27277 * r27280; double r27292 = r27290 - r27291; double r27293 = r27289 + r27292; double r27294 = r27288 * r27293; double r27295 = 1/2; double r27296 = -r27295; double r27297 = pow(r27275, r27296); double r27298 = r27297 - r27280; double r27299 = r27283 ? r27294 : r27298; return r27299; } 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 r27300, r27301, r27302, r27303, r27304, r27305, r27306, r27307; void setup_mpfr_f_im() { mpfr_set_default_prec(1104); mpfr_init_set_str(r27300, "1", 10, MPFR_RNDN); mpfr_init(r27301); mpfr_init(r27302); mpfr_init(r27303); mpfr_init(r27304); mpfr_init(r27305); mpfr_init(r27306); mpfr_init(r27307); } double f_im(double x) { ; mpfr_set_d(r27301, x, MPFR_RNDN); mpfr_sqrt(r27302, r27301, MPFR_RNDN); mpfr_div(r27303, r27300, r27302, MPFR_RNDN); mpfr_add(r27304, r27301, r27300, MPFR_RNDN); mpfr_sqrt(r27305, r27304, MPFR_RNDN); mpfr_div(r27306, r27300, r27305, MPFR_RNDN); mpfr_sub(r27307, r27303, r27306, MPFR_RNDN); return mpfr_get_d(r27307, MPFR_RNDN); } static mpfr_t r27308, r27309, r27310, r27311, r27312, r27313, r27314, r27315, r27316, r27317, r27318, r27319, r27320, r27321, r27322, r27323, r27324, r27325, r27326, r27327, r27328, r27329, r27330, r27331, r27332, r27333; void setup_mpfr_f_fm() { mpfr_set_default_prec(1104); mpfr_init_set_str(r27308, "1", 10, MPFR_RNDN); mpfr_init(r27309); mpfr_init(r27310); mpfr_init(r27311); mpfr_init(r27312); mpfr_init(r27313); mpfr_init(r27314); mpfr_init(r27315); mpfr_init_set_str(r27316, "0.0048865480386920505", 10, MPFR_RNDN); mpfr_init(r27317); mpfr_init(r27318); mpfr_init(r27319); mpfr_init(r27320); mpfr_init(r27321); mpfr_init(r27322); mpfr_init(r27323); mpfr_init(r27324); mpfr_init(r27325); mpfr_init(r27326); mpfr_init(r27327); mpfr_init(r27328); mpfr_init_set_str(r27329, "1/2", 10, MPFR_RNDN); mpfr_init(r27330); mpfr_init(r27331); mpfr_init(r27332); mpfr_init(r27333); } double f_fm(double x) { ; mpfr_set_d(r27309, x, MPFR_RNDN); mpfr_sqrt(r27310, r27309, MPFR_RNDN); mpfr_div(r27311, r27308, r27310, MPFR_RNDN); mpfr_add(r27312, r27309, r27308, MPFR_RNDN); mpfr_sqrt(r27313, r27312, MPFR_RNDN); mpfr_div(r27314, r27308, r27313, MPFR_RNDN); mpfr_sub(r27315, r27311, r27314, MPFR_RNDN); ; mpfr_set_si(r27317, mpfr_cmp(r27315, r27316) <= 0, MPFR_RNDN); mpfr_div(r27318, r27312, r27310, MPFR_RNDN); mpfr_div(r27319, r27318, r27309, MPFR_RNDN); mpfr_div(r27320, r27309, r27313, MPFR_RNDN); mpfr_fma(r27321, r27309, r27319, r27320, MPFR_RNDN); mpfr_div(r27322, r27308, r27321, MPFR_RNDN); mpfr_mul(r27323, r27311, r27311, MPFR_RNDN); mpfr_mul(r27324, r27314, r27314, MPFR_RNDN); mpfr_mul(r27325, r27311, r27314, MPFR_RNDN); mpfr_sub(r27326, r27324, r27325, MPFR_RNDN); mpfr_add(r27327, r27323, r27326, MPFR_RNDN); mpfr_mul(r27328, r27322, r27327, MPFR_RNDN); ; mpfr_neg(r27330, r27329, MPFR_RNDN); mpfr_pow(r27331, r27309, r27330, MPFR_RNDN); mpfr_sub(r27332, r27331, r27314, MPFR_RNDN); if (mpfr_get_si(r27317, MPFR_RNDN)) { mpfr_set(r27333, r27328, MPFR_RNDN); } else { mpfr_set(r27333, r27332, MPFR_RNDN); }; return mpfr_get_d(r27333, MPFR_RNDN); } static mpfr_t r27334, r27335, r27336, r27337, r27338, r27339, r27340, r27341, r27342, r27343, r27344, r27345, r27346, r27347, r27348, r27349, r27350, r27351, r27352, r27353, r27354, r27355, r27356, r27357, r27358, r27359; void setup_mpfr_f_dm() { mpfr_set_default_prec(1104); mpfr_init_set_str(r27334, "1", 10, MPFR_RNDN); mpfr_init(r27335); mpfr_init(r27336); mpfr_init(r27337); mpfr_init(r27338); mpfr_init(r27339); mpfr_init(r27340); mpfr_init(r27341); mpfr_init_set_str(r27342, "0.0048865480386920505", 10, MPFR_RNDN); mpfr_init(r27343); mpfr_init(r27344); mpfr_init(r27345); mpfr_init(r27346); mpfr_init(r27347); mpfr_init(r27348); mpfr_init(r27349); mpfr_init(r27350); mpfr_init(r27351); mpfr_init(r27352); mpfr_init(r27353); mpfr_init(r27354); mpfr_init_set_str(r27355, "1/2", 10, MPFR_RNDN); mpfr_init(r27356); mpfr_init(r27357); mpfr_init(r27358); mpfr_init(r27359); } double f_dm(double x) { ; mpfr_set_d(r27335, x, MPFR_RNDN); mpfr_sqrt(r27336, r27335, MPFR_RNDN); mpfr_div(r27337, r27334, r27336, MPFR_RNDN); mpfr_add(r27338, r27335, r27334, MPFR_RNDN); mpfr_sqrt(r27339, r27338, MPFR_RNDN); mpfr_div(r27340, r27334, r27339, MPFR_RNDN); mpfr_sub(r27341, r27337, r27340, MPFR_RNDN); ; mpfr_set_si(r27343, mpfr_cmp(r27341, r27342) <= 0, MPFR_RNDN); mpfr_div(r27344, r27338, r27336, MPFR_RNDN); mpfr_div(r27345, r27344, r27335, MPFR_RNDN); mpfr_div(r27346, r27335, r27339, MPFR_RNDN); mpfr_fma(r27347, r27335, r27345, r27346, MPFR_RNDN); mpfr_div(r27348, r27334, r27347, MPFR_RNDN); mpfr_mul(r27349, r27337, r27337, MPFR_RNDN); mpfr_mul(r27350, r27340, r27340, MPFR_RNDN); mpfr_mul(r27351, r27337, r27340, MPFR_RNDN); mpfr_sub(r27352, r27350, r27351, MPFR_RNDN); mpfr_add(r27353, r27349, r27352, MPFR_RNDN); mpfr_mul(r27354, r27348, r27353, MPFR_RNDN); ; mpfr_neg(r27356, r27355, MPFR_RNDN); mpfr_pow(r27357, r27335, r27356, MPFR_RNDN); mpfr_sub(r27358, r27357, r27340, MPFR_RNDN); if (mpfr_get_si(r27343, MPFR_RNDN)) { mpfr_set(r27359, r27354, MPFR_RNDN); } else { mpfr_set(r27359, r27358, MPFR_RNDN); }; return mpfr_get_d(r27359, MPFR_RNDN); }