#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Jmat.Real.erf";

double f_if(float x) {
        float r24762 = 1;
        float r24763 = 0.3275911;
        float r24764 = x;
        float r24765 = fabs(r24764);
        float r24766 = r24763 * r24765;
        float r24767 = r24762 + r24766;
        float r24768 = r24762 / r24767;
        float r24769 = 0.254829592;
        float r24770 = -0.284496736;
        float r24771 = 1.421413741;
        float r24772 = -1.453152027;
        float r24773 = 1.061405429;
        float r24774 = r24768 * r24773;
        float r24775 = r24772 + r24774;
        float r24776 = r24768 * r24775;
        float r24777 = r24771 + r24776;
        float r24778 = r24768 * r24777;
        float r24779 = r24770 + r24778;
        float r24780 = r24768 * r24779;
        float r24781 = r24769 + r24780;
        float r24782 = r24768 * r24781;
        float r24783 = r24765 * r24765;
        float r24784 = -r24783;
        float r24785 = exp(r24784);
        float r24786 = r24782 * r24785;
        float r24787 = r24762 - r24786;
        return r24787;
}

double f_id(double x) {
        double r24788 = 1;
        double r24789 = 0.3275911;
        double r24790 = x;
        double r24791 = fabs(r24790);
        double r24792 = r24789 * r24791;
        double r24793 = r24788 + r24792;
        double r24794 = r24788 / r24793;
        double r24795 = 0.254829592;
        double r24796 = -0.284496736;
        double r24797 = 1.421413741;
        double r24798 = -1.453152027;
        double r24799 = 1.061405429;
        double r24800 = r24794 * r24799;
        double r24801 = r24798 + r24800;
        double r24802 = r24794 * r24801;
        double r24803 = r24797 + r24802;
        double r24804 = r24794 * r24803;
        double r24805 = r24796 + r24804;
        double r24806 = r24794 * r24805;
        double r24807 = r24795 + r24806;
        double r24808 = r24794 * r24807;
        double r24809 = r24791 * r24791;
        double r24810 = -r24809;
        double r24811 = exp(r24810);
        double r24812 = r24808 * r24811;
        double r24813 = r24788 - r24812;
        return r24813;
}


double f_of(float x) {
        float r24814 = 1;
        float r24815 = -r24814;
        float r24816 = x;
        float r24817 = fabs(r24816);
        float r24818 = r24817 * r24817;
        float r24819 = exp(r24818);
        float r24820 = r24815 / r24819;
        float r24821 = 1.061405429;
        float r24822 = 0.3275911;
        float r24823 = fma(r24817, r24822, r24814);
        float r24824 = r24821 / r24823;
        float r24825 = -1.453152027;
        float r24826 = r24824 + r24825;
        float r24827 = r24814 / r24823;
        float r24828 = 1.421413741;
        float r24829 = fma(r24826, r24827, r24828);
        float r24830 = -0.284496736;
        float r24831 = fma(r24829, r24827, r24830);
        float r24832 = log1p(r24827);
        float r24833 = expm1(r24832);
        float r24834 = r24827 * r24833;
        float r24835 = 0.254829592;
        float r24836 = r24835 / r24823;
        float r24837 = fma(r24831, r24834, r24836);
        float r24838 = r24820 * r24837;
        float r24839 = r24838 + r24814;
        return r24839;
}

