#include #include #include #include #include char *name = "simple fma test"; double f_if(float x, float y, float z) { float r20591 = x; float r20592 = y; float r20593 = z; float r20594 = fma(r20591, r20592, r20593); float r20595 = 1; float r20596 = r20591 * r20592; float r20597 = r20596 + r20593; float r20598 = r20595 + r20597; float r20599 = r20594 - r20598; return r20599; } double f_id(double x, double y, double z) { double r20600 = x; double r20601 = y; double r20602 = z; double r20603 = fma(r20600, r20601, r20602); double r20604 = 1; double r20605 = r20600 * r20601; double r20606 = r20605 + r20602; double r20607 = r20604 + r20606; double r20608 = r20603 - r20607; return r20608; } double f_of(float x, float y, float z) { float r20609 = x; float r20610 = y; float r20611 = z; float r20612 = fma(r20609, r20610, r20611); float r20613 = 1; float r20614 = r20609 * r20610; float r20615 = r20614 + r20611; float r20616 = r20614 - r20611; float r20617 = cbrt(r20616); float r20618 = r20617 * r20617; float r20619 = r20615 / r20618; float r20620 = r20616 / r20617; float r20621 = r20619 * r20620; float r20622 = r20613 + r20621; float r20623 = r20612 - r20622; float r20624 = cbrt(r20623); float r20625 = r20624 * r20624; float r20626 = r20625 * r20624; return r20626; } double f_od(double x, double y, double z) { double r20627 = x; double r20628 = y; double r20629 = z; double r20630 = fma(r20627, r20628, r20629); double r20631 = 1; double r20632 = r20627 * r20628; double r20633 = r20632 + r20629; double r20634 = r20632 - r20629; double r20635 = cbrt(r20634); double r20636 = r20635 * r20635; double r20637 = r20633 / r20636; double r20638 = r20634 / r20635; double r20639 = r20637 * r20638; double r20640 = r20631 + r20639; double r20641 = r20630 - r20640; double r20642 = cbrt(r20641); double r20643 = r20642 * r20642; double r20644 = r20643 * r20642; return r20644; } 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 r20645, r20646, r20647, r20648, r20649, r20650, r20651, r20652, r20653; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r20645); mpfr_init(r20646); mpfr_init(r20647); mpfr_init(r20648); mpfr_init_set_str(r20649, "1", 10, MPFR_RNDN); mpfr_init(r20650); mpfr_init(r20651); mpfr_init(r20652); mpfr_init(r20653); } double f_im(double x, double y, double z) { mpfr_set_d(r20645, x, MPFR_RNDN); mpfr_set_d(r20646, y, MPFR_RNDN); mpfr_set_d(r20647, z, MPFR_RNDN); mpfr_fma(r20648, r20645, r20646, r20647, MPFR_RNDN); ; mpfr_mul(r20650, r20645, r20646, MPFR_RNDN); mpfr_add(r20651, r20650, r20647, MPFR_RNDN); mpfr_add(r20652, r20649, r20651, MPFR_RNDN); mpfr_sub(r20653, r20648, r20652, MPFR_RNDN); return mpfr_get_d(r20653, MPFR_RNDN); } static mpfr_t r20654, r20655, r20656, r20657, r20658, r20659, r20660, r20661, r20662, r20663, r20664, r20665, r20666, r20667, r20668, r20669, r20670, r20671; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r20654); mpfr_init(r20655); mpfr_init(r20656); mpfr_init(r20657); mpfr_init_set_str(r20658, "1", 10, MPFR_RNDN); mpfr_init(r20659); mpfr_init(r20660); mpfr_init(r20661); mpfr_init(r20662); mpfr_init(r20663); mpfr_init(r20664); mpfr_init(r20665); mpfr_init(r20666); mpfr_init(r20667); mpfr_init(r20668); mpfr_init(r20669); mpfr_init(r20670); mpfr_init(r20671); } double f_fm(double x, double y, double z) { mpfr_set_d(r20654, x, MPFR_RNDN); mpfr_set_d(r20655, y, MPFR_RNDN); mpfr_set_d(r20656, z, MPFR_RNDN); mpfr_fma(r20657, r20654, r20655, r20656, MPFR_RNDN); ; mpfr_mul(r20659, r20654, r20655, MPFR_RNDN); mpfr_add(r20660, r20659, r20656, MPFR_RNDN); mpfr_sub(r20661, r20659, r20656, MPFR_RNDN); mpfr_cbrt(r20662, r20661, MPFR_RNDN); mpfr_mul(r20663, r20662, r20662, MPFR_RNDN); mpfr_div(r20664, r20660, r20663, MPFR_RNDN); mpfr_div(r20665, r20661, r20662, MPFR_RNDN); mpfr_mul(r20666, r20664, r20665, MPFR_RNDN); mpfr_add(r20667, r20658, r20666, MPFR_RNDN); mpfr_sub(r20668, r20657, r20667, MPFR_RNDN); mpfr_cbrt(r20669, r20668, MPFR_RNDN); mpfr_mul(r20670, r20669, r20669, MPFR_RNDN); mpfr_mul(r20671, r20670, r20669, MPFR_RNDN); return mpfr_get_d(r20671, MPFR_RNDN); } static mpfr_t r20672, r20673, r20674, r20675, r20676, r20677, r20678, r20679, r20680, r20681, r20682, r20683, r20684, r20685, r20686, r20687, r20688, r20689; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r20672); mpfr_init(r20673); mpfr_init(r20674); mpfr_init(r20675); mpfr_init_set_str(r20676, "1", 10, MPFR_RNDN); mpfr_init(r20677); mpfr_init(r20678); mpfr_init(r20679); mpfr_init(r20680); mpfr_init(r20681); mpfr_init(r20682); mpfr_init(r20683); mpfr_init(r20684); mpfr_init(r20685); mpfr_init(r20686); mpfr_init(r20687); mpfr_init(r20688); mpfr_init(r20689); } double f_dm(double x, double y, double z) { mpfr_set_d(r20672, x, MPFR_RNDN); mpfr_set_d(r20673, y, MPFR_RNDN); mpfr_set_d(r20674, z, MPFR_RNDN); mpfr_fma(r20675, r20672, r20673, r20674, MPFR_RNDN); ; mpfr_mul(r20677, r20672, r20673, MPFR_RNDN); mpfr_add(r20678, r20677, r20674, MPFR_RNDN); mpfr_sub(r20679, r20677, r20674, MPFR_RNDN); mpfr_cbrt(r20680, r20679, MPFR_RNDN); mpfr_mul(r20681, r20680, r20680, MPFR_RNDN); mpfr_div(r20682, r20678, r20681, MPFR_RNDN); mpfr_div(r20683, r20679, r20680, MPFR_RNDN); mpfr_mul(r20684, r20682, r20683, MPFR_RNDN); mpfr_add(r20685, r20676, r20684, MPFR_RNDN); mpfr_sub(r20686, r20675, r20685, MPFR_RNDN); mpfr_cbrt(r20687, r20686, MPFR_RNDN); mpfr_mul(r20688, r20687, r20687, MPFR_RNDN); mpfr_mul(r20689, r20688, r20687, MPFR_RNDN); return mpfr_get_d(r20689, MPFR_RNDN); }