#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r25631 = x; float r25632 = y; float r25633 = z; float r25634 = fma(r25631, r25632, r25633); float r25635 = 1; float r25636 = r25631 * r25632; float r25637 = r25636 + r25633; float r25638 = r25635 + r25637; float r25639 = r25634 - r25638; return r25639; } double f_id(double x, double y, double z) { double r25640 = x; double r25641 = y; double r25642 = z; double r25643 = fma(r25640, r25641, r25642); double r25644 = 1; double r25645 = r25640 * r25641; double r25646 = r25645 + r25642; double r25647 = r25644 + r25646; double r25648 = r25643 - r25647; return r25648; } double f_of(float x, float y, float z) { float r25649 = x; float r25650 = y; float r25651 = z; float r25652 = fma(r25649, r25650, r25651); float r25653 = r25650 * r25649; float r25654 = r25652 - r25653; float r25655 = r25654 - r25651; float r25656 = 1; float r25657 = r25655 - r25656; float r25658 = exp(r25657); float r25659 = log(r25658); float r25660 = -1.015625; bool r25661 = r25659 <= r25660; float r25662 = r25652 - r25651; float r25663 = r25662 - r25653; float r25664 = r25663 - r25656; float r25665 = exp(r25664); float r25666 = log(r25665); float r25667 = exp(r25666); float r25668 = log(r25667); float r25669 = -0.9999999999567729; bool r25670 = r25659 <= r25669; float r25671 = r25670 ? r25659 : r25668; float r25672 = r25661 ? r25668 : r25671; return r25672; } double f_od(double x, double y, double z) { double r25673 = x; double r25674 = y; double r25675 = z; double r25676 = fma(r25673, r25674, r25675); double r25677 = r25674 * r25673; double r25678 = r25676 - r25677; double r25679 = r25678 - r25675; double r25680 = 1; double r25681 = r25679 - r25680; double r25682 = exp(r25681); double r25683 = log(r25682); double r25684 = -1.015625; bool r25685 = r25683 <= r25684; double r25686 = r25676 - r25675; double r25687 = r25686 - r25677; double r25688 = r25687 - r25680; double r25689 = exp(r25688); double r25690 = log(r25689); double r25691 = exp(r25690); double r25692 = log(r25691); double r25693 = -0.9999999999567729; bool r25694 = r25683 <= r25693; double r25695 = r25694 ? r25683 : r25692; double r25696 = r25685 ? r25692 : r25695; return r25696; } 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 r25697, r25698, r25699, r25700, r25701, r25702, r25703, r25704, r25705; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r25697); mpfr_init(r25698); mpfr_init(r25699); mpfr_init(r25700); mpfr_init_set_str(r25701, "1", 10, MPFR_RNDN); mpfr_init(r25702); mpfr_init(r25703); mpfr_init(r25704); mpfr_init(r25705); } double f_im(double x, double y, double z) { mpfr_set_d(r25697, x, MPFR_RNDN); mpfr_set_d(r25698, y, MPFR_RNDN); mpfr_set_d(r25699, z, MPFR_RNDN); mpfr_fma(r25700, r25697, r25698, r25699, MPFR_RNDN); ; mpfr_mul(r25702, r25697, r25698, MPFR_RNDN); mpfr_add(r25703, r25702, r25699, MPFR_RNDN); mpfr_add(r25704, r25701, r25703, MPFR_RNDN); mpfr_sub(r25705, r25700, r25704, MPFR_RNDN); return mpfr_get_d(r25705, MPFR_RNDN); } static mpfr_t r25706, r25707, r25708, r25709, r25710, r25711, r25712, r25713, r25714, r25715, r25716, r25717, r25718, r25719, r25720, r25721, r25722, r25723, r25724, r25725, r25726, r25727, r25728, r25729; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r25706); mpfr_init(r25707); mpfr_init(r25708); mpfr_init(r25709); mpfr_init(r25710); mpfr_init(r25711); mpfr_init(r25712); mpfr_init_set_str(r25713, "1", 10, MPFR_RNDN); mpfr_init(r25714); mpfr_init(r25715); mpfr_init(r25716); mpfr_init_set_str(r25717, "-1.015625", 10, MPFR_RNDN); mpfr_init(r25718); mpfr_init(r25719); mpfr_init(r25720); mpfr_init(r25721); mpfr_init(r25722); mpfr_init(r25723); mpfr_init(r25724); mpfr_init(r25725); mpfr_init_set_str(r25726, "-0.9999999999567729", 10, MPFR_RNDN); mpfr_init(r25727); mpfr_init(r25728); mpfr_init(r25729); } double f_fm(double x, double y, double z) { mpfr_set_d(r25706, x, MPFR_RNDN); mpfr_set_d(r25707, y, MPFR_RNDN); mpfr_set_d(r25708, z, MPFR_RNDN); mpfr_fma(r25709, r25706, r25707, r25708, MPFR_RNDN); mpfr_mul(r25710, r25707, r25706, MPFR_RNDN); mpfr_sub(r25711, r25709, r25710, MPFR_RNDN); mpfr_sub(r25712, r25711, r25708, MPFR_RNDN); ; mpfr_sub(r25714, r25712, r25713, MPFR_RNDN); mpfr_exp(r25715, r25714, MPFR_RNDN); mpfr_log(r25716, r25715, MPFR_RNDN); ; mpfr_set_si(r25718, mpfr_cmp(r25716, r25717) <= 0, MPFR_RNDN); mpfr_sub(r25719, r25709, r25708, MPFR_RNDN); mpfr_sub(r25720, r25719, r25710, MPFR_RNDN); mpfr_sub(r25721, r25720, r25713, MPFR_RNDN); mpfr_exp(r25722, r25721, MPFR_RNDN); mpfr_log(r25723, r25722, MPFR_RNDN); mpfr_exp(r25724, r25723, MPFR_RNDN); mpfr_log(r25725, r25724, MPFR_RNDN); ; mpfr_set_si(r25727, mpfr_cmp(r25716, r25726) <= 0, MPFR_RNDN); if (mpfr_get_si(r25727, MPFR_RNDN)) { mpfr_set(r25728, r25716, MPFR_RNDN); } else { mpfr_set(r25728, r25725, MPFR_RNDN); }; if (mpfr_get_si(r25718, MPFR_RNDN)) { mpfr_set(r25729, r25725, MPFR_RNDN); } else { mpfr_set(r25729, r25728, MPFR_RNDN); }; return mpfr_get_d(r25729, MPFR_RNDN); } static mpfr_t r25730, r25731, r25732, r25733, r25734, r25735, r25736, r25737, r25738, r25739, r25740, r25741, r25742, r25743, r25744, r25745, r25746, r25747, r25748, r25749, r25750, r25751, r25752, r25753; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r25730); mpfr_init(r25731); mpfr_init(r25732); mpfr_init(r25733); mpfr_init(r25734); mpfr_init(r25735); mpfr_init(r25736); mpfr_init_set_str(r25737, "1", 10, MPFR_RNDN); mpfr_init(r25738); mpfr_init(r25739); mpfr_init(r25740); mpfr_init_set_str(r25741, "-1.015625", 10, MPFR_RNDN); mpfr_init(r25742); mpfr_init(r25743); mpfr_init(r25744); mpfr_init(r25745); mpfr_init(r25746); mpfr_init(r25747); mpfr_init(r25748); mpfr_init(r25749); mpfr_init_set_str(r25750, "-0.9999999999567729", 10, MPFR_RNDN); mpfr_init(r25751); mpfr_init(r25752); mpfr_init(r25753); } double f_dm(double x, double y, double z) { mpfr_set_d(r25730, x, MPFR_RNDN); mpfr_set_d(r25731, y, MPFR_RNDN); mpfr_set_d(r25732, z, MPFR_RNDN); mpfr_fma(r25733, r25730, r25731, r25732, MPFR_RNDN); mpfr_mul(r25734, r25731, r25730, MPFR_RNDN); mpfr_sub(r25735, r25733, r25734, MPFR_RNDN); mpfr_sub(r25736, r25735, r25732, MPFR_RNDN); ; mpfr_sub(r25738, r25736, r25737, MPFR_RNDN); mpfr_exp(r25739, r25738, MPFR_RNDN); mpfr_log(r25740, r25739, MPFR_RNDN); ; mpfr_set_si(r25742, mpfr_cmp(r25740, r25741) <= 0, MPFR_RNDN); mpfr_sub(r25743, r25733, r25732, MPFR_RNDN); mpfr_sub(r25744, r25743, r25734, MPFR_RNDN); mpfr_sub(r25745, r25744, r25737, MPFR_RNDN); mpfr_exp(r25746, r25745, MPFR_RNDN); mpfr_log(r25747, r25746, MPFR_RNDN); mpfr_exp(r25748, r25747, MPFR_RNDN); mpfr_log(r25749, r25748, MPFR_RNDN); ; mpfr_set_si(r25751, mpfr_cmp(r25740, r25750) <= 0, MPFR_RNDN); if (mpfr_get_si(r25751, MPFR_RNDN)) { mpfr_set(r25752, r25740, MPFR_RNDN); } else { mpfr_set(r25752, r25749, MPFR_RNDN); }; if (mpfr_get_si(r25742, MPFR_RNDN)) { mpfr_set(r25753, r25749, MPFR_RNDN); } else { mpfr_set(r25753, r25752, MPFR_RNDN); }; return mpfr_get_d(r25753, MPFR_RNDN); }