#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 r24259 = 1; float r24260 = atan2(1.0, 0.0); float r24261 = sqrt(r24260); float r24262 = r24259 / r24261; float r24263 = 2; float r24264 = x; float r24265 = fabs(r24264); float r24266 = r24263 * r24265; float r24267 = 3; float r24268 = r24263 / r24267; float r24269 = r24265 * r24265; float r24270 = r24269 * r24265; float r24271 = r24268 * r24270; float r24272 = r24266 + r24271; float r24273 = 5; float r24274 = r24259 / r24273; float r24275 = r24270 * r24265; float r24276 = r24275 * r24265; float r24277 = r24274 * r24276; float r24278 = r24272 + r24277; float r24279 = 21; float r24280 = r24259 / r24279; float r24281 = r24276 * r24265; float r24282 = r24281 * r24265; float r24283 = r24280 * r24282; float r24284 = r24278 + r24283; float r24285 = r24262 * r24284; float r24286 = fabs(r24285); return r24286; } double f_id(double x) { double r24287 = 1; double r24288 = atan2(1.0, 0.0); double r24289 = sqrt(r24288); double r24290 = r24287 / r24289; double r24291 = 2; double r24292 = x; double r24293 = fabs(r24292); double r24294 = r24291 * r24293; double r24295 = 3; double r24296 = r24291 / r24295; double r24297 = r24293 * r24293; double r24298 = r24297 * r24293; double r24299 = r24296 * r24298; double r24300 = r24294 + r24299; double r24301 = 5; double r24302 = r24287 / r24301; double r24303 = r24298 * r24293; double r24304 = r24303 * r24293; double r24305 = r24302 * r24304; double r24306 = r24300 + r24305; double r24307 = 21; double r24308 = r24287 / r24307; double r24309 = r24304 * r24293; double r24310 = r24309 * r24293; double r24311 = r24308 * r24310; double r24312 = r24306 + r24311; double r24313 = r24290 * r24312; double r24314 = fabs(r24313); return r24314; } double f_of(float x) { float r24315 = 1; float r24316 = atan2(1.0, 0.0); float r24317 = sqrt(r24316); float r24318 = r24315 / r24317; float r24319 = 5; float r24320 = r24315 / r24319; float r24321 = x; float r24322 = fabs(r24321); float r24323 = r24320 * r24322; float r24324 = r24322 * r24322; float r24325 = r24323 * r24324; float r24326 = 2; float r24327 = 3; float r24328 = r24326 / r24327; float r24329 = r24322 * r24328; float r24330 = r24326 * r24322; float r24331 = fma(r24329, r24324, r24330); float r24332 = fma(r24325, r24324, r24331); float r24333 = pow(r24322, r24327); float r24334 = r24333 * r24333; float r24335 = 21; float r24336 = r24322 / r24335; float r24337 = pow(r24336, r24315); float r24338 = r24337 * r24318; float r24339 = r24334 * r24338; float r24340 = fma(r24318, r24332, r24339); float r24341 = fabs(r24340); return r24341; } double f_od(double x) { double r24342 = 1; double r24343 = atan2(1.0, 0.0); double r24344 = sqrt(r24343); double r24345 = r24342 / r24344; double r24346 = 5; double r24347 = r24342 / r24346; double r24348 = x; double r24349 = fabs(r24348); double r24350 = r24347 * r24349; double r24351 = r24349 * r24349; double r24352 = r24350 * r24351; double r24353 = 2; double r24354 = 3; double r24355 = r24353 / r24354; double r24356 = r24349 * r24355; double r24357 = r24353 * r24349; double r24358 = fma(r24356, r24351, r24357); double r24359 = fma(r24352, r24351, r24358); double r24360 = pow(r24349, r24354); double r24361 = r24360 * r24360; double r24362 = 21; double r24363 = r24349 / r24362; double r24364 = pow(r24363, r24342); double r24365 = r24364 * r24345; double r24366 = r24361 * r24365; double r24367 = fma(r24345, r24359, r24366); double r24368 = fabs(r24367); return r24368; } 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 r24369, r24370, r24371, r24372, r24373, r24374, r24375, r24376, r24377, r24378, r24379, r24380, r24381, r24382, r24383, r24384, r24385, r24386, r24387, r24388, r24389, r24390, r24391, r24392, r24393, r24394, r24395, r24396; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r24369, "1", 10, MPFR_RNDN); mpfr_init(r24370); mpfr_init(r24371); mpfr_init(r24372); mpfr_init_set_str(r24373, "2", 10, MPFR_RNDN); mpfr_init(r24374); mpfr_init(r24375); mpfr_init(r24376); mpfr_init_set_str(r24377, "3", 10, MPFR_RNDN); mpfr_init(r24378); mpfr_init(r24379); mpfr_init(r24380); mpfr_init(r24381); mpfr_init(r24382); mpfr_init_set_str(r24383, "5", 10, MPFR_RNDN); mpfr_init(r24384); mpfr_init(r24385); mpfr_init(r24386); mpfr_init(r24387); mpfr_init(r24388); mpfr_init_set_str(r24389, "21", 10, MPFR_RNDN); mpfr_init(r24390); mpfr_init(r24391); mpfr_init(r24392); mpfr_init(r24393); mpfr_init(r24394); mpfr_init(r24395); mpfr_init(r24396); } double f_im(double x) { ; mpfr_const_pi(r24370, MPFR_RNDN); mpfr_sqrt(r24371, r24370, MPFR_RNDN); mpfr_div(r24372, r24369, r24371, MPFR_RNDN); ; mpfr_set_d(r24374, x, MPFR_RNDN); mpfr_abs(r24375, r24374, MPFR_RNDN); mpfr_mul(r24376, r24373, r24375, MPFR_RNDN); ; mpfr_div(r24378, r24373, r24377, MPFR_RNDN); mpfr_mul(r24379, r24375, r24375, MPFR_RNDN); mpfr_mul(r24380, r24379, r24375, MPFR_RNDN); mpfr_mul(r24381, r24378, r24380, MPFR_RNDN); mpfr_add(r24382, r24376, r24381, MPFR_RNDN); ; mpfr_div(r24384, r24369, r24383, MPFR_RNDN); mpfr_mul(r24385, r24380, r24375, MPFR_RNDN); mpfr_mul(r24386, r24385, r24375, MPFR_RNDN); mpfr_mul(r24387, r24384, r24386, MPFR_RNDN); mpfr_add(r24388, r24382, r24387, MPFR_RNDN); ; mpfr_div(r24390, r24369, r24389, MPFR_RNDN); mpfr_mul(r24391, r24386, r24375, MPFR_RNDN); mpfr_mul(r24392, r24391, r24375, MPFR_RNDN); mpfr_mul(r24393, r24390, r24392, MPFR_RNDN); mpfr_add(r24394, r24388, r24393, MPFR_RNDN); mpfr_mul(r24395, r24372, r24394, MPFR_RNDN); mpfr_abs(r24396, r24395, MPFR_RNDN); return mpfr_get_d(r24396, MPFR_RNDN); } static mpfr_t r24397, r24398, r24399, r24400, r24401, r24402, r24403, r24404, r24405, r24406, r24407, r24408, r24409, r24410, r24411, r24412, r24413, r24414, r24415, r24416, r24417, r24418, r24419, r24420, r24421, r24422, r24423; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24397, "1", 10, MPFR_RNDN); mpfr_init(r24398); mpfr_init(r24399); mpfr_init(r24400); mpfr_init_set_str(r24401, "5", 10, MPFR_RNDN); mpfr_init(r24402); mpfr_init(r24403); mpfr_init(r24404); mpfr_init(r24405); mpfr_init(r24406); mpfr_init(r24407); mpfr_init_set_str(r24408, "2", 10, MPFR_RNDN); mpfr_init_set_str(r24409, "3", 10, MPFR_RNDN); mpfr_init(r24410); mpfr_init(r24411); mpfr_init(r24412); mpfr_init(r24413); mpfr_init(r24414); mpfr_init(r24415); mpfr_init(r24416); mpfr_init_set_str(r24417, "21", 10, MPFR_RNDN); mpfr_init(r24418); mpfr_init(r24419); mpfr_init(r24420); mpfr_init(r24421); mpfr_init(r24422); mpfr_init(r24423); } double f_fm(double x) { ; mpfr_const_pi(r24398, MPFR_RNDN); mpfr_sqrt(r24399, r24398, MPFR_RNDN); mpfr_div(r24400, r24397, r24399, MPFR_RNDN); ; mpfr_div(r24402, r24397, r24401, MPFR_RNDN); mpfr_set_d(r24403, x, MPFR_RNDN); mpfr_abs(r24404, r24403, MPFR_RNDN); mpfr_mul(r24405, r24402, r24404, MPFR_RNDN); mpfr_mul(r24406, r24404, r24404, MPFR_RNDN); mpfr_mul(r24407, r24405, r24406, MPFR_RNDN); ; ; mpfr_div(r24410, r24408, r24409, MPFR_RNDN); mpfr_mul(r24411, r24404, r24410, MPFR_RNDN); mpfr_mul(r24412, r24408, r24404, MPFR_RNDN); mpfr_fma(r24413, r24411, r24406, r24412, MPFR_RNDN); mpfr_fma(r24414, r24407, r24406, r24413, MPFR_RNDN); mpfr_pow(r24415, r24404, r24409, MPFR_RNDN); mpfr_mul(r24416, r24415, r24415, MPFR_RNDN); ; mpfr_div(r24418, r24404, r24417, MPFR_RNDN); mpfr_pow(r24419, r24418, r24397, MPFR_RNDN); mpfr_mul(r24420, r24419, r24400, MPFR_RNDN); mpfr_mul(r24421, r24416, r24420, MPFR_RNDN); mpfr_fma(r24422, r24400, r24414, r24421, MPFR_RNDN); mpfr_abs(r24423, r24422, MPFR_RNDN); return mpfr_get_d(r24423, MPFR_RNDN); } static mpfr_t r24424, r24425, r24426, r24427, r24428, r24429, r24430, r24431, r24432, r24433, r24434, r24435, r24436, r24437, r24438, r24439, r24440, r24441, r24442, r24443, r24444, r24445, r24446, r24447, r24448, r24449, r24450; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24424, "1", 10, MPFR_RNDN); mpfr_init(r24425); mpfr_init(r24426); mpfr_init(r24427); mpfr_init_set_str(r24428, "5", 10, MPFR_RNDN); mpfr_init(r24429); mpfr_init(r24430); mpfr_init(r24431); mpfr_init(r24432); mpfr_init(r24433); mpfr_init(r24434); mpfr_init_set_str(r24435, "2", 10, MPFR_RNDN); mpfr_init_set_str(r24436, "3", 10, MPFR_RNDN); mpfr_init(r24437); mpfr_init(r24438); mpfr_init(r24439); mpfr_init(r24440); mpfr_init(r24441); mpfr_init(r24442); mpfr_init(r24443); mpfr_init_set_str(r24444, "21", 10, MPFR_RNDN); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init(r24448); mpfr_init(r24449); mpfr_init(r24450); } double f_dm(double x) { ; mpfr_const_pi(r24425, MPFR_RNDN); mpfr_sqrt(r24426, r24425, MPFR_RNDN); mpfr_div(r24427, r24424, r24426, MPFR_RNDN); ; mpfr_div(r24429, r24424, r24428, MPFR_RNDN); mpfr_set_d(r24430, x, MPFR_RNDN); mpfr_abs(r24431, r24430, MPFR_RNDN); mpfr_mul(r24432, r24429, r24431, MPFR_RNDN); mpfr_mul(r24433, r24431, r24431, MPFR_RNDN); mpfr_mul(r24434, r24432, r24433, MPFR_RNDN); ; ; mpfr_div(r24437, r24435, r24436, MPFR_RNDN); mpfr_mul(r24438, r24431, r24437, MPFR_RNDN); mpfr_mul(r24439, r24435, r24431, MPFR_RNDN); mpfr_fma(r24440, r24438, r24433, r24439, MPFR_RNDN); mpfr_fma(r24441, r24434, r24433, r24440, MPFR_RNDN); mpfr_pow(r24442, r24431, r24436, MPFR_RNDN); mpfr_mul(r24443, r24442, r24442, MPFR_RNDN); ; mpfr_div(r24445, r24431, r24444, MPFR_RNDN); mpfr_pow(r24446, r24445, r24424, MPFR_RNDN); mpfr_mul(r24447, r24446, r24427, MPFR_RNDN); mpfr_mul(r24448, r24443, r24447, MPFR_RNDN); mpfr_fma(r24449, r24427, r24441, r24448, MPFR_RNDN); mpfr_abs(r24450, r24449, MPFR_RNDN); return mpfr_get_d(r24450, MPFR_RNDN); }