#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r18522 = a; float r18523 = r18522 * r18522; float r18524 = b; float r18525 = r18524 * r18524; float r18526 = r18523 + r18525; float r18527 = r18526 * r18526; float r18528 = 4.0f; float r18529 = 1.0f; float r18530 = r18529 + r18522; float r18531 = r18523 * r18530; float r18532 = 3.0f; float r18533 = r18532 * r18522; float r18534 = r18529 - r18533; float r18535 = r18525 * r18534; float r18536 = r18531 + r18535; float r18537 = r18528 * r18536; float r18538 = r18527 + r18537; float r18539 = r18538 - r18529; return r18539; } double f_id(double a, double b) { double r18540 = a; double r18541 = r18540 * r18540; double r18542 = b; double r18543 = r18542 * r18542; double r18544 = r18541 + r18543; double r18545 = r18544 * r18544; double r18546 = 4.0; double r18547 = 1.0; double r18548 = r18547 + r18540; double r18549 = r18541 * r18548; double r18550 = 3.0; double r18551 = r18550 * r18540; double r18552 = r18547 - r18551; double r18553 = r18543 * r18552; double r18554 = r18549 + r18553; double r18555 = r18546 * r18554; double r18556 = r18545 + r18555; double r18557 = r18556 - r18547; return r18557; } double f_of(float a, float b) { float r18558 = 1.0f; float r18559 = a; float r18560 = 3.0f; float r18561 = r18559 * r18560; float r18562 = r18558 - r18561; float r18563 = b; float r18564 = r18563 * r18563; float r18565 = fma(r18559, r18559, r18559); float r18566 = r18565 * r18559; float r18567 = fma(r18562, r18564, r18566); float r18568 = 4.0f; float r18569 = 2.0f; float r18570 = r18563 * r18569; float r18571 = r18570 * r18563; float r18572 = r18559 * r18559; float r18573 = pow(r18559, r18568); float r18574 = pow(r18563, r18568); float r18575 = r18573 + r18574; float r18576 = fma(r18571, r18572, r18575); float r18577 = fma(r18567, r18568, r18576); float r18578 = r18577 - r18558; return r18578; } double f_od(double a, double b) { double r18579 = 1.0; double r18580 = a; double r18581 = 3.0; double r18582 = r18580 * r18581; double r18583 = r18579 - r18582; double r18584 = b; double r18585 = r18584 * r18584; double r18586 = fma(r18580, r18580, r18580); double r18587 = r18586 * r18580; double r18588 = fma(r18583, r18585, r18587); double r18589 = 4.0; double r18590 = 2.0; double r18591 = r18584 * r18590; double r18592 = r18591 * r18584; double r18593 = r18580 * r18580; double r18594 = pow(r18580, r18589); double r18595 = pow(r18584, r18589); double r18596 = r18594 + r18595; double r18597 = fma(r18592, r18593, r18596); double r18598 = fma(r18588, r18589, r18597); double r18599 = r18598 - r18579; return r18599; } 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 r18600, r18601, r18602, r18603, r18604, r18605, r18606, r18607, r18608, r18609, r18610, r18611, r18612, r18613, r18614, r18615, r18616, r18617; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r18600); mpfr_init(r18601); mpfr_init(r18602); mpfr_init(r18603); mpfr_init(r18604); mpfr_init(r18605); mpfr_init_set_str(r18606, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18607, "1", 10, MPFR_RNDN); mpfr_init(r18608); mpfr_init(r18609); mpfr_init_set_str(r18610, "3", 10, MPFR_RNDN); mpfr_init(r18611); mpfr_init(r18612); mpfr_init(r18613); mpfr_init(r18614); mpfr_init(r18615); mpfr_init(r18616); mpfr_init(r18617); } double f_im(double a, double b) { mpfr_set_d(r18600, a, MPFR_RNDN); mpfr_sqr(r18601, r18600, MPFR_RNDN); mpfr_set_d(r18602, b, MPFR_RNDN); mpfr_sqr(r18603, r18602, MPFR_RNDN); mpfr_add(r18604, r18601, r18603, MPFR_RNDN); mpfr_sqr(r18605, r18604, MPFR_RNDN); ; ; mpfr_add(r18608, r18607, r18600, MPFR_RNDN); mpfr_mul(r18609, r18601, r18608, MPFR_RNDN); ; mpfr_mul(r18611, r18610, r18600, MPFR_RNDN); mpfr_sub(r18612, r18607, r18611, MPFR_RNDN); mpfr_mul(r18613, r18603, r18612, MPFR_RNDN); mpfr_add(r18614, r18609, r18613, MPFR_RNDN); mpfr_mul(r18615, r18606, r18614, MPFR_RNDN); mpfr_add(r18616, r18605, r18615, MPFR_RNDN); mpfr_sub(r18617, r18616, r18607, MPFR_RNDN); return mpfr_get_d(r18617, MPFR_RNDN); } static mpfr_t r18618, r18619, r18620, r18621, r18622, r18623, r18624, r18625, r18626, r18627, r18628, r18629, r18630, r18631, r18632, r18633, r18634, r18635, r18636, r18637, r18638; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18618, "1", 10, MPFR_RNDN); mpfr_init(r18619); mpfr_init_set_str(r18620, "3", 10, MPFR_RNDN); mpfr_init(r18621); mpfr_init(r18622); mpfr_init(r18623); mpfr_init(r18624); mpfr_init(r18625); mpfr_init(r18626); mpfr_init(r18627); mpfr_init_set_str(r18628, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18629, "2", 10, MPFR_RNDN); mpfr_init(r18630); mpfr_init(r18631); mpfr_init(r18632); mpfr_init(r18633); mpfr_init(r18634); mpfr_init(r18635); mpfr_init(r18636); mpfr_init(r18637); mpfr_init(r18638); } double f_fm(double a, double b) { ; mpfr_set_d(r18619, a, MPFR_RNDN); ; mpfr_mul(r18621, r18619, r18620, MPFR_RNDN); mpfr_sub(r18622, r18618, r18621, MPFR_RNDN); mpfr_set_d(r18623, b, MPFR_RNDN); mpfr_sqr(r18624, r18623, MPFR_RNDN); mpfr_fma(r18625, r18619, r18619, r18619, MPFR_RNDN); mpfr_mul(r18626, r18625, r18619, MPFR_RNDN); mpfr_fma(r18627, r18622, r18624, r18626, MPFR_RNDN); ; ; mpfr_mul(r18630, r18623, r18629, MPFR_RNDN); mpfr_mul(r18631, r18630, r18623, MPFR_RNDN); mpfr_mul(r18632, r18619, r18619, MPFR_RNDN); mpfr_pow(r18633, r18619, r18628, MPFR_RNDN); mpfr_pow(r18634, r18623, r18628, MPFR_RNDN); mpfr_add(r18635, r18633, r18634, MPFR_RNDN); mpfr_fma(r18636, r18631, r18632, r18635, MPFR_RNDN); mpfr_fma(r18637, r18627, r18628, r18636, MPFR_RNDN); mpfr_sub(r18638, r18637, r18618, MPFR_RNDN); return mpfr_get_d(r18638, MPFR_RNDN); } static mpfr_t r18639, r18640, r18641, r18642, r18643, r18644, r18645, r18646, r18647, r18648, r18649, r18650, r18651, r18652, r18653, r18654, r18655, r18656, r18657, r18658, r18659; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18639, "1", 10, MPFR_RNDN); mpfr_init(r18640); mpfr_init_set_str(r18641, "3", 10, MPFR_RNDN); mpfr_init(r18642); mpfr_init(r18643); mpfr_init(r18644); mpfr_init(r18645); mpfr_init(r18646); mpfr_init(r18647); mpfr_init(r18648); mpfr_init_set_str(r18649, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18650, "2", 10, MPFR_RNDN); mpfr_init(r18651); mpfr_init(r18652); mpfr_init(r18653); mpfr_init(r18654); mpfr_init(r18655); mpfr_init(r18656); mpfr_init(r18657); mpfr_init(r18658); mpfr_init(r18659); } double f_dm(double a, double b) { ; mpfr_set_d(r18640, a, MPFR_RNDN); ; mpfr_mul(r18642, r18640, r18641, MPFR_RNDN); mpfr_sub(r18643, r18639, r18642, MPFR_RNDN); mpfr_set_d(r18644, b, MPFR_RNDN); mpfr_sqr(r18645, r18644, MPFR_RNDN); mpfr_fma(r18646, r18640, r18640, r18640, MPFR_RNDN); mpfr_mul(r18647, r18646, r18640, MPFR_RNDN); mpfr_fma(r18648, r18643, r18645, r18647, MPFR_RNDN); ; ; mpfr_mul(r18651, r18644, r18650, MPFR_RNDN); mpfr_mul(r18652, r18651, r18644, MPFR_RNDN); mpfr_mul(r18653, r18640, r18640, MPFR_RNDN); mpfr_pow(r18654, r18640, r18649, MPFR_RNDN); mpfr_pow(r18655, r18644, r18649, MPFR_RNDN); mpfr_add(r18656, r18654, r18655, MPFR_RNDN); mpfr_fma(r18657, r18652, r18653, r18656, MPFR_RNDN); mpfr_fma(r18658, r18648, r18649, r18657, MPFR_RNDN); mpfr_sub(r18659, r18658, r18639, MPFR_RNDN); return mpfr_get_d(r18659, MPFR_RNDN); }