#include #include #include #include #include char *name = "Jmat.Real.erf"; double f_if(float x) { float r18177 = 1.0f; float r18178 = 0.32759109139442444f; float r18179 = x; float r18180 = fabs(r18179); float r18181 = r18178 * r18180; float r18182 = r18177 + r18181; float r18183 = r18177 / r18182; float r18184 = 0.2548295855522156f; float r18185 = -0.2844967246055603f; float r18186 = 1.421413779258728f; float r18187 = -1.453152060508728f; float r18188 = 1.0614054203033447f; float r18189 = r18183 * r18188; float r18190 = r18187 + r18189; float r18191 = r18183 * r18190; float r18192 = r18186 + r18191; float r18193 = r18183 * r18192; float r18194 = r18185 + r18193; float r18195 = r18183 * r18194; float r18196 = r18184 + r18195; float r18197 = r18183 * r18196; float r18198 = r18180 * r18180; float r18199 = -r18198; float r18200 = exp(r18199); float r18201 = r18197 * r18200; float r18202 = r18177 - r18201; return r18202; } double f_id(double x) { double r18203 = 1.0; double r18204 = 0.32759109139442444; double r18205 = x; double r18206 = fabs(r18205); double r18207 = r18204 * r18206; double r18208 = r18203 + r18207; double r18209 = r18203 / r18208; double r18210 = 0.2548295855522156; double r18211 = -0.2844967246055603; double r18212 = 1.421413779258728; double r18213 = -1.453152060508728; double r18214 = 1.0614054203033447; double r18215 = r18209 * r18214; double r18216 = r18213 + r18215; double r18217 = r18209 * r18216; double r18218 = r18212 + r18217; double r18219 = r18209 * r18218; double r18220 = r18211 + r18219; double r18221 = r18209 * r18220; double r18222 = r18210 + r18221; double r18223 = r18209 * r18222; double r18224 = r18206 * r18206; double r18225 = -r18224; double r18226 = exp(r18225); double r18227 = r18223 * r18226; double r18228 = r18203 - r18227; return r18228; } double f_of(float x) { float r18229 = 1.0f; float r18230 = -1.453152060508728f; float r18231 = 1.0614054203033447f; float r18232 = 0.32759109139442444f; float r18233 = x; float r18234 = fabs(r18233); float r18235 = fma(r18232, r18234, r18229); float r18236 = r18231 / r18235; float r18237 = r18230 + r18236; float r18238 = cbrt(r18237); float r18239 = cbrt(r18235); float r18240 = r18239 * r18239; float r18241 = r18238 / r18240; float r18242 = r18241 * (r18241 * r18241); float r18243 = 1.421413779258728f; float r18244 = r18243 / r18235; float r18245 = -0.2844967246055603f; float r18246 = r18244 + r18245; float r18247 = r18242 + r18246; float r18248 = r18229 / r18235; float r18249 = r18248 / r18235; float r18250 = 0.2548295855522156f; float r18251 = r18250 / r18235; float r18252 = fma(r18247, r18249, r18251); float r18253 = r18234 * r18234; float r18254 = exp(r18253); float r18255 = r18252 / r18254; float r18256 = r18229 - r18255; return r18256; } double f_od(double x) { double r18257 = 1.0; double r18258 = -1.453152060508728; double r18259 = 1.0614054203033447; double r18260 = 0.32759109139442444; double r18261 = x; double r18262 = fabs(r18261); double r18263 = fma(r18260, r18262, r18257); double r18264 = r18259 / r18263; double r18265 = r18258 + r18264; double r18266 = cbrt(r18265); double r18267 = cbrt(r18263); double r18268 = r18267 * r18267; double r18269 = r18266 / r18268; double r18270 = r18269 * (r18269 * r18269); double r18271 = 1.421413779258728; double r18272 = r18271 / r18263; double r18273 = -0.2844967246055603; double r18274 = r18272 + r18273; double r18275 = r18270 + r18274; double r18276 = r18257 / r18263; double r18277 = r18276 / r18263; double r18278 = 0.2548295855522156; double r18279 = r18278 / r18263; double r18280 = fma(r18275, r18277, r18279); double r18281 = r18262 * r18262; double r18282 = exp(r18281); double r18283 = r18280 / r18282; double r18284 = r18257 - r18283; return r18284; } 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 r18285, r18286, r18287, r18288, r18289, r18290, r18291, r18292, r18293, r18294, r18295, r18296, r18297, r18298, r18299, r18300, r18301, r18302, r18303, r18304, r18305, r18306, r18307, r18308, r18309, r18310; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18285, "1", 10, MPFR_RNDN); mpfr_init_set_str(r18286, "0.3275911", 10, MPFR_RNDN); mpfr_init(r18287); mpfr_init(r18288); mpfr_init(r18289); mpfr_init(r18290); mpfr_init(r18291); mpfr_init_set_str(r18292, "0.254829592", 10, MPFR_RNDN); mpfr_init_set_str(r18293, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r18294, "1.421413741", 10, MPFR_RNDN); mpfr_init_set_str(r18295, "-1.453152027", 10, MPFR_RNDN); mpfr_init_set_str(r18296, "1.061405429", 10, MPFR_RNDN); mpfr_init(r18297); mpfr_init(r18298); mpfr_init(r18299); mpfr_init(r18300); mpfr_init(r18301); mpfr_init(r18302); mpfr_init(r18303); mpfr_init(r18304); mpfr_init(r18305); mpfr_init(r18306); mpfr_init(r18307); mpfr_init(r18308); mpfr_init(r18309); mpfr_init(r18310); } double f_im(double x) { ; ; mpfr_set_d(r18287, x, MPFR_RNDN); mpfr_abs(r18288, r18287, MPFR_RNDN); mpfr_mul(r18289, r18286, r18288, MPFR_RNDN); mpfr_add(r18290, r18285, r18289, MPFR_RNDN); mpfr_div(r18291, r18285, r18290, MPFR_RNDN); ; ; ; ; ; mpfr_mul(r18297, r18291, r18296, MPFR_RNDN); mpfr_add(r18298, r18295, r18297, MPFR_RNDN); mpfr_mul(r18299, r18291, r18298, MPFR_RNDN); mpfr_add(r18300, r18294, r18299, MPFR_RNDN); mpfr_mul(r18301, r18291, r18300, MPFR_RNDN); mpfr_add(r18302, r18293, r18301, MPFR_RNDN); mpfr_mul(r18303, r18291, r18302, MPFR_RNDN); mpfr_add(r18304, r18292, r18303, MPFR_RNDN); mpfr_mul(r18305, r18291, r18304, MPFR_RNDN); mpfr_mul(r18306, r18288, r18288, MPFR_RNDN); mpfr_neg(r18307, r18306, MPFR_RNDN); mpfr_exp(r18308, r18307, MPFR_RNDN); mpfr_mul(r18309, r18305, r18308, MPFR_RNDN); mpfr_sub(r18310, r18285, r18309, MPFR_RNDN); return mpfr_get_d(r18310, MPFR_RNDN); } static mpfr_t r18311, r18312, r18313, r18314, r18315, r18316, r18317, r18318, r18319, r18320, r18321, r18322, r18323, r18324, r18325, r18326, r18327, r18328, r18329, r18330, r18331, r18332, r18333, r18334, r18335, r18336, r18337, r18338; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18311, "1", 10, MPFR_RNDN); mpfr_init_set_str(r18312, "-1.453152027", 10, MPFR_RNDN); mpfr_init_set_str(r18313, "1.061405429", 10, MPFR_RNDN); mpfr_init_set_str(r18314, "0.3275911", 10, MPFR_RNDN); mpfr_init(r18315); mpfr_init(r18316); mpfr_init(r18317); mpfr_init(r18318); mpfr_init(r18319); mpfr_init(r18320); mpfr_init(r18321); mpfr_init(r18322); mpfr_init(r18323); mpfr_init(r18324); mpfr_init_set_str(r18325, "1.421413741", 10, MPFR_RNDN); mpfr_init(r18326); mpfr_init_set_str(r18327, "-0.284496736", 10, MPFR_RNDN); mpfr_init(r18328); mpfr_init(r18329); mpfr_init(r18330); mpfr_init(r18331); mpfr_init_set_str(r18332, "0.254829592", 10, MPFR_RNDN); mpfr_init(r18333); mpfr_init(r18334); mpfr_init(r18335); mpfr_init(r18336); mpfr_init(r18337); mpfr_init(r18338); } double f_fm(double x) { ; ; ; ; mpfr_set_d(r18315, x, MPFR_RNDN); mpfr_abs(r18316, r18315, MPFR_RNDN); mpfr_fma(r18317, r18314, r18316, r18311, MPFR_RNDN); mpfr_div(r18318, r18313, r18317, MPFR_RNDN); mpfr_add(r18319, r18312, r18318, MPFR_RNDN); mpfr_cbrt(r18320, r18319, MPFR_RNDN); mpfr_cbrt(r18321, r18317, MPFR_RNDN); mpfr_mul(r18322, r18321, r18321, MPFR_RNDN); mpfr_div(r18323, r18320, r18322, MPFR_RNDN); mpfr_mul(r18324, r18323, r18323, MPFR_RNDN); mpfr_mul(r18324, r18324, r18323, MPFR_RNDN); ; mpfr_div(r18326, r18325, r18317, MPFR_RNDN); ; mpfr_add(r18328, r18326, r18327, MPFR_RNDN); mpfr_add(r18329, r18324, r18328, MPFR_RNDN); mpfr_div(r18330, r18311, r18317, MPFR_RNDN); mpfr_div(r18331, r18330, r18317, MPFR_RNDN); ; mpfr_div(r18333, r18332, r18317, MPFR_RNDN); mpfr_fma(r18334, r18329, r18331, r18333, MPFR_RNDN); mpfr_mul(r18335, r18316, r18316, MPFR_RNDN); mpfr_exp(r18336, r18335, MPFR_RNDN); mpfr_div(r18337, r18334, r18336, MPFR_RNDN); mpfr_sub(r18338, r18311, r18337, MPFR_RNDN); return mpfr_get_d(r18338, MPFR_RNDN); } static mpfr_t r18339, r18340, r18341, r18342, r18343, r18344, r18345, r18346, r18347, r18348, r18349, r18350, r18351, r18352, r18353, r18354, r18355, r18356, r18357, r18358, r18359, r18360, r18361, r18362, r18363, r18364, r18365, r18366; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18339, "1", 10, MPFR_RNDN); mpfr_init_set_str(r18340, "-1.453152027", 10, MPFR_RNDN); mpfr_init_set_str(r18341, "1.061405429", 10, MPFR_RNDN); mpfr_init_set_str(r18342, "0.3275911", 10, MPFR_RNDN); mpfr_init(r18343); mpfr_init(r18344); mpfr_init(r18345); mpfr_init(r18346); mpfr_init(r18347); mpfr_init(r18348); mpfr_init(r18349); mpfr_init(r18350); mpfr_init(r18351); mpfr_init(r18352); mpfr_init_set_str(r18353, "1.421413741", 10, MPFR_RNDN); mpfr_init(r18354); mpfr_init_set_str(r18355, "-0.284496736", 10, MPFR_RNDN); mpfr_init(r18356); mpfr_init(r18357); mpfr_init(r18358); mpfr_init(r18359); mpfr_init_set_str(r18360, "0.254829592", 10, MPFR_RNDN); mpfr_init(r18361); mpfr_init(r18362); mpfr_init(r18363); mpfr_init(r18364); mpfr_init(r18365); mpfr_init(r18366); } double f_dm(double x) { ; ; ; ; mpfr_set_d(r18343, x, MPFR_RNDN); mpfr_abs(r18344, r18343, MPFR_RNDN); mpfr_fma(r18345, r18342, r18344, r18339, MPFR_RNDN); mpfr_div(r18346, r18341, r18345, MPFR_RNDN); mpfr_add(r18347, r18340, r18346, MPFR_RNDN); mpfr_cbrt(r18348, r18347, MPFR_RNDN); mpfr_cbrt(r18349, r18345, MPFR_RNDN); mpfr_mul(r18350, r18349, r18349, MPFR_RNDN); mpfr_div(r18351, r18348, r18350, MPFR_RNDN); mpfr_mul(r18352, r18351, r18351, MPFR_RNDN); mpfr_mul(r18352, r18352, r18351, MPFR_RNDN); ; mpfr_div(r18354, r18353, r18345, MPFR_RNDN); ; mpfr_add(r18356, r18354, r18355, MPFR_RNDN); mpfr_add(r18357, r18352, r18356, MPFR_RNDN); mpfr_div(r18358, r18339, r18345, MPFR_RNDN); mpfr_div(r18359, r18358, r18345, MPFR_RNDN); ; mpfr_div(r18361, r18360, r18345, MPFR_RNDN); mpfr_fma(r18362, r18357, r18359, r18361, MPFR_RNDN); mpfr_mul(r18363, r18344, r18344, MPFR_RNDN); mpfr_exp(r18364, r18363, MPFR_RNDN); mpfr_div(r18365, r18362, r18364, MPFR_RNDN); mpfr_sub(r18366, r18339, r18365, MPFR_RNDN); return mpfr_get_d(r18366, MPFR_RNDN); }