#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 r22438 = x; float r22439 = y; float r22440 = r22438 * r22439; float r22441 = z; float r22442 = t; float r22443 = r22441 * r22442; float r22444 = r22440 + r22443; float r22445 = a; float r22446 = b; float r22447 = r22445 * r22446; float r22448 = r22444 + r22447; float r22449 = c; float r22450 = i; float r22451 = r22449 * r22450; float r22452 = r22448 + r22451; return r22452; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r22453 = x; double r22454 = y; double r22455 = r22453 * r22454; double r22456 = z; double r22457 = t; double r22458 = r22456 * r22457; double r22459 = r22455 + r22458; double r22460 = a; double r22461 = b; double r22462 = r22460 * r22461; double r22463 = r22459 + r22462; double r22464 = c; double r22465 = i; double r22466 = r22464 * r22465; double r22467 = r22463 + r22466; return r22467; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r22468 = z; float r22469 = t; float r22470 = i; float r22471 = c; float r22472 = b; float r22473 = a; float r22474 = r22472 * r22473; float r22475 = fma(r22470, r22471, r22474); float r22476 = fma(r22468, r22469, r22475); float r22477 = y; float r22478 = x; float r22479 = r22477 * r22478; float r22480 = r22476 + r22479; return r22480; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r22481 = z; double r22482 = t; double r22483 = i; double r22484 = c; double r22485 = b; double r22486 = a; double r22487 = r22485 * r22486; double r22488 = fma(r22483, r22484, r22487); double r22489 = fma(r22481, r22482, r22488); double r22490 = y; double r22491 = x; double r22492 = r22490 * r22491; double r22493 = r22489 + r22492; return r22493; } 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 r22494, r22495, r22496, r22497, r22498, r22499, r22500, r22501, r22502, r22503, r22504, r22505, r22506, r22507, r22508; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22494); mpfr_init(r22495); mpfr_init(r22496); mpfr_init(r22497); mpfr_init(r22498); mpfr_init(r22499); mpfr_init(r22500); mpfr_init(r22501); mpfr_init(r22502); mpfr_init(r22503); mpfr_init(r22504); mpfr_init(r22505); mpfr_init(r22506); mpfr_init(r22507); mpfr_init(r22508); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22494, x, MPFR_RNDN); mpfr_set_d(r22495, y, MPFR_RNDN); mpfr_mul(r22496, r22494, r22495, MPFR_RNDN); mpfr_set_d(r22497, z, MPFR_RNDN); mpfr_set_d(r22498, t, MPFR_RNDN); mpfr_mul(r22499, r22497, r22498, MPFR_RNDN); mpfr_add(r22500, r22496, r22499, MPFR_RNDN); mpfr_set_d(r22501, a, MPFR_RNDN); mpfr_set_d(r22502, b, MPFR_RNDN); mpfr_mul(r22503, r22501, r22502, MPFR_RNDN); mpfr_add(r22504, r22500, r22503, MPFR_RNDN); mpfr_set_d(r22505, c, MPFR_RNDN); mpfr_set_d(r22506, i, MPFR_RNDN); mpfr_mul(r22507, r22505, r22506, MPFR_RNDN); mpfr_add(r22508, r22504, r22507, MPFR_RNDN); return mpfr_get_d(r22508, MPFR_RNDN); } static mpfr_t r22509, r22510, r22511, r22512, r22513, r22514, r22515, r22516, r22517, r22518, r22519, r22520, r22521; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22509); mpfr_init(r22510); mpfr_init(r22511); mpfr_init(r22512); mpfr_init(r22513); mpfr_init(r22514); mpfr_init(r22515); mpfr_init(r22516); mpfr_init(r22517); mpfr_init(r22518); mpfr_init(r22519); mpfr_init(r22520); mpfr_init(r22521); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22509, z, MPFR_RNDN); mpfr_set_d(r22510, t, MPFR_RNDN); mpfr_set_d(r22511, i, MPFR_RNDN); mpfr_set_d(r22512, c, MPFR_RNDN); mpfr_set_d(r22513, b, MPFR_RNDN); mpfr_set_d(r22514, a, MPFR_RNDN); mpfr_mul(r22515, r22513, r22514, MPFR_RNDN); mpfr_fma(r22516, r22511, r22512, r22515, MPFR_RNDN); mpfr_fma(r22517, r22509, r22510, r22516, MPFR_RNDN); mpfr_set_d(r22518, y, MPFR_RNDN); mpfr_set_d(r22519, x, MPFR_RNDN); mpfr_mul(r22520, r22518, r22519, MPFR_RNDN); mpfr_add(r22521, r22517, r22520, MPFR_RNDN); return mpfr_get_d(r22521, MPFR_RNDN); } static mpfr_t r22522, r22523, r22524, r22525, r22526, r22527, r22528, r22529, r22530, r22531, r22532, r22533, r22534; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22522); mpfr_init(r22523); mpfr_init(r22524); mpfr_init(r22525); mpfr_init(r22526); mpfr_init(r22527); mpfr_init(r22528); mpfr_init(r22529); mpfr_init(r22530); mpfr_init(r22531); mpfr_init(r22532); mpfr_init(r22533); mpfr_init(r22534); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r22522, z, MPFR_RNDN); mpfr_set_d(r22523, t, MPFR_RNDN); mpfr_set_d(r22524, i, MPFR_RNDN); mpfr_set_d(r22525, c, MPFR_RNDN); mpfr_set_d(r22526, b, MPFR_RNDN); mpfr_set_d(r22527, a, MPFR_RNDN); mpfr_mul(r22528, r22526, r22527, MPFR_RNDN); mpfr_fma(r22529, r22524, r22525, r22528, MPFR_RNDN); mpfr_fma(r22530, r22522, r22523, r22529, MPFR_RNDN); mpfr_set_d(r22531, y, MPFR_RNDN); mpfr_set_d(r22532, x, MPFR_RNDN); mpfr_mul(r22533, r22531, r22532, MPFR_RNDN); mpfr_add(r22534, r22530, r22533, MPFR_RNDN); return mpfr_get_d(r22534, MPFR_RNDN); }