#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r22602 = x; float r22603 = y; float r22604 = z; float r22605 = fma(r22602, r22603, r22604); float r22606 = 1; float r22607 = r22602 * r22603; float r22608 = r22607 + r22604; float r22609 = r22606 + r22608; float r22610 = r22605 - r22609; return r22610; } double f_id(double x, double y, double z) { double r22611 = x; double r22612 = y; double r22613 = z; double r22614 = fma(r22611, r22612, r22613); double r22615 = 1; double r22616 = r22611 * r22612; double r22617 = r22616 + r22613; double r22618 = r22615 + r22617; double r22619 = r22614 - r22618; return r22619; } double f_of(float x, float y, float z) { float r22620 = x; float r22621 = y; float r22622 = z; float r22623 = fma(r22620, r22621, r22622); float r22624 = r22623 - r22622; float r22625 = r22621 * r22620; float r22626 = r22624 - r22625; float r22627 = 1; float r22628 = r22626 - r22627; float r22629 = exp(r22628); float r22630 = log(r22629); float r22631 = exp(r22630); float r22632 = log(r22631); float r22633 = -1.0000000000000462; bool r22634 = r22632 <= r22633; float r22635 = r22623 - r22625; float r22636 = r22635 - r22622; float r22637 = r22636 - r22627; float r22638 = exp(r22637); float r22639 = log(r22638); float r22640 = -1.0; bool r22641 = r22632 <= r22640; float r22642 = r22641 ? r22632 : r22639; float r22643 = r22634 ? r22639 : r22642; return r22643; } double f_od(double x, double y, double z) { double r22644 = x; double r22645 = y; double r22646 = z; double r22647 = fma(r22644, r22645, r22646); double r22648 = r22647 - r22646; double r22649 = r22645 * r22644; double r22650 = r22648 - r22649; double r22651 = 1; double r22652 = r22650 - r22651; double r22653 = exp(r22652); double r22654 = log(r22653); double r22655 = exp(r22654); double r22656 = log(r22655); double r22657 = -1.0000000000000462; bool r22658 = r22656 <= r22657; double r22659 = r22647 - r22649; double r22660 = r22659 - r22646; double r22661 = r22660 - r22651; double r22662 = exp(r22661); double r22663 = log(r22662); double r22664 = -1.0; bool r22665 = r22656 <= r22664; double r22666 = r22665 ? r22656 : r22663; double r22667 = r22658 ? r22663 : r22666; return r22667; } 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 r22668, r22669, r22670, r22671, r22672, r22673, r22674, r22675, r22676; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r22668); mpfr_init(r22669); mpfr_init(r22670); mpfr_init(r22671); mpfr_init_set_str(r22672, "1", 10, MPFR_RNDN); mpfr_init(r22673); mpfr_init(r22674); mpfr_init(r22675); mpfr_init(r22676); } double f_im(double x, double y, double z) { mpfr_set_d(r22668, x, MPFR_RNDN); mpfr_set_d(r22669, y, MPFR_RNDN); mpfr_set_d(r22670, z, MPFR_RNDN); mpfr_fma(r22671, r22668, r22669, r22670, MPFR_RNDN); ; mpfr_mul(r22673, r22668, r22669, MPFR_RNDN); mpfr_add(r22674, r22673, r22670, MPFR_RNDN); mpfr_add(r22675, r22672, r22674, MPFR_RNDN); mpfr_sub(r22676, r22671, r22675, MPFR_RNDN); return mpfr_get_d(r22676, MPFR_RNDN); } static mpfr_t r22677, r22678, r22679, r22680, r22681, r22682, r22683, r22684, r22685, r22686, r22687, r22688, r22689, r22690, r22691, r22692, r22693, r22694, r22695, r22696, r22697, r22698, r22699, r22700; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r22677); mpfr_init(r22678); mpfr_init(r22679); mpfr_init(r22680); mpfr_init(r22681); mpfr_init(r22682); mpfr_init(r22683); mpfr_init_set_str(r22684, "1", 10, MPFR_RNDN); mpfr_init(r22685); mpfr_init(r22686); mpfr_init(r22687); mpfr_init(r22688); mpfr_init(r22689); mpfr_init_set_str(r22690, "-1.0000000000000462", 10, MPFR_RNDN); mpfr_init(r22691); mpfr_init(r22692); mpfr_init(r22693); mpfr_init(r22694); mpfr_init(r22695); mpfr_init(r22696); mpfr_init_set_str(r22697, "-1.0", 10, MPFR_RNDN); mpfr_init(r22698); mpfr_init(r22699); mpfr_init(r22700); } double f_fm(double x, double y, double z) { mpfr_set_d(r22677, x, MPFR_RNDN); mpfr_set_d(r22678, y, MPFR_RNDN); mpfr_set_d(r22679, z, MPFR_RNDN); mpfr_fma(r22680, r22677, r22678, r22679, MPFR_RNDN); mpfr_sub(r22681, r22680, r22679, MPFR_RNDN); mpfr_mul(r22682, r22678, r22677, MPFR_RNDN); mpfr_sub(r22683, r22681, r22682, MPFR_RNDN); ; mpfr_sub(r22685, r22683, r22684, MPFR_RNDN); mpfr_exp(r22686, r22685, MPFR_RNDN); mpfr_log(r22687, r22686, MPFR_RNDN); mpfr_exp(r22688, r22687, MPFR_RNDN); mpfr_log(r22689, r22688, MPFR_RNDN); ; mpfr_set_si(r22691, mpfr_cmp(r22689, r22690) <= 0, MPFR_RNDN); mpfr_sub(r22692, r22680, r22682, MPFR_RNDN); mpfr_sub(r22693, r22692, r22679, MPFR_RNDN); mpfr_sub(r22694, r22693, r22684, MPFR_RNDN); mpfr_exp(r22695, r22694, MPFR_RNDN); mpfr_log(r22696, r22695, MPFR_RNDN); ; mpfr_set_si(r22698, mpfr_cmp(r22689, r22697) <= 0, MPFR_RNDN); if (mpfr_get_si(r22698, MPFR_RNDN)) { mpfr_set(r22699, r22689, MPFR_RNDN); } else { mpfr_set(r22699, r22696, MPFR_RNDN); }; if (mpfr_get_si(r22691, MPFR_RNDN)) { mpfr_set(r22700, r22696, MPFR_RNDN); } else { mpfr_set(r22700, r22699, MPFR_RNDN); }; return mpfr_get_d(r22700, MPFR_RNDN); } static mpfr_t r22701, r22702, r22703, r22704, r22705, r22706, r22707, r22708, r22709, r22710, r22711, r22712, r22713, r22714, r22715, r22716, r22717, r22718, r22719, r22720, r22721, r22722, r22723, r22724; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r22701); mpfr_init(r22702); mpfr_init(r22703); mpfr_init(r22704); mpfr_init(r22705); mpfr_init(r22706); mpfr_init(r22707); mpfr_init_set_str(r22708, "1", 10, MPFR_RNDN); mpfr_init(r22709); mpfr_init(r22710); mpfr_init(r22711); mpfr_init(r22712); mpfr_init(r22713); mpfr_init_set_str(r22714, "-1.0000000000000462", 10, MPFR_RNDN); mpfr_init(r22715); mpfr_init(r22716); mpfr_init(r22717); mpfr_init(r22718); mpfr_init(r22719); mpfr_init(r22720); mpfr_init_set_str(r22721, "-1.0", 10, MPFR_RNDN); mpfr_init(r22722); mpfr_init(r22723); mpfr_init(r22724); } double f_dm(double x, double y, double z) { mpfr_set_d(r22701, x, MPFR_RNDN); mpfr_set_d(r22702, y, MPFR_RNDN); mpfr_set_d(r22703, z, MPFR_RNDN); mpfr_fma(r22704, r22701, r22702, r22703, MPFR_RNDN); mpfr_sub(r22705, r22704, r22703, MPFR_RNDN); mpfr_mul(r22706, r22702, r22701, MPFR_RNDN); mpfr_sub(r22707, r22705, r22706, MPFR_RNDN); ; mpfr_sub(r22709, r22707, r22708, MPFR_RNDN); mpfr_exp(r22710, r22709, MPFR_RNDN); mpfr_log(r22711, r22710, MPFR_RNDN); mpfr_exp(r22712, r22711, MPFR_RNDN); mpfr_log(r22713, r22712, MPFR_RNDN); ; mpfr_set_si(r22715, mpfr_cmp(r22713, r22714) <= 0, MPFR_RNDN); mpfr_sub(r22716, r22704, r22706, MPFR_RNDN); mpfr_sub(r22717, r22716, r22703, MPFR_RNDN); mpfr_sub(r22718, r22717, r22708, MPFR_RNDN); mpfr_exp(r22719, r22718, MPFR_RNDN); mpfr_log(r22720, r22719, MPFR_RNDN); ; mpfr_set_si(r22722, mpfr_cmp(r22713, r22721) <= 0, MPFR_RNDN); if (mpfr_get_si(r22722, MPFR_RNDN)) { mpfr_set(r22723, r22713, MPFR_RNDN); } else { mpfr_set(r22723, r22720, MPFR_RNDN); }; if (mpfr_get_si(r22715, MPFR_RNDN)) { mpfr_set(r22724, r22720, MPFR_RNDN); } else { mpfr_set(r22724, r22723, MPFR_RNDN); }; return mpfr_get_d(r22724, MPFR_RNDN); }