#include #include #include #include #include char *name = "Complex division, imag part"; double f_if(float a, float b, float c, float d) { float r17007 = b; float r17008 = c; float r17009 = r17007 * r17008; float r17010 = a; float r17011 = d; float r17012 = r17010 * r17011; float r17013 = r17009 - r17012; float r17014 = r17008 * r17008; float r17015 = r17011 * r17011; float r17016 = r17014 + r17015; float r17017 = r17013 / r17016; return r17017; } double f_id(double a, double b, double c, double d) { double r17018 = b; double r17019 = c; double r17020 = r17018 * r17019; double r17021 = a; double r17022 = d; double r17023 = r17021 * r17022; double r17024 = r17020 - r17023; double r17025 = r17019 * r17019; double r17026 = r17022 * r17022; double r17027 = r17025 + r17026; double r17028 = r17024 / r17027; return r17028; } double f_of(float a, float b, float c, float d) { float r17029 = c; float r17030 = -4166742528.0f; bool r17031 = r17029 <= r17030; float r17032 = b; float r17033 = r17032 / r17029; float r17034 = d; float r17035 = r17034 / r17029; float r17036 = a; float r17037 = r17036 / r17029; float r17038 = r17035 * r17037; float r17039 = r17033 - r17038; float r17040 = 6144087082139648.0f; bool r17041 = r17029 <= r17040; float r17042 = r17032 * r17029; float r17043 = r17036 * r17034; float r17044 = r17042 - r17043; float r17045 = r17029 * r17029; float r17046 = r17034 * r17034; float r17047 = r17045 + r17046; float r17048 = sqrt(r17047); float r17049 = r17048 * r17048; float r17050 = r17044 / r17049; float r17051 = r17041 ? r17050 : r17039; float r17052 = r17031 ? r17039 : r17051; return r17052; } double f_od(double a, double b, double c, double d) { double r17053 = c; double r17054 = -4166742528.0; bool r17055 = r17053 <= r17054; double r17056 = b; double r17057 = r17056 / r17053; double r17058 = d; double r17059 = r17058 / r17053; double r17060 = a; double r17061 = r17060 / r17053; double r17062 = r17059 * r17061; double r17063 = r17057 - r17062; double r17064 = 6144087082139648.0; bool r17065 = r17053 <= r17064; double r17066 = r17056 * r17053; double r17067 = r17060 * r17058; double r17068 = r17066 - r17067; double r17069 = r17053 * r17053; double r17070 = r17058 * r17058; double r17071 = r17069 + r17070; double r17072 = sqrt(r17071); double r17073 = r17072 * r17072; double r17074 = r17068 / r17073; double r17075 = r17065 ? r17074 : r17063; double r17076 = r17055 ? r17063 : r17075; return r17076; } 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 r17077, r17078, r17079, r17080, r17081, r17082, r17083, r17084, r17085, r17086, r17087; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17077); mpfr_init(r17078); mpfr_init(r17079); mpfr_init(r17080); mpfr_init(r17081); mpfr_init(r17082); mpfr_init(r17083); mpfr_init(r17084); mpfr_init(r17085); mpfr_init(r17086); mpfr_init(r17087); } double f_im(double a, double b, double c, double d) { mpfr_set_d(r17077, b, MPFR_RNDN); mpfr_set_d(r17078, c, MPFR_RNDN); mpfr_mul(r17079, r17077, r17078, MPFR_RNDN); mpfr_set_d(r17080, a, MPFR_RNDN); mpfr_set_d(r17081, d, MPFR_RNDN); mpfr_mul(r17082, r17080, r17081, MPFR_RNDN); mpfr_sub(r17083, r17079, r17082, MPFR_RNDN); mpfr_sqr(r17084, r17078, MPFR_RNDN); mpfr_sqr(r17085, r17081, MPFR_RNDN); mpfr_add(r17086, r17084, r17085, MPFR_RNDN); mpfr_div(r17087, r17083, r17086, MPFR_RNDN); return mpfr_get_d(r17087, MPFR_RNDN); } static mpfr_t r17088, r17089, r17090, r17091, r17092, r17093, r17094, r17095, r17096, r17097, r17098, r17099, r17100, r17101, r17102, r17103, r17104, r17105, r17106, r17107, r17108, r17109, r17110, r17111; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17088); mpfr_init_set_str(r17089, "-4.1667425f+09", 10, MPFR_RNDN); mpfr_init(r17090); mpfr_init(r17091); mpfr_init(r17092); mpfr_init(r17093); mpfr_init(r17094); mpfr_init(r17095); mpfr_init(r17096); mpfr_init(r17097); mpfr_init(r17098); mpfr_init_set_str(r17099, "6.144087f+15", 10, MPFR_RNDN); mpfr_init(r17100); mpfr_init(r17101); mpfr_init(r17102); mpfr_init(r17103); mpfr_init(r17104); mpfr_init(r17105); mpfr_init(r17106); mpfr_init(r17107); mpfr_init(r17108); mpfr_init(r17109); mpfr_init(r17110); mpfr_init(r17111); } double f_fm(double a, double b, double c, double d) { mpfr_set_d(r17088, c, MPFR_RNDN); ; mpfr_set_si(r17090, mpfr_cmp(r17088, r17089) <= 0, MPFR_RNDN); mpfr_set_d(r17091, b, MPFR_RNDN); mpfr_div(r17092, r17091, r17088, MPFR_RNDN); mpfr_set_d(r17093, d, MPFR_RNDN); mpfr_div(r17094, r17093, r17088, MPFR_RNDN); mpfr_set_d(r17095, a, MPFR_RNDN); mpfr_div(r17096, r17095, r17088, MPFR_RNDN); mpfr_mul(r17097, r17094, r17096, MPFR_RNDN); mpfr_sub(r17098, r17092, r17097, MPFR_RNDN); ; mpfr_set_si(r17100, mpfr_cmp(r17088, r17099) <= 0, MPFR_RNDN); mpfr_mul(r17101, r17091, r17088, MPFR_RNDN); mpfr_mul(r17102, r17095, r17093, MPFR_RNDN); mpfr_sub(r17103, r17101, r17102, MPFR_RNDN); mpfr_sqr(r17104, r17088, MPFR_RNDN); mpfr_sqr(r17105, r17093, MPFR_RNDN); mpfr_add(r17106, r17104, r17105, MPFR_RNDN); mpfr_sqrt(r17107, r17106, MPFR_RNDN); mpfr_sqr(r17108, r17107, MPFR_RNDN); mpfr_div(r17109, r17103, r17108, MPFR_RNDN); if (mpfr_get_si(r17100, MPFR_RNDN)) { mpfr_set(r17110, r17109, MPFR_RNDN); } else { mpfr_set(r17110, r17098, MPFR_RNDN); }; if (mpfr_get_si(r17090, MPFR_RNDN)) { mpfr_set(r17111, r17098, MPFR_RNDN); } else { mpfr_set(r17111, r17110, MPFR_RNDN); }; return mpfr_get_d(r17111, MPFR_RNDN); } static mpfr_t r17112, r17113, r17114, r17115, r17116, r17117, r17118, r17119, r17120, r17121, r17122, r17123, r17124, r17125, r17126, r17127, r17128, r17129, r17130, r17131, r17132, r17133, r17134, r17135; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17112); mpfr_init_set_str(r17113, "-4.1667425f+09", 10, MPFR_RNDN); mpfr_init(r17114); mpfr_init(r17115); mpfr_init(r17116); mpfr_init(r17117); mpfr_init(r17118); mpfr_init(r17119); mpfr_init(r17120); mpfr_init(r17121); mpfr_init(r17122); mpfr_init_set_str(r17123, "6.144087f+15", 10, MPFR_RNDN); mpfr_init(r17124); mpfr_init(r17125); mpfr_init(r17126); mpfr_init(r17127); mpfr_init(r17128); mpfr_init(r17129); mpfr_init(r17130); mpfr_init(r17131); mpfr_init(r17132); mpfr_init(r17133); mpfr_init(r17134); mpfr_init(r17135); } double f_dm(double a, double b, double c, double d) { mpfr_set_d(r17112, c, MPFR_RNDN); ; mpfr_set_si(r17114, mpfr_cmp(r17112, r17113) <= 0, MPFR_RNDN); mpfr_set_d(r17115, b, MPFR_RNDN); mpfr_div(r17116, r17115, r17112, MPFR_RNDN); mpfr_set_d(r17117, d, MPFR_RNDN); mpfr_div(r17118, r17117, r17112, MPFR_RNDN); mpfr_set_d(r17119, a, MPFR_RNDN); mpfr_div(r17120, r17119, r17112, MPFR_RNDN); mpfr_mul(r17121, r17118, r17120, MPFR_RNDN); mpfr_sub(r17122, r17116, r17121, MPFR_RNDN); ; mpfr_set_si(r17124, mpfr_cmp(r17112, r17123) <= 0, MPFR_RNDN); mpfr_mul(r17125, r17115, r17112, MPFR_RNDN); mpfr_mul(r17126, r17119, r17117, MPFR_RNDN); mpfr_sub(r17127, r17125, r17126, MPFR_RNDN); mpfr_sqr(r17128, r17112, MPFR_RNDN); mpfr_sqr(r17129, r17117, MPFR_RNDN); mpfr_add(r17130, r17128, r17129, MPFR_RNDN); mpfr_sqrt(r17131, r17130, MPFR_RNDN); mpfr_sqr(r17132, r17131, MPFR_RNDN); mpfr_div(r17133, r17127, r17132, MPFR_RNDN); if (mpfr_get_si(r17124, MPFR_RNDN)) { mpfr_set(r17134, r17133, MPFR_RNDN); } else { mpfr_set(r17134, r17122, MPFR_RNDN); }; if (mpfr_get_si(r17114, MPFR_RNDN)) { mpfr_set(r17135, r17122, MPFR_RNDN); } else { mpfr_set(r17135, r17134, MPFR_RNDN); }; return mpfr_get_d(r17135, MPFR_RNDN); }