#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 r23894 = x_re; float r23895 = r23894 * r23894; float r23896 = x_im; float r23897 = r23896 * r23896; float r23898 = r23895 + r23897; float r23899 = sqrt(r23898); float r23900 = log(r23899); float r23901 = y_re; float r23902 = r23900 * r23901; float r23903 = atan2(r23896, r23894); float r23904 = y_im; float r23905 = r23903 * r23904; float r23906 = r23902 - r23905; float r23907 = exp(r23906); float r23908 = r23900 * r23904; float r23909 = r23903 * r23901; float r23910 = r23908 + r23909; float r23911 = cos(r23910); float r23912 = r23907 * r23911; return r23912; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r23913 = x_re; double r23914 = r23913 * r23913; double r23915 = x_im; double r23916 = r23915 * r23915; double r23917 = r23914 + r23916; double r23918 = sqrt(r23917); double r23919 = log(r23918); double r23920 = y_re; double r23921 = r23919 * r23920; double r23922 = atan2(r23915, r23913); double r23923 = y_im; double r23924 = r23922 * r23923; double r23925 = r23921 - r23924; double r23926 = exp(r23925); double r23927 = r23919 * r23923; double r23928 = r23922 * r23920; double r23929 = r23927 + r23928; double r23930 = cos(r23929); double r23931 = r23926 * r23930; return r23931; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r23932 = x_im; float r23933 = x_re; float r23934 = hypot(r23932, r23933); float r23935 = log(r23934); float r23936 = y_im; float r23937 = atan2(r23932, r23933); float r23938 = y_re; float r23939 = r23937 * r23938; float r23940 = fma(r23935, r23936, r23939); float r23941 = cos(r23940); float r23942 = expm1(r23941); float r23943 = exp(r23942); float r23944 = log(r23943); float r23945 = log1p(r23944); float r23946 = r23936 * r23937; float r23947 = r23935 * r23938; float r23948 = r23946 - r23947; float r23949 = exp(r23948); float r23950 = r23945 / r23949; return r23950; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r23951 = x_im; double r23952 = x_re; double r23953 = hypot(r23951, r23952); double r23954 = log(r23953); double r23955 = y_im; double r23956 = atan2(r23951, r23952); double r23957 = y_re; double r23958 = r23956 * r23957; double r23959 = fma(r23954, r23955, r23958); double r23960 = cos(r23959); double r23961 = expm1(r23960); double r23962 = exp(r23961); double r23963 = log(r23962); double r23964 = log1p(r23963); double r23965 = r23955 * r23956; double r23966 = r23954 * r23957; double r23967 = r23965 - r23966; double r23968 = exp(r23967); double r23969 = r23964 / r23968; return r23969; } 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 r23970, r23971, r23972, r23973, r23974, r23975, r23976, r23977, r23978, r23979, r23980, r23981, r23982, r23983, r23984, r23985, r23986, r23987, r23988; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r23970); mpfr_init(r23971); mpfr_init(r23972); mpfr_init(r23973); mpfr_init(r23974); mpfr_init(r23975); mpfr_init(r23976); mpfr_init(r23977); mpfr_init(r23978); mpfr_init(r23979); mpfr_init(r23980); mpfr_init(r23981); mpfr_init(r23982); mpfr_init(r23983); mpfr_init(r23984); mpfr_init(r23985); mpfr_init(r23986); mpfr_init(r23987); mpfr_init(r23988); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r23970, x_re, MPFR_RNDN); mpfr_mul(r23971, r23970, r23970, MPFR_RNDN); mpfr_set_d(r23972, x_im, MPFR_RNDN); mpfr_mul(r23973, r23972, r23972, MPFR_RNDN); mpfr_add(r23974, r23971, r23973, MPFR_RNDN); mpfr_sqrt(r23975, r23974, MPFR_RNDN); mpfr_log(r23976, r23975, MPFR_RNDN); mpfr_set_d(r23977, y_re, MPFR_RNDN); mpfr_mul(r23978, r23976, r23977, MPFR_RNDN); mpfr_atan2(r23979, r23972, r23970, MPFR_RNDN); mpfr_set_d(r23980, y_im, MPFR_RNDN); mpfr_mul(r23981, r23979, r23980, MPFR_RNDN); mpfr_sub(r23982, r23978, r23981, MPFR_RNDN); mpfr_exp(r23983, r23982, MPFR_RNDN); mpfr_mul(r23984, r23976, r23980, MPFR_RNDN); mpfr_mul(r23985, r23979, r23977, MPFR_RNDN); mpfr_add(r23986, r23984, r23985, MPFR_RNDN); mpfr_cos(r23987, r23986, MPFR_RNDN); mpfr_mul(r23988, r23983, r23987, MPFR_RNDN); return mpfr_get_d(r23988, MPFR_RNDN); } static mpfr_t r23989, r23990, r23991, r23992, r23993, r23994, r23995, r23996, r23997, r23998, r23999, r24000, r24001, r24002, r24003, r24004, r24005, r24006, r24007; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r23989); mpfr_init(r23990); mpfr_init(r23991); mpfr_init(r23992); mpfr_init(r23993); mpfr_init(r23994); mpfr_init(r23995); mpfr_init(r23996); mpfr_init(r23997); mpfr_init(r23998); mpfr_init(r23999); mpfr_init(r24000); mpfr_init(r24001); mpfr_init(r24002); mpfr_init(r24003); mpfr_init(r24004); mpfr_init(r24005); mpfr_init(r24006); mpfr_init(r24007); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r23989, x_im, MPFR_RNDN); mpfr_set_d(r23990, x_re, MPFR_RNDN); mpfr_hypot(r23991, r23989, r23990, MPFR_RNDN); mpfr_log(r23992, r23991, MPFR_RNDN); mpfr_set_d(r23993, y_im, MPFR_RNDN); mpfr_atan2(r23994, r23989, r23990, MPFR_RNDN); mpfr_set_d(r23995, y_re, MPFR_RNDN); mpfr_mul(r23996, r23994, r23995, MPFR_RNDN); mpfr_fma(r23997, r23992, r23993, r23996, MPFR_RNDN); mpfr_cos(r23998, r23997, MPFR_RNDN); mpfr_expm1(r23999, r23998, MPFR_RNDN); mpfr_exp(r24000, r23999, MPFR_RNDN); mpfr_log(r24001, r24000, MPFR_RNDN); mpfr_log1p(r24002, r24001, MPFR_RNDN); mpfr_mul(r24003, r23993, r23994, MPFR_RNDN); mpfr_mul(r24004, r23992, r23995, MPFR_RNDN); mpfr_sub(r24005, r24003, r24004, MPFR_RNDN); mpfr_exp(r24006, r24005, MPFR_RNDN); mpfr_div(r24007, r24002, r24006, MPFR_RNDN); return mpfr_get_d(r24007, MPFR_RNDN); } static mpfr_t r24008, r24009, r24010, r24011, r24012, r24013, r24014, r24015, r24016, r24017, r24018, r24019, r24020, r24021, r24022, r24023, r24024, r24025, r24026; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r24008); mpfr_init(r24009); mpfr_init(r24010); mpfr_init(r24011); mpfr_init(r24012); mpfr_init(r24013); mpfr_init(r24014); mpfr_init(r24015); mpfr_init(r24016); mpfr_init(r24017); mpfr_init(r24018); mpfr_init(r24019); mpfr_init(r24020); mpfr_init(r24021); mpfr_init(r24022); mpfr_init(r24023); mpfr_init(r24024); mpfr_init(r24025); mpfr_init(r24026); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r24008, x_im, MPFR_RNDN); mpfr_set_d(r24009, x_re, MPFR_RNDN); mpfr_hypot(r24010, r24008, r24009, MPFR_RNDN); mpfr_log(r24011, r24010, MPFR_RNDN); mpfr_set_d(r24012, y_im, MPFR_RNDN); mpfr_atan2(r24013, r24008, r24009, MPFR_RNDN); mpfr_set_d(r24014, y_re, MPFR_RNDN); mpfr_mul(r24015, r24013, r24014, MPFR_RNDN); mpfr_fma(r24016, r24011, r24012, r24015, MPFR_RNDN); mpfr_cos(r24017, r24016, MPFR_RNDN); mpfr_expm1(r24018, r24017, MPFR_RNDN); mpfr_exp(r24019, r24018, MPFR_RNDN); mpfr_log(r24020, r24019, MPFR_RNDN); mpfr_log1p(r24021, r24020, MPFR_RNDN); mpfr_mul(r24022, r24012, r24013, MPFR_RNDN); mpfr_mul(r24023, r24011, r24014, MPFR_RNDN); mpfr_sub(r24024, r24022, r24023, MPFR_RNDN); mpfr_exp(r24025, r24024, MPFR_RNDN); mpfr_div(r24026, r24021, r24025, MPFR_RNDN); return mpfr_get_d(r24026, MPFR_RNDN); }