#include #include #include #include #include char *name = "math.sqrt on complex, real part"; double f_if(float re, float im) { float r26004 = 0.5; float r26005 = 2.0; float r26006 = re; float r26007 = r26006 * r26006; float r26008 = im; float r26009 = r26008 * r26008; float r26010 = r26007 + r26009; float r26011 = sqrt(r26010); float r26012 = r26011 + r26006; float r26013 = r26005 * r26012; float r26014 = sqrt(r26013); float r26015 = r26004 * r26014; return r26015; } double f_id(double re, double im) { double r26016 = 0.5; double r26017 = 2.0; double r26018 = re; double r26019 = r26018 * r26018; double r26020 = im; double r26021 = r26020 * r26020; double r26022 = r26019 + r26021; double r26023 = sqrt(r26022); double r26024 = r26023 + r26018; double r26025 = r26017 * r26024; double r26026 = sqrt(r26025); double r26027 = r26016 * r26026; return r26027; } double f_of(float re, float im) { float r26028 = 2.0; float r26029 = im; float r26030 = re; float r26031 = r26029 + r26030; float r26032 = r26028 * r26031; float r26033 = -1.0827629684787965e-204; bool r26034 = r26032 <= r26033; float r26035 = 0.5; float r26036 = r26029 * r26029; float r26037 = r26030 * r26030; float r26038 = r26037 + r26036; float r26039 = sqrt(r26038); float r26040 = r26039 - r26030; float r26041 = r26036 / r26040; float r26042 = r26028 * r26041; float r26043 = sqrt(r26042); float r26044 = r26035 * r26043; float r26045 = 5.047375239993434e+166; bool r26046 = r26032 <= r26045; float r26047 = cbrt(r26039); float r26048 = r26047 * r26047; float r26049 = r26048 * r26047; float r26050 = r26049 + r26030; float r26051 = r26028 * r26050; float r26052 = sqrt(r26051); float r26053 = r26035 * r26052; float r26054 = 2.1298540120048136e+187; bool r26055 = r26032 <= r26054; float r26056 = sqrt(r26032); float r26057 = r26035 * r26056; float r26058 = r26030 + r26030; float r26059 = r26028 * r26058; float r26060 = sqrt(r26059); float r26061 = r26035 * r26060; float r26062 = r26055 ? r26057 : r26061; float r26063 = r26046 ? r26053 : r26062; float r26064 = r26034 ? r26044 : r26063; return r26064; } double f_od(double re, double im) { double r26065 = 2.0; double r26066 = im; double r26067 = re; double r26068 = r26066 + r26067; double r26069 = r26065 * r26068; double r26070 = -1.0827629684787965e-204; bool r26071 = r26069 <= r26070; double r26072 = 0.5; double r26073 = r26066 * r26066; double r26074 = r26067 * r26067; double r26075 = r26074 + r26073; double r26076 = sqrt(r26075); double r26077 = r26076 - r26067; double r26078 = r26073 / r26077; double r26079 = r26065 * r26078; double r26080 = sqrt(r26079); double r26081 = r26072 * r26080; double r26082 = 5.047375239993434e+166; bool r26083 = r26069 <= r26082; double r26084 = cbrt(r26076); double r26085 = r26084 * r26084; double r26086 = r26085 * r26084; double r26087 = r26086 + r26067; double r26088 = r26065 * r26087; double r26089 = sqrt(r26088); double r26090 = r26072 * r26089; double r26091 = 2.1298540120048136e+187; bool r26092 = r26069 <= r26091; double r26093 = sqrt(r26069); double r26094 = r26072 * r26093; double r26095 = r26067 + r26067; double r26096 = r26065 * r26095; double r26097 = sqrt(r26096); double r26098 = r26072 * r26097; double r26099 = r26092 ? r26094 : r26098; double r26100 = r26083 ? r26090 : r26099; double r26101 = r26071 ? r26081 : r26100; return r26101; } 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 r26102, r26103, r26104, r26105, r26106, r26107, r26108, r26109, r26110, r26111, r26112, r26113; void setup_mpfr_f_im() { mpfr_set_default_prec(3408); mpfr_init_set_str(r26102, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r26103, "2.0", 10, MPFR_RNDN); mpfr_init(r26104); mpfr_init(r26105); mpfr_init(r26106); mpfr_init(r26107); mpfr_init(r26108); mpfr_init(r26109); mpfr_init(r26110); mpfr_init(r26111); mpfr_init(r26112); mpfr_init(r26113); } double f_im(double re, double im) { ; ; mpfr_set_d(r26104, re, MPFR_RNDN); mpfr_mul(r26105, r26104, r26104, MPFR_RNDN); mpfr_set_d(r26106, im, MPFR_RNDN); mpfr_mul(r26107, r26106, r26106, MPFR_RNDN); mpfr_add(r26108, r26105, r26107, MPFR_RNDN); mpfr_sqrt(r26109, r26108, MPFR_RNDN); mpfr_add(r26110, r26109, r26104, MPFR_RNDN); mpfr_mul(r26111, r26103, r26110, MPFR_RNDN); mpfr_sqrt(r26112, r26111, MPFR_RNDN); mpfr_mul(r26113, r26102, r26112, MPFR_RNDN); return mpfr_get_d(r26113, MPFR_RNDN); } static mpfr_t r26114, r26115, r26116, r26117, r26118, r26119, r26120, r26121, r26122, r26123, r26124, r26125, r26126, r26127, r26128, r26129, r26130, r26131, r26132, r26133, r26134, r26135, r26136, r26137, r26138, r26139, r26140, r26141, r26142, r26143, r26144, r26145, r26146, r26147, r26148, r26149, r26150; void setup_mpfr_f_fm() { mpfr_set_default_prec(3408); mpfr_init_set_str(r26114, "2.0", 10, MPFR_RNDN); mpfr_init(r26115); mpfr_init(r26116); mpfr_init(r26117); mpfr_init(r26118); mpfr_init_set_str(r26119, "-1.0827629684787965e-204", 10, MPFR_RNDN); mpfr_init(r26120); mpfr_init_set_str(r26121, "0.5", 10, MPFR_RNDN); mpfr_init(r26122); mpfr_init(r26123); mpfr_init(r26124); mpfr_init(r26125); mpfr_init(r26126); mpfr_init(r26127); mpfr_init(r26128); mpfr_init(r26129); mpfr_init(r26130); mpfr_init_set_str(r26131, "5.047375239993434e+166", 10, MPFR_RNDN); mpfr_init(r26132); mpfr_init(r26133); mpfr_init(r26134); mpfr_init(r26135); mpfr_init(r26136); mpfr_init(r26137); mpfr_init(r26138); mpfr_init(r26139); mpfr_init_set_str(r26140, "2.1298540120048136e+187", 10, MPFR_RNDN); mpfr_init(r26141); mpfr_init(r26142); mpfr_init(r26143); mpfr_init(r26144); mpfr_init(r26145); mpfr_init(r26146); mpfr_init(r26147); mpfr_init(r26148); mpfr_init(r26149); mpfr_init(r26150); } double f_fm(double re, double im) { ; mpfr_set_d(r26115, im, MPFR_RNDN); mpfr_set_d(r26116, re, MPFR_RNDN); mpfr_add(r26117, r26115, r26116, MPFR_RNDN); mpfr_mul(r26118, r26114, r26117, MPFR_RNDN); ; mpfr_set_si(r26120, mpfr_cmp(r26118, r26119) <= 0, MPFR_RNDN); ; mpfr_mul(r26122, r26115, r26115, MPFR_RNDN); mpfr_mul(r26123, r26116, r26116, MPFR_RNDN); mpfr_add(r26124, r26123, r26122, MPFR_RNDN); mpfr_sqrt(r26125, r26124, MPFR_RNDN); mpfr_sub(r26126, r26125, r26116, MPFR_RNDN); mpfr_div(r26127, r26122, r26126, MPFR_RNDN); mpfr_mul(r26128, r26114, r26127, MPFR_RNDN); mpfr_sqrt(r26129, r26128, MPFR_RNDN); mpfr_mul(r26130, r26121, r26129, MPFR_RNDN); ; mpfr_set_si(r26132, mpfr_cmp(r26118, r26131) <= 0, MPFR_RNDN); mpfr_cbrt(r26133, r26125, MPFR_RNDN); mpfr_mul(r26134, r26133, r26133, MPFR_RNDN); mpfr_mul(r26135, r26134, r26133, MPFR_RNDN); mpfr_add(r26136, r26135, r26116, MPFR_RNDN); mpfr_mul(r26137, r26114, r26136, MPFR_RNDN); mpfr_sqrt(r26138, r26137, MPFR_RNDN); mpfr_mul(r26139, r26121, r26138, MPFR_RNDN); ; mpfr_set_si(r26141, mpfr_cmp(r26118, r26140) <= 0, MPFR_RNDN); mpfr_sqrt(r26142, r26118, MPFR_RNDN); mpfr_mul(r26143, r26121, r26142, MPFR_RNDN); mpfr_add(r26144, r26116, r26116, MPFR_RNDN); mpfr_mul(r26145, r26114, r26144, MPFR_RNDN); mpfr_sqrt(r26146, r26145, MPFR_RNDN); mpfr_mul(r26147, r26121, r26146, MPFR_RNDN); if (mpfr_get_si(r26141, MPFR_RNDN)) { mpfr_set(r26148, r26143, MPFR_RNDN); } else { mpfr_set(r26148, r26147, MPFR_RNDN); }; if (mpfr_get_si(r26132, MPFR_RNDN)) { mpfr_set(r26149, r26139, MPFR_RNDN); } else { mpfr_set(r26149, r26148, MPFR_RNDN); }; if (mpfr_get_si(r26120, MPFR_RNDN)) { mpfr_set(r26150, r26130, MPFR_RNDN); } else { mpfr_set(r26150, r26149, MPFR_RNDN); }; return mpfr_get_d(r26150, MPFR_RNDN); } static mpfr_t r26151, r26152, r26153, r26154, r26155, r26156, r26157, r26158, r26159, r26160, r26161, r26162, r26163, r26164, r26165, r26166, r26167, r26168, r26169, r26170, r26171, r26172, r26173, r26174, r26175, r26176, r26177, r26178, r26179, r26180, r26181, r26182, r26183, r26184, r26185, r26186, r26187; void setup_mpfr_f_dm() { mpfr_set_default_prec(3408); mpfr_init_set_str(r26151, "2.0", 10, MPFR_RNDN); mpfr_init(r26152); mpfr_init(r26153); mpfr_init(r26154); mpfr_init(r26155); mpfr_init_set_str(r26156, "-1.0827629684787965e-204", 10, MPFR_RNDN); mpfr_init(r26157); mpfr_init_set_str(r26158, "0.5", 10, MPFR_RNDN); mpfr_init(r26159); mpfr_init(r26160); mpfr_init(r26161); mpfr_init(r26162); mpfr_init(r26163); mpfr_init(r26164); mpfr_init(r26165); mpfr_init(r26166); mpfr_init(r26167); mpfr_init_set_str(r26168, "5.047375239993434e+166", 10, MPFR_RNDN); mpfr_init(r26169); mpfr_init(r26170); mpfr_init(r26171); mpfr_init(r26172); mpfr_init(r26173); mpfr_init(r26174); mpfr_init(r26175); mpfr_init(r26176); mpfr_init_set_str(r26177, "2.1298540120048136e+187", 10, MPFR_RNDN); mpfr_init(r26178); mpfr_init(r26179); mpfr_init(r26180); mpfr_init(r26181); mpfr_init(r26182); mpfr_init(r26183); mpfr_init(r26184); mpfr_init(r26185); mpfr_init(r26186); mpfr_init(r26187); } double f_dm(double re, double im) { ; mpfr_set_d(r26152, im, MPFR_RNDN); mpfr_set_d(r26153, re, MPFR_RNDN); mpfr_add(r26154, r26152, r26153, MPFR_RNDN); mpfr_mul(r26155, r26151, r26154, MPFR_RNDN); ; mpfr_set_si(r26157, mpfr_cmp(r26155, r26156) <= 0, MPFR_RNDN); ; mpfr_mul(r26159, r26152, r26152, MPFR_RNDN); mpfr_mul(r26160, r26153, r26153, MPFR_RNDN); mpfr_add(r26161, r26160, r26159, MPFR_RNDN); mpfr_sqrt(r26162, r26161, MPFR_RNDN); mpfr_sub(r26163, r26162, r26153, MPFR_RNDN); mpfr_div(r26164, r26159, r26163, MPFR_RNDN); mpfr_mul(r26165, r26151, r26164, MPFR_RNDN); mpfr_sqrt(r26166, r26165, MPFR_RNDN); mpfr_mul(r26167, r26158, r26166, MPFR_RNDN); ; mpfr_set_si(r26169, mpfr_cmp(r26155, r26168) <= 0, MPFR_RNDN); mpfr_cbrt(r26170, r26162, MPFR_RNDN); mpfr_mul(r26171, r26170, r26170, MPFR_RNDN); mpfr_mul(r26172, r26171, r26170, MPFR_RNDN); mpfr_add(r26173, r26172, r26153, MPFR_RNDN); mpfr_mul(r26174, r26151, r26173, MPFR_RNDN); mpfr_sqrt(r26175, r26174, MPFR_RNDN); mpfr_mul(r26176, r26158, r26175, MPFR_RNDN); ; mpfr_set_si(r26178, mpfr_cmp(r26155, r26177) <= 0, MPFR_RNDN); mpfr_sqrt(r26179, r26155, MPFR_RNDN); mpfr_mul(r26180, r26158, r26179, MPFR_RNDN); mpfr_add(r26181, r26153, r26153, MPFR_RNDN); mpfr_mul(r26182, r26151, r26181, MPFR_RNDN); mpfr_sqrt(r26183, r26182, MPFR_RNDN); mpfr_mul(r26184, r26158, r26183, MPFR_RNDN); if (mpfr_get_si(r26178, MPFR_RNDN)) { mpfr_set(r26185, r26180, MPFR_RNDN); } else { mpfr_set(r26185, r26184, MPFR_RNDN); }; if (mpfr_get_si(r26169, MPFR_RNDN)) { mpfr_set(r26186, r26176, MPFR_RNDN); } else { mpfr_set(r26186, r26185, MPFR_RNDN); }; if (mpfr_get_si(r26157, MPFR_RNDN)) { mpfr_set(r26187, r26167, MPFR_RNDN); } else { mpfr_set(r26187, r26186, MPFR_RNDN); }; return mpfr_get_d(r26187, MPFR_RNDN); }