#include #include #include #include #include char *name = "math.log/2 on complex, real part"; double f_if(float re, float im, float base) { float r22038 = re; float r22039 = r22038 * r22038; float r22040 = im; float r22041 = r22040 * r22040; float r22042 = r22039 + r22041; float r22043 = sqrt(r22042); float r22044 = log(r22043); float r22045 = base; float r22046 = log(r22045); float r22047 = r22044 * r22046; float r22048 = atan2(r22040, r22038); float r22049 = 0; float r22050 = r22048 * r22049; float r22051 = r22047 + r22050; float r22052 = r22046 * r22046; float r22053 = r22049 * r22049; float r22054 = r22052 + r22053; float r22055 = r22051 / r22054; return r22055; } double f_id(double re, double im, double base) { double r22056 = re; double r22057 = r22056 * r22056; double r22058 = im; double r22059 = r22058 * r22058; double r22060 = r22057 + r22059; double r22061 = sqrt(r22060); double r22062 = log(r22061); double r22063 = base; double r22064 = log(r22063); double r22065 = r22062 * r22064; double r22066 = atan2(r22058, r22056); double r22067 = 0; double r22068 = r22066 * r22067; double r22069 = r22065 + r22068; double r22070 = r22064 * r22064; double r22071 = r22067 * r22067; double r22072 = r22070 + r22071; double r22073 = r22069 / r22072; return r22073; } double f_of(float re, float im, float base) { float r22074 = re; float r22075 = im; float r22076 = hypot(r22074, r22075); float r22077 = cbrt(r22076); float r22078 = r22077 * r22077; float r22079 = log1p(r22077); float r22080 = expm1(r22079); float r22081 = r22078 * r22080; float r22082 = log(r22081); float r22083 = base; float r22084 = log(r22083); float r22085 = r22082 / r22084; return r22085; } double f_od(double re, double im, double base) { double r22086 = re; double r22087 = im; double r22088 = hypot(r22086, r22087); double r22089 = cbrt(r22088); double r22090 = r22089 * r22089; double r22091 = log1p(r22089); double r22092 = expm1(r22091); double r22093 = r22090 * r22092; double r22094 = log(r22093); double r22095 = base; double r22096 = log(r22095); double r22097 = r22094 / r22096; return r22097; } 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 r22098, r22099, r22100, r22101, r22102, r22103, r22104, r22105, r22106, r22107, r22108, r22109, r22110, r22111, r22112, r22113, r22114, r22115; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22098); mpfr_init(r22099); mpfr_init(r22100); mpfr_init(r22101); mpfr_init(r22102); mpfr_init(r22103); mpfr_init(r22104); mpfr_init(r22105); mpfr_init(r22106); mpfr_init(r22107); mpfr_init(r22108); mpfr_init_set_str(r22109, "0", 10, MPFR_RNDN); mpfr_init(r22110); mpfr_init(r22111); mpfr_init(r22112); mpfr_init(r22113); mpfr_init(r22114); mpfr_init(r22115); } double f_im(double re, double im, double base) { mpfr_set_d(r22098, re, MPFR_RNDN); mpfr_mul(r22099, r22098, r22098, MPFR_RNDN); mpfr_set_d(r22100, im, MPFR_RNDN); mpfr_mul(r22101, r22100, r22100, MPFR_RNDN); mpfr_add(r22102, r22099, r22101, MPFR_RNDN); mpfr_sqrt(r22103, r22102, MPFR_RNDN); mpfr_log(r22104, r22103, MPFR_RNDN); mpfr_set_d(r22105, base, MPFR_RNDN); mpfr_log(r22106, r22105, MPFR_RNDN); mpfr_mul(r22107, r22104, r22106, MPFR_RNDN); mpfr_atan2(r22108, r22100, r22098, MPFR_RNDN); ; mpfr_mul(r22110, r22108, r22109, MPFR_RNDN); mpfr_add(r22111, r22107, r22110, MPFR_RNDN); mpfr_mul(r22112, r22106, r22106, MPFR_RNDN); mpfr_mul(r22113, r22109, r22109, MPFR_RNDN); mpfr_add(r22114, r22112, r22113, MPFR_RNDN); mpfr_div(r22115, r22111, r22114, MPFR_RNDN); return mpfr_get_d(r22115, MPFR_RNDN); } static mpfr_t r22116, r22117, r22118, r22119, r22120, r22121, r22122, r22123, r22124, r22125, r22126, r22127; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22116); mpfr_init(r22117); mpfr_init(r22118); mpfr_init(r22119); mpfr_init(r22120); mpfr_init(r22121); mpfr_init(r22122); mpfr_init(r22123); mpfr_init(r22124); mpfr_init(r22125); mpfr_init(r22126); mpfr_init(r22127); } double f_fm(double re, double im, double base) { mpfr_set_d(r22116, re, MPFR_RNDN); mpfr_set_d(r22117, im, MPFR_RNDN); mpfr_hypot(r22118, r22116, r22117, MPFR_RNDN); mpfr_cbrt(r22119, r22118, MPFR_RNDN); mpfr_mul(r22120, r22119, r22119, MPFR_RNDN); mpfr_log1p(r22121, r22119, MPFR_RNDN); mpfr_expm1(r22122, r22121, MPFR_RNDN); mpfr_mul(r22123, r22120, r22122, MPFR_RNDN); mpfr_log(r22124, r22123, MPFR_RNDN); mpfr_set_d(r22125, base, MPFR_RNDN); mpfr_log(r22126, r22125, MPFR_RNDN); mpfr_div(r22127, r22124, r22126, MPFR_RNDN); return mpfr_get_d(r22127, MPFR_RNDN); } static mpfr_t r22128, r22129, r22130, r22131, r22132, r22133, r22134, r22135, r22136, r22137, r22138, r22139; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22128); mpfr_init(r22129); mpfr_init(r22130); mpfr_init(r22131); mpfr_init(r22132); mpfr_init(r22133); mpfr_init(r22134); mpfr_init(r22135); mpfr_init(r22136); mpfr_init(r22137); mpfr_init(r22138); mpfr_init(r22139); } double f_dm(double re, double im, double base) { mpfr_set_d(r22128, re, MPFR_RNDN); mpfr_set_d(r22129, im, MPFR_RNDN); mpfr_hypot(r22130, r22128, r22129, MPFR_RNDN); mpfr_cbrt(r22131, r22130, MPFR_RNDN); mpfr_mul(r22132, r22131, r22131, MPFR_RNDN); mpfr_log1p(r22133, r22131, MPFR_RNDN); mpfr_expm1(r22134, r22133, MPFR_RNDN); mpfr_mul(r22135, r22132, r22134, MPFR_RNDN); mpfr_log(r22136, r22135, MPFR_RNDN); mpfr_set_d(r22137, base, MPFR_RNDN); mpfr_log(r22138, r22137, MPFR_RNDN); mpfr_div(r22139, r22136, r22138, MPFR_RNDN); return mpfr_get_d(r22139, MPFR_RNDN); }