#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 r16292 = x; float r16293 = y; float r16294 = r16292 * r16293; float r16295 = z; float r16296 = t; float r16297 = r16295 * r16296; float r16298 = r16294 + r16297; float r16299 = a; float r16300 = b; float r16301 = r16299 * r16300; float r16302 = r16298 + r16301; float r16303 = c; float r16304 = i; float r16305 = r16303 * r16304; float r16306 = r16302 + r16305; return r16306; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r16307 = x; double r16308 = y; double r16309 = r16307 * r16308; double r16310 = z; double r16311 = t; double r16312 = r16310 * r16311; double r16313 = r16309 + r16312; double r16314 = a; double r16315 = b; double r16316 = r16314 * r16315; double r16317 = r16313 + r16316; double r16318 = c; double r16319 = i; double r16320 = r16318 * r16319; double r16321 = r16317 + r16320; return r16321; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r16322 = z; float r16323 = t; float r16324 = y; float r16325 = x; float r16326 = i; float r16327 = c; float r16328 = b; float r16329 = a; float r16330 = r16328 * r16329; float r16331 = fma(r16326, r16327, r16330); float r16332 = fma(r16324, r16325, r16331); float r16333 = fma(r16322, r16323, r16332); return r16333; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r16334 = z; double r16335 = t; double r16336 = y; double r16337 = x; double r16338 = i; double r16339 = c; double r16340 = b; double r16341 = a; double r16342 = r16340 * r16341; double r16343 = fma(r16338, r16339, r16342); double r16344 = fma(r16336, r16337, r16343); double r16345 = fma(r16334, r16335, r16344); return r16345; } 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 r16346, r16347, r16348, r16349, r16350, r16351, r16352, r16353, r16354, r16355, r16356, r16357, r16358, r16359, r16360; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16346); mpfr_init(r16347); mpfr_init(r16348); mpfr_init(r16349); mpfr_init(r16350); mpfr_init(r16351); mpfr_init(r16352); mpfr_init(r16353); mpfr_init(r16354); mpfr_init(r16355); mpfr_init(r16356); mpfr_init(r16357); mpfr_init(r16358); mpfr_init(r16359); mpfr_init(r16360); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r16346, x, MPFR_RNDN); mpfr_set_d(r16347, y, MPFR_RNDN); mpfr_mul(r16348, r16346, r16347, MPFR_RNDN); mpfr_set_d(r16349, z, MPFR_RNDN); mpfr_set_d(r16350, t, MPFR_RNDN); mpfr_mul(r16351, r16349, r16350, MPFR_RNDN); mpfr_add(r16352, r16348, r16351, MPFR_RNDN); mpfr_set_d(r16353, a, MPFR_RNDN); mpfr_set_d(r16354, b, MPFR_RNDN); mpfr_mul(r16355, r16353, r16354, MPFR_RNDN); mpfr_add(r16356, r16352, r16355, MPFR_RNDN); mpfr_set_d(r16357, c, MPFR_RNDN); mpfr_set_d(r16358, i, MPFR_RNDN); mpfr_mul(r16359, r16357, r16358, MPFR_RNDN); mpfr_add(r16360, r16356, r16359, MPFR_RNDN); return mpfr_get_d(r16360, MPFR_RNDN); } static mpfr_t r16361, r16362, r16363, r16364, r16365, r16366, r16367, r16368, r16369, r16370, r16371, r16372; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16361); mpfr_init(r16362); mpfr_init(r16363); mpfr_init(r16364); mpfr_init(r16365); mpfr_init(r16366); mpfr_init(r16367); mpfr_init(r16368); mpfr_init(r16369); mpfr_init(r16370); mpfr_init(r16371); mpfr_init(r16372); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r16361, z, MPFR_RNDN); mpfr_set_d(r16362, t, MPFR_RNDN); mpfr_set_d(r16363, y, MPFR_RNDN); mpfr_set_d(r16364, x, MPFR_RNDN); mpfr_set_d(r16365, i, MPFR_RNDN); mpfr_set_d(r16366, c, MPFR_RNDN); mpfr_set_d(r16367, b, MPFR_RNDN); mpfr_set_d(r16368, a, MPFR_RNDN); mpfr_mul(r16369, r16367, r16368, MPFR_RNDN); mpfr_fma(r16370, r16365, r16366, r16369, MPFR_RNDN); mpfr_fma(r16371, r16363, r16364, r16370, MPFR_RNDN); mpfr_fma(r16372, r16361, r16362, r16371, MPFR_RNDN); return mpfr_get_d(r16372, MPFR_RNDN); } static mpfr_t r16373, r16374, r16375, r16376, r16377, r16378, r16379, r16380, r16381, r16382, r16383, r16384; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16373); mpfr_init(r16374); mpfr_init(r16375); mpfr_init(r16376); mpfr_init(r16377); mpfr_init(r16378); mpfr_init(r16379); mpfr_init(r16380); mpfr_init(r16381); mpfr_init(r16382); mpfr_init(r16383); mpfr_init(r16384); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r16373, z, MPFR_RNDN); mpfr_set_d(r16374, t, MPFR_RNDN); mpfr_set_d(r16375, y, MPFR_RNDN); mpfr_set_d(r16376, x, MPFR_RNDN); mpfr_set_d(r16377, i, MPFR_RNDN); mpfr_set_d(r16378, c, MPFR_RNDN); mpfr_set_d(r16379, b, MPFR_RNDN); mpfr_set_d(r16380, a, MPFR_RNDN); mpfr_mul(r16381, r16379, r16380, MPFR_RNDN); mpfr_fma(r16382, r16377, r16378, r16381, MPFR_RNDN); mpfr_fma(r16383, r16375, r16376, r16382, MPFR_RNDN); mpfr_fma(r16384, r16373, r16374, r16383, MPFR_RNDN); return mpfr_get_d(r16384, MPFR_RNDN); }