#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 r23843 = x_re; float r23844 = r23843 * r23843; float r23845 = x_im; float r23846 = r23845 * r23845; float r23847 = r23844 + r23846; float r23848 = sqrt(r23847); float r23849 = log(r23848); float r23850 = y_re; float r23851 = r23849 * r23850; float r23852 = atan2(r23845, r23843); float r23853 = y_im; float r23854 = r23852 * r23853; float r23855 = r23851 - r23854; float r23856 = exp(r23855); float r23857 = r23849 * r23853; float r23858 = r23852 * r23850; float r23859 = r23857 + r23858; float r23860 = cos(r23859); float r23861 = r23856 * r23860; return r23861; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r23862 = x_re; double r23863 = r23862 * r23862; double r23864 = x_im; double r23865 = r23864 * r23864; double r23866 = r23863 + r23865; double r23867 = sqrt(r23866); double r23868 = log(r23867); double r23869 = y_re; double r23870 = r23868 * r23869; double r23871 = atan2(r23864, r23862); double r23872 = y_im; double r23873 = r23871 * r23872; double r23874 = r23870 - r23873; double r23875 = exp(r23874); double r23876 = r23868 * r23872; double r23877 = r23871 * r23869; double r23878 = r23876 + r23877; double r23879 = cos(r23878); double r23880 = r23875 * r23879; return r23880; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r23881 = y_im; float r23882 = x_im; float r23883 = x_re; float r23884 = hypot(r23882, r23883); float r23885 = log(r23884); float r23886 = y_re; float r23887 = atan2(r23882, r23883); float r23888 = r23886 * r23887; float r23889 = fma(r23881, r23885, r23888); float r23890 = cos(r23889); float r23891 = expm1(r23890); float r23892 = 3; float r23893 = pow(r23891, r23892); float r23894 = exp(r23893); float r23895 = log(r23894); float r23896 = cbrt(r23895); float r23897 = log1p(r23896); float r23898 = r23881 * r23887; float r23899 = r23885 * r23886; float r23900 = r23898 - r23899; float r23901 = exp(r23900); float r23902 = r23897 / r23901; return r23902; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r23903 = y_im; double r23904 = x_im; double r23905 = x_re; double r23906 = hypot(r23904, r23905); double r23907 = log(r23906); double r23908 = y_re; double r23909 = atan2(r23904, r23905); double r23910 = r23908 * r23909; double r23911 = fma(r23903, r23907, r23910); double r23912 = cos(r23911); double r23913 = expm1(r23912); double r23914 = 3; double r23915 = pow(r23913, r23914); double r23916 = exp(r23915); double r23917 = log(r23916); double r23918 = cbrt(r23917); double r23919 = log1p(r23918); double r23920 = r23903 * r23909; double r23921 = r23907 * r23908; double r23922 = r23920 - r23921; double r23923 = exp(r23922); double r23924 = r23919 / r23923; return r23924; } 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 r23925, r23926, r23927, r23928, r23929, r23930, r23931, r23932, r23933, r23934, r23935, r23936, r23937, r23938, r23939, r23940, r23941, r23942, r23943; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r23925); mpfr_init(r23926); mpfr_init(r23927); mpfr_init(r23928); mpfr_init(r23929); mpfr_init(r23930); mpfr_init(r23931); mpfr_init(r23932); mpfr_init(r23933); mpfr_init(r23934); mpfr_init(r23935); mpfr_init(r23936); mpfr_init(r23937); mpfr_init(r23938); mpfr_init(r23939); mpfr_init(r23940); mpfr_init(r23941); mpfr_init(r23942); mpfr_init(r23943); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r23925, x_re, MPFR_RNDN); mpfr_mul(r23926, r23925, r23925, MPFR_RNDN); mpfr_set_d(r23927, x_im, MPFR_RNDN); mpfr_mul(r23928, r23927, r23927, MPFR_RNDN); mpfr_add(r23929, r23926, r23928, MPFR_RNDN); mpfr_sqrt(r23930, r23929, MPFR_RNDN); mpfr_log(r23931, r23930, MPFR_RNDN); mpfr_set_d(r23932, y_re, MPFR_RNDN); mpfr_mul(r23933, r23931, r23932, MPFR_RNDN); mpfr_atan2(r23934, r23927, r23925, MPFR_RNDN); mpfr_set_d(r23935, y_im, MPFR_RNDN); mpfr_mul(r23936, r23934, r23935, MPFR_RNDN); mpfr_sub(r23937, r23933, r23936, MPFR_RNDN); mpfr_exp(r23938, r23937, MPFR_RNDN); mpfr_mul(r23939, r23931, r23935, MPFR_RNDN); mpfr_mul(r23940, r23934, r23932, MPFR_RNDN); mpfr_add(r23941, r23939, r23940, MPFR_RNDN); mpfr_cos(r23942, r23941, MPFR_RNDN); mpfr_mul(r23943, r23938, r23942, MPFR_RNDN); return mpfr_get_d(r23943, MPFR_RNDN); } static mpfr_t r23944, r23945, r23946, r23947, r23948, r23949, r23950, r23951, r23952, r23953, r23954, r23955, r23956, r23957, r23958, r23959, r23960, r23961, r23962, r23963, r23964, r23965; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r23944); mpfr_init(r23945); mpfr_init(r23946); mpfr_init(r23947); mpfr_init(r23948); mpfr_init(r23949); mpfr_init(r23950); mpfr_init(r23951); mpfr_init(r23952); mpfr_init(r23953); mpfr_init(r23954); mpfr_init_set_str(r23955, "3", 10, MPFR_RNDN); mpfr_init(r23956); mpfr_init(r23957); mpfr_init(r23958); mpfr_init(r23959); mpfr_init(r23960); mpfr_init(r23961); mpfr_init(r23962); mpfr_init(r23963); mpfr_init(r23964); mpfr_init(r23965); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r23944, y_im, MPFR_RNDN); mpfr_set_d(r23945, x_im, MPFR_RNDN); mpfr_set_d(r23946, x_re, MPFR_RNDN); mpfr_hypot(r23947, r23945, r23946, MPFR_RNDN); mpfr_log(r23948, r23947, MPFR_RNDN); mpfr_set_d(r23949, y_re, MPFR_RNDN); mpfr_atan2(r23950, r23945, r23946, MPFR_RNDN); mpfr_mul(r23951, r23949, r23950, MPFR_RNDN); mpfr_fma(r23952, r23944, r23948, r23951, MPFR_RNDN); mpfr_cos(r23953, r23952, MPFR_RNDN); mpfr_expm1(r23954, r23953, MPFR_RNDN); ; mpfr_pow(r23956, r23954, r23955, MPFR_RNDN); mpfr_exp(r23957, r23956, MPFR_RNDN); mpfr_log(r23958, r23957, MPFR_RNDN); mpfr_cbrt(r23959, r23958, MPFR_RNDN); mpfr_log1p(r23960, r23959, MPFR_RNDN); mpfr_mul(r23961, r23944, r23950, MPFR_RNDN); mpfr_mul(r23962, r23948, r23949, MPFR_RNDN); mpfr_sub(r23963, r23961, r23962, MPFR_RNDN); mpfr_exp(r23964, r23963, MPFR_RNDN); mpfr_div(r23965, r23960, r23964, MPFR_RNDN); return mpfr_get_d(r23965, MPFR_RNDN); } static mpfr_t r23966, r23967, r23968, r23969, r23970, r23971, r23972, r23973, r23974, r23975, r23976, r23977, r23978, r23979, r23980, r23981, r23982, r23983, r23984, r23985, r23986, r23987; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r23966); mpfr_init(r23967); mpfr_init(r23968); mpfr_init(r23969); mpfr_init(r23970); mpfr_init(r23971); mpfr_init(r23972); mpfr_init(r23973); mpfr_init(r23974); mpfr_init(r23975); mpfr_init(r23976); mpfr_init_set_str(r23977, "3", 10, MPFR_RNDN); 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); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r23966, y_im, MPFR_RNDN); mpfr_set_d(r23967, x_im, MPFR_RNDN); mpfr_set_d(r23968, x_re, MPFR_RNDN); mpfr_hypot(r23969, r23967, r23968, MPFR_RNDN); mpfr_log(r23970, r23969, MPFR_RNDN); mpfr_set_d(r23971, y_re, MPFR_RNDN); mpfr_atan2(r23972, r23967, r23968, MPFR_RNDN); mpfr_mul(r23973, r23971, r23972, MPFR_RNDN); mpfr_fma(r23974, r23966, r23970, r23973, MPFR_RNDN); mpfr_cos(r23975, r23974, MPFR_RNDN); mpfr_expm1(r23976, r23975, MPFR_RNDN); ; mpfr_pow(r23978, r23976, r23977, MPFR_RNDN); mpfr_exp(r23979, r23978, MPFR_RNDN); mpfr_log(r23980, r23979, MPFR_RNDN); mpfr_cbrt(r23981, r23980, MPFR_RNDN); mpfr_log1p(r23982, r23981, MPFR_RNDN); mpfr_mul(r23983, r23966, r23972, MPFR_RNDN); mpfr_mul(r23984, r23970, r23971, MPFR_RNDN); mpfr_sub(r23985, r23983, r23984, MPFR_RNDN); mpfr_exp(r23986, r23985, MPFR_RNDN); mpfr_div(r23987, r23982, r23986, MPFR_RNDN); return mpfr_get_d(r23987, MPFR_RNDN); }