#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 r27324 = 1; float r27325 = atan2(1.0, 0.0); float r27326 = sqrt(r27325); float r27327 = r27324 / r27326; float r27328 = 2; float r27329 = x; float r27330 = fabs(r27329); float r27331 = r27328 * r27330; float r27332 = 3; float r27333 = r27328 / r27332; float r27334 = r27330 * r27330; float r27335 = r27334 * r27330; float r27336 = r27333 * r27335; float r27337 = r27331 + r27336; float r27338 = 5; float r27339 = r27324 / r27338; float r27340 = r27335 * r27330; float r27341 = r27340 * r27330; float r27342 = r27339 * r27341; float r27343 = r27337 + r27342; float r27344 = 21; float r27345 = r27324 / r27344; float r27346 = r27341 * r27330; float r27347 = r27346 * r27330; float r27348 = r27345 * r27347; float r27349 = r27343 + r27348; float r27350 = r27327 * r27349; float r27351 = fabs(r27350); return r27351; } double f_id(double x) { double r27352 = 1; double r27353 = atan2(1.0, 0.0); double r27354 = sqrt(r27353); double r27355 = r27352 / r27354; double r27356 = 2; double r27357 = x; double r27358 = fabs(r27357); double r27359 = r27356 * r27358; double r27360 = 3; double r27361 = r27356 / r27360; double r27362 = r27358 * r27358; double r27363 = r27362 * r27358; double r27364 = r27361 * r27363; double r27365 = r27359 + r27364; double r27366 = 5; double r27367 = r27352 / r27366; double r27368 = r27363 * r27358; double r27369 = r27368 * r27358; double r27370 = r27367 * r27369; double r27371 = r27365 + r27370; double r27372 = 21; double r27373 = r27352 / r27372; double r27374 = r27369 * r27358; double r27375 = r27374 * r27358; double r27376 = r27373 * r27375; double r27377 = r27371 + r27376; double r27378 = r27355 * r27377; double r27379 = fabs(r27378); return r27379; } double f_of(float x) { float r27380 = 1; float r27381 = atan2(1.0, 0.0); float r27382 = sqrt(r27381); float r27383 = r27380 / r27382; float r27384 = 1/5; float r27385 = x; float r27386 = fabs(r27385); float r27387 = r27384 * r27386; float r27388 = r27386 * r27386; float r27389 = r27387 * r27388; float r27390 = 2; float r27391 = 3; float r27392 = r27390 / r27391; float r27393 = r27392 * r27386; float r27394 = r27390 * r27386; float r27395 = fma(r27393, r27388, r27394); float r27396 = fma(r27389, r27388, r27395); float r27397 = pow(r27386, r27391); float r27398 = sqrt(r27382); float r27399 = r27398 * r27398; float r27400 = r27397 / r27399; float r27401 = 21; float r27402 = r27401 / r27386; float r27403 = r27397 / r27402; float r27404 = r27400 * r27403; float r27405 = fma(r27383, r27396, r27404); float r27406 = fabs(r27405); return r27406; } double f_od(double x) { double r27407 = 1; double r27408 = atan2(1.0, 0.0); double r27409 = sqrt(r27408); double r27410 = r27407 / r27409; double r27411 = 1/5; double r27412 = x; double r27413 = fabs(r27412); double r27414 = r27411 * r27413; double r27415 = r27413 * r27413; double r27416 = r27414 * r27415; double r27417 = 2; double r27418 = 3; double r27419 = r27417 / r27418; double r27420 = r27419 * r27413; double r27421 = r27417 * r27413; double r27422 = fma(r27420, r27415, r27421); double r27423 = fma(r27416, r27415, r27422); double r27424 = pow(r27413, r27418); double r27425 = sqrt(r27409); double r27426 = r27425 * r27425; double r27427 = r27424 / r27426; double r27428 = 21; double r27429 = r27428 / r27413; double r27430 = r27424 / r27429; double r27431 = r27427 * r27430; double r27432 = fma(r27410, r27423, r27431); double r27433 = fabs(r27432); return r27433; } 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 r27434, r27435, r27436, r27437, r27438, r27439, r27440, r27441, r27442, r27443, r27444, r27445, r27446, r27447, r27448, r27449, r27450, r27451, r27452, r27453, r27454, r27455, r27456, r27457, r27458, r27459, r27460, r27461; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27434, "1", 10, MPFR_RNDN); mpfr_init(r27435); mpfr_init(r27436); mpfr_init(r27437); mpfr_init_set_str(r27438, "2", 10, MPFR_RNDN); mpfr_init(r27439); mpfr_init(r27440); mpfr_init(r27441); mpfr_init_set_str(r27442, "3", 10, MPFR_RNDN); mpfr_init(r27443); mpfr_init(r27444); mpfr_init(r27445); mpfr_init(r27446); mpfr_init(r27447); mpfr_init_set_str(r27448, "5", 10, MPFR_RNDN); mpfr_init(r27449); mpfr_init(r27450); mpfr_init(r27451); mpfr_init(r27452); mpfr_init(r27453); mpfr_init_set_str(r27454, "21", 10, MPFR_RNDN); mpfr_init(r27455); mpfr_init(r27456); mpfr_init(r27457); mpfr_init(r27458); mpfr_init(r27459); mpfr_init(r27460); mpfr_init(r27461); } double f_im(double x) { ; mpfr_const_pi(r27435, MPFR_RNDN); mpfr_sqrt(r27436, r27435, MPFR_RNDN); mpfr_div(r27437, r27434, r27436, MPFR_RNDN); ; mpfr_set_d(r27439, x, MPFR_RNDN); mpfr_abs(r27440, r27439, MPFR_RNDN); mpfr_mul(r27441, r27438, r27440, MPFR_RNDN); ; mpfr_div(r27443, r27438, r27442, MPFR_RNDN); mpfr_mul(r27444, r27440, r27440, MPFR_RNDN); mpfr_mul(r27445, r27444, r27440, MPFR_RNDN); mpfr_mul(r27446, r27443, r27445, MPFR_RNDN); mpfr_add(r27447, r27441, r27446, MPFR_RNDN); ; mpfr_div(r27449, r27434, r27448, MPFR_RNDN); mpfr_mul(r27450, r27445, r27440, MPFR_RNDN); mpfr_mul(r27451, r27450, r27440, MPFR_RNDN); mpfr_mul(r27452, r27449, r27451, MPFR_RNDN); mpfr_add(r27453, r27447, r27452, MPFR_RNDN); ; mpfr_div(r27455, r27434, r27454, MPFR_RNDN); mpfr_mul(r27456, r27451, r27440, MPFR_RNDN); mpfr_mul(r27457, r27456, r27440, MPFR_RNDN); mpfr_mul(r27458, r27455, r27457, MPFR_RNDN); mpfr_add(r27459, r27453, r27458, MPFR_RNDN); mpfr_mul(r27460, r27437, r27459, MPFR_RNDN); mpfr_abs(r27461, r27460, MPFR_RNDN); return mpfr_get_d(r27461, MPFR_RNDN); } static mpfr_t r27462, r27463, r27464, r27465, r27466, r27467, r27468, r27469, r27470, r27471, r27472, r27473, r27474, r27475, r27476, r27477, r27478, r27479, r27480, r27481, r27482, r27483, r27484, r27485, r27486, r27487, r27488; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27462, "1", 10, MPFR_RNDN); mpfr_init(r27463); mpfr_init(r27464); mpfr_init(r27465); mpfr_init_set_str(r27466, "1/5", 10, MPFR_RNDN); mpfr_init(r27467); mpfr_init(r27468); mpfr_init(r27469); mpfr_init(r27470); mpfr_init(r27471); mpfr_init_set_str(r27472, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27473, "3", 10, MPFR_RNDN); mpfr_init(r27474); mpfr_init(r27475); mpfr_init(r27476); mpfr_init(r27477); mpfr_init(r27478); mpfr_init(r27479); mpfr_init(r27480); mpfr_init(r27481); mpfr_init(r27482); mpfr_init_set_str(r27483, "21", 10, MPFR_RNDN); mpfr_init(r27484); mpfr_init(r27485); mpfr_init(r27486); mpfr_init(r27487); mpfr_init(r27488); } double f_fm(double x) { ; mpfr_const_pi(r27463, MPFR_RNDN); mpfr_sqrt(r27464, r27463, MPFR_RNDN); mpfr_div(r27465, r27462, r27464, MPFR_RNDN); ; mpfr_set_d(r27467, x, MPFR_RNDN); mpfr_abs(r27468, r27467, MPFR_RNDN); mpfr_mul(r27469, r27466, r27468, MPFR_RNDN); mpfr_mul(r27470, r27468, r27468, MPFR_RNDN); mpfr_mul(r27471, r27469, r27470, MPFR_RNDN); ; ; mpfr_div(r27474, r27472, r27473, MPFR_RNDN); mpfr_mul(r27475, r27474, r27468, MPFR_RNDN); mpfr_mul(r27476, r27472, r27468, MPFR_RNDN); mpfr_fma(r27477, r27475, r27470, r27476, MPFR_RNDN); mpfr_fma(r27478, r27471, r27470, r27477, MPFR_RNDN); mpfr_pow(r27479, r27468, r27473, MPFR_RNDN); mpfr_sqrt(r27480, r27464, MPFR_RNDN); mpfr_mul(r27481, r27480, r27480, MPFR_RNDN); mpfr_div(r27482, r27479, r27481, MPFR_RNDN); ; mpfr_div(r27484, r27483, r27468, MPFR_RNDN); mpfr_div(r27485, r27479, r27484, MPFR_RNDN); mpfr_mul(r27486, r27482, r27485, MPFR_RNDN); mpfr_fma(r27487, r27465, r27478, r27486, MPFR_RNDN); mpfr_abs(r27488, r27487, MPFR_RNDN); return mpfr_get_d(r27488, MPFR_RNDN); } static mpfr_t r27489, r27490, r27491, r27492, r27493, r27494, r27495, r27496, r27497, r27498, r27499, r27500, r27501, r27502, r27503, r27504, r27505, r27506, r27507, r27508, r27509, r27510, r27511, r27512, r27513, r27514, r27515; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27489, "1", 10, MPFR_RNDN); mpfr_init(r27490); mpfr_init(r27491); mpfr_init(r27492); mpfr_init_set_str(r27493, "1/5", 10, MPFR_RNDN); mpfr_init(r27494); mpfr_init(r27495); mpfr_init(r27496); mpfr_init(r27497); mpfr_init(r27498); mpfr_init_set_str(r27499, "2", 10, MPFR_RNDN); mpfr_init_set_str(r27500, "3", 10, MPFR_RNDN); mpfr_init(r27501); mpfr_init(r27502); mpfr_init(r27503); mpfr_init(r27504); mpfr_init(r27505); mpfr_init(r27506); mpfr_init(r27507); mpfr_init(r27508); mpfr_init(r27509); mpfr_init_set_str(r27510, "21", 10, MPFR_RNDN); mpfr_init(r27511); mpfr_init(r27512); mpfr_init(r27513); mpfr_init(r27514); mpfr_init(r27515); } double f_dm(double x) { ; mpfr_const_pi(r27490, MPFR_RNDN); mpfr_sqrt(r27491, r27490, MPFR_RNDN); mpfr_div(r27492, r27489, r27491, MPFR_RNDN); ; mpfr_set_d(r27494, x, MPFR_RNDN); mpfr_abs(r27495, r27494, MPFR_RNDN); mpfr_mul(r27496, r27493, r27495, MPFR_RNDN); mpfr_mul(r27497, r27495, r27495, MPFR_RNDN); mpfr_mul(r27498, r27496, r27497, MPFR_RNDN); ; ; mpfr_div(r27501, r27499, r27500, MPFR_RNDN); mpfr_mul(r27502, r27501, r27495, MPFR_RNDN); mpfr_mul(r27503, r27499, r27495, MPFR_RNDN); mpfr_fma(r27504, r27502, r27497, r27503, MPFR_RNDN); mpfr_fma(r27505, r27498, r27497, r27504, MPFR_RNDN); mpfr_pow(r27506, r27495, r27500, MPFR_RNDN); mpfr_sqrt(r27507, r27491, MPFR_RNDN); mpfr_mul(r27508, r27507, r27507, MPFR_RNDN); mpfr_div(r27509, r27506, r27508, MPFR_RNDN); ; mpfr_div(r27511, r27510, r27495, MPFR_RNDN); mpfr_div(r27512, r27506, r27511, MPFR_RNDN); mpfr_mul(r27513, r27509, r27512, MPFR_RNDN); mpfr_fma(r27514, r27492, r27505, r27513, MPFR_RNDN); mpfr_abs(r27515, r27514, MPFR_RNDN); return mpfr_get_d(r27515, MPFR_RNDN); }