#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 r21062 = x_re; float r21063 = r21062 * r21062; float r21064 = x_im; float r21065 = r21064 * r21064; float r21066 = r21063 + r21065; float r21067 = sqrt(r21066); float r21068 = log(r21067); float r21069 = y_re; float r21070 = r21068 * r21069; float r21071 = atan2(r21064, r21062); float r21072 = y_im; float r21073 = r21071 * r21072; float r21074 = r21070 - r21073; float r21075 = exp(r21074); float r21076 = r21068 * r21072; float r21077 = r21071 * r21069; float r21078 = r21076 + r21077; float r21079 = sin(r21078); float r21080 = r21075 * r21079; return r21080; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r21081 = x_re; double r21082 = r21081 * r21081; double r21083 = x_im; double r21084 = r21083 * r21083; double r21085 = r21082 + r21084; double r21086 = sqrt(r21085); double r21087 = log(r21086); double r21088 = y_re; double r21089 = r21087 * r21088; double r21090 = atan2(r21083, r21081); double r21091 = y_im; double r21092 = r21090 * r21091; double r21093 = r21089 - r21092; double r21094 = exp(r21093); double r21095 = r21087 * r21091; double r21096 = r21090 * r21088; double r21097 = r21095 + r21096; double r21098 = sin(r21097); double r21099 = r21094 * r21098; return r21099; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r21100 = x_im; float r21101 = x_re; float r21102 = hypot(r21100, r21101); float r21103 = log(r21102); float r21104 = y_re; float r21105 = r21103 * r21104; float r21106 = y_im; float r21107 = atan2(r21100, r21101); float r21108 = r21106 * r21107; float r21109 = r21105 - r21108; float r21110 = exp(r21109); float r21111 = r21107 * r21104; float r21112 = fma(r21103, r21106, r21111); float r21113 = sin(r21112); float r21114 = expm1(r21113); float r21115 = log1p(r21114); float r21116 = r21110 * r21115; return r21116; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r21117 = x_im; double r21118 = x_re; double r21119 = hypot(r21117, r21118); double r21120 = log(r21119); double r21121 = y_re; double r21122 = r21120 * r21121; double r21123 = y_im; double r21124 = atan2(r21117, r21118); double r21125 = r21123 * r21124; double r21126 = r21122 - r21125; double r21127 = exp(r21126); double r21128 = r21124 * r21121; double r21129 = fma(r21120, r21123, r21128); double r21130 = sin(r21129); double r21131 = expm1(r21130); double r21132 = log1p(r21131); double r21133 = r21127 * r21132; return r21133; } 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 r21134, r21135, r21136, r21137, r21138, r21139, r21140, r21141, r21142, r21143, r21144, r21145, r21146, r21147, r21148, r21149, r21150, r21151, r21152; void setup_mpfr_f_im() { mpfr_set_default_prec(2192); mpfr_init(r21134); mpfr_init(r21135); mpfr_init(r21136); mpfr_init(r21137); mpfr_init(r21138); mpfr_init(r21139); mpfr_init(r21140); mpfr_init(r21141); mpfr_init(r21142); mpfr_init(r21143); mpfr_init(r21144); mpfr_init(r21145); mpfr_init(r21146); mpfr_init(r21147); mpfr_init(r21148); mpfr_init(r21149); mpfr_init(r21150); mpfr_init(r21151); mpfr_init(r21152); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21134, x_re, MPFR_RNDN); mpfr_mul(r21135, r21134, r21134, MPFR_RNDN); mpfr_set_d(r21136, x_im, MPFR_RNDN); mpfr_mul(r21137, r21136, r21136, MPFR_RNDN); mpfr_add(r21138, r21135, r21137, MPFR_RNDN); mpfr_sqrt(r21139, r21138, MPFR_RNDN); mpfr_log(r21140, r21139, MPFR_RNDN); mpfr_set_d(r21141, y_re, MPFR_RNDN); mpfr_mul(r21142, r21140, r21141, MPFR_RNDN); mpfr_atan2(r21143, r21136, r21134, MPFR_RNDN); mpfr_set_d(r21144, y_im, MPFR_RNDN); mpfr_mul(r21145, r21143, r21144, MPFR_RNDN); mpfr_sub(r21146, r21142, r21145, MPFR_RNDN); mpfr_exp(r21147, r21146, MPFR_RNDN); mpfr_mul(r21148, r21140, r21144, MPFR_RNDN); mpfr_mul(r21149, r21143, r21141, MPFR_RNDN); mpfr_add(r21150, r21148, r21149, MPFR_RNDN); mpfr_sin(r21151, r21150, MPFR_RNDN); mpfr_mul(r21152, r21147, r21151, MPFR_RNDN); return mpfr_get_d(r21152, MPFR_RNDN); } static mpfr_t r21153, r21154, r21155, r21156, r21157, r21158, r21159, r21160, r21161, r21162, r21163, r21164, r21165, r21166, r21167, r21168, r21169; void setup_mpfr_f_fm() { mpfr_set_default_prec(2192); mpfr_init(r21153); mpfr_init(r21154); mpfr_init(r21155); mpfr_init(r21156); mpfr_init(r21157); mpfr_init(r21158); mpfr_init(r21159); mpfr_init(r21160); mpfr_init(r21161); mpfr_init(r21162); mpfr_init(r21163); mpfr_init(r21164); mpfr_init(r21165); mpfr_init(r21166); mpfr_init(r21167); mpfr_init(r21168); mpfr_init(r21169); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21153, x_im, MPFR_RNDN); mpfr_set_d(r21154, x_re, MPFR_RNDN); mpfr_hypot(r21155, r21153, r21154, MPFR_RNDN); mpfr_log(r21156, r21155, MPFR_RNDN); mpfr_set_d(r21157, y_re, MPFR_RNDN); mpfr_mul(r21158, r21156, r21157, MPFR_RNDN); mpfr_set_d(r21159, y_im, MPFR_RNDN); mpfr_atan2(r21160, r21153, r21154, MPFR_RNDN); mpfr_mul(r21161, r21159, r21160, MPFR_RNDN); mpfr_sub(r21162, r21158, r21161, MPFR_RNDN); mpfr_exp(r21163, r21162, MPFR_RNDN); mpfr_mul(r21164, r21160, r21157, MPFR_RNDN); mpfr_fma(r21165, r21156, r21159, r21164, MPFR_RNDN); mpfr_sin(r21166, r21165, MPFR_RNDN); mpfr_expm1(r21167, r21166, MPFR_RNDN); mpfr_log1p(r21168, r21167, MPFR_RNDN); mpfr_mul(r21169, r21163, r21168, MPFR_RNDN); return mpfr_get_d(r21169, MPFR_RNDN); } static mpfr_t r21170, r21171, r21172, r21173, r21174, r21175, r21176, r21177, r21178, r21179, r21180, r21181, r21182, r21183, r21184, r21185, r21186; void setup_mpfr_f_dm() { mpfr_set_default_prec(2192); mpfr_init(r21170); mpfr_init(r21171); mpfr_init(r21172); mpfr_init(r21173); mpfr_init(r21174); mpfr_init(r21175); mpfr_init(r21176); mpfr_init(r21177); mpfr_init(r21178); mpfr_init(r21179); mpfr_init(r21180); mpfr_init(r21181); mpfr_init(r21182); mpfr_init(r21183); mpfr_init(r21184); mpfr_init(r21185); mpfr_init(r21186); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21170, x_im, MPFR_RNDN); mpfr_set_d(r21171, x_re, MPFR_RNDN); mpfr_hypot(r21172, r21170, r21171, MPFR_RNDN); mpfr_log(r21173, r21172, MPFR_RNDN); mpfr_set_d(r21174, y_re, MPFR_RNDN); mpfr_mul(r21175, r21173, r21174, MPFR_RNDN); mpfr_set_d(r21176, y_im, MPFR_RNDN); mpfr_atan2(r21177, r21170, r21171, MPFR_RNDN); mpfr_mul(r21178, r21176, r21177, MPFR_RNDN); mpfr_sub(r21179, r21175, r21178, MPFR_RNDN); mpfr_exp(r21180, r21179, MPFR_RNDN); mpfr_mul(r21181, r21177, r21174, MPFR_RNDN); mpfr_fma(r21182, r21173, r21176, r21181, MPFR_RNDN); mpfr_sin(r21183, r21182, MPFR_RNDN); mpfr_expm1(r21184, r21183, MPFR_RNDN); mpfr_log1p(r21185, r21184, MPFR_RNDN); mpfr_mul(r21186, r21180, r21185, MPFR_RNDN); return mpfr_get_d(r21186, MPFR_RNDN); }