#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 r56099 = x; float r56100 = 0.5; float r56101 = r56099 * r56100; float r56102 = y; float r56103 = r56101 - r56102; float r56104 = z; float r56105 = 2.0; float r56106 = r56104 * r56105; float r56107 = sqrt(r56106); float r56108 = r56103 * r56107; float r56109 = t; float r56110 = r56109 * r56109; float r56111 = r56110 / r56105; float r56112 = exp(r56111); float r56113 = r56108 * r56112; return r56113; } double f_id(double x, double y, double z, double t) { double r56114 = x; double r56115 = 0.5; double r56116 = r56114 * r56115; double r56117 = y; double r56118 = r56116 - r56117; double r56119 = z; double r56120 = 2.0; double r56121 = r56119 * r56120; double r56122 = sqrt(r56121); double r56123 = r56118 * r56122; double r56124 = t; double r56125 = r56124 * r56124; double r56126 = r56125 / r56120; double r56127 = exp(r56126); double r56128 = r56123 * r56127; return r56128; } double f_of(float x, float y, float z, float t) { float r56129 = x; float r56130 = 0.5; float r56131 = r56129 * r56130; float r56132 = y; float r56133 = r56131 - r56132; float r56134 = z; float r56135 = 2.0; float r56136 = r56134 * r56135; float r56137 = sqrt(r56136); float r56138 = r56133 * r56137; float r56139 = t; float r56140 = r56139 * r56139; float r56141 = r56140 / r56135; float r56142 = exp(r56141); float r56143 = r56138 * r56142; return r56143; } double f_od(double x, double y, double z, double t) { double r56144 = x; double r56145 = 0.5; double r56146 = r56144 * r56145; double r56147 = y; double r56148 = r56146 - r56147; double r56149 = z; double r56150 = 2.0; double r56151 = r56149 * r56150; double r56152 = sqrt(r56151); double r56153 = r56148 * r56152; double r56154 = t; double r56155 = r56154 * r56154; double r56156 = r56155 / r56150; double r56157 = exp(r56156); double r56158 = r56153 * r56157; return r56158; } 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 r56159, r56160, r56161, r56162, r56163, r56164, r56165, r56166, r56167, r56168, r56169, r56170, r56171, r56172, r56173; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r56159); mpfr_init_set_str(r56160, "0.5", 10, MPFR_RNDN); mpfr_init(r56161); mpfr_init(r56162); mpfr_init(r56163); mpfr_init(r56164); mpfr_init_set_str(r56165, "2.0", 10, MPFR_RNDN); mpfr_init(r56166); mpfr_init(r56167); mpfr_init(r56168); mpfr_init(r56169); mpfr_init(r56170); mpfr_init(r56171); mpfr_init(r56172); mpfr_init(r56173); } double f_im(double x, double y, double z, double t) { mpfr_set_d(r56159, x, MPFR_RNDN); ; mpfr_mul(r56161, r56159, r56160, MPFR_RNDN); mpfr_set_d(r56162, y, MPFR_RNDN); mpfr_sub(r56163, r56161, r56162, MPFR_RNDN); mpfr_set_d(r56164, z, MPFR_RNDN); ; mpfr_mul(r56166, r56164, r56165, MPFR_RNDN); mpfr_sqrt(r56167, r56166, MPFR_RNDN); mpfr_mul(r56168, r56163, r56167, MPFR_RNDN); mpfr_set_d(r56169, t, MPFR_RNDN); mpfr_mul(r56170, r56169, r56169, MPFR_RNDN); mpfr_div(r56171, r56170, r56165, MPFR_RNDN); mpfr_exp(r56172, r56171, MPFR_RNDN); mpfr_mul(r56173, r56168, r56172, MPFR_RNDN); return mpfr_get_d(r56173, MPFR_RNDN); } static mpfr_t r56174, r56175, r56176, r56177, r56178, r56179, r56180, r56181, r56182, r56183, r56184, r56185, r56186, r56187, r56188; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r56174); mpfr_init_set_str(r56175, "0.5", 10, MPFR_RNDN); mpfr_init(r56176); mpfr_init(r56177); mpfr_init(r56178); mpfr_init(r56179); mpfr_init_set_str(r56180, "2.0", 10, MPFR_RNDN); mpfr_init(r56181); mpfr_init(r56182); mpfr_init(r56183); mpfr_init(r56184); mpfr_init(r56185); mpfr_init(r56186); mpfr_init(r56187); mpfr_init(r56188); } double f_fm(double x, double y, double z, double t) { mpfr_set_d(r56174, x, MPFR_RNDN); ; mpfr_mul(r56176, r56174, r56175, MPFR_RNDN); mpfr_set_d(r56177, y, MPFR_RNDN); mpfr_sub(r56178, r56176, r56177, MPFR_RNDN); mpfr_set_d(r56179, z, MPFR_RNDN); ; mpfr_mul(r56181, r56179, r56180, MPFR_RNDN); mpfr_sqrt(r56182, r56181, MPFR_RNDN); mpfr_mul(r56183, r56178, r56182, MPFR_RNDN); mpfr_set_d(r56184, t, MPFR_RNDN); mpfr_mul(r56185, r56184, r56184, MPFR_RNDN); mpfr_div(r56186, r56185, r56180, MPFR_RNDN); mpfr_exp(r56187, r56186, MPFR_RNDN); mpfr_mul(r56188, r56183, r56187, MPFR_RNDN); return mpfr_get_d(r56188, MPFR_RNDN); } static mpfr_t r56189, r56190, r56191, r56192, r56193, r56194, r56195, r56196, r56197, r56198, r56199, r56200, r56201, r56202, r56203; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r56189); mpfr_init_set_str(r56190, "0.5", 10, MPFR_RNDN); mpfr_init(r56191); mpfr_init(r56192); mpfr_init(r56193); mpfr_init(r56194); mpfr_init_set_str(r56195, "2.0", 10, MPFR_RNDN); mpfr_init(r56196); mpfr_init(r56197); mpfr_init(r56198); mpfr_init(r56199); mpfr_init(r56200); mpfr_init(r56201); mpfr_init(r56202); mpfr_init(r56203); } double f_dm(double x, double y, double z, double t) { mpfr_set_d(r56189, x, MPFR_RNDN); ; mpfr_mul(r56191, r56189, r56190, MPFR_RNDN); mpfr_set_d(r56192, y, MPFR_RNDN); mpfr_sub(r56193, r56191, r56192, MPFR_RNDN); mpfr_set_d(r56194, z, MPFR_RNDN); ; mpfr_mul(r56196, r56194, r56195, MPFR_RNDN); mpfr_sqrt(r56197, r56196, MPFR_RNDN); mpfr_mul(r56198, r56193, r56197, MPFR_RNDN); mpfr_set_d(r56199, t, MPFR_RNDN); mpfr_mul(r56200, r56199, r56199, MPFR_RNDN); mpfr_div(r56201, r56200, r56195, MPFR_RNDN); mpfr_exp(r56202, r56201, MPFR_RNDN); mpfr_mul(r56203, r56198, r56202, MPFR_RNDN); return mpfr_get_d(r56203, MPFR_RNDN); }