#include #include #include #include #include char *name = "Linear.V4:$cdot from linear-1.19.1.3"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i) { float r23689 = x; float r23690 = y; float r23691 = r23689 * r23690; float r23692 = z; float r23693 = t; float r23694 = r23692 * r23693; float r23695 = r23691 + r23694; float r23696 = a; float r23697 = b; float r23698 = r23696 * r23697; float r23699 = r23695 + r23698; float r23700 = c; float r23701 = i; float r23702 = r23700 * r23701; float r23703 = r23699 + r23702; return r23703; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r23704 = x; double r23705 = y; double r23706 = r23704 * r23705; double r23707 = z; double r23708 = t; double r23709 = r23707 * r23708; double r23710 = r23706 + r23709; double r23711 = a; double r23712 = b; double r23713 = r23711 * r23712; double r23714 = r23710 + r23713; double r23715 = c; double r23716 = i; double r23717 = r23715 * r23716; double r23718 = r23714 + r23717; return r23718; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r23719 = z; float r23720 = t; float r23721 = y; float r23722 = x; float r23723 = r23721 * r23722; float r23724 = fma(r23719, r23720, r23723); float r23725 = i; float r23726 = c; float r23727 = b; float r23728 = a; float r23729 = r23727 * r23728; float r23730 = fma(r23725, r23726, r23729); float r23731 = r23724 + r23730; return r23731; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r23732 = z; double r23733 = t; double r23734 = y; double r23735 = x; double r23736 = r23734 * r23735; double r23737 = fma(r23732, r23733, r23736); double r23738 = i; double r23739 = c; double r23740 = b; double r23741 = a; double r23742 = r23740 * r23741; double r23743 = fma(r23738, r23739, r23742); double r23744 = r23737 + r23743; return r23744; } 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 r23745, r23746, r23747, r23748, r23749, r23750, r23751, r23752, r23753, r23754, r23755, r23756, r23757, r23758, r23759; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r23745); mpfr_init(r23746); mpfr_init(r23747); mpfr_init(r23748); mpfr_init(r23749); mpfr_init(r23750); mpfr_init(r23751); mpfr_init(r23752); mpfr_init(r23753); mpfr_init(r23754); mpfr_init(r23755); mpfr_init(r23756); mpfr_init(r23757); mpfr_init(r23758); mpfr_init(r23759); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r23745, x, MPFR_RNDN); mpfr_set_d(r23746, y, MPFR_RNDN); mpfr_mul(r23747, r23745, r23746, MPFR_RNDN); mpfr_set_d(r23748, z, MPFR_RNDN); mpfr_set_d(r23749, t, MPFR_RNDN); mpfr_mul(r23750, r23748, r23749, MPFR_RNDN); mpfr_add(r23751, r23747, r23750, MPFR_RNDN); mpfr_set_d(r23752, a, MPFR_RNDN); mpfr_set_d(r23753, b, MPFR_RNDN); mpfr_mul(r23754, r23752, r23753, MPFR_RNDN); mpfr_add(r23755, r23751, r23754, MPFR_RNDN); mpfr_set_d(r23756, c, MPFR_RNDN); mpfr_set_d(r23757, i, MPFR_RNDN); mpfr_mul(r23758, r23756, r23757, MPFR_RNDN); mpfr_add(r23759, r23755, r23758, MPFR_RNDN); return mpfr_get_d(r23759, MPFR_RNDN); } static mpfr_t r23760, r23761, r23762, r23763, r23764, r23765, r23766, r23767, r23768, r23769, r23770, r23771, r23772; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r23760); mpfr_init(r23761); mpfr_init(r23762); mpfr_init(r23763); mpfr_init(r23764); mpfr_init(r23765); mpfr_init(r23766); mpfr_init(r23767); mpfr_init(r23768); mpfr_init(r23769); mpfr_init(r23770); mpfr_init(r23771); mpfr_init(r23772); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r23760, z, MPFR_RNDN); mpfr_set_d(r23761, t, MPFR_RNDN); mpfr_set_d(r23762, y, MPFR_RNDN); mpfr_set_d(r23763, x, MPFR_RNDN); mpfr_mul(r23764, r23762, r23763, MPFR_RNDN); mpfr_fma(r23765, r23760, r23761, r23764, MPFR_RNDN); mpfr_set_d(r23766, i, MPFR_RNDN); mpfr_set_d(r23767, c, MPFR_RNDN); mpfr_set_d(r23768, b, MPFR_RNDN); mpfr_set_d(r23769, a, MPFR_RNDN); mpfr_mul(r23770, r23768, r23769, MPFR_RNDN); mpfr_fma(r23771, r23766, r23767, r23770, MPFR_RNDN); mpfr_add(r23772, r23765, r23771, MPFR_RNDN); return mpfr_get_d(r23772, MPFR_RNDN); } static mpfr_t r23773, r23774, r23775, r23776, r23777, r23778, r23779, r23780, r23781, r23782, r23783, r23784, r23785; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r23773); mpfr_init(r23774); mpfr_init(r23775); mpfr_init(r23776); mpfr_init(r23777); mpfr_init(r23778); mpfr_init(r23779); mpfr_init(r23780); mpfr_init(r23781); mpfr_init(r23782); mpfr_init(r23783); mpfr_init(r23784); mpfr_init(r23785); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r23773, z, MPFR_RNDN); mpfr_set_d(r23774, t, MPFR_RNDN); mpfr_set_d(r23775, y, MPFR_RNDN); mpfr_set_d(r23776, x, MPFR_RNDN); mpfr_mul(r23777, r23775, r23776, MPFR_RNDN); mpfr_fma(r23778, r23773, r23774, r23777, MPFR_RNDN); mpfr_set_d(r23779, i, MPFR_RNDN); mpfr_set_d(r23780, c, MPFR_RNDN); mpfr_set_d(r23781, b, MPFR_RNDN); mpfr_set_d(r23782, a, MPFR_RNDN); mpfr_mul(r23783, r23781, r23782, MPFR_RNDN); mpfr_fma(r23784, r23779, r23780, r23783, MPFR_RNDN); mpfr_add(r23785, r23778, r23784, MPFR_RNDN); return mpfr_get_d(r23785, MPFR_RNDN); }