#include #include #include #include #include char *name = "fabs fraction 1"; double f_if(float x, float y, float z) { float r8613 = x; float r8614 = 4; float r8615 = r8613 + r8614; float r8616 = y; float r8617 = r8615 / r8616; float r8618 = r8613 / r8616; float r8619 = z; float r8620 = r8618 * r8619; float r8621 = r8617 - r8620; float r8622 = fabs(r8621); return r8622; } double f_id(double x, double y, double z) { double r8623 = x; double r8624 = 4; double r8625 = r8623 + r8624; double r8626 = y; double r8627 = r8625 / r8626; double r8628 = r8623 / r8626; double r8629 = z; double r8630 = r8628 * r8629; double r8631 = r8627 - r8630; double r8632 = fabs(r8631); return r8632; } double f_of(float x, float y, float z) { float r8633 = x; float r8634 = 4; float r8635 = r8633 + r8634; float r8636 = y; float r8637 = r8635 / r8636; float r8638 = z; float r8639 = r8638 / r8636; float r8640 = r8633 * r8639; float r8641 = r8637 - r8640; float r8642 = fabs(r8641); float r8643 = 1.7391963550747397e+215; bool r8644 = r8642 <= r8643; float r8645 = 1; float r8646 = r8645 / r8636; float r8647 = r8633 / r8636; float r8648 = r8647 * r8638; float r8649 = -r8648; float r8650 = fma(r8635, r8646, r8649); float r8651 = fabs(r8650); float r8652 = r8644 ? r8642 : r8651; return r8652; } double f_od(double x, double y, double z) { double r8653 = x; double r8654 = 4; double r8655 = r8653 + r8654; double r8656 = y; double r8657 = r8655 / r8656; double r8658 = z; double r8659 = r8658 / r8656; double r8660 = r8653 * r8659; double r8661 = r8657 - r8660; double r8662 = fabs(r8661); double r8663 = 1.7391963550747397e+215; bool r8664 = r8662 <= r8663; double r8665 = 1; double r8666 = r8665 / r8656; double r8667 = r8653 / r8656; double r8668 = r8667 * r8658; double r8669 = -r8668; double r8670 = fma(r8655, r8666, r8669); double r8671 = fabs(r8670); double r8672 = r8664 ? r8662 : r8671; return r8672; } 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 r8673, r8674, r8675, r8676, r8677, r8678, r8679, r8680, r8681, r8682; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r8673); mpfr_init_set_str(r8674, "4", 10, MPFR_RNDN); mpfr_init(r8675); mpfr_init(r8676); mpfr_init(r8677); mpfr_init(r8678); mpfr_init(r8679); mpfr_init(r8680); mpfr_init(r8681); mpfr_init(r8682); } double f_im(double x, double y, double z) { mpfr_set_d(r8673, x, MPFR_RNDN); ; mpfr_add(r8675, r8673, r8674, MPFR_RNDN); mpfr_set_d(r8676, y, MPFR_RNDN); mpfr_div(r8677, r8675, r8676, MPFR_RNDN); mpfr_div(r8678, r8673, r8676, MPFR_RNDN); mpfr_set_d(r8679, z, MPFR_RNDN); mpfr_mul(r8680, r8678, r8679, MPFR_RNDN); mpfr_sub(r8681, r8677, r8680, MPFR_RNDN); mpfr_abs(r8682, r8681, MPFR_RNDN); return mpfr_get_d(r8682, MPFR_RNDN); } static mpfr_t r8683, r8684, r8685, r8686, r8687, r8688, r8689, r8690, r8691, r8692, r8693, r8694, r8695, r8696, r8697, r8698, r8699, r8700, r8701, r8702; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r8683); mpfr_init_set_str(r8684, "4", 10, MPFR_RNDN); mpfr_init(r8685); mpfr_init(r8686); mpfr_init(r8687); mpfr_init(r8688); mpfr_init(r8689); mpfr_init(r8690); mpfr_init(r8691); mpfr_init(r8692); mpfr_init_set_str(r8693, "1.7391963550747397e+215", 10, MPFR_RNDN); mpfr_init(r8694); mpfr_init_set_str(r8695, "1", 10, MPFR_RNDN); mpfr_init(r8696); mpfr_init(r8697); mpfr_init(r8698); mpfr_init(r8699); mpfr_init(r8700); mpfr_init(r8701); mpfr_init(r8702); } double f_fm(double x, double y, double z) { mpfr_set_d(r8683, x, MPFR_RNDN); ; mpfr_add(r8685, r8683, r8684, MPFR_RNDN); mpfr_set_d(r8686, y, MPFR_RNDN); mpfr_div(r8687, r8685, r8686, MPFR_RNDN); mpfr_set_d(r8688, z, MPFR_RNDN); mpfr_div(r8689, r8688, r8686, MPFR_RNDN); mpfr_mul(r8690, r8683, r8689, MPFR_RNDN); mpfr_sub(r8691, r8687, r8690, MPFR_RNDN); mpfr_abs(r8692, r8691, MPFR_RNDN); ; mpfr_set_si(r8694, mpfr_cmp(r8692, r8693) <= 0, MPFR_RNDN); ; mpfr_div(r8696, r8695, r8686, MPFR_RNDN); mpfr_div(r8697, r8683, r8686, MPFR_RNDN); mpfr_mul(r8698, r8697, r8688, MPFR_RNDN); mpfr_neg(r8699, r8698, MPFR_RNDN); mpfr_fma(r8700, r8685, r8696, r8699, MPFR_RNDN); mpfr_abs(r8701, r8700, MPFR_RNDN); if (mpfr_get_si(r8694, MPFR_RNDN)) { mpfr_set(r8702, r8692, MPFR_RNDN); } else { mpfr_set(r8702, r8701, MPFR_RNDN); }; return mpfr_get_d(r8702, MPFR_RNDN); } static mpfr_t r8703, r8704, r8705, r8706, r8707, r8708, r8709, r8710, r8711, r8712, r8713, r8714, r8715, r8716, r8717, r8718, r8719, r8720, r8721, r8722; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r8703); mpfr_init_set_str(r8704, "4", 10, MPFR_RNDN); mpfr_init(r8705); mpfr_init(r8706); mpfr_init(r8707); mpfr_init(r8708); mpfr_init(r8709); mpfr_init(r8710); mpfr_init(r8711); mpfr_init(r8712); mpfr_init_set_str(r8713, "1.7391963550747397e+215", 10, MPFR_RNDN); mpfr_init(r8714); mpfr_init_set_str(r8715, "1", 10, MPFR_RNDN); mpfr_init(r8716); mpfr_init(r8717); mpfr_init(r8718); mpfr_init(r8719); mpfr_init(r8720); mpfr_init(r8721); mpfr_init(r8722); } double f_dm(double x, double y, double z) { mpfr_set_d(r8703, x, MPFR_RNDN); ; mpfr_add(r8705, r8703, r8704, MPFR_RNDN); mpfr_set_d(r8706, y, MPFR_RNDN); mpfr_div(r8707, r8705, r8706, MPFR_RNDN); mpfr_set_d(r8708, z, MPFR_RNDN); mpfr_div(r8709, r8708, r8706, MPFR_RNDN); mpfr_mul(r8710, r8703, r8709, MPFR_RNDN); mpfr_sub(r8711, r8707, r8710, MPFR_RNDN); mpfr_abs(r8712, r8711, MPFR_RNDN); ; mpfr_set_si(r8714, mpfr_cmp(r8712, r8713) <= 0, MPFR_RNDN); ; mpfr_div(r8716, r8715, r8706, MPFR_RNDN); mpfr_div(r8717, r8703, r8706, MPFR_RNDN); mpfr_mul(r8718, r8717, r8708, MPFR_RNDN); mpfr_neg(r8719, r8718, MPFR_RNDN); mpfr_fma(r8720, r8705, r8716, r8719, MPFR_RNDN); mpfr_abs(r8721, r8720, MPFR_RNDN); if (mpfr_get_si(r8714, MPFR_RNDN)) { mpfr_set(r8722, r8712, MPFR_RNDN); } else { mpfr_set(r8722, r8721, MPFR_RNDN); }; return mpfr_get_d(r8722, MPFR_RNDN); }