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