#include #include #include #include #include char *name = "Jmat.Real.erfi, branch x less than or equal to 0.5"; double f_if(float x) { float r26766 = 1; float r26767 = atan2(1.0, 0.0); float r26768 = sqrt(r26767); float r26769 = r26766 / r26768; float r26770 = 2; float r26771 = x; float r26772 = fabs(r26771); float r26773 = r26770 * r26772; float r26774 = 3; float r26775 = r26770 / r26774; float r26776 = r26772 * r26772; float r26777 = r26776 * r26772; float r26778 = r26775 * r26777; float r26779 = r26773 + r26778; float r26780 = 5; float r26781 = r26766 / r26780; float r26782 = r26777 * r26772; float r26783 = r26782 * r26772; float r26784 = r26781 * r26783; float r26785 = r26779 + r26784; float r26786 = 21; float r26787 = r26766 / r26786; float r26788 = r26783 * r26772; float r26789 = r26788 * r26772; float r26790 = r26787 * r26789; float r26791 = r26785 + r26790; float r26792 = r26769 * r26791; float r26793 = fabs(r26792); return r26793; } double f_id(double x) { double r26794 = 1; double r26795 = atan2(1.0, 0.0); double r26796 = sqrt(r26795); double r26797 = r26794 / r26796; double r26798 = 2; double r26799 = x; double r26800 = fabs(r26799); double r26801 = r26798 * r26800; double r26802 = 3; double r26803 = r26798 / r26802; double r26804 = r26800 * r26800; double r26805 = r26804 * r26800; double r26806 = r26803 * r26805; double r26807 = r26801 + r26806; double r26808 = 5; double r26809 = r26794 / r26808; double r26810 = r26805 * r26800; double r26811 = r26810 * r26800; double r26812 = r26809 * r26811; double r26813 = r26807 + r26812; double r26814 = 21; double r26815 = r26794 / r26814; double r26816 = r26811 * r26800; double r26817 = r26816 * r26800; double r26818 = r26815 * r26817; double r26819 = r26813 + r26818; double r26820 = r26797 * r26819; double r26821 = fabs(r26820); return r26821; } double f_of(float x) { float r26822 = 1; float r26823 = atan2(1.0, 0.0); float r26824 = sqrt(r26823); float r26825 = r26822 / r26824; float r26826 = 5; float r26827 = r26822 / r26826; float r26828 = x; float r26829 = fabs(r26828); float r26830 = r26827 * r26829; float r26831 = r26829 * r26829; float r26832 = r26830 * r26831; float r26833 = 2; float r26834 = 3; float r26835 = r26833 / r26834; float r26836 = r26829 * r26835; float r26837 = r26833 * r26829; float r26838 = fma(r26836, r26831, r26837); float r26839 = fma(r26832, r26831, r26838); float r26840 = pow(r26829, r26834); float r26841 = r26840 * r26840; float r26842 = 21; float r26843 = r26829 / r26842; float r26844 = pow(r26843, r26822); float r26845 = r26844 * r26825; float r26846 = r26841 * r26845; float r26847 = fma(r26825, r26839, r26846); float r26848 = fabs(r26847); return r26848; } double f_od(double x) { double r26849 = 1; double r26850 = atan2(1.0, 0.0); double r26851 = sqrt(r26850); double r26852 = r26849 / r26851; double r26853 = 5; double r26854 = r26849 / r26853; double r26855 = x; double r26856 = fabs(r26855); double r26857 = r26854 * r26856; double r26858 = r26856 * r26856; double r26859 = r26857 * r26858; double r26860 = 2; double r26861 = 3; double r26862 = r26860 / r26861; double r26863 = r26856 * r26862; double r26864 = r26860 * r26856; double r26865 = fma(r26863, r26858, r26864); double r26866 = fma(r26859, r26858, r26865); double r26867 = pow(r26856, r26861); double r26868 = r26867 * r26867; double r26869 = 21; double r26870 = r26856 / r26869; double r26871 = pow(r26870, r26849); double r26872 = r26871 * r26852; double r26873 = r26868 * r26872; double r26874 = fma(r26852, r26866, r26873); double r26875 = fabs(r26874); return r26875; } 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 r26876, r26877, r26878, r26879, r26880, r26881, r26882, r26883, r26884, r26885, r26886, r26887, r26888, r26889, r26890, r26891, r26892, r26893, r26894, r26895, r26896, r26897, r26898, r26899, r26900, r26901, r26902, r26903; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r26876, "1", 10, MPFR_RNDN); mpfr_init(r26877); mpfr_init(r26878); mpfr_init(r26879); mpfr_init_set_str(r26880, "2", 10, MPFR_RNDN); mpfr_init(r26881); mpfr_init(r26882); mpfr_init(r26883); mpfr_init_set_str(r26884, "3", 10, MPFR_RNDN); mpfr_init(r26885); mpfr_init(r26886); mpfr_init(r26887); mpfr_init(r26888); mpfr_init(r26889); mpfr_init_set_str(r26890, "5", 10, MPFR_RNDN); mpfr_init(r26891); mpfr_init(r26892); mpfr_init(r26893); mpfr_init(r26894); mpfr_init(r26895); mpfr_init_set_str(r26896, "21", 10, MPFR_RNDN); mpfr_init(r26897); mpfr_init(r26898); mpfr_init(r26899); mpfr_init(r26900); mpfr_init(r26901); mpfr_init(r26902); mpfr_init(r26903); } double f_im(double x) { ; mpfr_const_pi(r26877, MPFR_RNDN); mpfr_sqrt(r26878, r26877, MPFR_RNDN); mpfr_div(r26879, r26876, r26878, MPFR_RNDN); ; mpfr_set_d(r26881, x, MPFR_RNDN); mpfr_abs(r26882, r26881, MPFR_RNDN); mpfr_mul(r26883, r26880, r26882, MPFR_RNDN); ; mpfr_div(r26885, r26880, r26884, MPFR_RNDN); mpfr_mul(r26886, r26882, r26882, MPFR_RNDN); mpfr_mul(r26887, r26886, r26882, MPFR_RNDN); mpfr_mul(r26888, r26885, r26887, MPFR_RNDN); mpfr_add(r26889, r26883, r26888, MPFR_RNDN); ; mpfr_div(r26891, r26876, r26890, MPFR_RNDN); mpfr_mul(r26892, r26887, r26882, MPFR_RNDN); mpfr_mul(r26893, r26892, r26882, MPFR_RNDN); mpfr_mul(r26894, r26891, r26893, MPFR_RNDN); mpfr_add(r26895, r26889, r26894, MPFR_RNDN); ; mpfr_div(r26897, r26876, r26896, MPFR_RNDN); mpfr_mul(r26898, r26893, r26882, MPFR_RNDN); mpfr_mul(r26899, r26898, r26882, MPFR_RNDN); mpfr_mul(r26900, r26897, r26899, MPFR_RNDN); mpfr_add(r26901, r26895, r26900, MPFR_RNDN); mpfr_mul(r26902, r26879, r26901, MPFR_RNDN); mpfr_abs(r26903, r26902, MPFR_RNDN); return mpfr_get_d(r26903, MPFR_RNDN); } static mpfr_t r26904, r26905, r26906, r26907, r26908, r26909, r26910, r26911, r26912, r26913, r26914, r26915, r26916, r26917, r26918, r26919, r26920, r26921, r26922, r26923, r26924, r26925, r26926, r26927, r26928, r26929, r26930; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26904, "1", 10, MPFR_RNDN); mpfr_init(r26905); mpfr_init(r26906); mpfr_init(r26907); mpfr_init_set_str(r26908, "5", 10, MPFR_RNDN); mpfr_init(r26909); mpfr_init(r26910); mpfr_init(r26911); mpfr_init(r26912); mpfr_init(r26913); mpfr_init(r26914); mpfr_init_set_str(r26915, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26916, "3", 10, MPFR_RNDN); mpfr_init(r26917); mpfr_init(r26918); mpfr_init(r26919); mpfr_init(r26920); mpfr_init(r26921); mpfr_init(r26922); mpfr_init(r26923); mpfr_init_set_str(r26924, "21", 10, MPFR_RNDN); mpfr_init(r26925); mpfr_init(r26926); mpfr_init(r26927); mpfr_init(r26928); mpfr_init(r26929); mpfr_init(r26930); } double f_fm(double x) { ; mpfr_const_pi(r26905, MPFR_RNDN); mpfr_sqrt(r26906, r26905, MPFR_RNDN); mpfr_div(r26907, r26904, r26906, MPFR_RNDN); ; mpfr_div(r26909, r26904, r26908, MPFR_RNDN); mpfr_set_d(r26910, x, MPFR_RNDN); mpfr_abs(r26911, r26910, MPFR_RNDN); mpfr_mul(r26912, r26909, r26911, MPFR_RNDN); mpfr_mul(r26913, r26911, r26911, MPFR_RNDN); mpfr_mul(r26914, r26912, r26913, MPFR_RNDN); ; ; mpfr_div(r26917, r26915, r26916, MPFR_RNDN); mpfr_mul(r26918, r26911, r26917, MPFR_RNDN); mpfr_mul(r26919, r26915, r26911, MPFR_RNDN); mpfr_fma(r26920, r26918, r26913, r26919, MPFR_RNDN); mpfr_fma(r26921, r26914, r26913, r26920, MPFR_RNDN); mpfr_pow(r26922, r26911, r26916, MPFR_RNDN); mpfr_mul(r26923, r26922, r26922, MPFR_RNDN); ; mpfr_div(r26925, r26911, r26924, MPFR_RNDN); mpfr_pow(r26926, r26925, r26904, MPFR_RNDN); mpfr_mul(r26927, r26926, r26907, MPFR_RNDN); mpfr_mul(r26928, r26923, r26927, MPFR_RNDN); mpfr_fma(r26929, r26907, r26921, r26928, MPFR_RNDN); mpfr_abs(r26930, r26929, MPFR_RNDN); return mpfr_get_d(r26930, MPFR_RNDN); } static mpfr_t r26931, r26932, r26933, r26934, r26935, r26936, r26937, r26938, r26939, r26940, r26941, r26942, r26943, r26944, r26945, r26946, r26947, r26948, r26949, r26950, r26951, r26952, r26953, r26954, r26955, r26956, r26957; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26931, "1", 10, MPFR_RNDN); mpfr_init(r26932); mpfr_init(r26933); mpfr_init(r26934); mpfr_init_set_str(r26935, "5", 10, MPFR_RNDN); mpfr_init(r26936); mpfr_init(r26937); mpfr_init(r26938); mpfr_init(r26939); mpfr_init(r26940); mpfr_init(r26941); mpfr_init_set_str(r26942, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26943, "3", 10, MPFR_RNDN); mpfr_init(r26944); mpfr_init(r26945); mpfr_init(r26946); mpfr_init(r26947); mpfr_init(r26948); mpfr_init(r26949); mpfr_init(r26950); mpfr_init_set_str(r26951, "21", 10, MPFR_RNDN); mpfr_init(r26952); mpfr_init(r26953); mpfr_init(r26954); mpfr_init(r26955); mpfr_init(r26956); mpfr_init(r26957); } double f_dm(double x) { ; mpfr_const_pi(r26932, MPFR_RNDN); mpfr_sqrt(r26933, r26932, MPFR_RNDN); mpfr_div(r26934, r26931, r26933, MPFR_RNDN); ; mpfr_div(r26936, r26931, r26935, MPFR_RNDN); mpfr_set_d(r26937, x, MPFR_RNDN); mpfr_abs(r26938, r26937, MPFR_RNDN); mpfr_mul(r26939, r26936, r26938, MPFR_RNDN); mpfr_mul(r26940, r26938, r26938, MPFR_RNDN); mpfr_mul(r26941, r26939, r26940, MPFR_RNDN); ; ; mpfr_div(r26944, r26942, r26943, MPFR_RNDN); mpfr_mul(r26945, r26938, r26944, MPFR_RNDN); mpfr_mul(r26946, r26942, r26938, MPFR_RNDN); mpfr_fma(r26947, r26945, r26940, r26946, MPFR_RNDN); mpfr_fma(r26948, r26941, r26940, r26947, MPFR_RNDN); mpfr_pow(r26949, r26938, r26943, MPFR_RNDN); mpfr_mul(r26950, r26949, r26949, MPFR_RNDN); ; mpfr_div(r26952, r26938, r26951, MPFR_RNDN); mpfr_pow(r26953, r26952, r26931, MPFR_RNDN); mpfr_mul(r26954, r26953, r26934, MPFR_RNDN); mpfr_mul(r26955, r26950, r26954, MPFR_RNDN); mpfr_fma(r26956, r26934, r26948, r26955, MPFR_RNDN); mpfr_abs(r26957, r26956, MPFR_RNDN); return mpfr_get_d(r26957, MPFR_RNDN); }