#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 r21411 = x_re; float r21412 = r21411 * r21411; float r21413 = x_im; float r21414 = r21413 * r21413; float r21415 = r21412 + r21414; float r21416 = sqrt(r21415); float r21417 = log(r21416); float r21418 = y_re; float r21419 = r21417 * r21418; float r21420 = atan2(r21413, r21411); float r21421 = y_im; float r21422 = r21420 * r21421; float r21423 = r21419 - r21422; float r21424 = exp(r21423); float r21425 = r21417 * r21421; float r21426 = r21420 * r21418; float r21427 = r21425 + r21426; float r21428 = sin(r21427); float r21429 = r21424 * r21428; return r21429; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r21430 = x_re; double r21431 = r21430 * r21430; double r21432 = x_im; double r21433 = r21432 * r21432; double r21434 = r21431 + r21433; double r21435 = sqrt(r21434); double r21436 = log(r21435); double r21437 = y_re; double r21438 = r21436 * r21437; double r21439 = atan2(r21432, r21430); double r21440 = y_im; double r21441 = r21439 * r21440; double r21442 = r21438 - r21441; double r21443 = exp(r21442); double r21444 = r21436 * r21440; double r21445 = r21439 * r21437; double r21446 = r21444 + r21445; double r21447 = sin(r21446); double r21448 = r21443 * r21447; return r21448; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r21449 = x_im; float r21450 = x_re; float r21451 = hypot(r21449, r21450); float r21452 = log(r21451); float r21453 = y_im; float r21454 = atan2(r21449, r21450); float r21455 = y_re; float r21456 = r21454 * r21455; float r21457 = fma(r21452, r21453, r21456); float r21458 = sin(r21457); float r21459 = cbrt(r21458); float r21460 = r21459 * r21459; float r21461 = r21460 * r21459; float r21462 = r21453 * r21454; float r21463 = r21452 * r21455; float r21464 = r21462 - r21463; float r21465 = exp(r21464); float r21466 = r21461 / r21465; return r21466; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r21467 = x_im; double r21468 = x_re; double r21469 = hypot(r21467, r21468); double r21470 = log(r21469); double r21471 = y_im; double r21472 = atan2(r21467, r21468); double r21473 = y_re; double r21474 = r21472 * r21473; double r21475 = fma(r21470, r21471, r21474); double r21476 = sin(r21475); double r21477 = cbrt(r21476); double r21478 = r21477 * r21477; double r21479 = r21478 * r21477; double r21480 = r21471 * r21472; double r21481 = r21470 * r21473; double r21482 = r21480 - r21481; double r21483 = exp(r21482); double r21484 = r21479 / r21483; return r21484; } 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 r21485, r21486, r21487, r21488, r21489, r21490, r21491, r21492, r21493, r21494, r21495, r21496, r21497, r21498, r21499, r21500, r21501, r21502, r21503; void setup_mpfr_f_im() { mpfr_set_default_prec(2128); 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); mpfr_init(r21495); mpfr_init(r21496); mpfr_init(r21497); mpfr_init(r21498); mpfr_init(r21499); mpfr_init(r21500); mpfr_init(r21501); mpfr_init(r21502); mpfr_init(r21503); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21485, x_re, MPFR_RNDN); mpfr_mul(r21486, r21485, r21485, MPFR_RNDN); mpfr_set_d(r21487, x_im, MPFR_RNDN); mpfr_mul(r21488, r21487, r21487, MPFR_RNDN); mpfr_add(r21489, r21486, r21488, MPFR_RNDN); mpfr_sqrt(r21490, r21489, MPFR_RNDN); mpfr_log(r21491, r21490, MPFR_RNDN); mpfr_set_d(r21492, y_re, MPFR_RNDN); mpfr_mul(r21493, r21491, r21492, MPFR_RNDN); mpfr_atan2(r21494, r21487, r21485, MPFR_RNDN); mpfr_set_d(r21495, y_im, MPFR_RNDN); mpfr_mul(r21496, r21494, r21495, MPFR_RNDN); mpfr_sub(r21497, r21493, r21496, MPFR_RNDN); mpfr_exp(r21498, r21497, MPFR_RNDN); mpfr_mul(r21499, r21491, r21495, MPFR_RNDN); mpfr_mul(r21500, r21494, r21492, MPFR_RNDN); mpfr_add(r21501, r21499, r21500, MPFR_RNDN); mpfr_sin(r21502, r21501, MPFR_RNDN); mpfr_mul(r21503, r21498, r21502, MPFR_RNDN); return mpfr_get_d(r21503, MPFR_RNDN); } static mpfr_t r21504, r21505, r21506, r21507, r21508, r21509, r21510, r21511, r21512, r21513, r21514, r21515, r21516, r21517, r21518, r21519, r21520, r21521; void setup_mpfr_f_fm() { mpfr_set_default_prec(2128); mpfr_init(r21504); mpfr_init(r21505); mpfr_init(r21506); mpfr_init(r21507); mpfr_init(r21508); mpfr_init(r21509); mpfr_init(r21510); mpfr_init(r21511); mpfr_init(r21512); mpfr_init(r21513); mpfr_init(r21514); mpfr_init(r21515); mpfr_init(r21516); mpfr_init(r21517); mpfr_init(r21518); mpfr_init(r21519); mpfr_init(r21520); mpfr_init(r21521); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21504, x_im, MPFR_RNDN); mpfr_set_d(r21505, x_re, MPFR_RNDN); mpfr_hypot(r21506, r21504, r21505, MPFR_RNDN); mpfr_log(r21507, r21506, MPFR_RNDN); mpfr_set_d(r21508, y_im, MPFR_RNDN); mpfr_atan2(r21509, r21504, r21505, MPFR_RNDN); mpfr_set_d(r21510, y_re, MPFR_RNDN); mpfr_mul(r21511, r21509, r21510, MPFR_RNDN); mpfr_fma(r21512, r21507, r21508, r21511, MPFR_RNDN); mpfr_sin(r21513, r21512, MPFR_RNDN); mpfr_cbrt(r21514, r21513, MPFR_RNDN); mpfr_mul(r21515, r21514, r21514, MPFR_RNDN); mpfr_mul(r21516, r21515, r21514, MPFR_RNDN); mpfr_mul(r21517, r21508, r21509, MPFR_RNDN); mpfr_mul(r21518, r21507, r21510, MPFR_RNDN); mpfr_sub(r21519, r21517, r21518, MPFR_RNDN); mpfr_exp(r21520, r21519, MPFR_RNDN); mpfr_div(r21521, r21516, r21520, MPFR_RNDN); return mpfr_get_d(r21521, MPFR_RNDN); } static mpfr_t r21522, r21523, r21524, r21525, r21526, r21527, r21528, r21529, r21530, r21531, r21532, r21533, r21534, r21535, r21536, r21537, r21538, r21539; void setup_mpfr_f_dm() { mpfr_set_default_prec(2128); mpfr_init(r21522); mpfr_init(r21523); mpfr_init(r21524); mpfr_init(r21525); mpfr_init(r21526); mpfr_init(r21527); mpfr_init(r21528); mpfr_init(r21529); mpfr_init(r21530); mpfr_init(r21531); mpfr_init(r21532); mpfr_init(r21533); mpfr_init(r21534); mpfr_init(r21535); mpfr_init(r21536); mpfr_init(r21537); mpfr_init(r21538); mpfr_init(r21539); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21522, x_im, MPFR_RNDN); mpfr_set_d(r21523, x_re, MPFR_RNDN); mpfr_hypot(r21524, r21522, r21523, MPFR_RNDN); mpfr_log(r21525, r21524, MPFR_RNDN); mpfr_set_d(r21526, y_im, MPFR_RNDN); mpfr_atan2(r21527, r21522, r21523, MPFR_RNDN); mpfr_set_d(r21528, y_re, MPFR_RNDN); mpfr_mul(r21529, r21527, r21528, MPFR_RNDN); mpfr_fma(r21530, r21525, r21526, r21529, MPFR_RNDN); mpfr_sin(r21531, r21530, MPFR_RNDN); mpfr_cbrt(r21532, r21531, MPFR_RNDN); mpfr_mul(r21533, r21532, r21532, MPFR_RNDN); mpfr_mul(r21534, r21533, r21532, MPFR_RNDN); mpfr_mul(r21535, r21526, r21527, MPFR_RNDN); mpfr_mul(r21536, r21525, r21528, MPFR_RNDN); mpfr_sub(r21537, r21535, r21536, MPFR_RNDN); mpfr_exp(r21538, r21537, MPFR_RNDN); mpfr_div(r21539, r21534, r21538, MPFR_RNDN); return mpfr_get_d(r21539, MPFR_RNDN); }