#include #include #include #include #include char *name = "math.sqrt on complex, real part"; double f_if(float re, float im) { float r18819 = 0.5f; float r18820 = 2.0f; float r18821 = re; float r18822 = r18821 * r18821; float r18823 = im; float r18824 = r18823 * r18823; float r18825 = r18822 + r18824; float r18826 = sqrt(r18825); float r18827 = r18826 + r18821; float r18828 = r18820 * r18827; float r18829 = sqrt(r18828); float r18830 = r18819 * r18829; return r18830; } double f_id(double re, double im) { double r18831 = 0.5; double r18832 = 2.0; double r18833 = re; double r18834 = r18833 * r18833; double r18835 = im; double r18836 = r18835 * r18835; double r18837 = r18834 + r18836; double r18838 = sqrt(r18837); double r18839 = r18838 + r18833; double r18840 = r18832 * r18839; double r18841 = sqrt(r18840); double r18842 = r18831 * r18841; return r18842; } double f_of(float re, float im) { float r18843 = 0.5f; float r18844 = 2.0f; float r18845 = re; float r18846 = r18845 * r18845; float r18847 = im; float r18848 = r18847 * r18847; float r18849 = r18846 + r18848; float r18850 = sqrt(r18849); float r18851 = sqrt(r18850); float r18852 = r18851 * r18851; float r18853 = r18852 + r18845; float r18854 = r18844 * r18853; float r18855 = sqrt(r18854); float r18856 = r18855 * (r18855 * r18855); float r18857 = cbrt(r18856); float r18858 = r18843 * r18857; return r18858; } double f_od(double re, double im) { double r18859 = 0.5; double r18860 = 2.0; double r18861 = re; double r18862 = r18861 * r18861; double r18863 = im; double r18864 = r18863 * r18863; double r18865 = r18862 + r18864; double r18866 = sqrt(r18865); double r18867 = sqrt(r18866); double r18868 = r18867 * r18867; double r18869 = r18868 + r18861; double r18870 = r18860 * r18869; double r18871 = sqrt(r18870); double r18872 = r18871 * (r18871 * r18871); double r18873 = cbrt(r18872); double r18874 = r18859 * r18873; return r18874; } 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 r18875, r18876, r18877, r18878, r18879, r18880, r18881, r18882, r18883, r18884, r18885, r18886; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18875, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r18876, "2.0", 10, MPFR_RNDN); mpfr_init(r18877); mpfr_init(r18878); mpfr_init(r18879); mpfr_init(r18880); mpfr_init(r18881); mpfr_init(r18882); mpfr_init(r18883); mpfr_init(r18884); mpfr_init(r18885); mpfr_init(r18886); } double f_im(double re, double im) { ; ; mpfr_set_d(r18877, re, MPFR_RNDN); mpfr_mul(r18878, r18877, r18877, MPFR_RNDN); mpfr_set_d(r18879, im, MPFR_RNDN); mpfr_mul(r18880, r18879, r18879, MPFR_RNDN); mpfr_add(r18881, r18878, r18880, MPFR_RNDN); mpfr_sqrt(r18882, r18881, MPFR_RNDN); mpfr_add(r18883, r18882, r18877, MPFR_RNDN); mpfr_mul(r18884, r18876, r18883, MPFR_RNDN); mpfr_sqrt(r18885, r18884, MPFR_RNDN); mpfr_mul(r18886, r18875, r18885, MPFR_RNDN); return mpfr_get_d(r18886, MPFR_RNDN); } static mpfr_t r18887, r18888, r18889, r18890, r18891, r18892, r18893, r18894, r18895, r18896, r18897, r18898, r18899, r18900, r18901, r18902; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18887, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r18888, "2.0", 10, MPFR_RNDN); mpfr_init(r18889); mpfr_init(r18890); mpfr_init(r18891); mpfr_init(r18892); mpfr_init(r18893); mpfr_init(r18894); mpfr_init(r18895); mpfr_init(r18896); mpfr_init(r18897); mpfr_init(r18898); mpfr_init(r18899); mpfr_init(r18900); mpfr_init(r18901); mpfr_init(r18902); } double f_fm(double re, double im) { ; ; mpfr_set_d(r18889, re, MPFR_RNDN); mpfr_sqr(r18890, r18889, MPFR_RNDN); mpfr_set_d(r18891, im, MPFR_RNDN); mpfr_mul(r18892, r18891, r18891, MPFR_RNDN); mpfr_add(r18893, r18890, r18892, MPFR_RNDN); mpfr_sqrt(r18894, r18893, MPFR_RNDN); mpfr_sqrt(r18895, r18894, MPFR_RNDN); mpfr_sqr(r18896, r18895, MPFR_RNDN); mpfr_add(r18897, r18896, r18889, MPFR_RNDN); mpfr_mul(r18898, r18888, r18897, MPFR_RNDN); mpfr_sqrt(r18899, r18898, MPFR_RNDN); mpfr_mul(r18900, r18899, r18899, MPFR_RNDN); mpfr_mul(r18900, r18900, r18899, MPFR_RNDN); mpfr_cbrt(r18901, r18900, MPFR_RNDN); mpfr_mul(r18902, r18887, r18901, MPFR_RNDN); return mpfr_get_d(r18902, MPFR_RNDN); } static mpfr_t r18903, r18904, r18905, r18906, r18907, r18908, r18909, r18910, r18911, r18912, r18913, r18914, r18915, r18916, r18917, r18918; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18903, "0.5", 10, MPFR_RNDN); mpfr_init_set_str(r18904, "2.0", 10, MPFR_RNDN); mpfr_init(r18905); mpfr_init(r18906); mpfr_init(r18907); mpfr_init(r18908); mpfr_init(r18909); mpfr_init(r18910); mpfr_init(r18911); mpfr_init(r18912); mpfr_init(r18913); mpfr_init(r18914); mpfr_init(r18915); mpfr_init(r18916); mpfr_init(r18917); mpfr_init(r18918); } double f_dm(double re, double im) { ; ; mpfr_set_d(r18905, re, MPFR_RNDN); mpfr_sqr(r18906, r18905, MPFR_RNDN); mpfr_set_d(r18907, im, MPFR_RNDN); mpfr_mul(r18908, r18907, r18907, MPFR_RNDN); mpfr_add(r18909, r18906, r18908, MPFR_RNDN); mpfr_sqrt(r18910, r18909, MPFR_RNDN); mpfr_sqrt(r18911, r18910, MPFR_RNDN); mpfr_sqr(r18912, r18911, MPFR_RNDN); mpfr_add(r18913, r18912, r18905, MPFR_RNDN); mpfr_mul(r18914, r18904, r18913, MPFR_RNDN); mpfr_sqrt(r18915, r18914, MPFR_RNDN); mpfr_mul(r18916, r18915, r18915, MPFR_RNDN); mpfr_mul(r18916, r18916, r18915, MPFR_RNDN); mpfr_cbrt(r18917, r18916, MPFR_RNDN); mpfr_mul(r18918, r18903, r18917, MPFR_RNDN); return mpfr_get_d(r18918, MPFR_RNDN); }