#include #include #include #include #include char *name = "powComplex, real part"; double f_if(float x_re, float x_im, float y_re, float y_im) { float r21237 = x_re; float r21238 = r21237 * r21237; float r21239 = x_im; float r21240 = r21239 * r21239; float r21241 = r21238 + r21240; float r21242 = sqrt(r21241); float r21243 = log(r21242); float r21244 = y_re; float r21245 = r21243 * r21244; float r21246 = atan2(r21239, r21237); float r21247 = y_im; float r21248 = r21246 * r21247; float r21249 = r21245 - r21248; float r21250 = exp(r21249); float r21251 = r21243 * r21247; float r21252 = r21246 * r21244; float r21253 = r21251 + r21252; float r21254 = cos(r21253); float r21255 = r21250 * r21254; return r21255; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r21256 = x_re; double r21257 = r21256 * r21256; double r21258 = x_im; double r21259 = r21258 * r21258; double r21260 = r21257 + r21259; double r21261 = sqrt(r21260); double r21262 = log(r21261); double r21263 = y_re; double r21264 = r21262 * r21263; double r21265 = atan2(r21258, r21256); double r21266 = y_im; double r21267 = r21265 * r21266; double r21268 = r21264 - r21267; double r21269 = exp(r21268); double r21270 = r21262 * r21266; double r21271 = r21265 * r21263; double r21272 = r21270 + r21271; double r21273 = cos(r21272); double r21274 = r21269 * r21273; return r21274; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r21275 = x_im; float r21276 = x_re; float r21277 = hypot(r21275, r21276); float r21278 = log(r21277); float r21279 = y_im; float r21280 = atan2(r21275, r21276); float r21281 = y_re; float r21282 = r21280 * r21281; float r21283 = fma(r21278, r21279, r21282); float r21284 = cos(r21283); float r21285 = expm1(r21284); float r21286 = exp(r21285); float r21287 = log(r21286); float r21288 = log1p(r21287); float r21289 = r21279 * r21280; float r21290 = r21278 * r21281; float r21291 = r21289 - r21290; float r21292 = exp(r21291); float r21293 = r21288 / r21292; return r21293; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r21294 = x_im; double r21295 = x_re; double r21296 = hypot(r21294, r21295); double r21297 = log(r21296); double r21298 = y_im; double r21299 = atan2(r21294, r21295); double r21300 = y_re; double r21301 = r21299 * r21300; double r21302 = fma(r21297, r21298, r21301); double r21303 = cos(r21302); double r21304 = expm1(r21303); double r21305 = exp(r21304); double r21306 = log(r21305); double r21307 = log1p(r21306); double r21308 = r21298 * r21299; double r21309 = r21297 * r21300; double r21310 = r21308 - r21309; double r21311 = exp(r21310); double r21312 = r21307 / r21311; return r21312; } 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 r21313, r21314, r21315, r21316, r21317, r21318, r21319, r21320, r21321, r21322, r21323, r21324, r21325, r21326, r21327, r21328, r21329, r21330, r21331; void setup_mpfr_f_im() { mpfr_set_default_prec(4112); mpfr_init(r21313); mpfr_init(r21314); mpfr_init(r21315); mpfr_init(r21316); mpfr_init(r21317); mpfr_init(r21318); mpfr_init(r21319); mpfr_init(r21320); mpfr_init(r21321); mpfr_init(r21322); mpfr_init(r21323); mpfr_init(r21324); mpfr_init(r21325); mpfr_init(r21326); mpfr_init(r21327); mpfr_init(r21328); mpfr_init(r21329); mpfr_init(r21330); mpfr_init(r21331); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21313, x_re, MPFR_RNDN); mpfr_mul(r21314, r21313, r21313, MPFR_RNDN); mpfr_set_d(r21315, x_im, MPFR_RNDN); mpfr_mul(r21316, r21315, r21315, MPFR_RNDN); mpfr_add(r21317, r21314, r21316, MPFR_RNDN); mpfr_sqrt(r21318, r21317, MPFR_RNDN); mpfr_log(r21319, r21318, MPFR_RNDN); mpfr_set_d(r21320, y_re, MPFR_RNDN); mpfr_mul(r21321, r21319, r21320, MPFR_RNDN); mpfr_atan2(r21322, r21315, r21313, MPFR_RNDN); mpfr_set_d(r21323, y_im, MPFR_RNDN); mpfr_mul(r21324, r21322, r21323, MPFR_RNDN); mpfr_sub(r21325, r21321, r21324, MPFR_RNDN); mpfr_exp(r21326, r21325, MPFR_RNDN); mpfr_mul(r21327, r21319, r21323, MPFR_RNDN); mpfr_mul(r21328, r21322, r21320, MPFR_RNDN); mpfr_add(r21329, r21327, r21328, MPFR_RNDN); mpfr_cos(r21330, r21329, MPFR_RNDN); mpfr_mul(r21331, r21326, r21330, MPFR_RNDN); return mpfr_get_d(r21331, MPFR_RNDN); } static mpfr_t r21332, r21333, r21334, r21335, r21336, r21337, r21338, r21339, r21340, r21341, r21342, r21343, r21344, r21345, r21346, r21347, r21348, r21349, r21350; void setup_mpfr_f_fm() { mpfr_set_default_prec(4112); mpfr_init(r21332); mpfr_init(r21333); mpfr_init(r21334); mpfr_init(r21335); mpfr_init(r21336); mpfr_init(r21337); mpfr_init(r21338); mpfr_init(r21339); mpfr_init(r21340); mpfr_init(r21341); mpfr_init(r21342); mpfr_init(r21343); mpfr_init(r21344); mpfr_init(r21345); mpfr_init(r21346); mpfr_init(r21347); mpfr_init(r21348); mpfr_init(r21349); mpfr_init(r21350); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21332, x_im, MPFR_RNDN); mpfr_set_d(r21333, x_re, MPFR_RNDN); mpfr_hypot(r21334, r21332, r21333, MPFR_RNDN); mpfr_log(r21335, r21334, MPFR_RNDN); mpfr_set_d(r21336, y_im, MPFR_RNDN); mpfr_atan2(r21337, r21332, r21333, MPFR_RNDN); mpfr_set_d(r21338, y_re, MPFR_RNDN); mpfr_mul(r21339, r21337, r21338, MPFR_RNDN); mpfr_fma(r21340, r21335, r21336, r21339, MPFR_RNDN); mpfr_cos(r21341, r21340, MPFR_RNDN); mpfr_expm1(r21342, r21341, MPFR_RNDN); mpfr_exp(r21343, r21342, MPFR_RNDN); mpfr_log(r21344, r21343, MPFR_RNDN); mpfr_log1p(r21345, r21344, MPFR_RNDN); mpfr_mul(r21346, r21336, r21337, MPFR_RNDN); mpfr_mul(r21347, r21335, r21338, MPFR_RNDN); mpfr_sub(r21348, r21346, r21347, MPFR_RNDN); mpfr_exp(r21349, r21348, MPFR_RNDN); mpfr_div(r21350, r21345, r21349, MPFR_RNDN); return mpfr_get_d(r21350, MPFR_RNDN); } static mpfr_t r21351, r21352, r21353, r21354, r21355, r21356, r21357, r21358, r21359, r21360, r21361, r21362, r21363, r21364, r21365, r21366, r21367, r21368, r21369; void setup_mpfr_f_dm() { mpfr_set_default_prec(4112); mpfr_init(r21351); mpfr_init(r21352); mpfr_init(r21353); mpfr_init(r21354); mpfr_init(r21355); mpfr_init(r21356); mpfr_init(r21357); mpfr_init(r21358); mpfr_init(r21359); mpfr_init(r21360); mpfr_init(r21361); mpfr_init(r21362); mpfr_init(r21363); mpfr_init(r21364); mpfr_init(r21365); mpfr_init(r21366); mpfr_init(r21367); mpfr_init(r21368); mpfr_init(r21369); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21351, x_im, MPFR_RNDN); mpfr_set_d(r21352, x_re, MPFR_RNDN); mpfr_hypot(r21353, r21351, r21352, MPFR_RNDN); mpfr_log(r21354, r21353, MPFR_RNDN); mpfr_set_d(r21355, y_im, MPFR_RNDN); mpfr_atan2(r21356, r21351, r21352, MPFR_RNDN); mpfr_set_d(r21357, y_re, MPFR_RNDN); mpfr_mul(r21358, r21356, r21357, MPFR_RNDN); mpfr_fma(r21359, r21354, r21355, r21358, MPFR_RNDN); mpfr_cos(r21360, r21359, MPFR_RNDN); mpfr_expm1(r21361, r21360, MPFR_RNDN); mpfr_exp(r21362, r21361, MPFR_RNDN); mpfr_log(r21363, r21362, MPFR_RNDN); mpfr_log1p(r21364, r21363, MPFR_RNDN); mpfr_mul(r21365, r21355, r21356, MPFR_RNDN); mpfr_mul(r21366, r21354, r21357, MPFR_RNDN); mpfr_sub(r21367, r21365, r21366, MPFR_RNDN); mpfr_exp(r21368, r21367, MPFR_RNDN); mpfr_div(r21369, r21364, r21368, MPFR_RNDN); return mpfr_get_d(r21369, MPFR_RNDN); }