#include #include #include #include #include char *name = "Data.Number.Erf:$cinvnormcdf from erf-2.0.0.0, A"; double f_if(float x, float y, float z, float t) { float r54880 = x; float r54881 = 0.5f; float r54882 = r54880 * r54881; float r54883 = y; float r54884 = r54882 - r54883; float r54885 = z; float r54886 = 2.0f; float r54887 = r54885 * r54886; float r54888 = sqrt(r54887); float r54889 = r54884 * r54888; float r54890 = t; float r54891 = r54890 * r54890; float r54892 = r54891 / r54886; float r54893 = exp(r54892); float r54894 = r54889 * r54893; return r54894; } double f_id(double x, double y, double z, double t) { double r54895 = x; double r54896 = 0.5; double r54897 = r54895 * r54896; double r54898 = y; double r54899 = r54897 - r54898; double r54900 = z; double r54901 = 2.0; double r54902 = r54900 * r54901; double r54903 = sqrt(r54902); double r54904 = r54899 * r54903; double r54905 = t; double r54906 = r54905 * r54905; double r54907 = r54906 / r54901; double r54908 = exp(r54907); double r54909 = r54904 * r54908; return r54909; } double f_of(float x, float y, float z, float t) { float r54910 = x; float r54911 = 0.5f; float r54912 = r54910 * r54911; float r54913 = y; float r54914 = r54912 - r54913; float r54915 = z; float r54916 = 2.0f; float r54917 = r54915 * r54916; float r54918 = sqrt(r54917); float r54919 = r54914 * r54918; float r54920 = t; float r54921 = r54920 * r54920; float r54922 = r54921 / r54916; float r54923 = exp(r54922); float r54924 = r54919 * r54923; return r54924; } double f_od(double x, double y, double z, double t) { double r54925 = x; double r54926 = 0.5; double r54927 = r54925 * r54926; double r54928 = y; double r54929 = r54927 - r54928; double r54930 = z; double r54931 = 2.0; double r54932 = r54930 * r54931; double r54933 = sqrt(r54932); double r54934 = r54929 * r54933; double r54935 = t; double r54936 = r54935 * r54935; double r54937 = r54936 / r54931; double r54938 = exp(r54937); double r54939 = r54934 * r54938; return r54939; } 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 r54940, r54941, r54942, r54943, r54944, r54945, r54946, r54947, r54948, r54949, r54950, r54951, r54952, r54953, r54954; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r54940); mpfr_init_set_str(r54941, "0.5", 10, MPFR_RNDN); mpfr_init(r54942); mpfr_init(r54943); mpfr_init(r54944); mpfr_init(r54945); mpfr_init_set_str(r54946, "2.0", 10, MPFR_RNDN); mpfr_init(r54947); mpfr_init(r54948); mpfr_init(r54949); mpfr_init(r54950); mpfr_init(r54951); mpfr_init(r54952); mpfr_init(r54953); mpfr_init(r54954); } double f_im(double x, double y, double z, double t) { mpfr_set_d(r54940, x, MPFR_RNDN); ; mpfr_mul(r54942, r54940, r54941, MPFR_RNDN); mpfr_set_d(r54943, y, MPFR_RNDN); mpfr_sub(r54944, r54942, r54943, MPFR_RNDN); mpfr_set_d(r54945, z, MPFR_RNDN); ; mpfr_mul(r54947, r54945, r54946, MPFR_RNDN); mpfr_sqrt(r54948, r54947, MPFR_RNDN); mpfr_mul(r54949, r54944, r54948, MPFR_RNDN); mpfr_set_d(r54950, t, MPFR_RNDN); mpfr_mul(r54951, r54950, r54950, MPFR_RNDN); mpfr_div(r54952, r54951, r54946, MPFR_RNDN); mpfr_exp(r54953, r54952, MPFR_RNDN); mpfr_mul(r54954, r54949, r54953, MPFR_RNDN); return mpfr_get_d(r54954, MPFR_RNDN); } static mpfr_t r54955, r54956, r54957, r54958, r54959, r54960, r54961, r54962, r54963, r54964, r54965, r54966, r54967, r54968, r54969; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r54955); mpfr_init_set_str(r54956, "0.5", 10, MPFR_RNDN); mpfr_init(r54957); mpfr_init(r54958); mpfr_init(r54959); mpfr_init(r54960); mpfr_init_set_str(r54961, "2.0", 10, MPFR_RNDN); mpfr_init(r54962); mpfr_init(r54963); mpfr_init(r54964); mpfr_init(r54965); mpfr_init(r54966); mpfr_init(r54967); mpfr_init(r54968); mpfr_init(r54969); } double f_fm(double x, double y, double z, double t) { mpfr_set_d(r54955, x, MPFR_RNDN); ; mpfr_mul(r54957, r54955, r54956, MPFR_RNDN); mpfr_set_d(r54958, y, MPFR_RNDN); mpfr_sub(r54959, r54957, r54958, MPFR_RNDN); mpfr_set_d(r54960, z, MPFR_RNDN); ; mpfr_mul(r54962, r54960, r54961, MPFR_RNDN); mpfr_sqrt(r54963, r54962, MPFR_RNDN); mpfr_mul(r54964, r54959, r54963, MPFR_RNDN); mpfr_set_d(r54965, t, MPFR_RNDN); mpfr_mul(r54966, r54965, r54965, MPFR_RNDN); mpfr_div(r54967, r54966, r54961, MPFR_RNDN); mpfr_exp(r54968, r54967, MPFR_RNDN); mpfr_mul(r54969, r54964, r54968, MPFR_RNDN); return mpfr_get_d(r54969, MPFR_RNDN); } static mpfr_t r54970, r54971, r54972, r54973, r54974, r54975, r54976, r54977, r54978, r54979, r54980, r54981, r54982, r54983, r54984; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r54970); mpfr_init_set_str(r54971, "0.5", 10, MPFR_RNDN); mpfr_init(r54972); mpfr_init(r54973); mpfr_init(r54974); mpfr_init(r54975); mpfr_init_set_str(r54976, "2.0", 10, MPFR_RNDN); mpfr_init(r54977); mpfr_init(r54978); mpfr_init(r54979); mpfr_init(r54980); mpfr_init(r54981); mpfr_init(r54982); mpfr_init(r54983); mpfr_init(r54984); } double f_dm(double x, double y, double z, double t) { mpfr_set_d(r54970, x, MPFR_RNDN); ; mpfr_mul(r54972, r54970, r54971, MPFR_RNDN); mpfr_set_d(r54973, y, MPFR_RNDN); mpfr_sub(r54974, r54972, r54973, MPFR_RNDN); mpfr_set_d(r54975, z, MPFR_RNDN); ; mpfr_mul(r54977, r54975, r54976, MPFR_RNDN); mpfr_sqrt(r54978, r54977, MPFR_RNDN); mpfr_mul(r54979, r54974, r54978, MPFR_RNDN); mpfr_set_d(r54980, t, MPFR_RNDN); mpfr_mul(r54981, r54980, r54980, MPFR_RNDN); mpfr_div(r54982, r54981, r54976, MPFR_RNDN); mpfr_exp(r54983, r54982, MPFR_RNDN); mpfr_mul(r54984, r54979, r54983, MPFR_RNDN); return mpfr_get_d(r54984, MPFR_RNDN); }