#include #include #include #include #include char *name = "powComplex, imaginary part"; double f_if(float x_re, float x_im, float y_re, float y_im) { float r14868 = x_re; float r14869 = r14868 * r14868; float r14870 = x_im; float r14871 = r14870 * r14870; float r14872 = r14869 + r14871; float r14873 = sqrt(r14872); float r14874 = log(r14873); float r14875 = y_re; float r14876 = r14874 * r14875; float r14877 = atan2(r14870, r14868); float r14878 = y_im; float r14879 = r14877 * r14878; float r14880 = r14876 - r14879; float r14881 = exp(r14880); float r14882 = r14874 * r14878; float r14883 = r14877 * r14875; float r14884 = r14882 + r14883; float r14885 = sin(r14884); float r14886 = r14881 * r14885; return r14886; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r14887 = x_re; double r14888 = r14887 * r14887; double r14889 = x_im; double r14890 = r14889 * r14889; double r14891 = r14888 + r14890; double r14892 = sqrt(r14891); double r14893 = log(r14892); double r14894 = y_re; double r14895 = r14893 * r14894; double r14896 = atan2(r14889, r14887); double r14897 = y_im; double r14898 = r14896 * r14897; double r14899 = r14895 - r14898; double r14900 = exp(r14899); double r14901 = r14893 * r14897; double r14902 = r14896 * r14894; double r14903 = r14901 + r14902; double r14904 = sin(r14903); double r14905 = r14900 * r14904; return r14905; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r14906 = y_re; float r14907 = x_im; float r14908 = x_re; float r14909 = atan2(r14907, r14908); float r14910 = r14906 * r14909; float r14911 = r14908 * r14908; float r14912 = r14907 * r14907; float r14913 = r14911 + r14912; float r14914 = sqrt(r14913); float r14915 = log(r14914); float r14916 = y_im; float r14917 = r14915 * r14916; float r14918 = r14910 + r14917; float r14919 = sin(r14918); float r14920 = cbrt(r14919); float r14921 = r14920 * (r14920 * r14920); float r14922 = exp(r14916); float r14923 = pow(r14922, r14909); float r14924 = pow(r14914, r14906); float r14925 = r14923 / r14924; float r14926 = r14921 / r14925; return r14926; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r14927 = y_re; double r14928 = x_im; double r14929 = x_re; double r14930 = atan2(r14928, r14929); double r14931 = r14927 * r14930; double r14932 = r14929 * r14929; double r14933 = r14928 * r14928; double r14934 = r14932 + r14933; double r14935 = sqrt(r14934); double r14936 = log(r14935); double r14937 = y_im; double r14938 = r14936 * r14937; double r14939 = r14931 + r14938; double r14940 = sin(r14939); double r14941 = cbrt(r14940); double r14942 = r14941 * (r14941 * r14941); double r14943 = exp(r14937); double r14944 = pow(r14943, r14930); double r14945 = pow(r14935, r14927); double r14946 = r14944 / r14945; double r14947 = r14942 / r14946; return r14947; } 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 r14948, r14949, r14950, r14951, r14952, r14953, r14954, r14955, r14956, r14957, r14958, r14959, r14960, r14961, r14962, r14963, r14964, r14965, r14966; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r14948); mpfr_init(r14949); mpfr_init(r14950); mpfr_init(r14951); mpfr_init(r14952); mpfr_init(r14953); mpfr_init(r14954); mpfr_init(r14955); mpfr_init(r14956); mpfr_init(r14957); mpfr_init(r14958); mpfr_init(r14959); mpfr_init(r14960); mpfr_init(r14961); mpfr_init(r14962); mpfr_init(r14963); mpfr_init(r14964); mpfr_init(r14965); mpfr_init(r14966); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r14948, x_re, MPFR_RNDN); mpfr_mul(r14949, r14948, r14948, MPFR_RNDN); mpfr_set_d(r14950, x_im, MPFR_RNDN); mpfr_mul(r14951, r14950, r14950, MPFR_RNDN); mpfr_add(r14952, r14949, r14951, MPFR_RNDN); mpfr_sqrt(r14953, r14952, MPFR_RNDN); mpfr_log(r14954, r14953, MPFR_RNDN); mpfr_set_d(r14955, y_re, MPFR_RNDN); mpfr_mul(r14956, r14954, r14955, MPFR_RNDN); mpfr_atan2(r14957, r14950, r14948, MPFR_RNDN); mpfr_set_d(r14958, y_im, MPFR_RNDN); mpfr_mul(r14959, r14957, r14958, MPFR_RNDN); mpfr_sub(r14960, r14956, r14959, MPFR_RNDN); mpfr_exp(r14961, r14960, MPFR_RNDN); mpfr_mul(r14962, r14954, r14958, MPFR_RNDN); mpfr_mul(r14963, r14957, r14955, MPFR_RNDN); mpfr_add(r14964, r14962, r14963, MPFR_RNDN); mpfr_sin(r14965, r14964, MPFR_RNDN); mpfr_mul(r14966, r14961, r14965, MPFR_RNDN); return mpfr_get_d(r14966, MPFR_RNDN); } static mpfr_t r14967, r14968, r14969, r14970, r14971, r14972, r14973, r14974, r14975, r14976, r14977, r14978, r14979, r14980, r14981, r14982, r14983, r14984, r14985, r14986, r14987; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r14967); mpfr_init(r14968); mpfr_init(r14969); mpfr_init(r14970); mpfr_init(r14971); mpfr_init(r14972); mpfr_init(r14973); mpfr_init(r14974); mpfr_init(r14975); mpfr_init(r14976); mpfr_init(r14977); mpfr_init(r14978); mpfr_init(r14979); mpfr_init(r14980); mpfr_init(r14981); mpfr_init(r14982); mpfr_init(r14983); mpfr_init(r14984); mpfr_init(r14985); mpfr_init(r14986); mpfr_init(r14987); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r14967, y_re, MPFR_RNDN); mpfr_set_d(r14968, x_im, MPFR_RNDN); mpfr_set_d(r14969, x_re, MPFR_RNDN); mpfr_atan2(r14970, r14968, r14969, MPFR_RNDN); mpfr_mul(r14971, r14967, r14970, MPFR_RNDN); mpfr_sqr(r14972, r14969, MPFR_RNDN); mpfr_mul(r14973, r14968, r14968, MPFR_RNDN); mpfr_add(r14974, r14972, r14973, MPFR_RNDN); mpfr_sqrt(r14975, r14974, MPFR_RNDN); mpfr_log(r14976, r14975, MPFR_RNDN); mpfr_set_d(r14977, y_im, MPFR_RNDN); mpfr_mul(r14978, r14976, r14977, MPFR_RNDN); mpfr_add(r14979, r14971, r14978, MPFR_RNDN); mpfr_sin(r14980, r14979, MPFR_RNDN); mpfr_cbrt(r14981, r14980, MPFR_RNDN); mpfr_mul(r14982, r14981, r14981, MPFR_RNDN); mpfr_mul(r14982, r14982, r14981, MPFR_RNDN); mpfr_exp(r14983, r14977, MPFR_RNDN); mpfr_pow(r14984, r14983, r14970, MPFR_RNDN); mpfr_pow(r14985, r14975, r14967, MPFR_RNDN); mpfr_div(r14986, r14984, r14985, MPFR_RNDN); mpfr_div(r14987, r14982, r14986, MPFR_RNDN); return mpfr_get_d(r14987, MPFR_RNDN); } static mpfr_t r14988, r14989, r14990, r14991, r14992, r14993, r14994, r14995, r14996, r14997, r14998, r14999, r15000, r15001, r15002, r15003, r15004, r15005, r15006, r15007, r15008; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r14988); mpfr_init(r14989); mpfr_init(r14990); mpfr_init(r14991); mpfr_init(r14992); mpfr_init(r14993); mpfr_init(r14994); mpfr_init(r14995); mpfr_init(r14996); mpfr_init(r14997); mpfr_init(r14998); mpfr_init(r14999); mpfr_init(r15000); mpfr_init(r15001); mpfr_init(r15002); mpfr_init(r15003); mpfr_init(r15004); mpfr_init(r15005); mpfr_init(r15006); mpfr_init(r15007); mpfr_init(r15008); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r14988, y_re, MPFR_RNDN); mpfr_set_d(r14989, x_im, MPFR_RNDN); mpfr_set_d(r14990, x_re, MPFR_RNDN); mpfr_atan2(r14991, r14989, r14990, MPFR_RNDN); mpfr_mul(r14992, r14988, r14991, MPFR_RNDN); mpfr_sqr(r14993, r14990, MPFR_RNDN); mpfr_mul(r14994, r14989, r14989, MPFR_RNDN); mpfr_add(r14995, r14993, r14994, MPFR_RNDN); mpfr_sqrt(r14996, r14995, MPFR_RNDN); mpfr_log(r14997, r14996, MPFR_RNDN); mpfr_set_d(r14998, y_im, MPFR_RNDN); mpfr_mul(r14999, r14997, r14998, MPFR_RNDN); mpfr_add(r15000, r14992, r14999, MPFR_RNDN); mpfr_sin(r15001, r15000, MPFR_RNDN); mpfr_cbrt(r15002, r15001, MPFR_RNDN); mpfr_mul(r15003, r15002, r15002, MPFR_RNDN); mpfr_mul(r15003, r15003, r15002, MPFR_RNDN); mpfr_exp(r15004, r14998, MPFR_RNDN); mpfr_pow(r15005, r15004, r14991, MPFR_RNDN); mpfr_pow(r15006, r14996, r14988, MPFR_RNDN); mpfr_div(r15007, r15005, r15006, MPFR_RNDN); mpfr_div(r15008, r15003, r15007, MPFR_RNDN); return mpfr_get_d(r15008, MPFR_RNDN); }