#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 r16546 = x; float r16547 = y; float r16548 = r16546 * r16547; float r16549 = z; float r16550 = t; float r16551 = r16549 * r16550; float r16552 = r16548 + r16551; float r16553 = a; float r16554 = b; float r16555 = r16553 * r16554; float r16556 = r16552 + r16555; float r16557 = c; float r16558 = i; float r16559 = r16557 * r16558; float r16560 = r16556 + r16559; return r16560; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r16561 = x; double r16562 = y; double r16563 = r16561 * r16562; double r16564 = z; double r16565 = t; double r16566 = r16564 * r16565; double r16567 = r16563 + r16566; double r16568 = a; double r16569 = b; double r16570 = r16568 * r16569; double r16571 = r16567 + r16570; double r16572 = c; double r16573 = i; double r16574 = r16572 * r16573; double r16575 = r16571 + r16574; return r16575; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r16576 = z; float r16577 = t; float r16578 = y; float r16579 = x; float r16580 = i; float r16581 = c; float r16582 = b; float r16583 = a; float r16584 = r16582 * r16583; float r16585 = fma(r16580, r16581, r16584); float r16586 = fma(r16578, r16579, r16585); float r16587 = fma(r16576, r16577, r16586); return r16587; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r16588 = z; double r16589 = t; double r16590 = y; double r16591 = x; double r16592 = i; double r16593 = c; double r16594 = b; double r16595 = a; double r16596 = r16594 * r16595; double r16597 = fma(r16592, r16593, r16596); double r16598 = fma(r16590, r16591, r16597); double r16599 = fma(r16588, r16589, r16598); return r16599; } 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 r16600, r16601, r16602, r16603, r16604, r16605, r16606, r16607, r16608, r16609, r16610, r16611, r16612, r16613, r16614; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16600); mpfr_init(r16601); mpfr_init(r16602); mpfr_init(r16603); mpfr_init(r16604); mpfr_init(r16605); mpfr_init(r16606); mpfr_init(r16607); mpfr_init(r16608); mpfr_init(r16609); mpfr_init(r16610); mpfr_init(r16611); mpfr_init(r16612); mpfr_init(r16613); mpfr_init(r16614); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r16600, x, MPFR_RNDN); mpfr_set_d(r16601, y, MPFR_RNDN); mpfr_mul(r16602, r16600, r16601, MPFR_RNDN); mpfr_set_d(r16603, z, MPFR_RNDN); mpfr_set_d(r16604, t, MPFR_RNDN); mpfr_mul(r16605, r16603, r16604, MPFR_RNDN); mpfr_add(r16606, r16602, r16605, MPFR_RNDN); mpfr_set_d(r16607, a, MPFR_RNDN); mpfr_set_d(r16608, b, MPFR_RNDN); mpfr_mul(r16609, r16607, r16608, MPFR_RNDN); mpfr_add(r16610, r16606, r16609, MPFR_RNDN); mpfr_set_d(r16611, c, MPFR_RNDN); mpfr_set_d(r16612, i, MPFR_RNDN); mpfr_mul(r16613, r16611, r16612, MPFR_RNDN); mpfr_add(r16614, r16610, r16613, MPFR_RNDN); return mpfr_get_d(r16614, MPFR_RNDN); } static mpfr_t r16615, r16616, r16617, r16618, r16619, r16620, r16621, r16622, r16623, r16624, r16625, r16626; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16615); mpfr_init(r16616); mpfr_init(r16617); mpfr_init(r16618); mpfr_init(r16619); mpfr_init(r16620); mpfr_init(r16621); mpfr_init(r16622); mpfr_init(r16623); mpfr_init(r16624); mpfr_init(r16625); mpfr_init(r16626); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r16615, z, MPFR_RNDN); mpfr_set_d(r16616, t, MPFR_RNDN); mpfr_set_d(r16617, y, MPFR_RNDN); mpfr_set_d(r16618, x, MPFR_RNDN); mpfr_set_d(r16619, i, MPFR_RNDN); mpfr_set_d(r16620, c, MPFR_RNDN); mpfr_set_d(r16621, b, MPFR_RNDN); mpfr_set_d(r16622, a, MPFR_RNDN); mpfr_mul(r16623, r16621, r16622, MPFR_RNDN); mpfr_fma(r16624, r16619, r16620, r16623, MPFR_RNDN); mpfr_fma(r16625, r16617, r16618, r16624, MPFR_RNDN); mpfr_fma(r16626, r16615, r16616, r16625, MPFR_RNDN); return mpfr_get_d(r16626, MPFR_RNDN); } static mpfr_t r16627, r16628, r16629, r16630, r16631, r16632, r16633, r16634, r16635, r16636, r16637, r16638; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16627); mpfr_init(r16628); mpfr_init(r16629); mpfr_init(r16630); mpfr_init(r16631); mpfr_init(r16632); mpfr_init(r16633); mpfr_init(r16634); mpfr_init(r16635); mpfr_init(r16636); mpfr_init(r16637); mpfr_init(r16638); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r16627, z, MPFR_RNDN); mpfr_set_d(r16628, t, MPFR_RNDN); mpfr_set_d(r16629, y, MPFR_RNDN); mpfr_set_d(r16630, x, MPFR_RNDN); mpfr_set_d(r16631, i, MPFR_RNDN); mpfr_set_d(r16632, c, MPFR_RNDN); mpfr_set_d(r16633, b, MPFR_RNDN); mpfr_set_d(r16634, a, MPFR_RNDN); mpfr_mul(r16635, r16633, r16634, MPFR_RNDN); mpfr_fma(r16636, r16631, r16632, r16635, MPFR_RNDN); mpfr_fma(r16637, r16629, r16630, r16636, MPFR_RNDN); mpfr_fma(r16638, r16627, r16628, r16637, MPFR_RNDN); return mpfr_get_d(r16638, MPFR_RNDN); }