#include #include #include #include #include char *name = "Complex division, real part"; double f_if(float a, float b, float c, float d) { float r26182 = a; float r26183 = c; float r26184 = r26182 * r26183; float r26185 = b; float r26186 = d; float r26187 = r26185 * r26186; float r26188 = r26184 + r26187; float r26189 = r26183 * r26183; float r26190 = r26186 * r26186; float r26191 = r26189 + r26190; float r26192 = r26188 / r26191; return r26192; } double f_id(double a, double b, double c, double d) { double r26193 = a; double r26194 = c; double r26195 = r26193 * r26194; double r26196 = b; double r26197 = d; double r26198 = r26196 * r26197; double r26199 = r26195 + r26198; double r26200 = r26194 * r26194; double r26201 = r26197 * r26197; double r26202 = r26200 + r26201; double r26203 = r26199 / r26202; return r26203; } double f_of(float a, float b, float c, float d) { float r26204 = 1; float r26205 = c; float r26206 = d; float r26207 = hypot(r26205, r26206); float r26208 = r26204 / r26207; float r26209 = b; float r26210 = a; float r26211 = r26205 * r26210; float r26212 = fma(r26209, r26206, r26211); float r26213 = r26212 / r26207; float r26214 = cbrt(r26213); float r26215 = r26214 * r26214; float r26216 = r26215 * r26214; float r26217 = r26208 * r26216; float r26218 = -inf.0; bool r26219 = r26217 <= r26218; float r26220 = -r26209; float r26221 = r26220 / r26207; float r26222 = 1.7976931348623157e+308; bool r26223 = r26217 <= r26222; float r26224 = r26210 * r26205; float r26225 = fma(r26206, r26209, r26224); float r26226 = r26225 / r26207; float r26227 = r26226 / r26207; float r26228 = r26210 / r26207; float r26229 = r26223 ? r26227 : r26228; float r26230 = r26219 ? r26221 : r26229; return r26230; } double f_od(double a, double b, double c, double d) { double r26231 = 1; double r26232 = c; double r26233 = d; double r26234 = hypot(r26232, r26233); double r26235 = r26231 / r26234; double r26236 = b; double r26237 = a; double r26238 = r26232 * r26237; double r26239 = fma(r26236, r26233, r26238); double r26240 = r26239 / r26234; double r26241 = cbrt(r26240); double r26242 = r26241 * r26241; double r26243 = r26242 * r26241; double r26244 = r26235 * r26243; double r26245 = -inf.0; bool r26246 = r26244 <= r26245; double r26247 = -r26236; double r26248 = r26247 / r26234; double r26249 = 1.7976931348623157e+308; bool r26250 = r26244 <= r26249; double r26251 = r26237 * r26232; double r26252 = fma(r26233, r26236, r26251); double r26253 = r26252 / r26234; double r26254 = r26253 / r26234; double r26255 = r26237 / r26234; double r26256 = r26250 ? r26254 : r26255; double r26257 = r26246 ? r26248 : r26256; return r26257; } 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 r26258, r26259, r26260, r26261, r26262, r26263, r26264, r26265, r26266, r26267, r26268; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26258); mpfr_init(r26259); mpfr_init(r26260); mpfr_init(r26261); mpfr_init(r26262); mpfr_init(r26263); mpfr_init(r26264); mpfr_init(r26265); mpfr_init(r26266); mpfr_init(r26267); mpfr_init(r26268); } double f_im(double a, double b, double c, double d) { mpfr_set_d(r26258, a, MPFR_RNDN); mpfr_set_d(r26259, c, MPFR_RNDN); mpfr_mul(r26260, r26258, r26259, MPFR_RNDN); mpfr_set_d(r26261, b, MPFR_RNDN); mpfr_set_d(r26262, d, MPFR_RNDN); mpfr_mul(r26263, r26261, r26262, MPFR_RNDN); mpfr_add(r26264, r26260, r26263, MPFR_RNDN); mpfr_mul(r26265, r26259, r26259, MPFR_RNDN); mpfr_mul(r26266, r26262, r26262, MPFR_RNDN); mpfr_add(r26267, r26265, r26266, MPFR_RNDN); mpfr_div(r26268, r26264, r26267, MPFR_RNDN); return mpfr_get_d(r26268, MPFR_RNDN); } static mpfr_t r26269, r26270, r26271, r26272, r26273, r26274, r26275, r26276, r26277, r26278, r26279, r26280, r26281, r26282, r26283, r26284, r26285, r26286, r26287, r26288, r26289, r26290, r26291, r26292, r26293, r26294, r26295; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26269, "1", 10, MPFR_RNDN); mpfr_init(r26270); mpfr_init(r26271); mpfr_init(r26272); mpfr_init(r26273); mpfr_init(r26274); mpfr_init(r26275); mpfr_init(r26276); mpfr_init(r26277); mpfr_init(r26278); mpfr_init(r26279); mpfr_init(r26280); mpfr_init(r26281); mpfr_init(r26282); mpfr_init_set_str(r26283, "-inf.0", 10, MPFR_RNDN); mpfr_init(r26284); mpfr_init(r26285); mpfr_init(r26286); mpfr_init_set_str(r26287, "1.7976931348623157e+308", 10, MPFR_RNDN); mpfr_init(r26288); mpfr_init(r26289); mpfr_init(r26290); mpfr_init(r26291); mpfr_init(r26292); mpfr_init(r26293); mpfr_init(r26294); mpfr_init(r26295); } double f_fm(double a, double b, double c, double d) { ; mpfr_set_d(r26270, c, MPFR_RNDN); mpfr_set_d(r26271, d, MPFR_RNDN); mpfr_hypot(r26272, r26270, r26271, MPFR_RNDN); mpfr_div(r26273, r26269, r26272, MPFR_RNDN); mpfr_set_d(r26274, b, MPFR_RNDN); mpfr_set_d(r26275, a, MPFR_RNDN); mpfr_mul(r26276, r26270, r26275, MPFR_RNDN); mpfr_fma(r26277, r26274, r26271, r26276, MPFR_RNDN); mpfr_div(r26278, r26277, r26272, MPFR_RNDN); mpfr_cbrt(r26279, r26278, MPFR_RNDN); mpfr_mul(r26280, r26279, r26279, MPFR_RNDN); mpfr_mul(r26281, r26280, r26279, MPFR_RNDN); mpfr_mul(r26282, r26273, r26281, MPFR_RNDN); ; mpfr_set_si(r26284, mpfr_cmp(r26282, r26283) <= 0, MPFR_RNDN); mpfr_neg(r26285, r26274, MPFR_RNDN); mpfr_div(r26286, r26285, r26272, MPFR_RNDN); ; mpfr_set_si(r26288, mpfr_cmp(r26282, r26287) <= 0, MPFR_RNDN); mpfr_mul(r26289, r26275, r26270, MPFR_RNDN); mpfr_fma(r26290, r26271, r26274, r26289, MPFR_RNDN); mpfr_div(r26291, r26290, r26272, MPFR_RNDN); mpfr_div(r26292, r26291, r26272, MPFR_RNDN); mpfr_div(r26293, r26275, r26272, MPFR_RNDN); if (mpfr_get_si(r26288, MPFR_RNDN)) { mpfr_set(r26294, r26292, MPFR_RNDN); } else { mpfr_set(r26294, r26293, MPFR_RNDN); }; if (mpfr_get_si(r26284, MPFR_RNDN)) { mpfr_set(r26295, r26286, MPFR_RNDN); } else { mpfr_set(r26295, r26294, MPFR_RNDN); }; return mpfr_get_d(r26295, MPFR_RNDN); } static mpfr_t r26296, r26297, r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321, r26322; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26296, "1", 10, MPFR_RNDN); mpfr_init(r26297); mpfr_init(r26298); mpfr_init(r26299); mpfr_init(r26300); mpfr_init(r26301); mpfr_init(r26302); mpfr_init(r26303); mpfr_init(r26304); mpfr_init(r26305); mpfr_init(r26306); mpfr_init(r26307); mpfr_init(r26308); mpfr_init(r26309); mpfr_init_set_str(r26310, "-inf.0", 10, MPFR_RNDN); mpfr_init(r26311); mpfr_init(r26312); mpfr_init(r26313); mpfr_init_set_str(r26314, "1.7976931348623157e+308", 10, MPFR_RNDN); mpfr_init(r26315); mpfr_init(r26316); mpfr_init(r26317); mpfr_init(r26318); mpfr_init(r26319); mpfr_init(r26320); mpfr_init(r26321); mpfr_init(r26322); } double f_dm(double a, double b, double c, double d) { ; mpfr_set_d(r26297, c, MPFR_RNDN); mpfr_set_d(r26298, d, MPFR_RNDN); mpfr_hypot(r26299, r26297, r26298, MPFR_RNDN); mpfr_div(r26300, r26296, r26299, MPFR_RNDN); mpfr_set_d(r26301, b, MPFR_RNDN); mpfr_set_d(r26302, a, MPFR_RNDN); mpfr_mul(r26303, r26297, r26302, MPFR_RNDN); mpfr_fma(r26304, r26301, r26298, r26303, MPFR_RNDN); mpfr_div(r26305, r26304, r26299, MPFR_RNDN); mpfr_cbrt(r26306, r26305, MPFR_RNDN); mpfr_mul(r26307, r26306, r26306, MPFR_RNDN); mpfr_mul(r26308, r26307, r26306, MPFR_RNDN); mpfr_mul(r26309, r26300, r26308, MPFR_RNDN); ; mpfr_set_si(r26311, mpfr_cmp(r26309, r26310) <= 0, MPFR_RNDN); mpfr_neg(r26312, r26301, MPFR_RNDN); mpfr_div(r26313, r26312, r26299, MPFR_RNDN); ; mpfr_set_si(r26315, mpfr_cmp(r26309, r26314) <= 0, MPFR_RNDN); mpfr_mul(r26316, r26302, r26297, MPFR_RNDN); mpfr_fma(r26317, r26298, r26301, r26316, MPFR_RNDN); mpfr_div(r26318, r26317, r26299, MPFR_RNDN); mpfr_div(r26319, r26318, r26299, MPFR_RNDN); mpfr_div(r26320, r26302, r26299, MPFR_RNDN); if (mpfr_get_si(r26315, MPFR_RNDN)) { mpfr_set(r26321, r26319, MPFR_RNDN); } else { mpfr_set(r26321, r26320, MPFR_RNDN); }; if (mpfr_get_si(r26311, MPFR_RNDN)) { mpfr_set(r26322, r26313, MPFR_RNDN); } else { mpfr_set(r26322, r26321, MPFR_RNDN); }; return mpfr_get_d(r26322, MPFR_RNDN); }