#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 r14792 = x_re; float r14793 = r14792 * r14792; float r14794 = x_im; float r14795 = r14794 * r14794; float r14796 = r14793 + r14795; float r14797 = sqrt(r14796); float r14798 = log(r14797); float r14799 = y_re; float r14800 = r14798 * r14799; float r14801 = atan2(r14794, r14792); float r14802 = y_im; float r14803 = r14801 * r14802; float r14804 = r14800 - r14803; float r14805 = exp(r14804); float r14806 = r14798 * r14802; float r14807 = r14801 * r14799; float r14808 = r14806 + r14807; float r14809 = sin(r14808); float r14810 = r14805 * r14809; return r14810; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r14811 = x_re; double r14812 = r14811 * r14811; double r14813 = x_im; double r14814 = r14813 * r14813; double r14815 = r14812 + r14814; double r14816 = sqrt(r14815); double r14817 = log(r14816); double r14818 = y_re; double r14819 = r14817 * r14818; double r14820 = atan2(r14813, r14811); double r14821 = y_im; double r14822 = r14820 * r14821; double r14823 = r14819 - r14822; double r14824 = exp(r14823); double r14825 = r14817 * r14821; double r14826 = r14820 * r14818; double r14827 = r14825 + r14826; double r14828 = sin(r14827); double r14829 = r14824 * r14828; return r14829; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r14830 = y_re; float r14831 = x_im; float r14832 = x_re; float r14833 = atan2(r14831, r14832); float r14834 = r14830 * r14833; float r14835 = r14832 * r14832; float r14836 = r14831 * r14831; float r14837 = r14835 + r14836; float r14838 = sqrt(r14837); float r14839 = log(r14838); float r14840 = y_im; float r14841 = r14839 * r14840; float r14842 = r14834 + r14841; float r14843 = sin(r14842); float r14844 = exp(r14840); float r14845 = pow(r14844, r14833); float r14846 = pow(r14838, r14830); float r14847 = r14845 / r14846; float r14848 = r14843 / r14847; return r14848; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r14849 = y_re; double r14850 = x_im; double r14851 = x_re; double r14852 = atan2(r14850, r14851); double r14853 = r14849 * r14852; double r14854 = r14851 * r14851; double r14855 = r14850 * r14850; double r14856 = r14854 + r14855; double r14857 = sqrt(r14856); double r14858 = log(r14857); double r14859 = y_im; double r14860 = r14858 * r14859; double r14861 = r14853 + r14860; double r14862 = sin(r14861); double r14863 = exp(r14859); double r14864 = pow(r14863, r14852); double r14865 = pow(r14857, r14849); double r14866 = r14864 / r14865; double r14867 = r14862 / r14866; return r14867; } 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 r14868, r14869, r14870, r14871, r14872, r14873, r14874, r14875, r14876, r14877, r14878, r14879, r14880, r14881, r14882, r14883, r14884, r14885, r14886; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r14868); mpfr_init(r14869); mpfr_init(r14870); mpfr_init(r14871); mpfr_init(r14872); mpfr_init(r14873); mpfr_init(r14874); mpfr_init(r14875); mpfr_init(r14876); mpfr_init(r14877); mpfr_init(r14878); mpfr_init(r14879); mpfr_init(r14880); mpfr_init(r14881); mpfr_init(r14882); mpfr_init(r14883); mpfr_init(r14884); mpfr_init(r14885); mpfr_init(r14886); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r14868, x_re, MPFR_RNDN); mpfr_mul(r14869, r14868, r14868, MPFR_RNDN); mpfr_set_d(r14870, x_im, MPFR_RNDN); mpfr_mul(r14871, r14870, r14870, MPFR_RNDN); mpfr_add(r14872, r14869, r14871, MPFR_RNDN); mpfr_sqrt(r14873, r14872, MPFR_RNDN); mpfr_log(r14874, r14873, MPFR_RNDN); mpfr_set_d(r14875, y_re, MPFR_RNDN); mpfr_mul(r14876, r14874, r14875, MPFR_RNDN); mpfr_atan2(r14877, r14870, r14868, MPFR_RNDN); mpfr_set_d(r14878, y_im, MPFR_RNDN); mpfr_mul(r14879, r14877, r14878, MPFR_RNDN); mpfr_sub(r14880, r14876, r14879, MPFR_RNDN); mpfr_exp(r14881, r14880, MPFR_RNDN); mpfr_mul(r14882, r14874, r14878, MPFR_RNDN); mpfr_mul(r14883, r14877, r14875, MPFR_RNDN); mpfr_add(r14884, r14882, r14883, MPFR_RNDN); mpfr_sin(r14885, r14884, MPFR_RNDN); mpfr_mul(r14886, r14881, r14885, MPFR_RNDN); return mpfr_get_d(r14886, MPFR_RNDN); } static mpfr_t r14887, r14888, r14889, r14890, r14891, r14892, r14893, r14894, r14895, r14896, r14897, r14898, r14899, r14900, r14901, r14902, r14903, r14904, r14905; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r14887); mpfr_init(r14888); mpfr_init(r14889); mpfr_init(r14890); mpfr_init(r14891); mpfr_init(r14892); mpfr_init(r14893); mpfr_init(r14894); mpfr_init(r14895); mpfr_init(r14896); mpfr_init(r14897); mpfr_init(r14898); mpfr_init(r14899); mpfr_init(r14900); mpfr_init(r14901); mpfr_init(r14902); mpfr_init(r14903); mpfr_init(r14904); mpfr_init(r14905); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r14887, y_re, MPFR_RNDN); mpfr_set_d(r14888, x_im, MPFR_RNDN); mpfr_set_d(r14889, x_re, MPFR_RNDN); mpfr_atan2(r14890, r14888, r14889, MPFR_RNDN); mpfr_mul(r14891, r14887, r14890, MPFR_RNDN); mpfr_sqr(r14892, r14889, MPFR_RNDN); mpfr_mul(r14893, r14888, r14888, MPFR_RNDN); mpfr_add(r14894, r14892, r14893, MPFR_RNDN); mpfr_sqrt(r14895, r14894, MPFR_RNDN); mpfr_log(r14896, r14895, MPFR_RNDN); mpfr_set_d(r14897, y_im, MPFR_RNDN); mpfr_mul(r14898, r14896, r14897, MPFR_RNDN); mpfr_add(r14899, r14891, r14898, MPFR_RNDN); mpfr_sin(r14900, r14899, MPFR_RNDN); mpfr_exp(r14901, r14897, MPFR_RNDN); mpfr_pow(r14902, r14901, r14890, MPFR_RNDN); mpfr_pow(r14903, r14895, r14887, MPFR_RNDN); mpfr_div(r14904, r14902, r14903, MPFR_RNDN); mpfr_div(r14905, r14900, r14904, MPFR_RNDN); return mpfr_get_d(r14905, MPFR_RNDN); } static mpfr_t r14906, r14907, r14908, r14909, r14910, r14911, r14912, r14913, r14914, r14915, r14916, r14917, r14918, r14919, r14920, r14921, r14922, r14923, r14924; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r14906); mpfr_init(r14907); mpfr_init(r14908); mpfr_init(r14909); mpfr_init(r14910); mpfr_init(r14911); mpfr_init(r14912); mpfr_init(r14913); mpfr_init(r14914); mpfr_init(r14915); mpfr_init(r14916); mpfr_init(r14917); mpfr_init(r14918); mpfr_init(r14919); mpfr_init(r14920); mpfr_init(r14921); mpfr_init(r14922); mpfr_init(r14923); mpfr_init(r14924); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r14906, y_re, MPFR_RNDN); mpfr_set_d(r14907, x_im, MPFR_RNDN); mpfr_set_d(r14908, x_re, MPFR_RNDN); mpfr_atan2(r14909, r14907, r14908, MPFR_RNDN); mpfr_mul(r14910, r14906, r14909, MPFR_RNDN); mpfr_sqr(r14911, r14908, MPFR_RNDN); mpfr_mul(r14912, r14907, r14907, MPFR_RNDN); mpfr_add(r14913, r14911, r14912, MPFR_RNDN); mpfr_sqrt(r14914, r14913, MPFR_RNDN); mpfr_log(r14915, r14914, MPFR_RNDN); mpfr_set_d(r14916, y_im, MPFR_RNDN); mpfr_mul(r14917, r14915, r14916, MPFR_RNDN); mpfr_add(r14918, r14910, r14917, MPFR_RNDN); mpfr_sin(r14919, r14918, MPFR_RNDN); mpfr_exp(r14920, r14916, MPFR_RNDN); mpfr_pow(r14921, r14920, r14909, MPFR_RNDN); mpfr_pow(r14922, r14914, r14906, MPFR_RNDN); mpfr_div(r14923, r14921, r14922, MPFR_RNDN); mpfr_div(r14924, r14919, r14923, MPFR_RNDN); return mpfr_get_d(r14924, MPFR_RNDN); }