#include #include #include #include #include char *name = "_divideComplex, imaginary part"; double f_if(float x_re, float x_im, float y_re, float y_im) { float r22153 = x_im; float r22154 = y_re; float r22155 = r22153 * r22154; float r22156 = x_re; float r22157 = y_im; float r22158 = r22156 * r22157; float r22159 = r22155 - r22158; float r22160 = r22154 * r22154; float r22161 = r22157 * r22157; float r22162 = r22160 + r22161; float r22163 = r22159 / r22162; return r22163; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r22164 = x_im; double r22165 = y_re; double r22166 = r22164 * r22165; double r22167 = x_re; double r22168 = y_im; double r22169 = r22167 * r22168; double r22170 = r22166 - r22169; double r22171 = r22165 * r22165; double r22172 = r22168 * r22168; double r22173 = r22171 + r22172; double r22174 = r22170 / r22173; return r22174; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r22175 = y_re; float r22176 = -1.694155349417558e+107; bool r22177 = r22175 <= r22176; float r22178 = x_im; float r22179 = -r22178; float r22180 = y_im; float r22181 = hypot(r22175, r22180); float r22182 = r22179 / r22181; float r22183 = r22178 * r22175; float r22184 = x_re; float r22185 = r22184 * r22180; float r22186 = r22183 - r22185; float r22187 = r22186 / r22181; float r22188 = r22187 / r22181; float r22189 = r22177 ? r22182 : r22188; return r22189; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r22190 = y_re; double r22191 = -1.694155349417558e+107; bool r22192 = r22190 <= r22191; double r22193 = x_im; double r22194 = -r22193; double r22195 = y_im; double r22196 = hypot(r22190, r22195); double r22197 = r22194 / r22196; double r22198 = r22193 * r22190; double r22199 = x_re; double r22200 = r22199 * r22195; double r22201 = r22198 - r22200; double r22202 = r22201 / r22196; double r22203 = r22202 / r22196; double r22204 = r22192 ? r22197 : r22203; return r22204; } 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 r22205, r22206, r22207, r22208, r22209, r22210, r22211, r22212, r22213, r22214, r22215; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22205); mpfr_init(r22206); mpfr_init(r22207); mpfr_init(r22208); mpfr_init(r22209); mpfr_init(r22210); mpfr_init(r22211); mpfr_init(r22212); mpfr_init(r22213); mpfr_init(r22214); mpfr_init(r22215); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r22205, x_im, MPFR_RNDN); mpfr_set_d(r22206, y_re, MPFR_RNDN); mpfr_mul(r22207, r22205, r22206, MPFR_RNDN); mpfr_set_d(r22208, x_re, MPFR_RNDN); mpfr_set_d(r22209, y_im, MPFR_RNDN); mpfr_mul(r22210, r22208, r22209, MPFR_RNDN); mpfr_sub(r22211, r22207, r22210, MPFR_RNDN); mpfr_mul(r22212, r22206, r22206, MPFR_RNDN); mpfr_mul(r22213, r22209, r22209, MPFR_RNDN); mpfr_add(r22214, r22212, r22213, MPFR_RNDN); mpfr_div(r22215, r22211, r22214, MPFR_RNDN); return mpfr_get_d(r22215, MPFR_RNDN); } static mpfr_t r22216, r22217, r22218, r22219, r22220, r22221, r22222, r22223, r22224, r22225, r22226, r22227, r22228, r22229, r22230; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22216); mpfr_init_set_str(r22217, "-1.694155349417558e+107", 10, MPFR_RNDN); mpfr_init(r22218); mpfr_init(r22219); mpfr_init(r22220); mpfr_init(r22221); mpfr_init(r22222); mpfr_init(r22223); mpfr_init(r22224); mpfr_init(r22225); mpfr_init(r22226); mpfr_init(r22227); mpfr_init(r22228); mpfr_init(r22229); mpfr_init(r22230); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r22216, y_re, MPFR_RNDN); ; mpfr_set_si(r22218, mpfr_cmp(r22216, r22217) <= 0, MPFR_RNDN); mpfr_set_d(r22219, x_im, MPFR_RNDN); mpfr_neg(r22220, r22219, MPFR_RNDN); mpfr_set_d(r22221, y_im, MPFR_RNDN); mpfr_hypot(r22222, r22216, r22221, MPFR_RNDN); mpfr_div(r22223, r22220, r22222, MPFR_RNDN); mpfr_mul(r22224, r22219, r22216, MPFR_RNDN); mpfr_set_d(r22225, x_re, MPFR_RNDN); mpfr_mul(r22226, r22225, r22221, MPFR_RNDN); mpfr_sub(r22227, r22224, r22226, MPFR_RNDN); mpfr_div(r22228, r22227, r22222, MPFR_RNDN); mpfr_div(r22229, r22228, r22222, MPFR_RNDN); if (mpfr_get_si(r22218, MPFR_RNDN)) { mpfr_set(r22230, r22223, MPFR_RNDN); } else { mpfr_set(r22230, r22229, MPFR_RNDN); }; return mpfr_get_d(r22230, MPFR_RNDN); } static mpfr_t r22231, r22232, r22233, r22234, r22235, r22236, r22237, r22238, r22239, r22240, r22241, r22242, r22243, r22244, r22245; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22231); mpfr_init_set_str(r22232, "-1.694155349417558e+107", 10, MPFR_RNDN); mpfr_init(r22233); mpfr_init(r22234); mpfr_init(r22235); mpfr_init(r22236); mpfr_init(r22237); mpfr_init(r22238); mpfr_init(r22239); mpfr_init(r22240); mpfr_init(r22241); mpfr_init(r22242); mpfr_init(r22243); mpfr_init(r22244); mpfr_init(r22245); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r22231, y_re, MPFR_RNDN); ; mpfr_set_si(r22233, mpfr_cmp(r22231, r22232) <= 0, MPFR_RNDN); mpfr_set_d(r22234, x_im, MPFR_RNDN); mpfr_neg(r22235, r22234, MPFR_RNDN); mpfr_set_d(r22236, y_im, MPFR_RNDN); mpfr_hypot(r22237, r22231, r22236, MPFR_RNDN); mpfr_div(r22238, r22235, r22237, MPFR_RNDN); mpfr_mul(r22239, r22234, r22231, MPFR_RNDN); mpfr_set_d(r22240, x_re, MPFR_RNDN); mpfr_mul(r22241, r22240, r22236, MPFR_RNDN); mpfr_sub(r22242, r22239, r22241, MPFR_RNDN); mpfr_div(r22243, r22242, r22237, MPFR_RNDN); mpfr_div(r22244, r22243, r22237, MPFR_RNDN); if (mpfr_get_si(r22233, MPFR_RNDN)) { mpfr_set(r22245, r22238, MPFR_RNDN); } else { mpfr_set(r22245, r22244, MPFR_RNDN); }; return mpfr_get_d(r22245, MPFR_RNDN); }