double f_od(double x) {
        double r24840 = 1;
        double r24841 = -r24840;
        double r24842 = x;
        double r24843 = fabs(r24842);
        double r24844 = r24843 * r24843;
        double r24845 = exp(r24844);
        double r24846 = r24841 / r24845;
        double r24847 = 1.061405429;
        double r24848 = 0.3275911;
        double r24849 = fma(r24843, r24848, r24840);
        double r24850 = r24847 / r24849;
        double r24851 = -1.453152027;
        double r24852 = r24850 + r24851;
        double r24853 = r24840 / r24849;
        double r24854 = 1.421413741;
        double r24855 = fma(r24852, r24853, r24854);
        double r24856 = -0.284496736;
        double r24857 = fma(r24855, r24853, r24856);
        double r24858 = log1p(r24853);
        double r24859 = expm1(r24858);
        double r24860 = r24853 * r24859;
        double r24861 = 0.254829592;
        double r24862 = r24861 / r24849;
        double r24863 = fma(r24857, r24860, r24862);
        double r24864 = r24846 * r24863;
        double r24865 = r24864 + r24840;
        return r24865;
}

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 r24866, r24867, r24868, r24869, r24870, r24871, r24872, r24873, r24874, r24875, r24876, r24877, r24878, r24879, r24880, r24881, r24882, r24883, r24884, r24885, r24886, r24887, r24888, r24889, r24890, r24891;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r24866, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r24867, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r24868);
        mpfr_init(r24869);
        mpfr_init(r24870);
        mpfr_init(r24871);
        mpfr_init(r24872);
        mpfr_init_set_str(r24873, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r24874, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r24875, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r24876, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r24877, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r24878);
        mpfr_init(r24879);
        mpfr_init(r24880);
        mpfr_init(r24881);
        mpfr_init(r24882);
        mpfr_init(r24883);
        mpfr_init(r24884);
        mpfr_init(r24885);
        mpfr_init(r24886);
        mpfr_init(r24887);
        mpfr_init(r24888);
        mpfr_init(r24889);
        mpfr_init(r24890);
        mpfr_init(r24891);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r24868, x, MPFR_RNDN);
        mpfr_abs(r24869, r24868, MPFR_RNDN);
        mpfr_mul(r24870, r24867, r24869, MPFR_RNDN);
        mpfr_add(r24871, r24866, r24870, MPFR_RNDN);
        mpfr_div(r24872, r24866, r24871, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r24878, r24872, r24877, MPFR_RNDN);
        mpfr_add(r24879, r24876, r24878, MPFR_RNDN);
        mpfr_mul(r24880, r24872, r24879, MPFR_RNDN);
        mpfr_add(r24881, r24875, r24880, MPFR_RNDN);
        mpfr_mul(r24882, r24872, r24881, MPFR_RNDN);
        mpfr_add(r24883, r24874, r24882, MPFR_RNDN);
        mpfr_mul(r24884, r24872, r24883, MPFR_RNDN);
        mpfr_add(r24885, r24873, r24884, MPFR_RNDN);
        mpfr_mul(r24886, r24872, r24885, MPFR_RNDN);
        mpfr_mul(r24887, r24869, r24869, MPFR_RNDN);
        mpfr_neg(r24888, r24887, MPFR_RNDN);
        mpfr_exp(r24889, r24888, MPFR_RNDN);
        mpfr_mul(r24890, r24886, r24889, MPFR_RNDN);
        mpfr_sub(r24891, r24866, r24890, MPFR_RNDN);
        return mpfr_get_d(r24891, MPFR_RNDN);
}

static mpfr_t r24892, r24893, r24894, r24895, r24896, r24897, r24898, r24899, r24900, r24901, r24902, r24903, r24904, r24905, r24906, r24907, r24908, r24909, r24910, r24911, r24912, r24913, r24914, r24915, r24916, r24917;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r24892, "1", 10, MPFR_RNDN);
        mpfr_init(r24893);
        mpfr_init(r24894);
        mpfr_init(r24895);
        mpfr_init(r24896);
        mpfr_init(r24897);
        mpfr_init(r24898);
        mpfr_init_set_str(r24899, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r24900, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r24901);
        mpfr_init(r24902);
        mpfr_init_set_str(r24903, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r24904);
        mpfr_init(r24905);
        mpfr_init_set_str(r24906, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r24907);
        mpfr_init_set_str(r24908, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r24909);
        mpfr_init(r24910);
        mpfr_init(r24911);
        mpfr_init(r24912);
        mpfr_init_set_str(r24913, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r24914);
        mpfr_init(r24915);
        mpfr_init(r24916);
        mpfr_init(r24917);
}

