#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 r22615 = x; float r22616 = y; float r22617 = r22615 * r22616; float r22618 = z; float r22619 = t; float r22620 = r22618 * r22619; float r22621 = r22617 + r22620; float r22622 = a; float r22623 = b; float r22624 = r22622 * r22623; float r22625 = r22621 + r22624; float r22626 = c; float r22627 = i; float r22628 = r22626 * r22627; float r22629 = r22625 + r22628; return r22629; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r22630 = x; double r22631 = y; double r22632 = r22630 * r22631; double r22633 = z; double r22634 = t; double r22635 = r22633 * r22634; double r22636 = r22632 + r22635; double r22637 = a; double r22638 = b; double r22639 = r22637 * r22638; double r22640 = r22636 + r22639; double r22641 = c; double r22642 = i; double r22643 = r22641 * r22642; double r22644 = r22640 + r22643; return r22644; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r22645 = x; float r22646 = y; float r22647 = r22645 * r22646; float r22648 = z; float r22649 = t; float r22650 = r22648 * r22649; float r22651 = r22647 + r22650; float r22652 = a; float r22653 = b; float r22654 = r22652 * r22653; float r22655 = r22651 + r22654; float r22656 = c; float r22657 = i; float r22658 = r22656 * r22657; float r22659 = r22655 + r22658; return r22659; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r22660 = x; double r22661 = y; double r22662 = r22660 * r22661; double r22663 = z; double r22664 = t; double r22665 = r22663 * r22664; double r22666 = r22662 + r22665; double r22667 = a; double r22668 = b; double r22669 = r22667 * r22668; double r22670 = r22666 + r22669; double r22671 = c; double r22672 = i; double r22673 = r22671 * r22672; double r22674 = r22670 + r22673; return r22674; } 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 r22675, r22676, r22677, r22678, r22679, r22680, r22681, r22682, r22683, r22684, r22685, r22686, r22687, r22688, r22689; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22675); mpfr_init(r22676); mpfr_init(r22677); mpfr_init(r22678); mpfr_init(r22679); mpfr_init(r22680); mpfr_init(r22681); mpfr_init(r22682); mpfr_init(r22683); mpfr_init(r22684); mpfr_init(r22685); mpfr_init(r22686); mpfr_init(r22687); mpfr_init(r22688); mpfr_init(r22689); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22675, x, MPFR_RNDN); mpfr_set_d(r22676, y, MPFR_RNDN); mpfr_mul(r22677, r22675, r22676, MPFR_RNDN); mpfr_set_d(r22678, z, MPFR_RNDN); mpfr_set_d(r22679, t, MPFR_RNDN); mpfr_mul(r22680, r22678, r22679, MPFR_RNDN); mpfr_add(r22681, r22677, r22680, MPFR_RNDN); mpfr_set_d(r22682, a, MPFR_RNDN); mpfr_set_d(r22683, b, MPFR_RNDN); mpfr_mul(r22684, r22682, r22683, MPFR_RNDN); mpfr_add(r22685, r22681, r22684, MPFR_RNDN); mpfr_set_d(r22686, c, MPFR_RNDN); mpfr_set_d(r22687, i, MPFR_RNDN); mpfr_mul(r22688, r22686, r22687, MPFR_RNDN); mpfr_add(r22689, r22685, r22688, MPFR_RNDN); return mpfr_get_d(r22689, MPFR_RNDN); } static mpfr_t r22690, r22691, r22692, r22693, r22694, r22695, r22696, r22697, r22698, r22699, r22700, r22701, r22702, r22703, r22704; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22690); mpfr_init(r22691); mpfr_init(r22692); mpfr_init(r22693); mpfr_init(r22694); mpfr_init(r22695); mpfr_init(r22696); mpfr_init(r22697); mpfr_init(r22698); mpfr_init(r22699); mpfr_init(r22700); mpfr_init(r22701); mpfr_init(r22702); mpfr_init(r22703); mpfr_init(r22704); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22690, x, MPFR_RNDN); mpfr_set_d(r22691, y, MPFR_RNDN); mpfr_mul(r22692, r22690, r22691, MPFR_RNDN); mpfr_set_d(r22693, z, MPFR_RNDN); mpfr_set_d(r22694, t, MPFR_RNDN); mpfr_mul(r22695, r22693, r22694, MPFR_RNDN); mpfr_add(r22696, r22692, r22695, MPFR_RNDN); mpfr_set_d(r22697, a, MPFR_RNDN); mpfr_set_d(r22698, b, MPFR_RNDN); mpfr_mul(r22699, r22697, r22698, MPFR_RNDN); mpfr_add(r22700, r22696, r22699, MPFR_RNDN); mpfr_set_d(r22701, c, MPFR_RNDN); mpfr_set_d(r22702, i, MPFR_RNDN); mpfr_mul(r22703, r22701, r22702, MPFR_RNDN); mpfr_add(r22704, r22700, r22703, MPFR_RNDN); return mpfr_get_d(r22704, MPFR_RNDN); } static mpfr_t r22705, r22706, r22707, r22708, r22709, r22710, r22711, r22712, r22713, r22714, r22715, r22716, r22717, r22718, r22719; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22705); mpfr_init(r22706); mpfr_init(r22707); mpfr_init(r22708); mpfr_init(r22709); mpfr_init(r22710); mpfr_init(r22711); mpfr_init(r22712); mpfr_init(r22713); mpfr_init(r22714); mpfr_init(r22715); mpfr_init(r22716); mpfr_init(r22717); mpfr_init(r22718); mpfr_init(r22719); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22705, x, MPFR_RNDN); mpfr_set_d(r22706, y, MPFR_RNDN); mpfr_mul(r22707, r22705, r22706, MPFR_RNDN); mpfr_set_d(r22708, z, MPFR_RNDN); mpfr_set_d(r22709, t, MPFR_RNDN); mpfr_mul(r22710, r22708, r22709, MPFR_RNDN); mpfr_add(r22711, r22707, r22710, MPFR_RNDN); mpfr_set_d(r22712, a, MPFR_RNDN); mpfr_set_d(r22713, b, MPFR_RNDN); mpfr_mul(r22714, r22712, r22713, MPFR_RNDN); mpfr_add(r22715, r22711, r22714, MPFR_RNDN); mpfr_set_d(r22716, c, MPFR_RNDN); mpfr_set_d(r22717, i, MPFR_RNDN); mpfr_mul(r22718, r22716, r22717, MPFR_RNDN); mpfr_add(r22719, r22715, r22718, MPFR_RNDN); return mpfr_get_d(r22719, MPFR_RNDN); }