#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r22646 = x; float r22647 = y; float r22648 = z; float r22649 = fma(r22646, r22647, r22648); float r22650 = 1; float r22651 = r22646 * r22647; float r22652 = r22651 + r22648; float r22653 = r22650 + r22652; float r22654 = r22649 - r22653; return r22654; } double f_id(double x, double y, double z) { double r22655 = x; double r22656 = y; double r22657 = z; double r22658 = fma(r22655, r22656, r22657); double r22659 = 1; double r22660 = r22655 * r22656; double r22661 = r22660 + r22657; double r22662 = r22659 + r22661; double r22663 = r22658 - r22662; return r22663; } double f_of(float x, float y, float z) { float r22664 = x; float r22665 = y; float r22666 = z; float r22667 = fma(r22664, r22665, r22666); float r22668 = r22665 * r22664; float r22669 = r22667 - r22668; float r22670 = r22669 - r22666; float r22671 = 1; float r22672 = r22670 - r22671; float r22673 = exp(r22672); float r22674 = log(r22673); float r22675 = -1.0036083895155425; bool r22676 = r22674 <= r22675; float r22677 = r22667 - r22666; float r22678 = r22677 - r22668; float r22679 = r22678 - r22671; float r22680 = exp(r22679); float r22681 = log(r22680); float r22682 = exp(r22681); float r22683 = log(r22682); float r22684 = -1.0; bool r22685 = r22674 <= r22684; float r22686 = r22685 ? r22674 : r22683; float r22687 = r22676 ? r22683 : r22686; return r22687; } double f_od(double x, double y, double z) { double r22688 = x; double r22689 = y; double r22690 = z; double r22691 = fma(r22688, r22689, r22690); double r22692 = r22689 * r22688; double r22693 = r22691 - r22692; double r22694 = r22693 - r22690; double r22695 = 1; double r22696 = r22694 - r22695; double r22697 = exp(r22696); double r22698 = log(r22697); double r22699 = -1.0036083895155425; bool r22700 = r22698 <= r22699; double r22701 = r22691 - r22690; double r22702 = r22701 - r22692; double r22703 = r22702 - r22695; double r22704 = exp(r22703); double r22705 = log(r22704); double r22706 = exp(r22705); double r22707 = log(r22706); double r22708 = -1.0; bool r22709 = r22698 <= r22708; double r22710 = r22709 ? r22698 : r22707; double r22711 = r22700 ? r22707 : r22710; return r22711; } 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 r22712, r22713, r22714, r22715, r22716, r22717, r22718, r22719, r22720; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r22712); mpfr_init(r22713); mpfr_init(r22714); mpfr_init(r22715); mpfr_init_set_str(r22716, "1", 10, MPFR_RNDN); mpfr_init(r22717); mpfr_init(r22718); mpfr_init(r22719); mpfr_init(r22720); } double f_im(double x, double y, double z) { mpfr_set_d(r22712, x, MPFR_RNDN); mpfr_set_d(r22713, y, MPFR_RNDN); mpfr_set_d(r22714, z, MPFR_RNDN); mpfr_fma(r22715, r22712, r22713, r22714, MPFR_RNDN); ; mpfr_mul(r22717, r22712, r22713, MPFR_RNDN); mpfr_add(r22718, r22717, r22714, MPFR_RNDN); mpfr_add(r22719, r22716, r22718, MPFR_RNDN); mpfr_sub(r22720, r22715, r22719, MPFR_RNDN); return mpfr_get_d(r22720, MPFR_RNDN); } static mpfr_t r22721, r22722, r22723, r22724, r22725, r22726, r22727, r22728, r22729, r22730, r22731, r22732, r22733, r22734, r22735, r22736, r22737, r22738, r22739, r22740, r22741, r22742, r22743, r22744; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r22721); mpfr_init(r22722); mpfr_init(r22723); mpfr_init(r22724); mpfr_init(r22725); mpfr_init(r22726); mpfr_init(r22727); mpfr_init_set_str(r22728, "1", 10, MPFR_RNDN); mpfr_init(r22729); mpfr_init(r22730); mpfr_init(r22731); mpfr_init_set_str(r22732, "-1.0036083895155425", 10, MPFR_RNDN); mpfr_init(r22733); mpfr_init(r22734); mpfr_init(r22735); mpfr_init(r22736); mpfr_init(r22737); mpfr_init(r22738); mpfr_init(r22739); mpfr_init(r22740); mpfr_init_set_str(r22741, "-1.0", 10, MPFR_RNDN); mpfr_init(r22742); mpfr_init(r22743); mpfr_init(r22744); } double f_fm(double x, double y, double z) { mpfr_set_d(r22721, x, MPFR_RNDN); mpfr_set_d(r22722, y, MPFR_RNDN); mpfr_set_d(r22723, z, MPFR_RNDN); mpfr_fma(r22724, r22721, r22722, r22723, MPFR_RNDN); mpfr_mul(r22725, r22722, r22721, MPFR_RNDN); mpfr_sub(r22726, r22724, r22725, MPFR_RNDN); mpfr_sub(r22727, r22726, r22723, MPFR_RNDN); ; mpfr_sub(r22729, r22727, r22728, MPFR_RNDN); mpfr_exp(r22730, r22729, MPFR_RNDN); mpfr_log(r22731, r22730, MPFR_RNDN); ; mpfr_set_si(r22733, mpfr_cmp(r22731, r22732) <= 0, MPFR_RNDN); mpfr_sub(r22734, r22724, r22723, MPFR_RNDN); mpfr_sub(r22735, r22734, r22725, MPFR_RNDN); mpfr_sub(r22736, r22735, r22728, MPFR_RNDN); mpfr_exp(r22737, r22736, MPFR_RNDN); mpfr_log(r22738, r22737, MPFR_RNDN); mpfr_exp(r22739, r22738, MPFR_RNDN); mpfr_log(r22740, r22739, MPFR_RNDN); ; mpfr_set_si(r22742, mpfr_cmp(r22731, r22741) <= 0, MPFR_RNDN); if (mpfr_get_si(r22742, MPFR_RNDN)) { mpfr_set(r22743, r22731, MPFR_RNDN); } else { mpfr_set(r22743, r22740, MPFR_RNDN); }; if (mpfr_get_si(r22733, MPFR_RNDN)) { mpfr_set(r22744, r22740, MPFR_RNDN); } else { mpfr_set(r22744, r22743, MPFR_RNDN); }; return mpfr_get_d(r22744, MPFR_RNDN); } static mpfr_t r22745, r22746, r22747, r22748, r22749, r22750, r22751, r22752, r22753, r22754, r22755, r22756, r22757, r22758, r22759, r22760, r22761, r22762, r22763, r22764, r22765, r22766, r22767, r22768; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r22745); mpfr_init(r22746); mpfr_init(r22747); mpfr_init(r22748); mpfr_init(r22749); mpfr_init(r22750); mpfr_init(r22751); mpfr_init_set_str(r22752, "1", 10, MPFR_RNDN); mpfr_init(r22753); mpfr_init(r22754); mpfr_init(r22755); mpfr_init_set_str(r22756, "-1.0036083895155425", 10, MPFR_RNDN); mpfr_init(r22757); mpfr_init(r22758); mpfr_init(r22759); mpfr_init(r22760); mpfr_init(r22761); mpfr_init(r22762); mpfr_init(r22763); mpfr_init(r22764); mpfr_init_set_str(r22765, "-1.0", 10, MPFR_RNDN); mpfr_init(r22766); mpfr_init(r22767); mpfr_init(r22768); } double f_dm(double x, double y, double z) { mpfr_set_d(r22745, x, MPFR_RNDN); mpfr_set_d(r22746, y, MPFR_RNDN); mpfr_set_d(r22747, z, MPFR_RNDN); mpfr_fma(r22748, r22745, r22746, r22747, MPFR_RNDN); mpfr_mul(r22749, r22746, r22745, MPFR_RNDN); mpfr_sub(r22750, r22748, r22749, MPFR_RNDN); mpfr_sub(r22751, r22750, r22747, MPFR_RNDN); ; mpfr_sub(r22753, r22751, r22752, MPFR_RNDN); mpfr_exp(r22754, r22753, MPFR_RNDN); mpfr_log(r22755, r22754, MPFR_RNDN); ; mpfr_set_si(r22757, mpfr_cmp(r22755, r22756) <= 0, MPFR_RNDN); mpfr_sub(r22758, r22748, r22747, MPFR_RNDN); mpfr_sub(r22759, r22758, r22749, MPFR_RNDN); mpfr_sub(r22760, r22759, r22752, MPFR_RNDN); mpfr_exp(r22761, r22760, MPFR_RNDN); mpfr_log(r22762, r22761, MPFR_RNDN); mpfr_exp(r22763, r22762, MPFR_RNDN); mpfr_log(r22764, r22763, MPFR_RNDN); ; mpfr_set_si(r22766, mpfr_cmp(r22755, r22765) <= 0, MPFR_RNDN); if (mpfr_get_si(r22766, MPFR_RNDN)) { mpfr_set(r22767, r22755, MPFR_RNDN); } else { mpfr_set(r22767, r22764, MPFR_RNDN); }; if (mpfr_get_si(r22757, MPFR_RNDN)) { mpfr_set(r22768, r22764, MPFR_RNDN); } else { mpfr_set(r22768, r22767, MPFR_RNDN); }; return mpfr_get_d(r22768, MPFR_RNDN); }