double f_fm(double x) {
        ;
        mpfr_neg(r24893, r24892, MPFR_RNDN);
        mpfr_set_d(r24894, x, MPFR_RNDN);
        mpfr_abs(r24895, r24894, MPFR_RNDN);
        mpfr_mul(r24896, r24895, r24895, MPFR_RNDN);
        mpfr_exp(r24897, r24896, MPFR_RNDN);
        mpfr_div(r24898, r24893, r24897, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r24901, r24895, r24900, r24892, MPFR_RNDN);
        mpfr_div(r24902, r24899, r24901, MPFR_RNDN);
        ;
        mpfr_add(r24904, r24902, r24903, MPFR_RNDN);
        mpfr_div(r24905, r24892, r24901, MPFR_RNDN);
        ;
        mpfr_fma(r24907, r24904, r24905, r24906, MPFR_RNDN);
        ;
        mpfr_fma(r24909, r24907, r24905, r24908, MPFR_RNDN);
        mpfr_log1p(r24910, r24905, MPFR_RNDN);
        mpfr_expm1(r24911, r24910, MPFR_RNDN);
        mpfr_mul(r24912, r24905, r24911, MPFR_RNDN);
        ;
        mpfr_div(r24914, r24913, r24901, MPFR_RNDN);
        mpfr_fma(r24915, r24909, r24912, r24914, MPFR_RNDN);
        mpfr_mul(r24916, r24898, r24915, MPFR_RNDN);
        mpfr_add(r24917, r24916, r24892, MPFR_RNDN);
        return mpfr_get_d(r24917, MPFR_RNDN);
}

static mpfr_t r24918, r24919, r24920, r24921, r24922, r24923, r24924, r24925, r24926, r24927, r24928, r24929, r24930, r24931, r24932, r24933, r24934, r24935, r24936, r24937, r24938, r24939, r24940, r24941, r24942, r24943;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r24918, "1", 10, MPFR_RNDN);
        mpfr_init(r24919);
        mpfr_init(r24920);
        mpfr_init(r24921);
        mpfr_init(r24922);
        mpfr_init(r24923);
        mpfr_init(r24924);
        mpfr_init_set_str(r24925, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r24926, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r24927);
        mpfr_init(r24928);
        mpfr_init_set_str(r24929, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r24930);
        mpfr_init(r24931);
        mpfr_init_set_str(r24932, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r24933);
        mpfr_init_set_str(r24934, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r24935);
        mpfr_init(r24936);
        mpfr_init(r24937);
        mpfr_init(r24938);
        mpfr_init_set_str(r24939, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r24940);
        mpfr_init(r24941);
        mpfr_init(r24942);
        mpfr_init(r24943);
}

double f_dm(double x) {
        ;
        mpfr_neg(r24919, r24918, MPFR_RNDN);
        mpfr_set_d(r24920, x, MPFR_RNDN);
        mpfr_abs(r24921, r24920, MPFR_RNDN);
        mpfr_mul(r24922, r24921, r24921, MPFR_RNDN);
        mpfr_exp(r24923, r24922, MPFR_RNDN);
        mpfr_div(r24924, r24919, r24923, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r24927, r24921, r24926, r24918, MPFR_RNDN);
        mpfr_div(r24928, r24925, r24927, MPFR_RNDN);
        ;
        mpfr_add(r24930, r24928, r24929, MPFR_RNDN);
        mpfr_div(r24931, r24918, r24927, MPFR_RNDN);
        ;
        mpfr_fma(r24933, r24930, r24931, r24932, MPFR_RNDN);
        ;
        mpfr_fma(r24935, r24933, r24931, r24934, MPFR_RNDN);
        mpfr_log1p(r24936, r24931, MPFR_RNDN);
        mpfr_expm1(r24937, r24936, MPFR_RNDN);
        mpfr_mul(r24938, r24931, r24937, MPFR_RNDN);
        ;
        mpfr_div(r24940, r24939, r24927, MPFR_RNDN);
        mpfr_fma(r24941, r24935, r24938, r24940, MPFR_RNDN);
        mpfr_mul(r24942, r24924, r24941, MPFR_RNDN);
        mpfr_add(r24943, r24942, r24918, MPFR_RNDN);
        return mpfr_get_d(r24943, MPFR_RNDN);
}