#include #include #include #include #include char *name = "2isqrt (example 3.6)"; double f_if(float x) { float r23154 = 1; float r23155 = x; float r23156 = sqrt(r23155); float r23157 = r23154 / r23156; float r23158 = r23155 + r23154; float r23159 = sqrt(r23158); float r23160 = r23154 / r23159; float r23161 = r23157 - r23160; return r23161; } double f_id(double x) { double r23162 = 1; double r23163 = x; double r23164 = sqrt(r23163); double r23165 = r23162 / r23164; double r23166 = r23163 + r23162; double r23167 = sqrt(r23166); double r23168 = r23162 / r23167; double r23169 = r23165 - r23168; return r23169; } double f_of(float x) { float r23170 = x; float r23171 = 1/2; float r23172 = -r23171; float r23173 = pow(r23170, r23172); float r23174 = r23173 * r23173; float r23175 = 1; float r23176 = r23170 + r23175; float r23177 = r23175 / r23176; float r23178 = r23174 - r23177; float r23179 = sqrt(r23176); float r23180 = r23175 / r23179; float r23181 = r23173 + r23180; float r23182 = r23178 / r23181; float r23183 = 4.381860511644953e-78; bool r23184 = r23182 <= r23183; float r23185 = sqrt(r23170); float r23186 = r23175 / r23185; float r23187 = sqrt(r23186); float r23188 = sqrt(r23180); float r23189 = r23187 + r23188; float r23190 = r23187 - r23188; float r23191 = log(r23190); float r23192 = exp(r23191); float r23193 = r23189 * r23192; float r23194 = r23184 ? r23193 : r23182; return r23194; } double f_od(double x) { double r23195 = x; double r23196 = 1/2; double r23197 = -r23196; double r23198 = pow(r23195, r23197); double r23199 = r23198 * r23198; double r23200 = 1; double r23201 = r23195 + r23200; double r23202 = r23200 / r23201; double r23203 = r23199 - r23202; double r23204 = sqrt(r23201); double r23205 = r23200 / r23204; double r23206 = r23198 + r23205; double r23207 = r23203 / r23206; double r23208 = 4.381860511644953e-78; bool r23209 = r23207 <= r23208; double r23210 = sqrt(r23195); double r23211 = r23200 / r23210; double r23212 = sqrt(r23211); double r23213 = sqrt(r23205); double r23214 = r23212 + r23213; double r23215 = r23212 - r23213; double r23216 = log(r23215); double r23217 = exp(r23216); double r23218 = r23214 * r23217; double r23219 = r23209 ? r23218 : r23207; return r23219; } 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 r23220, r23221, r23222, r23223, r23224, r23225, r23226, r23227; void setup_mpfr_f_im() { mpfr_set_default_prec(1168); mpfr_init_set_str(r23220, "1", 10, MPFR_RNDN); mpfr_init(r23221); mpfr_init(r23222); mpfr_init(r23223); mpfr_init(r23224); mpfr_init(r23225); mpfr_init(r23226); mpfr_init(r23227); } double f_im(double x) { ; mpfr_set_d(r23221, x, MPFR_RNDN); mpfr_sqrt(r23222, r23221, MPFR_RNDN); mpfr_div(r23223, r23220, r23222, MPFR_RNDN); mpfr_add(r23224, r23221, r23220, MPFR_RNDN); mpfr_sqrt(r23225, r23224, MPFR_RNDN); mpfr_div(r23226, r23220, r23225, MPFR_RNDN); mpfr_sub(r23227, r23223, r23226, MPFR_RNDN); return mpfr_get_d(r23227, MPFR_RNDN); } static mpfr_t r23228, r23229, r23230, r23231, r23232, r23233, r23234, r23235, r23236, r23237, r23238, r23239, r23240, r23241, r23242, r23243, r23244, r23245, r23246, r23247, r23248, r23249, r23250, r23251, r23252; void setup_mpfr_f_fm() { mpfr_set_default_prec(1168); mpfr_init(r23228); mpfr_init_set_str(r23229, "1/2", 10, MPFR_RNDN); mpfr_init(r23230); mpfr_init(r23231); mpfr_init(r23232); mpfr_init_set_str(r23233, "1", 10, MPFR_RNDN); mpfr_init(r23234); mpfr_init(r23235); mpfr_init(r23236); mpfr_init(r23237); mpfr_init(r23238); mpfr_init(r23239); mpfr_init(r23240); mpfr_init_set_str(r23241, "4.381860511644953e-78", 10, MPFR_RNDN); mpfr_init(r23242); mpfr_init(r23243); mpfr_init(r23244); mpfr_init(r23245); mpfr_init(r23246); mpfr_init(r23247); mpfr_init(r23248); mpfr_init(r23249); mpfr_init(r23250); mpfr_init(r23251); mpfr_init(r23252); } double f_fm(double x) { mpfr_set_d(r23228, x, MPFR_RNDN); ; mpfr_neg(r23230, r23229, MPFR_RNDN); mpfr_pow(r23231, r23228, r23230, MPFR_RNDN); mpfr_mul(r23232, r23231, r23231, MPFR_RNDN); ; mpfr_add(r23234, r23228, r23233, MPFR_RNDN); mpfr_div(r23235, r23233, r23234, MPFR_RNDN); mpfr_sub(r23236, r23232, r23235, MPFR_RNDN); mpfr_sqrt(r23237, r23234, MPFR_RNDN); mpfr_div(r23238, r23233, r23237, MPFR_RNDN); mpfr_add(r23239, r23231, r23238, MPFR_RNDN); mpfr_div(r23240, r23236, r23239, MPFR_RNDN); ; mpfr_set_si(r23242, mpfr_cmp(r23240, r23241) <= 0, MPFR_RNDN); mpfr_sqrt(r23243, r23228, MPFR_RNDN); mpfr_div(r23244, r23233, r23243, MPFR_RNDN); mpfr_sqrt(r23245, r23244, MPFR_RNDN); mpfr_sqrt(r23246, r23238, MPFR_RNDN); mpfr_add(r23247, r23245, r23246, MPFR_RNDN); mpfr_sub(r23248, r23245, r23246, MPFR_RNDN); mpfr_log(r23249, r23248, MPFR_RNDN); mpfr_exp(r23250, r23249, MPFR_RNDN); mpfr_mul(r23251, r23247, r23250, MPFR_RNDN); if (mpfr_get_si(r23242, MPFR_RNDN)) { mpfr_set(r23252, r23251, MPFR_RNDN); } else { mpfr_set(r23252, r23240, MPFR_RNDN); }; return mpfr_get_d(r23252, MPFR_RNDN); } static mpfr_t r23253, r23254, r23255, r23256, r23257, r23258, r23259, r23260, r23261, r23262, r23263, r23264, r23265, r23266, r23267, r23268, r23269, r23270, r23271, r23272, r23273, r23274, r23275, r23276, r23277; void setup_mpfr_f_dm() { mpfr_set_default_prec(1168); mpfr_init(r23253); mpfr_init_set_str(r23254, "1/2", 10, MPFR_RNDN); mpfr_init(r23255); mpfr_init(r23256); mpfr_init(r23257); mpfr_init_set_str(r23258, "1", 10, MPFR_RNDN); mpfr_init(r23259); mpfr_init(r23260); mpfr_init(r23261); mpfr_init(r23262); mpfr_init(r23263); mpfr_init(r23264); mpfr_init(r23265); mpfr_init_set_str(r23266, "4.381860511644953e-78", 10, MPFR_RNDN); mpfr_init(r23267); mpfr_init(r23268); mpfr_init(r23269); mpfr_init(r23270); mpfr_init(r23271); mpfr_init(r23272); mpfr_init(r23273); mpfr_init(r23274); mpfr_init(r23275); mpfr_init(r23276); mpfr_init(r23277); } double f_dm(double x) { mpfr_set_d(r23253, x, MPFR_RNDN); ; mpfr_neg(r23255, r23254, MPFR_RNDN); mpfr_pow(r23256, r23253, r23255, MPFR_RNDN); mpfr_mul(r23257, r23256, r23256, MPFR_RNDN); ; mpfr_add(r23259, r23253, r23258, MPFR_RNDN); mpfr_div(r23260, r23258, r23259, MPFR_RNDN); mpfr_sub(r23261, r23257, r23260, MPFR_RNDN); mpfr_sqrt(r23262, r23259, MPFR_RNDN); mpfr_div(r23263, r23258, r23262, MPFR_RNDN); mpfr_add(r23264, r23256, r23263, MPFR_RNDN); mpfr_div(r23265, r23261, r23264, MPFR_RNDN); ; mpfr_set_si(r23267, mpfr_cmp(r23265, r23266) <= 0, MPFR_RNDN); mpfr_sqrt(r23268, r23253, MPFR_RNDN); mpfr_div(r23269, r23258, r23268, MPFR_RNDN); mpfr_sqrt(r23270, r23269, MPFR_RNDN); mpfr_sqrt(r23271, r23263, MPFR_RNDN); mpfr_add(r23272, r23270, r23271, MPFR_RNDN); mpfr_sub(r23273, r23270, r23271, MPFR_RNDN); mpfr_log(r23274, r23273, MPFR_RNDN); mpfr_exp(r23275, r23274, MPFR_RNDN); mpfr_mul(r23276, r23272, r23275, MPFR_RNDN); if (mpfr_get_si(r23267, MPFR_RNDN)) { mpfr_set(r23277, r23276, MPFR_RNDN); } else { mpfr_set(r23277, r23265, MPFR_RNDN); }; return mpfr_get_d(r23277, MPFR_RNDN); }