#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 r20925 = x_re; float r20926 = r20925 * r20925; float r20927 = x_im; float r20928 = r20927 * r20927; float r20929 = r20926 + r20928; float r20930 = sqrt(r20929); float r20931 = log(r20930); float r20932 = y_re; float r20933 = r20931 * r20932; float r20934 = atan2(r20927, r20925); float r20935 = y_im; float r20936 = r20934 * r20935; float r20937 = r20933 - r20936; float r20938 = exp(r20937); float r20939 = r20931 * r20935; float r20940 = r20934 * r20932; float r20941 = r20939 + r20940; float r20942 = cos(r20941); float r20943 = r20938 * r20942; return r20943; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r20944 = x_re; double r20945 = r20944 * r20944; double r20946 = x_im; double r20947 = r20946 * r20946; double r20948 = r20945 + r20947; double r20949 = sqrt(r20948); double r20950 = log(r20949); double r20951 = y_re; double r20952 = r20950 * r20951; double r20953 = atan2(r20946, r20944); double r20954 = y_im; double r20955 = r20953 * r20954; double r20956 = r20952 - r20955; double r20957 = exp(r20956); double r20958 = r20950 * r20954; double r20959 = r20953 * r20951; double r20960 = r20958 + r20959; double r20961 = cos(r20960); double r20962 = r20957 * r20961; return r20962; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r20963 = x_im; float r20964 = x_re; float r20965 = hypot(r20963, r20964); float r20966 = log(r20965); float r20967 = y_re; float r20968 = r20966 * r20967; float r20969 = y_im; float r20970 = atan2(r20963, r20964); float r20971 = r20969 * r20970; float r20972 = cbrt(r20971); float r20973 = r20972 * r20972; float r20974 = r20973 * r20972; float r20975 = r20968 - r20974; float r20976 = exp(r20975); float r20977 = r20970 * r20967; float r20978 = fma(r20966, r20969, r20977); float r20979 = cos(r20978); float r20980 = expm1(r20979); float r20981 = log1p(r20980); float r20982 = r20976 * r20981; return r20982; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r20983 = x_im; double r20984 = x_re; double r20985 = hypot(r20983, r20984); double r20986 = log(r20985); double r20987 = y_re; double r20988 = r20986 * r20987; double r20989 = y_im; double r20990 = atan2(r20983, r20984); double r20991 = r20989 * r20990; double r20992 = cbrt(r20991); double r20993 = r20992 * r20992; double r20994 = r20993 * r20992; double r20995 = r20988 - r20994; double r20996 = exp(r20995); double r20997 = r20990 * r20987; double r20998 = fma(r20986, r20989, r20997); double r20999 = cos(r20998); double r21000 = expm1(r20999); double r21001 = log1p(r21000); double r21002 = r20996 * r21001; return r21002; } 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 r21003, r21004, r21005, r21006, r21007, r21008, r21009, r21010, r21011, r21012, r21013, r21014, r21015, r21016, r21017, r21018, r21019, r21020, r21021; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r21003); mpfr_init(r21004); mpfr_init(r21005); mpfr_init(r21006); mpfr_init(r21007); mpfr_init(r21008); mpfr_init(r21009); mpfr_init(r21010); mpfr_init(r21011); mpfr_init(r21012); mpfr_init(r21013); mpfr_init(r21014); mpfr_init(r21015); mpfr_init(r21016); mpfr_init(r21017); mpfr_init(r21018); mpfr_init(r21019); mpfr_init(r21020); mpfr_init(r21021); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21003, x_re, MPFR_RNDN); mpfr_mul(r21004, r21003, r21003, MPFR_RNDN); mpfr_set_d(r21005, x_im, MPFR_RNDN); mpfr_mul(r21006, r21005, r21005, MPFR_RNDN); mpfr_add(r21007, r21004, r21006, MPFR_RNDN); mpfr_sqrt(r21008, r21007, MPFR_RNDN); mpfr_log(r21009, r21008, MPFR_RNDN); mpfr_set_d(r21010, y_re, MPFR_RNDN); mpfr_mul(r21011, r21009, r21010, MPFR_RNDN); mpfr_atan2(r21012, r21005, r21003, MPFR_RNDN); mpfr_set_d(r21013, y_im, MPFR_RNDN); mpfr_mul(r21014, r21012, r21013, MPFR_RNDN); mpfr_sub(r21015, r21011, r21014, MPFR_RNDN); mpfr_exp(r21016, r21015, MPFR_RNDN); mpfr_mul(r21017, r21009, r21013, MPFR_RNDN); mpfr_mul(r21018, r21012, r21010, MPFR_RNDN); mpfr_add(r21019, r21017, r21018, MPFR_RNDN); mpfr_cos(r21020, r21019, MPFR_RNDN); mpfr_mul(r21021, r21016, r21020, MPFR_RNDN); return mpfr_get_d(r21021, MPFR_RNDN); } static mpfr_t r21022, r21023, r21024, r21025, r21026, r21027, r21028, r21029, r21030, r21031, r21032, r21033, r21034, r21035, r21036, r21037, r21038, r21039, r21040, r21041; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r21022); mpfr_init(r21023); mpfr_init(r21024); mpfr_init(r21025); mpfr_init(r21026); mpfr_init(r21027); mpfr_init(r21028); mpfr_init(r21029); mpfr_init(r21030); mpfr_init(r21031); mpfr_init(r21032); mpfr_init(r21033); mpfr_init(r21034); mpfr_init(r21035); mpfr_init(r21036); mpfr_init(r21037); mpfr_init(r21038); mpfr_init(r21039); mpfr_init(r21040); mpfr_init(r21041); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21022, x_im, MPFR_RNDN); mpfr_set_d(r21023, x_re, MPFR_RNDN); mpfr_hypot(r21024, r21022, r21023, MPFR_RNDN); mpfr_log(r21025, r21024, MPFR_RNDN); mpfr_set_d(r21026, y_re, MPFR_RNDN); mpfr_mul(r21027, r21025, r21026, MPFR_RNDN); mpfr_set_d(r21028, y_im, MPFR_RNDN); mpfr_atan2(r21029, r21022, r21023, MPFR_RNDN); mpfr_mul(r21030, r21028, r21029, MPFR_RNDN); mpfr_cbrt(r21031, r21030, MPFR_RNDN); mpfr_mul(r21032, r21031, r21031, MPFR_RNDN); mpfr_mul(r21033, r21032, r21031, MPFR_RNDN); mpfr_sub(r21034, r21027, r21033, MPFR_RNDN); mpfr_exp(r21035, r21034, MPFR_RNDN); mpfr_mul(r21036, r21029, r21026, MPFR_RNDN); mpfr_fma(r21037, r21025, r21028, r21036, MPFR_RNDN); mpfr_cos(r21038, r21037, MPFR_RNDN); mpfr_expm1(r21039, r21038, MPFR_RNDN); mpfr_log1p(r21040, r21039, MPFR_RNDN); mpfr_mul(r21041, r21035, r21040, MPFR_RNDN); return mpfr_get_d(r21041, MPFR_RNDN); } static mpfr_t r21042, r21043, r21044, r21045, r21046, r21047, r21048, r21049, r21050, r21051, r21052, r21053, r21054, r21055, r21056, r21057, r21058, r21059, r21060, r21061; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r21042); mpfr_init(r21043); mpfr_init(r21044); mpfr_init(r21045); mpfr_init(r21046); mpfr_init(r21047); mpfr_init(r21048); mpfr_init(r21049); mpfr_init(r21050); mpfr_init(r21051); mpfr_init(r21052); mpfr_init(r21053); mpfr_init(r21054); mpfr_init(r21055); mpfr_init(r21056); mpfr_init(r21057); mpfr_init(r21058); mpfr_init(r21059); mpfr_init(r21060); mpfr_init(r21061); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r21042, x_im, MPFR_RNDN); mpfr_set_d(r21043, x_re, MPFR_RNDN); mpfr_hypot(r21044, r21042, r21043, MPFR_RNDN); mpfr_log(r21045, r21044, MPFR_RNDN); mpfr_set_d(r21046, y_re, MPFR_RNDN); mpfr_mul(r21047, r21045, r21046, MPFR_RNDN); mpfr_set_d(r21048, y_im, MPFR_RNDN); mpfr_atan2(r21049, r21042, r21043, MPFR_RNDN); mpfr_mul(r21050, r21048, r21049, MPFR_RNDN); mpfr_cbrt(r21051, r21050, MPFR_RNDN); mpfr_mul(r21052, r21051, r21051, MPFR_RNDN); mpfr_mul(r21053, r21052, r21051, MPFR_RNDN); mpfr_sub(r21054, r21047, r21053, MPFR_RNDN); mpfr_exp(r21055, r21054, MPFR_RNDN); mpfr_mul(r21056, r21049, r21046, MPFR_RNDN); mpfr_fma(r21057, r21045, r21048, r21056, MPFR_RNDN); mpfr_cos(r21058, r21057, MPFR_RNDN); mpfr_expm1(r21059, r21058, MPFR_RNDN); mpfr_log1p(r21060, r21059, MPFR_RNDN); mpfr_mul(r21061, r21055, r21060, MPFR_RNDN); return mpfr_get_d(r21061, MPFR_RNDN); }