#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r22548 = x; float r22549 = y; float r22550 = z; float r22551 = fma(r22548, r22549, r22550); float r22552 = 1; float r22553 = r22548 * r22549; float r22554 = r22553 + r22550; float r22555 = r22552 + r22554; float r22556 = r22551 - r22555; return r22556; } double f_id(double x, double y, double z) { double r22557 = x; double r22558 = y; double r22559 = z; double r22560 = fma(r22557, r22558, r22559); double r22561 = 1; double r22562 = r22557 * r22558; double r22563 = r22562 + r22559; double r22564 = r22561 + r22563; double r22565 = r22560 - r22564; return r22565; } double f_of(float x, float y, float z) { float r22566 = x; float r22567 = y; float r22568 = z; float r22569 = fma(r22566, r22567, r22568); float r22570 = r22569 - r22568; float r22571 = r22567 * r22566; float r22572 = r22570 - r22571; float r22573 = 1; float r22574 = r22572 - r22573; float r22575 = exp(r22574); float r22576 = log(r22575); float r22577 = exp(r22576); float r22578 = log(r22577); float r22579 = -1.0000002541637818; bool r22580 = r22578 <= r22579; float r22581 = r22569 - r22571; float r22582 = r22581 - r22568; float r22583 = r22582 - r22573; float r22584 = exp(r22583); float r22585 = log(r22584); float r22586 = -1.0; bool r22587 = r22578 <= r22586; float r22588 = r22587 ? r22578 : r22585; float r22589 = r22580 ? r22585 : r22588; return r22589; } double f_od(double x, double y, double z) { double r22590 = x; double r22591 = y; double r22592 = z; double r22593 = fma(r22590, r22591, r22592); double r22594 = r22593 - r22592; double r22595 = r22591 * r22590; double r22596 = r22594 - r22595; double r22597 = 1; double r22598 = r22596 - r22597; double r22599 = exp(r22598); double r22600 = log(r22599); double r22601 = exp(r22600); double r22602 = log(r22601); double r22603 = -1.0000002541637818; bool r22604 = r22602 <= r22603; double r22605 = r22593 - r22595; double r22606 = r22605 - r22592; double r22607 = r22606 - r22597; double r22608 = exp(r22607); double r22609 = log(r22608); double r22610 = -1.0; bool r22611 = r22602 <= r22610; double r22612 = r22611 ? r22602 : r22609; double r22613 = r22604 ? r22609 : r22612; return r22613; } 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 r22614, r22615, r22616, r22617, r22618, r22619, r22620, r22621, r22622; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r22614); mpfr_init(r22615); mpfr_init(r22616); mpfr_init(r22617); mpfr_init_set_str(r22618, "1", 10, MPFR_RNDN); mpfr_init(r22619); mpfr_init(r22620); mpfr_init(r22621); mpfr_init(r22622); } double f_im(double x, double y, double z) { mpfr_set_d(r22614, x, MPFR_RNDN); mpfr_set_d(r22615, y, MPFR_RNDN); mpfr_set_d(r22616, z, MPFR_RNDN); mpfr_fma(r22617, r22614, r22615, r22616, MPFR_RNDN); ; mpfr_mul(r22619, r22614, r22615, MPFR_RNDN); mpfr_add(r22620, r22619, r22616, MPFR_RNDN); mpfr_add(r22621, r22618, r22620, MPFR_RNDN); mpfr_sub(r22622, r22617, r22621, MPFR_RNDN); return mpfr_get_d(r22622, MPFR_RNDN); } static mpfr_t r22623, r22624, r22625, r22626, r22627, r22628, r22629, r22630, r22631, r22632, r22633, r22634, r22635, r22636, r22637, r22638, r22639, r22640, r22641, r22642, r22643, r22644, r22645, r22646; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r22623); mpfr_init(r22624); mpfr_init(r22625); mpfr_init(r22626); mpfr_init(r22627); mpfr_init(r22628); mpfr_init(r22629); mpfr_init_set_str(r22630, "1", 10, MPFR_RNDN); mpfr_init(r22631); mpfr_init(r22632); mpfr_init(r22633); mpfr_init(r22634); mpfr_init(r22635); mpfr_init_set_str(r22636, "-1.0000002541637818", 10, MPFR_RNDN); mpfr_init(r22637); mpfr_init(r22638); mpfr_init(r22639); mpfr_init(r22640); mpfr_init(r22641); mpfr_init(r22642); mpfr_init_set_str(r22643, "-1.0", 10, MPFR_RNDN); mpfr_init(r22644); mpfr_init(r22645); mpfr_init(r22646); } double f_fm(double x, double y, double z) { mpfr_set_d(r22623, x, MPFR_RNDN); mpfr_set_d(r22624, y, MPFR_RNDN); mpfr_set_d(r22625, z, MPFR_RNDN); mpfr_fma(r22626, r22623, r22624, r22625, MPFR_RNDN); mpfr_sub(r22627, r22626, r22625, MPFR_RNDN); mpfr_mul(r22628, r22624, r22623, MPFR_RNDN); mpfr_sub(r22629, r22627, r22628, MPFR_RNDN); ; mpfr_sub(r22631, r22629, r22630, MPFR_RNDN); mpfr_exp(r22632, r22631, MPFR_RNDN); mpfr_log(r22633, r22632, MPFR_RNDN); mpfr_exp(r22634, r22633, MPFR_RNDN); mpfr_log(r22635, r22634, MPFR_RNDN); ; mpfr_set_si(r22637, mpfr_cmp(r22635, r22636) <= 0, MPFR_RNDN); mpfr_sub(r22638, r22626, r22628, MPFR_RNDN); mpfr_sub(r22639, r22638, r22625, MPFR_RNDN); mpfr_sub(r22640, r22639, r22630, MPFR_RNDN); mpfr_exp(r22641, r22640, MPFR_RNDN); mpfr_log(r22642, r22641, MPFR_RNDN); ; mpfr_set_si(r22644, mpfr_cmp(r22635, r22643) <= 0, MPFR_RNDN); if (mpfr_get_si(r22644, MPFR_RNDN)) { mpfr_set(r22645, r22635, MPFR_RNDN); } else { mpfr_set(r22645, r22642, MPFR_RNDN); }; if (mpfr_get_si(r22637, MPFR_RNDN)) { mpfr_set(r22646, r22642, MPFR_RNDN); } else { mpfr_set(r22646, r22645, MPFR_RNDN); }; return mpfr_get_d(r22646, MPFR_RNDN); } static mpfr_t r22647, r22648, r22649, r22650, r22651, r22652, r22653, r22654, r22655, r22656, r22657, r22658, r22659, r22660, r22661, r22662, r22663, r22664, r22665, r22666, r22667, r22668, r22669, r22670; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r22647); mpfr_init(r22648); mpfr_init(r22649); mpfr_init(r22650); mpfr_init(r22651); mpfr_init(r22652); mpfr_init(r22653); mpfr_init_set_str(r22654, "1", 10, MPFR_RNDN); mpfr_init(r22655); mpfr_init(r22656); mpfr_init(r22657); mpfr_init(r22658); mpfr_init(r22659); mpfr_init_set_str(r22660, "-1.0000002541637818", 10, MPFR_RNDN); mpfr_init(r22661); mpfr_init(r22662); mpfr_init(r22663); mpfr_init(r22664); mpfr_init(r22665); mpfr_init(r22666); mpfr_init_set_str(r22667, "-1.0", 10, MPFR_RNDN); mpfr_init(r22668); mpfr_init(r22669); mpfr_init(r22670); } double f_dm(double x, double y, double z) { mpfr_set_d(r22647, x, MPFR_RNDN); mpfr_set_d(r22648, y, MPFR_RNDN); mpfr_set_d(r22649, z, MPFR_RNDN); mpfr_fma(r22650, r22647, r22648, r22649, MPFR_RNDN); mpfr_sub(r22651, r22650, r22649, MPFR_RNDN); mpfr_mul(r22652, r22648, r22647, MPFR_RNDN); mpfr_sub(r22653, r22651, r22652, MPFR_RNDN); ; mpfr_sub(r22655, r22653, r22654, MPFR_RNDN); mpfr_exp(r22656, r22655, MPFR_RNDN); mpfr_log(r22657, r22656, MPFR_RNDN); mpfr_exp(r22658, r22657, MPFR_RNDN); mpfr_log(r22659, r22658, MPFR_RNDN); ; mpfr_set_si(r22661, mpfr_cmp(r22659, r22660) <= 0, MPFR_RNDN); mpfr_sub(r22662, r22650, r22652, MPFR_RNDN); mpfr_sub(r22663, r22662, r22649, MPFR_RNDN); mpfr_sub(r22664, r22663, r22654, MPFR_RNDN); mpfr_exp(r22665, r22664, MPFR_RNDN); mpfr_log(r22666, r22665, MPFR_RNDN); ; mpfr_set_si(r22668, mpfr_cmp(r22659, r22667) <= 0, MPFR_RNDN); if (mpfr_get_si(r22668, MPFR_RNDN)) { mpfr_set(r22669, r22659, MPFR_RNDN); } else { mpfr_set(r22669, r22666, MPFR_RNDN); }; if (mpfr_get_si(r22661, MPFR_RNDN)) { mpfr_set(r22670, r22666, MPFR_RNDN); } else { mpfr_set(r22670, r22669, MPFR_RNDN); }; return mpfr_get_d(r22670, MPFR_RNDN); }