#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r22577 = x; float r22578 = y; float r22579 = z; float r22580 = fma(r22577, r22578, r22579); float r22581 = 1; float r22582 = r22577 * r22578; float r22583 = r22582 + r22579; float r22584 = r22581 + r22583; float r22585 = r22580 - r22584; return r22585; } double f_id(double x, double y, double z) { double r22586 = x; double r22587 = y; double r22588 = z; double r22589 = fma(r22586, r22587, r22588); double r22590 = 1; double r22591 = r22586 * r22587; double r22592 = r22591 + r22588; double r22593 = r22590 + r22592; double r22594 = r22589 - r22593; return r22594; } double f_of(float x, float y, float z) { float r22595 = x; float r22596 = y; float r22597 = z; float r22598 = fma(r22595, r22596, r22597); float r22599 = r22596 * r22595; float r22600 = r22598 - r22599; float r22601 = r22600 - r22597; float r22602 = 1; float r22603 = r22601 - r22602; float r22604 = exp(r22603); float r22605 = log(r22604); float r22606 = -1.0000000000000036; bool r22607 = r22605 <= r22606; float r22608 = r22598 - r22597; float r22609 = r22608 - r22599; float r22610 = r22609 - r22602; float r22611 = exp(r22610); float r22612 = log(r22611); float r22613 = exp(r22612); float r22614 = log(r22613); float r22615 = -1.0; bool r22616 = r22605 <= r22615; float r22617 = r22616 ? r22605 : r22614; float r22618 = r22607 ? r22614 : r22617; return r22618; } double f_od(double x, double y, double z) { double r22619 = x; double r22620 = y; double r22621 = z; double r22622 = fma(r22619, r22620, r22621); double r22623 = r22620 * r22619; double r22624 = r22622 - r22623; double r22625 = r22624 - r22621; double r22626 = 1; double r22627 = r22625 - r22626; double r22628 = exp(r22627); double r22629 = log(r22628); double r22630 = -1.0000000000000036; bool r22631 = r22629 <= r22630; double r22632 = r22622 - r22621; double r22633 = r22632 - r22623; double r22634 = r22633 - r22626; double r22635 = exp(r22634); double r22636 = log(r22635); double r22637 = exp(r22636); double r22638 = log(r22637); double r22639 = -1.0; bool r22640 = r22629 <= r22639; double r22641 = r22640 ? r22629 : r22638; double r22642 = r22631 ? r22638 : r22641; return r22642; } 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 r22643, r22644, r22645, r22646, r22647, r22648, r22649, r22650, r22651; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r22643); mpfr_init(r22644); mpfr_init(r22645); mpfr_init(r22646); mpfr_init_set_str(r22647, "1", 10, MPFR_RNDN); mpfr_init(r22648); mpfr_init(r22649); mpfr_init(r22650); mpfr_init(r22651); } double f_im(double x, double y, double z) { mpfr_set_d(r22643, x, MPFR_RNDN); mpfr_set_d(r22644, y, MPFR_RNDN); mpfr_set_d(r22645, z, MPFR_RNDN); mpfr_fma(r22646, r22643, r22644, r22645, MPFR_RNDN); ; mpfr_mul(r22648, r22643, r22644, MPFR_RNDN); mpfr_add(r22649, r22648, r22645, MPFR_RNDN); mpfr_add(r22650, r22647, r22649, MPFR_RNDN); mpfr_sub(r22651, r22646, r22650, MPFR_RNDN); return mpfr_get_d(r22651, MPFR_RNDN); } static mpfr_t r22652, r22653, r22654, r22655, r22656, r22657, r22658, r22659, r22660, r22661, r22662, r22663, r22664, r22665, r22666, r22667, r22668, r22669, r22670, r22671, r22672, r22673, r22674, r22675; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r22652); mpfr_init(r22653); mpfr_init(r22654); mpfr_init(r22655); mpfr_init(r22656); mpfr_init(r22657); mpfr_init(r22658); mpfr_init_set_str(r22659, "1", 10, MPFR_RNDN); mpfr_init(r22660); mpfr_init(r22661); mpfr_init(r22662); mpfr_init_set_str(r22663, "-1.0000000000000036", 10, MPFR_RNDN); mpfr_init(r22664); mpfr_init(r22665); mpfr_init(r22666); mpfr_init(r22667); mpfr_init(r22668); mpfr_init(r22669); mpfr_init(r22670); mpfr_init(r22671); mpfr_init_set_str(r22672, "-1.0", 10, MPFR_RNDN); mpfr_init(r22673); mpfr_init(r22674); mpfr_init(r22675); } double f_fm(double x, double y, double z) { mpfr_set_d(r22652, x, MPFR_RNDN); mpfr_set_d(r22653, y, MPFR_RNDN); mpfr_set_d(r22654, z, MPFR_RNDN); mpfr_fma(r22655, r22652, r22653, r22654, MPFR_RNDN); mpfr_mul(r22656, r22653, r22652, MPFR_RNDN); mpfr_sub(r22657, r22655, r22656, MPFR_RNDN); mpfr_sub(r22658, r22657, r22654, MPFR_RNDN); ; mpfr_sub(r22660, r22658, r22659, MPFR_RNDN); mpfr_exp(r22661, r22660, MPFR_RNDN); mpfr_log(r22662, r22661, MPFR_RNDN); ; mpfr_set_si(r22664, mpfr_cmp(r22662, r22663) <= 0, MPFR_RNDN); mpfr_sub(r22665, r22655, r22654, MPFR_RNDN); mpfr_sub(r22666, r22665, r22656, MPFR_RNDN); mpfr_sub(r22667, r22666, r22659, MPFR_RNDN); mpfr_exp(r22668, r22667, MPFR_RNDN); mpfr_log(r22669, r22668, MPFR_RNDN); mpfr_exp(r22670, r22669, MPFR_RNDN); mpfr_log(r22671, r22670, MPFR_RNDN); ; mpfr_set_si(r22673, mpfr_cmp(r22662, r22672) <= 0, MPFR_RNDN); if (mpfr_get_si(r22673, MPFR_RNDN)) { mpfr_set(r22674, r22662, MPFR_RNDN); } else { mpfr_set(r22674, r22671, MPFR_RNDN); }; if (mpfr_get_si(r22664, MPFR_RNDN)) { mpfr_set(r22675, r22671, MPFR_RNDN); } else { mpfr_set(r22675, r22674, MPFR_RNDN); }; return mpfr_get_d(r22675, MPFR_RNDN); } static mpfr_t r22676, r22677, r22678, r22679, r22680, r22681, r22682, r22683, r22684, r22685, r22686, r22687, r22688, r22689, r22690, r22691, r22692, r22693, r22694, r22695, r22696, r22697, r22698, r22699; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r22676); mpfr_init(r22677); mpfr_init(r22678); mpfr_init(r22679); mpfr_init(r22680); mpfr_init(r22681); mpfr_init(r22682); mpfr_init_set_str(r22683, "1", 10, MPFR_RNDN); mpfr_init(r22684); mpfr_init(r22685); mpfr_init(r22686); mpfr_init_set_str(r22687, "-1.0000000000000036", 10, MPFR_RNDN); mpfr_init(r22688); mpfr_init(r22689); mpfr_init(r22690); mpfr_init(r22691); mpfr_init(r22692); mpfr_init(r22693); mpfr_init(r22694); mpfr_init(r22695); mpfr_init_set_str(r22696, "-1.0", 10, MPFR_RNDN); mpfr_init(r22697); mpfr_init(r22698); mpfr_init(r22699); } double f_dm(double x, double y, double z) { mpfr_set_d(r22676, x, MPFR_RNDN); mpfr_set_d(r22677, y, MPFR_RNDN); mpfr_set_d(r22678, z, MPFR_RNDN); mpfr_fma(r22679, r22676, r22677, r22678, MPFR_RNDN); mpfr_mul(r22680, r22677, r22676, MPFR_RNDN); mpfr_sub(r22681, r22679, r22680, MPFR_RNDN); mpfr_sub(r22682, r22681, r22678, MPFR_RNDN); ; mpfr_sub(r22684, r22682, r22683, MPFR_RNDN); mpfr_exp(r22685, r22684, MPFR_RNDN); mpfr_log(r22686, r22685, MPFR_RNDN); ; mpfr_set_si(r22688, mpfr_cmp(r22686, r22687) <= 0, MPFR_RNDN); mpfr_sub(r22689, r22679, r22678, MPFR_RNDN); mpfr_sub(r22690, r22689, r22680, MPFR_RNDN); mpfr_sub(r22691, r22690, r22683, MPFR_RNDN); mpfr_exp(r22692, r22691, MPFR_RNDN); mpfr_log(r22693, r22692, MPFR_RNDN); mpfr_exp(r22694, r22693, MPFR_RNDN); mpfr_log(r22695, r22694, MPFR_RNDN); ; mpfr_set_si(r22697, mpfr_cmp(r22686, r22696) <= 0, MPFR_RNDN); if (mpfr_get_si(r22697, MPFR_RNDN)) { mpfr_set(r22698, r22686, MPFR_RNDN); } else { mpfr_set(r22698, r22695, MPFR_RNDN); }; if (mpfr_get_si(r22688, MPFR_RNDN)) { mpfr_set(r22699, r22695, MPFR_RNDN); } else { mpfr_set(r22699, r22698, MPFR_RNDN); }; return mpfr_get_d(r22699, MPFR_RNDN); }