#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 r24464 = 1; float r24465 = atan2(1.0, 0.0); float r24466 = sqrt(r24465); float r24467 = r24464 / r24466; float r24468 = 2; float r24469 = x; float r24470 = fabs(r24469); float r24471 = r24468 * r24470; float r24472 = 3; float r24473 = r24468 / r24472; float r24474 = r24470 * r24470; float r24475 = r24474 * r24470; float r24476 = r24473 * r24475; float r24477 = r24471 + r24476; float r24478 = 5; float r24479 = r24464 / r24478; float r24480 = r24475 * r24470; float r24481 = r24480 * r24470; float r24482 = r24479 * r24481; float r24483 = r24477 + r24482; float r24484 = 21; float r24485 = r24464 / r24484; float r24486 = r24481 * r24470; float r24487 = r24486 * r24470; float r24488 = r24485 * r24487; float r24489 = r24483 + r24488; float r24490 = r24467 * r24489; float r24491 = fabs(r24490); return r24491; } double f_id(double x) { double r24492 = 1; double r24493 = atan2(1.0, 0.0); double r24494 = sqrt(r24493); double r24495 = r24492 / r24494; double r24496 = 2; double r24497 = x; double r24498 = fabs(r24497); double r24499 = r24496 * r24498; double r24500 = 3; double r24501 = r24496 / r24500; double r24502 = r24498 * r24498; double r24503 = r24502 * r24498; double r24504 = r24501 * r24503; double r24505 = r24499 + r24504; double r24506 = 5; double r24507 = r24492 / r24506; double r24508 = r24503 * r24498; double r24509 = r24508 * r24498; double r24510 = r24507 * r24509; double r24511 = r24505 + r24510; double r24512 = 21; double r24513 = r24492 / r24512; double r24514 = r24509 * r24498; double r24515 = r24514 * r24498; double r24516 = r24513 * r24515; double r24517 = r24511 + r24516; double r24518 = r24495 * r24517; double r24519 = fabs(r24518); return r24519; } double f_of(float x) { float r24520 = 1; float r24521 = atan2(1.0, 0.0); float r24522 = sqrt(r24521); float r24523 = r24520 / r24522; float r24524 = 5; float r24525 = r24520 / r24524; float r24526 = x; float r24527 = fabs(r24526); float r24528 = 3; float r24529 = pow(r24527, r24528); float r24530 = r24527 * r24527; float r24531 = r24529 * r24530; float r24532 = 2/3; float r24533 = 2; float r24534 = r24533 * r24527; float r24535 = fma(r24532, r24529, r24534); float r24536 = fma(r24525, r24531, r24535); float r24537 = 1/21; float r24538 = r24537 / r24522; float r24539 = 7; float r24540 = pow(r24527, r24539); float r24541 = r24538 * r24540; float r24542 = fma(r24523, r24536, r24541); float r24543 = fabs(r24542); return r24543; } double f_od(double x) { double r24544 = 1; double r24545 = atan2(1.0, 0.0); double r24546 = sqrt(r24545); double r24547 = r24544 / r24546; double r24548 = 5; double r24549 = r24544 / r24548; double r24550 = x; double r24551 = fabs(r24550); double r24552 = 3; double r24553 = pow(r24551, r24552); double r24554 = r24551 * r24551; double r24555 = r24553 * r24554; double r24556 = 2/3; double r24557 = 2; double r24558 = r24557 * r24551; double r24559 = fma(r24556, r24553, r24558); double r24560 = fma(r24549, r24555, r24559); double r24561 = 1/21; double r24562 = r24561 / r24546; double r24563 = 7; double r24564 = pow(r24551, r24563); double r24565 = r24562 * r24564; double r24566 = fma(r24547, r24560, r24565); double r24567 = fabs(r24566); return r24567; } 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 r24568, r24569, r24570, r24571, r24572, r24573, r24574, r24575, r24576, r24577, r24578, r24579, r24580, r24581, r24582, r24583, r24584, r24585, r24586, r24587, r24588, r24589, r24590, r24591, r24592, r24593, r24594, r24595; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r24568, "1", 10, MPFR_RNDN); mpfr_init(r24569); mpfr_init(r24570); mpfr_init(r24571); mpfr_init_set_str(r24572, "2", 10, MPFR_RNDN); mpfr_init(r24573); mpfr_init(r24574); mpfr_init(r24575); mpfr_init_set_str(r24576, "3", 10, MPFR_RNDN); mpfr_init(r24577); mpfr_init(r24578); mpfr_init(r24579); mpfr_init(r24580); mpfr_init(r24581); mpfr_init_set_str(r24582, "5", 10, MPFR_RNDN); mpfr_init(r24583); mpfr_init(r24584); mpfr_init(r24585); mpfr_init(r24586); mpfr_init(r24587); mpfr_init_set_str(r24588, "21", 10, MPFR_RNDN); mpfr_init(r24589); mpfr_init(r24590); mpfr_init(r24591); mpfr_init(r24592); mpfr_init(r24593); mpfr_init(r24594); mpfr_init(r24595); } double f_im(double x) { ; mpfr_const_pi(r24569, MPFR_RNDN); mpfr_sqrt(r24570, r24569, MPFR_RNDN); mpfr_div(r24571, r24568, r24570, MPFR_RNDN); ; mpfr_set_d(r24573, x, MPFR_RNDN); mpfr_abs(r24574, r24573, MPFR_RNDN); mpfr_mul(r24575, r24572, r24574, MPFR_RNDN); ; mpfr_div(r24577, r24572, r24576, MPFR_RNDN); mpfr_mul(r24578, r24574, r24574, MPFR_RNDN); mpfr_mul(r24579, r24578, r24574, MPFR_RNDN); mpfr_mul(r24580, r24577, r24579, MPFR_RNDN); mpfr_add(r24581, r24575, r24580, MPFR_RNDN); ; mpfr_div(r24583, r24568, r24582, MPFR_RNDN); mpfr_mul(r24584, r24579, r24574, MPFR_RNDN); mpfr_mul(r24585, r24584, r24574, MPFR_RNDN); mpfr_mul(r24586, r24583, r24585, MPFR_RNDN); mpfr_add(r24587, r24581, r24586, MPFR_RNDN); ; mpfr_div(r24589, r24568, r24588, MPFR_RNDN); mpfr_mul(r24590, r24585, r24574, MPFR_RNDN); mpfr_mul(r24591, r24590, r24574, MPFR_RNDN); mpfr_mul(r24592, r24589, r24591, MPFR_RNDN); mpfr_add(r24593, r24587, r24592, MPFR_RNDN); mpfr_mul(r24594, r24571, r24593, MPFR_RNDN); mpfr_abs(r24595, r24594, MPFR_RNDN); return mpfr_get_d(r24595, MPFR_RNDN); } static mpfr_t r24596, r24597, r24598, r24599, r24600, r24601, r24602, r24603, r24604, r24605, r24606, r24607, r24608, r24609, r24610, r24611, r24612, r24613, r24614, r24615, r24616, r24617, r24618, r24619; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24596, "1", 10, MPFR_RNDN); mpfr_init(r24597); mpfr_init(r24598); mpfr_init(r24599); mpfr_init_set_str(r24600, "5", 10, MPFR_RNDN); mpfr_init(r24601); mpfr_init(r24602); mpfr_init(r24603); mpfr_init_set_str(r24604, "3", 10, MPFR_RNDN); mpfr_init(r24605); mpfr_init(r24606); mpfr_init(r24607); mpfr_init_set_str(r24608, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r24609, "2", 10, MPFR_RNDN); mpfr_init(r24610); mpfr_init(r24611); mpfr_init(r24612); mpfr_init_set_str(r24613, "1/21", 10, MPFR_RNDN); mpfr_init(r24614); mpfr_init_set_str(r24615, "7", 10, MPFR_RNDN); mpfr_init(r24616); mpfr_init(r24617); mpfr_init(r24618); mpfr_init(r24619); } double f_fm(double x) { ; mpfr_const_pi(r24597, MPFR_RNDN); mpfr_sqrt(r24598, r24597, MPFR_RNDN); mpfr_div(r24599, r24596, r24598, MPFR_RNDN); ; mpfr_div(r24601, r24596, r24600, MPFR_RNDN); mpfr_set_d(r24602, x, MPFR_RNDN); mpfr_abs(r24603, r24602, MPFR_RNDN); ; mpfr_pow(r24605, r24603, r24604, MPFR_RNDN); mpfr_mul(r24606, r24603, r24603, MPFR_RNDN); mpfr_mul(r24607, r24605, r24606, MPFR_RNDN); ; ; mpfr_mul(r24610, r24609, r24603, MPFR_RNDN); mpfr_fma(r24611, r24608, r24605, r24610, MPFR_RNDN); mpfr_fma(r24612, r24601, r24607, r24611, MPFR_RNDN); ; mpfr_div(r24614, r24613, r24598, MPFR_RNDN); ; mpfr_pow(r24616, r24603, r24615, MPFR_RNDN); mpfr_mul(r24617, r24614, r24616, MPFR_RNDN); mpfr_fma(r24618, r24599, r24612, r24617, MPFR_RNDN); mpfr_abs(r24619, r24618, MPFR_RNDN); return mpfr_get_d(r24619, MPFR_RNDN); } static mpfr_t r24620, r24621, r24622, r24623, r24624, r24625, r24626, r24627, r24628, r24629, r24630, r24631, r24632, r24633, r24634, r24635, r24636, r24637, r24638, r24639, r24640, r24641, r24642, r24643; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24620, "1", 10, MPFR_RNDN); mpfr_init(r24621); mpfr_init(r24622); mpfr_init(r24623); mpfr_init_set_str(r24624, "5", 10, MPFR_RNDN); mpfr_init(r24625); mpfr_init(r24626); mpfr_init(r24627); mpfr_init_set_str(r24628, "3", 10, MPFR_RNDN); mpfr_init(r24629); mpfr_init(r24630); mpfr_init(r24631); mpfr_init_set_str(r24632, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r24633, "2", 10, MPFR_RNDN); mpfr_init(r24634); mpfr_init(r24635); mpfr_init(r24636); mpfr_init_set_str(r24637, "1/21", 10, MPFR_RNDN); mpfr_init(r24638); mpfr_init_set_str(r24639, "7", 10, MPFR_RNDN); mpfr_init(r24640); mpfr_init(r24641); mpfr_init(r24642); mpfr_init(r24643); } double f_dm(double x) { ; mpfr_const_pi(r24621, MPFR_RNDN); mpfr_sqrt(r24622, r24621, MPFR_RNDN); mpfr_div(r24623, r24620, r24622, MPFR_RNDN); ; mpfr_div(r24625, r24620, r24624, MPFR_RNDN); mpfr_set_d(r24626, x, MPFR_RNDN); mpfr_abs(r24627, r24626, MPFR_RNDN); ; mpfr_pow(r24629, r24627, r24628, MPFR_RNDN); mpfr_mul(r24630, r24627, r24627, MPFR_RNDN); mpfr_mul(r24631, r24629, r24630, MPFR_RNDN); ; ; mpfr_mul(r24634, r24633, r24627, MPFR_RNDN); mpfr_fma(r24635, r24632, r24629, r24634, MPFR_RNDN); mpfr_fma(r24636, r24625, r24631, r24635, MPFR_RNDN); ; mpfr_div(r24638, r24637, r24622, MPFR_RNDN); ; mpfr_pow(r24640, r24627, r24639, MPFR_RNDN); mpfr_mul(r24641, r24638, r24640, MPFR_RNDN); mpfr_fma(r24642, r24623, r24636, r24641, MPFR_RNDN); mpfr_abs(r24643, r24642, MPFR_RNDN); return mpfr_get_d(r24643, MPFR_RNDN); }