#include #include #include #include #include char *name = "powComplex, imaginary part"; double f_if(float x_re, float x_im, float y_re, float y_im) { float r21370 = x_re; float r21371 = r21370 * r21370; float r21372 = x_im; float r21373 = r21372 * r21372; float r21374 = r21371 + r21373; float r21375 = sqrt(r21374); float r21376 = log(r21375); float r21377 = y_re; float r21378 = r21376 * r21377; float r21379 = atan2(r21372, r21370); float r21380 = y_im; float r21381 = r21379 * r21380; float r21382 = r21378 - r21381; float r21383 = exp(r21382); float r21384 = r21376 * r21380; float r21385 = r21379 * r21377; float r21386 = r21384 + r21385; float r21387 = sin(r21386); float r21388 = r21383 * r21387; return r21388; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r21389 = x_re; double r21390 = r21389 * r21389; double r21391 = x_im; double r21392 = r21391 * r21391; double r21393 = r21390 + r21392; double r21394 = sqrt(r21393); double r21395 = log(r21394); double r21396 = y_re; double r21397 = r21395 * r21396; double r21398 = atan2(r21391, r21389); double r21399 = y_im; double r21400 = r21398 * r21399; double r21401 = r21397 - r21400; double r21402 = exp(r21401); double r21403 = r21395 * r21399; double r21404 = r21398 * r21396; double r21405 = r21403 + r21404; double r21406 = sin(r21405); double r21407 = r21402 * r21406; return r21407; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r21408 = x_im; float r21409 = x_re; float r21410 = hypot(r21408, r21409); float r21411 = log(r21410); float r21412 = y_im; float r21413 = atan2(r21408, r21409); float r21414 = y_re; float r21415 = r21413 * r21414; float r21416 = fma(r21411, r21412, r21415); float r21417 = sin(r21416); float r21418 = expm1(r21417); float r21419 = log1p(r21418); float r21420 = r21412 * r21413; float r21421 = r21411 * r21414; float r21422 = r21420 - r21421; float r21423 = exp(r21422); float r21424 = r21419 / r21423; return r21424; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r21425 = x_im; double r21426 = x_re; double r21427 = hypot(r21425, r21426); double r21428 = log(r21427); double r21429 = y_im; double r21430 = atan2(r21425, r21426); double r21431 = y_re; double r21432 = r21430 * r21431; double r21433 = fma(r21428, r21429, r21432); double r21434 = sin(r21433); double r21435 = expm1(r21434); double r21436 = log1p(r21435); double r21437 = r21429 * r21430; double r21438 = r21428 * r21431; double r21439 = r21437 - r21438; double r21440 = exp(r21439); double r21441 = r21436 / r21440; return r21441; } 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 r21442, r21443, r21444, r21445, r21446, r21447, r21448, r21449, r21450, r21451, r21452, r21453, r21454, r21455, r21456, r21457, r21458, r21459, r21460; void setup_mpfr_f_im() { mpfr_set_default_prec(4112); mpfr_init(r21442); mpfr_init(r21443); mpfr_init(r21444); mpfr_init(r21445); mpfr_init(r21446); mpfr_init(r21447); mpfr_init(r21448); mpfr_init(r21449); mpfr_init(r21450); mpfr_init(r21451); mpfr_init(r21452); mpfr_init(r21453); mpfr_init(r21454); mpfr_init(r21455); mpfr_init(r21456); mpfr_init(r21457); mpfr_init(r21458); mpfr_init(r21459); mpfr_init(r21460); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21442, x_re, MPFR_RNDN); mpfr_mul(r21443, r21442, r21442, MPFR_RNDN); mpfr_set_d(r21444, x_im, MPFR_RNDN); mpfr_mul(r21445, r21444, r21444, MPFR_RNDN); mpfr_add(r21446, r21443, r21445, MPFR_RNDN); mpfr_sqrt(r21447, r21446, MPFR_RNDN); mpfr_log(r21448, r21447, MPFR_RNDN); mpfr_set_d(r21449, y_re, MPFR_RNDN); mpfr_mul(r21450, r21448, r21449, MPFR_RNDN); mpfr_atan2(r21451, r21444, r21442, MPFR_RNDN); mpfr_set_d(r21452, y_im, MPFR_RNDN); mpfr_mul(r21453, r21451, r21452, MPFR_RNDN); mpfr_sub(r21454, r21450, r21453, MPFR_RNDN); mpfr_exp(r21455, r21454, MPFR_RNDN); mpfr_mul(r21456, r21448, r21452, MPFR_RNDN); mpfr_mul(r21457, r21451, r21449, MPFR_RNDN); mpfr_add(r21458, r21456, r21457, MPFR_RNDN); mpfr_sin(r21459, r21458, MPFR_RNDN); mpfr_mul(r21460, r21455, r21459, MPFR_RNDN); return mpfr_get_d(r21460, MPFR_RNDN); } static mpfr_t r21461, r21462, r21463, r21464, r21465, r21466, r21467, r21468, r21469, r21470, r21471, r21472, r21473, r21474, r21475, r21476, r21477; void setup_mpfr_f_fm() { mpfr_set_default_prec(4112); mpfr_init(r21461); mpfr_init(r21462); mpfr_init(r21463); mpfr_init(r21464); mpfr_init(r21465); mpfr_init(r21466); mpfr_init(r21467); mpfr_init(r21468); mpfr_init(r21469); mpfr_init(r21470); mpfr_init(r21471); mpfr_init(r21472); mpfr_init(r21473); mpfr_init(r21474); mpfr_init(r21475); mpfr_init(r21476); mpfr_init(r21477); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21461, x_im, MPFR_RNDN); mpfr_set_d(r21462, x_re, MPFR_RNDN); mpfr_hypot(r21463, r21461, r21462, MPFR_RNDN); mpfr_log(r21464, r21463, MPFR_RNDN); mpfr_set_d(r21465, y_im, MPFR_RNDN); mpfr_atan2(r21466, r21461, r21462, MPFR_RNDN); mpfr_set_d(r21467, y_re, MPFR_RNDN); mpfr_mul(r21468, r21466, r21467, MPFR_RNDN); mpfr_fma(r21469, r21464, r21465, r21468, MPFR_RNDN); mpfr_sin(r21470, r21469, MPFR_RNDN); mpfr_expm1(r21471, r21470, MPFR_RNDN); mpfr_log1p(r21472, r21471, MPFR_RNDN); mpfr_mul(r21473, r21465, r21466, MPFR_RNDN); mpfr_mul(r21474, r21464, r21467, MPFR_RNDN); mpfr_sub(r21475, r21473, r21474, MPFR_RNDN); mpfr_exp(r21476, r21475, MPFR_RNDN); mpfr_div(r21477, r21472, r21476, MPFR_RNDN); return mpfr_get_d(r21477, MPFR_RNDN); } static mpfr_t r21478, r21479, r21480, r21481, r21482, r21483, r21484, r21485, r21486, r21487, r21488, r21489, r21490, r21491, r21492, r21493, r21494; void setup_mpfr_f_dm() { mpfr_set_default_prec(4112); mpfr_init(r21478); mpfr_init(r21479); mpfr_init(r21480); mpfr_init(r21481); mpfr_init(r21482); mpfr_init(r21483); mpfr_init(r21484); mpfr_init(r21485); mpfr_init(r21486); mpfr_init(r21487); mpfr_init(r21488); mpfr_init(r21489); mpfr_init(r21490); mpfr_init(r21491); mpfr_init(r21492); mpfr_init(r21493); mpfr_init(r21494); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21478, x_im, MPFR_RNDN); mpfr_set_d(r21479, x_re, MPFR_RNDN); mpfr_hypot(r21480, r21478, r21479, MPFR_RNDN); mpfr_log(r21481, r21480, MPFR_RNDN); mpfr_set_d(r21482, y_im, MPFR_RNDN); mpfr_atan2(r21483, r21478, r21479, MPFR_RNDN); mpfr_set_d(r21484, y_re, MPFR_RNDN); mpfr_mul(r21485, r21483, r21484, MPFR_RNDN); mpfr_fma(r21486, r21481, r21482, r21485, MPFR_RNDN); mpfr_sin(r21487, r21486, MPFR_RNDN); mpfr_expm1(r21488, r21487, MPFR_RNDN); mpfr_log1p(r21489, r21488, MPFR_RNDN); mpfr_mul(r21490, r21482, r21483, MPFR_RNDN); mpfr_mul(r21491, r21481, r21484, MPFR_RNDN); mpfr_sub(r21492, r21490, r21491, MPFR_RNDN); mpfr_exp(r21493, r21492, MPFR_RNDN); mpfr_div(r21494, r21489, r21493, MPFR_RNDN); return mpfr_get_d(r21494, MPFR_RNDN); }