#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 r20392 = x_re; float r20393 = r20392 * r20392; float r20394 = x_im; float r20395 = r20394 * r20394; float r20396 = r20393 + r20395; float r20397 = sqrt(r20396); float r20398 = log(r20397); float r20399 = y_re; float r20400 = r20398 * r20399; float r20401 = atan2(r20394, r20392); float r20402 = y_im; float r20403 = r20401 * r20402; float r20404 = r20400 - r20403; float r20405 = exp(r20404); float r20406 = r20398 * r20402; float r20407 = r20401 * r20399; float r20408 = r20406 + r20407; float r20409 = cos(r20408); float r20410 = r20405 * r20409; return r20410; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r20411 = x_re; double r20412 = r20411 * r20411; double r20413 = x_im; double r20414 = r20413 * r20413; double r20415 = r20412 + r20414; double r20416 = sqrt(r20415); double r20417 = log(r20416); double r20418 = y_re; double r20419 = r20417 * r20418; double r20420 = atan2(r20413, r20411); double r20421 = y_im; double r20422 = r20420 * r20421; double r20423 = r20419 - r20422; double r20424 = exp(r20423); double r20425 = r20417 * r20421; double r20426 = r20420 * r20418; double r20427 = r20425 + r20426; double r20428 = cos(r20427); double r20429 = r20424 * r20428; return r20429; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r20430 = x_im; float r20431 = x_re; float r20432 = hypot(r20430, r20431); float r20433 = log(r20432); float r20434 = y_im; float r20435 = atan2(r20430, r20431); float r20436 = y_re; float r20437 = r20435 * r20436; float r20438 = fma(r20433, r20434, r20437); float r20439 = cos(r20438); float r20440 = expm1(r20439); float r20441 = log1p(r20440); float r20442 = r20434 * r20435; float r20443 = cbrt(r20442); float r20444 = r20443 * r20443; float r20445 = r20444 * r20443; float r20446 = r20433 * r20436; float r20447 = r20445 - r20446; float r20448 = exp(r20447); float r20449 = r20441 / r20448; return r20449; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r20450 = x_im; double r20451 = x_re; double r20452 = hypot(r20450, r20451); double r20453 = log(r20452); double r20454 = y_im; double r20455 = atan2(r20450, r20451); double r20456 = y_re; double r20457 = r20455 * r20456; double r20458 = fma(r20453, r20454, r20457); double r20459 = cos(r20458); double r20460 = expm1(r20459); double r20461 = log1p(r20460); double r20462 = r20454 * r20455; double r20463 = cbrt(r20462); double r20464 = r20463 * r20463; double r20465 = r20464 * r20463; double r20466 = r20453 * r20456; double r20467 = r20465 - r20466; double r20468 = exp(r20467); double r20469 = r20461 / r20468; return r20469; } 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 r20470, r20471, r20472, r20473, r20474, r20475, r20476, r20477, r20478, r20479, r20480, r20481, r20482, r20483, r20484, r20485, r20486, r20487, r20488; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r20470); mpfr_init(r20471); mpfr_init(r20472); mpfr_init(r20473); mpfr_init(r20474); mpfr_init(r20475); mpfr_init(r20476); mpfr_init(r20477); mpfr_init(r20478); mpfr_init(r20479); mpfr_init(r20480); mpfr_init(r20481); mpfr_init(r20482); mpfr_init(r20483); mpfr_init(r20484); mpfr_init(r20485); mpfr_init(r20486); mpfr_init(r20487); mpfr_init(r20488); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r20470, x_re, MPFR_RNDN); mpfr_mul(r20471, r20470, r20470, MPFR_RNDN); mpfr_set_d(r20472, x_im, MPFR_RNDN); mpfr_mul(r20473, r20472, r20472, MPFR_RNDN); mpfr_add(r20474, r20471, r20473, MPFR_RNDN); mpfr_sqrt(r20475, r20474, MPFR_RNDN); mpfr_log(r20476, r20475, MPFR_RNDN); mpfr_set_d(r20477, y_re, MPFR_RNDN); mpfr_mul(r20478, r20476, r20477, MPFR_RNDN); mpfr_atan2(r20479, r20472, r20470, MPFR_RNDN); mpfr_set_d(r20480, y_im, MPFR_RNDN); mpfr_mul(r20481, r20479, r20480, MPFR_RNDN); mpfr_sub(r20482, r20478, r20481, MPFR_RNDN); mpfr_exp(r20483, r20482, MPFR_RNDN); mpfr_mul(r20484, r20476, r20480, MPFR_RNDN); mpfr_mul(r20485, r20479, r20477, MPFR_RNDN); mpfr_add(r20486, r20484, r20485, MPFR_RNDN); mpfr_cos(r20487, r20486, MPFR_RNDN); mpfr_mul(r20488, r20483, r20487, MPFR_RNDN); return mpfr_get_d(r20488, MPFR_RNDN); } static mpfr_t r20489, r20490, r20491, r20492, r20493, r20494, r20495, r20496, r20497, r20498, r20499, r20500, r20501, r20502, r20503, r20504, r20505, r20506, r20507, r20508; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r20489); mpfr_init(r20490); mpfr_init(r20491); mpfr_init(r20492); mpfr_init(r20493); mpfr_init(r20494); mpfr_init(r20495); mpfr_init(r20496); mpfr_init(r20497); mpfr_init(r20498); mpfr_init(r20499); mpfr_init(r20500); mpfr_init(r20501); mpfr_init(r20502); mpfr_init(r20503); mpfr_init(r20504); mpfr_init(r20505); mpfr_init(r20506); mpfr_init(r20507); mpfr_init(r20508); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r20489, x_im, MPFR_RNDN); mpfr_set_d(r20490, x_re, MPFR_RNDN); mpfr_hypot(r20491, r20489, r20490, MPFR_RNDN); mpfr_log(r20492, r20491, MPFR_RNDN); mpfr_set_d(r20493, y_im, MPFR_RNDN); mpfr_atan2(r20494, r20489, r20490, MPFR_RNDN); mpfr_set_d(r20495, y_re, MPFR_RNDN); mpfr_mul(r20496, r20494, r20495, MPFR_RNDN); mpfr_fma(r20497, r20492, r20493, r20496, MPFR_RNDN); mpfr_cos(r20498, r20497, MPFR_RNDN); mpfr_expm1(r20499, r20498, MPFR_RNDN); mpfr_log1p(r20500, r20499, MPFR_RNDN); mpfr_mul(r20501, r20493, r20494, MPFR_RNDN); mpfr_cbrt(r20502, r20501, MPFR_RNDN); mpfr_mul(r20503, r20502, r20502, MPFR_RNDN); mpfr_mul(r20504, r20503, r20502, MPFR_RNDN); mpfr_mul(r20505, r20492, r20495, MPFR_RNDN); mpfr_sub(r20506, r20504, r20505, MPFR_RNDN); mpfr_exp(r20507, r20506, MPFR_RNDN); mpfr_div(r20508, r20500, r20507, MPFR_RNDN); return mpfr_get_d(r20508, MPFR_RNDN); } static mpfr_t r20509, r20510, r20511, r20512, r20513, r20514, r20515, r20516, r20517, r20518, r20519, r20520, r20521, r20522, r20523, r20524, r20525, r20526, r20527, r20528; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r20509); mpfr_init(r20510); mpfr_init(r20511); mpfr_init(r20512); mpfr_init(r20513); mpfr_init(r20514); mpfr_init(r20515); mpfr_init(r20516); mpfr_init(r20517); mpfr_init(r20518); mpfr_init(r20519); mpfr_init(r20520); mpfr_init(r20521); mpfr_init(r20522); mpfr_init(r20523); mpfr_init(r20524); mpfr_init(r20525); mpfr_init(r20526); mpfr_init(r20527); mpfr_init(r20528); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r20509, x_im, MPFR_RNDN); mpfr_set_d(r20510, x_re, MPFR_RNDN); mpfr_hypot(r20511, r20509, r20510, MPFR_RNDN); mpfr_log(r20512, r20511, MPFR_RNDN); mpfr_set_d(r20513, y_im, MPFR_RNDN); mpfr_atan2(r20514, r20509, r20510, MPFR_RNDN); mpfr_set_d(r20515, y_re, MPFR_RNDN); mpfr_mul(r20516, r20514, r20515, MPFR_RNDN); mpfr_fma(r20517, r20512, r20513, r20516, MPFR_RNDN); mpfr_cos(r20518, r20517, MPFR_RNDN); mpfr_expm1(r20519, r20518, MPFR_RNDN); mpfr_log1p(r20520, r20519, MPFR_RNDN); mpfr_mul(r20521, r20513, r20514, MPFR_RNDN); mpfr_cbrt(r20522, r20521, MPFR_RNDN); mpfr_mul(r20523, r20522, r20522, MPFR_RNDN); mpfr_mul(r20524, r20523, r20522, MPFR_RNDN); mpfr_mul(r20525, r20512, r20515, MPFR_RNDN); mpfr_sub(r20526, r20524, r20525, MPFR_RNDN); mpfr_exp(r20527, r20526, MPFR_RNDN); mpfr_div(r20528, r20520, r20527, MPFR_RNDN); return mpfr_get_d(r20528, MPFR_RNDN); }