#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 r24419 = 1; float r24420 = atan2(1.0, 0.0); float r24421 = sqrt(r24420); float r24422 = r24419 / r24421; float r24423 = 2; float r24424 = x; float r24425 = fabs(r24424); float r24426 = r24423 * r24425; float r24427 = 3; float r24428 = r24423 / r24427; float r24429 = r24425 * r24425; float r24430 = r24429 * r24425; float r24431 = r24428 * r24430; float r24432 = r24426 + r24431; float r24433 = 5; float r24434 = r24419 / r24433; float r24435 = r24430 * r24425; float r24436 = r24435 * r24425; float r24437 = r24434 * r24436; float r24438 = r24432 + r24437; float r24439 = 21; float r24440 = r24419 / r24439; float r24441 = r24436 * r24425; float r24442 = r24441 * r24425; float r24443 = r24440 * r24442; float r24444 = r24438 + r24443; float r24445 = r24422 * r24444; float r24446 = fabs(r24445); return r24446; } double f_id(double x) { double r24447 = 1; double r24448 = atan2(1.0, 0.0); double r24449 = sqrt(r24448); double r24450 = r24447 / r24449; double r24451 = 2; double r24452 = x; double r24453 = fabs(r24452); double r24454 = r24451 * r24453; double r24455 = 3; double r24456 = r24451 / r24455; double r24457 = r24453 * r24453; double r24458 = r24457 * r24453; double r24459 = r24456 * r24458; double r24460 = r24454 + r24459; double r24461 = 5; double r24462 = r24447 / r24461; double r24463 = r24458 * r24453; double r24464 = r24463 * r24453; double r24465 = r24462 * r24464; double r24466 = r24460 + r24465; double r24467 = 21; double r24468 = r24447 / r24467; double r24469 = r24464 * r24453; double r24470 = r24469 * r24453; double r24471 = r24468 * r24470; double r24472 = r24466 + r24471; double r24473 = r24450 * r24472; double r24474 = fabs(r24473); return r24474; } double f_of(float x) { float r24475 = 1; float r24476 = atan2(1.0, 0.0); float r24477 = sqrt(r24476); float r24478 = r24475 / r24477; float r24479 = x; float r24480 = fabs(r24479); float r24481 = 5; float r24482 = r24480 / r24481; float r24483 = r24480 * r24480; float r24484 = r24483 * r24483; float r24485 = 3; float r24486 = 2; float r24487 = r24485 / r24486; float r24488 = r24480 / r24487; float r24489 = r24486 * r24480; float r24490 = fma(r24483, r24488, r24489); float r24491 = fma(r24482, r24484, r24490); float r24492 = r24480 / r24477; float r24493 = 21; float r24494 = r24492 / r24493; float r24495 = pow(r24480, r24485); float r24496 = r24495 * r24495; float r24497 = r24494 * r24496; float r24498 = fma(r24478, r24491, r24497); float r24499 = fabs(r24498); return r24499; } double f_od(double x) { double r24500 = 1; double r24501 = atan2(1.0, 0.0); double r24502 = sqrt(r24501); double r24503 = r24500 / r24502; double r24504 = x; double r24505 = fabs(r24504); double r24506 = 5; double r24507 = r24505 / r24506; double r24508 = r24505 * r24505; double r24509 = r24508 * r24508; double r24510 = 3; double r24511 = 2; double r24512 = r24510 / r24511; double r24513 = r24505 / r24512; double r24514 = r24511 * r24505; double r24515 = fma(r24508, r24513, r24514); double r24516 = fma(r24507, r24509, r24515); double r24517 = r24505 / r24502; double r24518 = 21; double r24519 = r24517 / r24518; double r24520 = pow(r24505, r24510); double r24521 = r24520 * r24520; double r24522 = r24519 * r24521; double r24523 = fma(r24503, r24516, r24522); double r24524 = fabs(r24523); return r24524; } 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 r24525, r24526, r24527, r24528, r24529, r24530, r24531, r24532, r24533, r24534, r24535, r24536, r24537, r24538, r24539, r24540, r24541, r24542, r24543, r24544, r24545, r24546, r24547, r24548, r24549, r24550, r24551, r24552; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r24525, "1", 10, MPFR_RNDN); mpfr_init(r24526); mpfr_init(r24527); mpfr_init(r24528); mpfr_init_set_str(r24529, "2", 10, MPFR_RNDN); mpfr_init(r24530); mpfr_init(r24531); mpfr_init(r24532); mpfr_init_set_str(r24533, "3", 10, MPFR_RNDN); mpfr_init(r24534); mpfr_init(r24535); mpfr_init(r24536); mpfr_init(r24537); mpfr_init(r24538); mpfr_init_set_str(r24539, "5", 10, MPFR_RNDN); mpfr_init(r24540); mpfr_init(r24541); mpfr_init(r24542); mpfr_init(r24543); mpfr_init(r24544); mpfr_init_set_str(r24545, "21", 10, MPFR_RNDN); mpfr_init(r24546); mpfr_init(r24547); mpfr_init(r24548); mpfr_init(r24549); mpfr_init(r24550); mpfr_init(r24551); mpfr_init(r24552); } double f_im(double x) { ; mpfr_const_pi(r24526, MPFR_RNDN); mpfr_sqrt(r24527, r24526, MPFR_RNDN); mpfr_div(r24528, r24525, r24527, MPFR_RNDN); ; mpfr_set_d(r24530, x, MPFR_RNDN); mpfr_abs(r24531, r24530, MPFR_RNDN); mpfr_mul(r24532, r24529, r24531, MPFR_RNDN); ; mpfr_div(r24534, r24529, r24533, MPFR_RNDN); mpfr_mul(r24535, r24531, r24531, MPFR_RNDN); mpfr_mul(r24536, r24535, r24531, MPFR_RNDN); mpfr_mul(r24537, r24534, r24536, MPFR_RNDN); mpfr_add(r24538, r24532, r24537, MPFR_RNDN); ; mpfr_div(r24540, r24525, r24539, MPFR_RNDN); mpfr_mul(r24541, r24536, r24531, MPFR_RNDN); mpfr_mul(r24542, r24541, r24531, MPFR_RNDN); mpfr_mul(r24543, r24540, r24542, MPFR_RNDN); mpfr_add(r24544, r24538, r24543, MPFR_RNDN); ; mpfr_div(r24546, r24525, r24545, MPFR_RNDN); mpfr_mul(r24547, r24542, r24531, MPFR_RNDN); mpfr_mul(r24548, r24547, r24531, MPFR_RNDN); mpfr_mul(r24549, r24546, r24548, MPFR_RNDN); mpfr_add(r24550, r24544, r24549, MPFR_RNDN); mpfr_mul(r24551, r24528, r24550, MPFR_RNDN); mpfr_abs(r24552, r24551, MPFR_RNDN); return mpfr_get_d(r24552, MPFR_RNDN); } static mpfr_t r24553, r24554, r24555, r24556, r24557, r24558, r24559, r24560, r24561, r24562, r24563, r24564, r24565, r24566, r24567, r24568, r24569, r24570, r24571, r24572, r24573, r24574, r24575, r24576, r24577; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24553, "1", 10, MPFR_RNDN); mpfr_init(r24554); mpfr_init(r24555); mpfr_init(r24556); mpfr_init(r24557); mpfr_init(r24558); mpfr_init_set_str(r24559, "5", 10, MPFR_RNDN); mpfr_init(r24560); mpfr_init(r24561); mpfr_init(r24562); mpfr_init_set_str(r24563, "3", 10, MPFR_RNDN); mpfr_init_set_str(r24564, "2", 10, MPFR_RNDN); mpfr_init(r24565); mpfr_init(r24566); mpfr_init(r24567); mpfr_init(r24568); mpfr_init(r24569); mpfr_init(r24570); mpfr_init_set_str(r24571, "21", 10, MPFR_RNDN); mpfr_init(r24572); mpfr_init(r24573); mpfr_init(r24574); mpfr_init(r24575); mpfr_init(r24576); mpfr_init(r24577); } double f_fm(double x) { ; mpfr_const_pi(r24554, MPFR_RNDN); mpfr_sqrt(r24555, r24554, MPFR_RNDN); mpfr_div(r24556, r24553, r24555, MPFR_RNDN); mpfr_set_d(r24557, x, MPFR_RNDN); mpfr_abs(r24558, r24557, MPFR_RNDN); ; mpfr_div(r24560, r24558, r24559, MPFR_RNDN); mpfr_mul(r24561, r24558, r24558, MPFR_RNDN); mpfr_mul(r24562, r24561, r24561, MPFR_RNDN); ; ; mpfr_div(r24565, r24563, r24564, MPFR_RNDN); mpfr_div(r24566, r24558, r24565, MPFR_RNDN); mpfr_mul(r24567, r24564, r24558, MPFR_RNDN); mpfr_fma(r24568, r24561, r24566, r24567, MPFR_RNDN); mpfr_fma(r24569, r24560, r24562, r24568, MPFR_RNDN); mpfr_div(r24570, r24558, r24555, MPFR_RNDN); ; mpfr_div(r24572, r24570, r24571, MPFR_RNDN); mpfr_pow(r24573, r24558, r24563, MPFR_RNDN); mpfr_mul(r24574, r24573, r24573, MPFR_RNDN); mpfr_mul(r24575, r24572, r24574, MPFR_RNDN); mpfr_fma(r24576, r24556, r24569, r24575, MPFR_RNDN); mpfr_abs(r24577, r24576, MPFR_RNDN); return mpfr_get_d(r24577, MPFR_RNDN); } static mpfr_t r24578, r24579, r24580, r24581, r24582, r24583, r24584, r24585, r24586, r24587, r24588, r24589, r24590, r24591, r24592, r24593, r24594, r24595, r24596, r24597, r24598, r24599, r24600, r24601, r24602; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24578, "1", 10, MPFR_RNDN); mpfr_init(r24579); mpfr_init(r24580); mpfr_init(r24581); mpfr_init(r24582); mpfr_init(r24583); mpfr_init_set_str(r24584, "5", 10, MPFR_RNDN); mpfr_init(r24585); mpfr_init(r24586); mpfr_init(r24587); mpfr_init_set_str(r24588, "3", 10, MPFR_RNDN); mpfr_init_set_str(r24589, "2", 10, MPFR_RNDN); mpfr_init(r24590); mpfr_init(r24591); mpfr_init(r24592); mpfr_init(r24593); mpfr_init(r24594); mpfr_init(r24595); mpfr_init_set_str(r24596, "21", 10, MPFR_RNDN); mpfr_init(r24597); mpfr_init(r24598); mpfr_init(r24599); mpfr_init(r24600); mpfr_init(r24601); mpfr_init(r24602); } double f_dm(double x) { ; mpfr_const_pi(r24579, MPFR_RNDN); mpfr_sqrt(r24580, r24579, MPFR_RNDN); mpfr_div(r24581, r24578, r24580, MPFR_RNDN); mpfr_set_d(r24582, x, MPFR_RNDN); mpfr_abs(r24583, r24582, MPFR_RNDN); ; mpfr_div(r24585, r24583, r24584, MPFR_RNDN); mpfr_mul(r24586, r24583, r24583, MPFR_RNDN); mpfr_mul(r24587, r24586, r24586, MPFR_RNDN); ; ; mpfr_div(r24590, r24588, r24589, MPFR_RNDN); mpfr_div(r24591, r24583, r24590, MPFR_RNDN); mpfr_mul(r24592, r24589, r24583, MPFR_RNDN); mpfr_fma(r24593, r24586, r24591, r24592, MPFR_RNDN); mpfr_fma(r24594, r24585, r24587, r24593, MPFR_RNDN); mpfr_div(r24595, r24583, r24580, MPFR_RNDN); ; mpfr_div(r24597, r24595, r24596, MPFR_RNDN); mpfr_pow(r24598, r24583, r24588, MPFR_RNDN); mpfr_mul(r24599, r24598, r24598, MPFR_RNDN); mpfr_mul(r24600, r24597, r24599, MPFR_RNDN); mpfr_fma(r24601, r24581, r24594, r24600, MPFR_RNDN); mpfr_abs(r24602, r24601, MPFR_RNDN); return mpfr_get_d(r24602, MPFR_RNDN); }