#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 r24789 = 1; float r24790 = atan2(1.0, 0.0); float r24791 = sqrt(r24790); float r24792 = r24789 / r24791; float r24793 = 2; float r24794 = x; float r24795 = fabs(r24794); float r24796 = r24793 * r24795; float r24797 = 3; float r24798 = r24793 / r24797; float r24799 = r24795 * r24795; float r24800 = r24799 * r24795; float r24801 = r24798 * r24800; float r24802 = r24796 + r24801; float r24803 = 5; float r24804 = r24789 / r24803; float r24805 = r24800 * r24795; float r24806 = r24805 * r24795; float r24807 = r24804 * r24806; float r24808 = r24802 + r24807; float r24809 = 21; float r24810 = r24789 / r24809; float r24811 = r24806 * r24795; float r24812 = r24811 * r24795; float r24813 = r24810 * r24812; float r24814 = r24808 + r24813; float r24815 = r24792 * r24814; float r24816 = fabs(r24815); return r24816; } double f_id(double x) { double r24817 = 1; double r24818 = atan2(1.0, 0.0); double r24819 = sqrt(r24818); double r24820 = r24817 / r24819; double r24821 = 2; double r24822 = x; double r24823 = fabs(r24822); double r24824 = r24821 * r24823; double r24825 = 3; double r24826 = r24821 / r24825; double r24827 = r24823 * r24823; double r24828 = r24827 * r24823; double r24829 = r24826 * r24828; double r24830 = r24824 + r24829; double r24831 = 5; double r24832 = r24817 / r24831; double r24833 = r24828 * r24823; double r24834 = r24833 * r24823; double r24835 = r24832 * r24834; double r24836 = r24830 + r24835; double r24837 = 21; double r24838 = r24817 / r24837; double r24839 = r24834 * r24823; double r24840 = r24839 * r24823; double r24841 = r24838 * r24840; double r24842 = r24836 + r24841; double r24843 = r24820 * r24842; double r24844 = fabs(r24843); return r24844; } double f_of(float x) { float r24845 = 1; float r24846 = atan2(1.0, 0.0); float r24847 = sqrt(r24846); float r24848 = r24845 / r24847; float r24849 = 2; float r24850 = x; float r24851 = fabs(r24850); float r24852 = r24849 * r24851; float r24853 = 3; float r24854 = r24849 / r24853; float r24855 = r24851 * r24851; float r24856 = r24855 * r24851; float r24857 = r24854 * r24856; float r24858 = r24852 + r24857; float r24859 = 5; float r24860 = r24845 / r24859; float r24861 = r24856 * r24851; float r24862 = r24861 * r24851; float r24863 = r24860 * r24862; float r24864 = r24858 + r24863; float r24865 = 21; float r24866 = r24845 / r24865; float r24867 = r24862 * r24851; float r24868 = r24867 * r24851; float r24869 = r24866 * r24868; float r24870 = r24864 + r24869; float r24871 = r24848 * r24870; float r24872 = fabs(r24871); return r24872; } double f_od(double x) { double r24873 = 1; double r24874 = atan2(1.0, 0.0); double r24875 = sqrt(r24874); double r24876 = r24873 / r24875; double r24877 = 2; double r24878 = x; double r24879 = fabs(r24878); double r24880 = r24877 * r24879; double r24881 = 3; double r24882 = r24877 / r24881; double r24883 = r24879 * r24879; double r24884 = r24883 * r24879; double r24885 = r24882 * r24884; double r24886 = r24880 + r24885; double r24887 = 5; double r24888 = r24873 / r24887; double r24889 = r24884 * r24879; double r24890 = r24889 * r24879; double r24891 = r24888 * r24890; double r24892 = r24886 + r24891; double r24893 = 21; double r24894 = r24873 / r24893; double r24895 = r24890 * r24879; double r24896 = r24895 * r24879; double r24897 = r24894 * r24896; double r24898 = r24892 + r24897; double r24899 = r24876 * r24898; double r24900 = fabs(r24899); return r24900; } 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 r24901, r24902, r24903, r24904, r24905, r24906, r24907, r24908, r24909, r24910, r24911, r24912, r24913, r24914, r24915, r24916, r24917, r24918, r24919, r24920, r24921, r24922, r24923, r24924, r24925, r24926, r24927, r24928; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r24901, "1", 10, MPFR_RNDN); mpfr_init(r24902); mpfr_init(r24903); mpfr_init(r24904); mpfr_init_set_str(r24905, "2", 10, MPFR_RNDN); mpfr_init(r24906); mpfr_init(r24907); mpfr_init(r24908); mpfr_init_set_str(r24909, "3", 10, MPFR_RNDN); mpfr_init(r24910); mpfr_init(r24911); mpfr_init(r24912); mpfr_init(r24913); mpfr_init(r24914); mpfr_init_set_str(r24915, "5", 10, MPFR_RNDN); mpfr_init(r24916); mpfr_init(r24917); mpfr_init(r24918); mpfr_init(r24919); mpfr_init(r24920); mpfr_init_set_str(r24921, "21", 10, MPFR_RNDN); mpfr_init(r24922); mpfr_init(r24923); mpfr_init(r24924); mpfr_init(r24925); mpfr_init(r24926); mpfr_init(r24927); mpfr_init(r24928); } double f_im(double x) { ; mpfr_const_pi(r24902, MPFR_RNDN); mpfr_sqrt(r24903, r24902, MPFR_RNDN); mpfr_div(r24904, r24901, r24903, MPFR_RNDN); ; mpfr_set_d(r24906, x, MPFR_RNDN); mpfr_abs(r24907, r24906, MPFR_RNDN); mpfr_mul(r24908, r24905, r24907, MPFR_RNDN); ; mpfr_div(r24910, r24905, r24909, MPFR_RNDN); mpfr_mul(r24911, r24907, r24907, MPFR_RNDN); mpfr_mul(r24912, r24911, r24907, MPFR_RNDN); mpfr_mul(r24913, r24910, r24912, MPFR_RNDN); mpfr_add(r24914, r24908, r24913, MPFR_RNDN); ; mpfr_div(r24916, r24901, r24915, MPFR_RNDN); mpfr_mul(r24917, r24912, r24907, MPFR_RNDN); mpfr_mul(r24918, r24917, r24907, MPFR_RNDN); mpfr_mul(r24919, r24916, r24918, MPFR_RNDN); mpfr_add(r24920, r24914, r24919, MPFR_RNDN); ; mpfr_div(r24922, r24901, r24921, MPFR_RNDN); mpfr_mul(r24923, r24918, r24907, MPFR_RNDN); mpfr_mul(r24924, r24923, r24907, MPFR_RNDN); mpfr_mul(r24925, r24922, r24924, MPFR_RNDN); mpfr_add(r24926, r24920, r24925, MPFR_RNDN); mpfr_mul(r24927, r24904, r24926, MPFR_RNDN); mpfr_abs(r24928, r24927, MPFR_RNDN); return mpfr_get_d(r24928, MPFR_RNDN); } static mpfr_t r24929, r24930, r24931, r24932, r24933, r24934, r24935, r24936, r24937, r24938, r24939, r24940, r24941, r24942, r24943, r24944, r24945, r24946, r24947, r24948, r24949, r24950, r24951, r24952, r24953, r24954, r24955, r24956; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24929, "1", 10, MPFR_RNDN); mpfr_init(r24930); mpfr_init(r24931); mpfr_init(r24932); mpfr_init_set_str(r24933, "2", 10, MPFR_RNDN); mpfr_init(r24934); mpfr_init(r24935); mpfr_init(r24936); mpfr_init_set_str(r24937, "3", 10, MPFR_RNDN); mpfr_init(r24938); mpfr_init(r24939); mpfr_init(r24940); mpfr_init(r24941); mpfr_init(r24942); mpfr_init_set_str(r24943, "5", 10, MPFR_RNDN); mpfr_init(r24944); mpfr_init(r24945); mpfr_init(r24946); mpfr_init(r24947); mpfr_init(r24948); mpfr_init_set_str(r24949, "21", 10, MPFR_RNDN); mpfr_init(r24950); mpfr_init(r24951); mpfr_init(r24952); mpfr_init(r24953); mpfr_init(r24954); mpfr_init(r24955); mpfr_init(r24956); } double f_fm(double x) { ; mpfr_const_pi(r24930, MPFR_RNDN); mpfr_sqrt(r24931, r24930, MPFR_RNDN); mpfr_div(r24932, r24929, r24931, MPFR_RNDN); ; mpfr_set_d(r24934, x, MPFR_RNDN); mpfr_abs(r24935, r24934, MPFR_RNDN); mpfr_mul(r24936, r24933, r24935, MPFR_RNDN); ; mpfr_div(r24938, r24933, r24937, MPFR_RNDN); mpfr_mul(r24939, r24935, r24935, MPFR_RNDN); mpfr_mul(r24940, r24939, r24935, MPFR_RNDN); mpfr_mul(r24941, r24938, r24940, MPFR_RNDN); mpfr_add(r24942, r24936, r24941, MPFR_RNDN); ; mpfr_div(r24944, r24929, r24943, MPFR_RNDN); mpfr_mul(r24945, r24940, r24935, MPFR_RNDN); mpfr_mul(r24946, r24945, r24935, MPFR_RNDN); mpfr_mul(r24947, r24944, r24946, MPFR_RNDN); mpfr_add(r24948, r24942, r24947, MPFR_RNDN); ; mpfr_div(r24950, r24929, r24949, MPFR_RNDN); mpfr_mul(r24951, r24946, r24935, MPFR_RNDN); mpfr_mul(r24952, r24951, r24935, MPFR_RNDN); mpfr_mul(r24953, r24950, r24952, MPFR_RNDN); mpfr_add(r24954, r24948, r24953, MPFR_RNDN); mpfr_mul(r24955, r24932, r24954, MPFR_RNDN); mpfr_abs(r24956, r24955, MPFR_RNDN); return mpfr_get_d(r24956, MPFR_RNDN); } static mpfr_t r24957, r24958, r24959, r24960, r24961, r24962, r24963, r24964, r24965, r24966, r24967, r24968, r24969, r24970, r24971, r24972, r24973, r24974, r24975, r24976, r24977, r24978, r24979, r24980, r24981, r24982, r24983, r24984; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24957, "1", 10, MPFR_RNDN); mpfr_init(r24958); mpfr_init(r24959); mpfr_init(r24960); mpfr_init_set_str(r24961, "2", 10, MPFR_RNDN); mpfr_init(r24962); mpfr_init(r24963); mpfr_init(r24964); mpfr_init_set_str(r24965, "3", 10, MPFR_RNDN); mpfr_init(r24966); mpfr_init(r24967); mpfr_init(r24968); mpfr_init(r24969); mpfr_init(r24970); mpfr_init_set_str(r24971, "5", 10, MPFR_RNDN); mpfr_init(r24972); mpfr_init(r24973); mpfr_init(r24974); mpfr_init(r24975); mpfr_init(r24976); mpfr_init_set_str(r24977, "21", 10, MPFR_RNDN); mpfr_init(r24978); mpfr_init(r24979); mpfr_init(r24980); mpfr_init(r24981); mpfr_init(r24982); mpfr_init(r24983); mpfr_init(r24984); } double f_dm(double x) { ; mpfr_const_pi(r24958, MPFR_RNDN); mpfr_sqrt(r24959, r24958, MPFR_RNDN); mpfr_div(r24960, r24957, r24959, MPFR_RNDN); ; mpfr_set_d(r24962, x, MPFR_RNDN); mpfr_abs(r24963, r24962, MPFR_RNDN); mpfr_mul(r24964, r24961, r24963, MPFR_RNDN); ; mpfr_div(r24966, r24961, r24965, MPFR_RNDN); mpfr_mul(r24967, r24963, r24963, MPFR_RNDN); mpfr_mul(r24968, r24967, r24963, MPFR_RNDN); mpfr_mul(r24969, r24966, r24968, MPFR_RNDN); mpfr_add(r24970, r24964, r24969, MPFR_RNDN); ; mpfr_div(r24972, r24957, r24971, MPFR_RNDN); mpfr_mul(r24973, r24968, r24963, MPFR_RNDN); mpfr_mul(r24974, r24973, r24963, MPFR_RNDN); mpfr_mul(r24975, r24972, r24974, MPFR_RNDN); mpfr_add(r24976, r24970, r24975, MPFR_RNDN); ; mpfr_div(r24978, r24957, r24977, MPFR_RNDN); mpfr_mul(r24979, r24974, r24963, MPFR_RNDN); mpfr_mul(r24980, r24979, r24963, MPFR_RNDN); mpfr_mul(r24981, r24978, r24980, MPFR_RNDN); mpfr_add(r24982, r24976, r24981, MPFR_RNDN); mpfr_mul(r24983, r24960, r24982, MPFR_RNDN); mpfr_abs(r24984, r24983, MPFR_RNDN); return mpfr_get_d(r24984, MPFR_RNDN); }