#include #include #include #include #include char *name = "_divideComplex, imaginary part"; double f_if(float x_re, float x_im, float y_re, float y_im) { float r25789 = x_im; float r25790 = y_re; float r25791 = r25789 * r25790; float r25792 = x_re; float r25793 = y_im; float r25794 = r25792 * r25793; float r25795 = r25791 - r25794; float r25796 = r25790 * r25790; float r25797 = r25793 * r25793; float r25798 = r25796 + r25797; float r25799 = r25795 / r25798; return r25799; } double f_id(double x_re, double x_im, double y_re, double y_im) { double r25800 = x_im; double r25801 = y_re; double r25802 = r25800 * r25801; double r25803 = x_re; double r25804 = y_im; double r25805 = r25803 * r25804; double r25806 = r25802 - r25805; double r25807 = r25801 * r25801; double r25808 = r25804 * r25804; double r25809 = r25807 + r25808; double r25810 = r25806 / r25809; return r25810; } double f_of(float x_re, float x_im, float y_re, float y_im) { float r25811 = x_im; float r25812 = y_re; float r25813 = r25811 * r25812; float r25814 = x_re; float r25815 = y_im; float r25816 = r25814 * r25815; float r25817 = r25813 - r25816; float r25818 = r25812 * r25812; float r25819 = r25815 * r25815; float r25820 = r25818 + r25819; float r25821 = sqrt(r25820); float r25822 = r25817 / r25821; float r25823 = r25819 + r25818; float r25824 = sqrt(r25823); float r25825 = r25823 * r25824; float r25826 = cbrt(r25825); float r25827 = r25822 / r25826; float r25828 = 1.7828154783041087e+308; bool r25829 = r25827 <= r25828; float r25830 = r25822 / r25821; float r25831 = r25812 / r25814; float r25832 = r25815 / r25831; float r25833 = r25811 - r25832; float r25834 = r25833 / r25824; float r25835 = r25829 ? r25830 : r25834; return r25835; } double f_od(double x_re, double x_im, double y_re, double y_im) { double r25836 = x_im; double r25837 = y_re; double r25838 = r25836 * r25837; double r25839 = x_re; double r25840 = y_im; double r25841 = r25839 * r25840; double r25842 = r25838 - r25841; double r25843 = r25837 * r25837; double r25844 = r25840 * r25840; double r25845 = r25843 + r25844; double r25846 = sqrt(r25845); double r25847 = r25842 / r25846; double r25848 = r25844 + r25843; double r25849 = sqrt(r25848); double r25850 = r25848 * r25849; double r25851 = cbrt(r25850); double r25852 = r25847 / r25851; double r25853 = 1.7828154783041087e+308; bool r25854 = r25852 <= r25853; double r25855 = r25847 / r25846; double r25856 = r25837 / r25839; double r25857 = r25840 / r25856; double r25858 = r25836 - r25857; double r25859 = r25858 / r25849; double r25860 = r25854 ? r25855 : r25859; return r25860; } 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 r25861, r25862, r25863, r25864, r25865, r25866, r25867, r25868, r25869, r25870, r25871; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25861); mpfr_init(r25862); mpfr_init(r25863); mpfr_init(r25864); mpfr_init(r25865); mpfr_init(r25866); mpfr_init(r25867); mpfr_init(r25868); mpfr_init(r25869); mpfr_init(r25870); mpfr_init(r25871); } double f_im(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r25861, x_im, MPFR_RNDN); mpfr_set_d(r25862, y_re, MPFR_RNDN); mpfr_mul(r25863, r25861, r25862, MPFR_RNDN); mpfr_set_d(r25864, x_re, MPFR_RNDN); mpfr_set_d(r25865, y_im, MPFR_RNDN); mpfr_mul(r25866, r25864, r25865, MPFR_RNDN); mpfr_sub(r25867, r25863, r25866, MPFR_RNDN); mpfr_mul(r25868, r25862, r25862, MPFR_RNDN); mpfr_mul(r25869, r25865, r25865, MPFR_RNDN); mpfr_add(r25870, r25868, r25869, MPFR_RNDN); mpfr_div(r25871, r25867, r25870, MPFR_RNDN); return mpfr_get_d(r25871, MPFR_RNDN); } static mpfr_t r25872, r25873, r25874, r25875, r25876, r25877, r25878, r25879, r25880, r25881, r25882, r25883, r25884, r25885, r25886, r25887, r25888, r25889, r25890, r25891, r25892, r25893, r25894, r25895, r25896; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25872); mpfr_init(r25873); mpfr_init(r25874); mpfr_init(r25875); mpfr_init(r25876); mpfr_init(r25877); mpfr_init(r25878); mpfr_init(r25879); mpfr_init(r25880); mpfr_init(r25881); mpfr_init(r25882); mpfr_init(r25883); mpfr_init(r25884); mpfr_init(r25885); mpfr_init(r25886); mpfr_init(r25887); mpfr_init(r25888); mpfr_init_set_str(r25889, "1.7828154783041087e+308", 10, MPFR_RNDN); mpfr_init(r25890); mpfr_init(r25891); mpfr_init(r25892); mpfr_init(r25893); mpfr_init(r25894); mpfr_init(r25895); mpfr_init(r25896); } double f_fm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r25872, x_im, MPFR_RNDN); mpfr_set_d(r25873, y_re, MPFR_RNDN); mpfr_mul(r25874, r25872, r25873, MPFR_RNDN); mpfr_set_d(r25875, x_re, MPFR_RNDN); mpfr_set_d(r25876, y_im, MPFR_RNDN); mpfr_mul(r25877, r25875, r25876, MPFR_RNDN); mpfr_sub(r25878, r25874, r25877, MPFR_RNDN); mpfr_mul(r25879, r25873, r25873, MPFR_RNDN); mpfr_mul(r25880, r25876, r25876, MPFR_RNDN); mpfr_add(r25881, r25879, r25880, MPFR_RNDN); mpfr_sqrt(r25882, r25881, MPFR_RNDN); mpfr_div(r25883, r25878, r25882, MPFR_RNDN); mpfr_add(r25884, r25880, r25879, MPFR_RNDN); mpfr_sqrt(r25885, r25884, MPFR_RNDN); mpfr_mul(r25886, r25884, r25885, MPFR_RNDN); mpfr_cbrt(r25887, r25886, MPFR_RNDN); mpfr_div(r25888, r25883, r25887, MPFR_RNDN); ; mpfr_set_si(r25890, mpfr_cmp(r25888, r25889) <= 0, MPFR_RNDN); mpfr_div(r25891, r25883, r25882, MPFR_RNDN); mpfr_div(r25892, r25873, r25875, MPFR_RNDN); mpfr_div(r25893, r25876, r25892, MPFR_RNDN); mpfr_sub(r25894, r25872, r25893, MPFR_RNDN); mpfr_div(r25895, r25894, r25885, MPFR_RNDN); if (mpfr_get_si(r25890, MPFR_RNDN)) { mpfr_set(r25896, r25891, MPFR_RNDN); } else { mpfr_set(r25896, r25895, MPFR_RNDN); }; return mpfr_get_d(r25896, MPFR_RNDN); } static mpfr_t r25897, r25898, r25899, r25900, r25901, r25902, r25903, r25904, r25905, r25906, r25907, r25908, r25909, r25910, r25911, r25912, r25913, r25914, r25915, r25916, r25917, r25918, r25919, r25920, r25921; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25897); mpfr_init(r25898); mpfr_init(r25899); mpfr_init(r25900); mpfr_init(r25901); mpfr_init(r25902); mpfr_init(r25903); mpfr_init(r25904); mpfr_init(r25905); mpfr_init(r25906); mpfr_init(r25907); mpfr_init(r25908); mpfr_init(r25909); mpfr_init(r25910); mpfr_init(r25911); mpfr_init(r25912); mpfr_init(r25913); mpfr_init_set_str(r25914, "1.7828154783041087e+308", 10, MPFR_RNDN); mpfr_init(r25915); mpfr_init(r25916); mpfr_init(r25917); mpfr_init(r25918); mpfr_init(r25919); mpfr_init(r25920); mpfr_init(r25921); } double f_dm(double x_re, double x_im, double y_re, double y_im) { mpfr_set_d(r25897, x_im, MPFR_RNDN); mpfr_set_d(r25898, y_re, MPFR_RNDN); mpfr_mul(r25899, r25897, r25898, MPFR_RNDN); mpfr_set_d(r25900, x_re, MPFR_RNDN); mpfr_set_d(r25901, y_im, MPFR_RNDN); mpfr_mul(r25902, r25900, r25901, MPFR_RNDN); mpfr_sub(r25903, r25899, r25902, MPFR_RNDN); mpfr_mul(r25904, r25898, r25898, MPFR_RNDN); mpfr_mul(r25905, r25901, r25901, MPFR_RNDN); mpfr_add(r25906, r25904, r25905, MPFR_RNDN); mpfr_sqrt(r25907, r25906, MPFR_RNDN); mpfr_div(r25908, r25903, r25907, MPFR_RNDN); mpfr_add(r25909, r25905, r25904, MPFR_RNDN); mpfr_sqrt(r25910, r25909, MPFR_RNDN); mpfr_mul(r25911, r25909, r25910, MPFR_RNDN); mpfr_cbrt(r25912, r25911, MPFR_RNDN); mpfr_div(r25913, r25908, r25912, MPFR_RNDN); ; mpfr_set_si(r25915, mpfr_cmp(r25913, r25914) <= 0, MPFR_RNDN); mpfr_div(r25916, r25908, r25907, MPFR_RNDN); mpfr_div(r25917, r25898, r25900, MPFR_RNDN); mpfr_div(r25918, r25901, r25917, MPFR_RNDN); mpfr_sub(r25919, r25897, r25918, MPFR_RNDN); mpfr_div(r25920, r25919, r25910, MPFR_RNDN); if (mpfr_get_si(r25915, MPFR_RNDN)) { mpfr_set(r25921, r25916, MPFR_RNDN); } else { mpfr_set(r25921, r25920, MPFR_RNDN); }; return mpfr_get_d(r25921, MPFR_RNDN); }