#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 r23118 = x; float r23119 = y; float r23120 = r23118 * r23119; float r23121 = z; float r23122 = t; float r23123 = r23121 * r23122; float r23124 = r23120 + r23123; float r23125 = a; float r23126 = b; float r23127 = r23125 * r23126; float r23128 = r23124 + r23127; float r23129 = c; float r23130 = i; float r23131 = r23129 * r23130; float r23132 = r23128 + r23131; return r23132; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r23133 = x; double r23134 = y; double r23135 = r23133 * r23134; double r23136 = z; double r23137 = t; double r23138 = r23136 * r23137; double r23139 = r23135 + r23138; double r23140 = a; double r23141 = b; double r23142 = r23140 * r23141; double r23143 = r23139 + r23142; double r23144 = c; double r23145 = i; double r23146 = r23144 * r23145; double r23147 = r23143 + r23146; return r23147; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r23148 = x; float r23149 = y; float r23150 = r23148 * r23149; float r23151 = z; float r23152 = t; float r23153 = r23151 * r23152; float r23154 = r23150 + r23153; float r23155 = a; float r23156 = b; float r23157 = r23155 * r23156; float r23158 = r23154 + r23157; float r23159 = c; float r23160 = i; float r23161 = r23159 * r23160; float r23162 = r23158 + r23161; return r23162; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r23163 = x; double r23164 = y; double r23165 = r23163 * r23164; double r23166 = z; double r23167 = t; double r23168 = r23166 * r23167; double r23169 = r23165 + r23168; double r23170 = a; double r23171 = b; double r23172 = r23170 * r23171; double r23173 = r23169 + r23172; double r23174 = c; double r23175 = i; double r23176 = r23174 * r23175; double r23177 = r23173 + r23176; return r23177; } 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 r23178, r23179, r23180, r23181, r23182, r23183, r23184, r23185, r23186, r23187, r23188, r23189, r23190, r23191, r23192; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23178); mpfr_init(r23179); mpfr_init(r23180); mpfr_init(r23181); mpfr_init(r23182); mpfr_init(r23183); mpfr_init(r23184); mpfr_init(r23185); mpfr_init(r23186); mpfr_init(r23187); mpfr_init(r23188); mpfr_init(r23189); mpfr_init(r23190); mpfr_init(r23191); mpfr_init(r23192); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r23178, x, MPFR_RNDN); mpfr_set_d(r23179, y, MPFR_RNDN); mpfr_mul(r23180, r23178, r23179, MPFR_RNDN); mpfr_set_d(r23181, z, MPFR_RNDN); mpfr_set_d(r23182, t, MPFR_RNDN); mpfr_mul(r23183, r23181, r23182, MPFR_RNDN); mpfr_add(r23184, r23180, r23183, MPFR_RNDN); mpfr_set_d(r23185, a, MPFR_RNDN); mpfr_set_d(r23186, b, MPFR_RNDN); mpfr_mul(r23187, r23185, r23186, MPFR_RNDN); mpfr_add(r23188, r23184, r23187, MPFR_RNDN); mpfr_set_d(r23189, c, MPFR_RNDN); mpfr_set_d(r23190, i, MPFR_RNDN); mpfr_mul(r23191, r23189, r23190, MPFR_RNDN); mpfr_add(r23192, r23188, r23191, MPFR_RNDN); return mpfr_get_d(r23192, MPFR_RNDN); } static mpfr_t r23193, r23194, r23195, r23196, r23197, r23198, r23199, r23200, r23201, r23202, r23203, r23204, r23205, r23206, r23207; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23193); mpfr_init(r23194); mpfr_init(r23195); mpfr_init(r23196); mpfr_init(r23197); mpfr_init(r23198); mpfr_init(r23199); mpfr_init(r23200); mpfr_init(r23201); mpfr_init(r23202); mpfr_init(r23203); mpfr_init(r23204); mpfr_init(r23205); mpfr_init(r23206); mpfr_init(r23207); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r23193, x, MPFR_RNDN); mpfr_set_d(r23194, y, MPFR_RNDN); mpfr_mul(r23195, r23193, r23194, MPFR_RNDN); mpfr_set_d(r23196, z, MPFR_RNDN); mpfr_set_d(r23197, t, MPFR_RNDN); mpfr_mul(r23198, r23196, r23197, MPFR_RNDN); mpfr_add(r23199, r23195, r23198, MPFR_RNDN); mpfr_set_d(r23200, a, MPFR_RNDN); mpfr_set_d(r23201, b, MPFR_RNDN); mpfr_mul(r23202, r23200, r23201, MPFR_RNDN); mpfr_add(r23203, r23199, r23202, MPFR_RNDN); mpfr_set_d(r23204, c, MPFR_RNDN); mpfr_set_d(r23205, i, MPFR_RNDN); mpfr_mul(r23206, r23204, r23205, MPFR_RNDN); mpfr_add(r23207, r23203, r23206, MPFR_RNDN); return mpfr_get_d(r23207, MPFR_RNDN); } static mpfr_t r23208, r23209, r23210, r23211, r23212, r23213, r23214, r23215, r23216, r23217, r23218, r23219, r23220, r23221, r23222; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23208); mpfr_init(r23209); mpfr_init(r23210); mpfr_init(r23211); mpfr_init(r23212); mpfr_init(r23213); mpfr_init(r23214); mpfr_init(r23215); mpfr_init(r23216); mpfr_init(r23217); mpfr_init(r23218); mpfr_init(r23219); mpfr_init(r23220); mpfr_init(r23221); mpfr_init(r23222); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r23208, x, MPFR_RNDN); mpfr_set_d(r23209, y, MPFR_RNDN); mpfr_mul(r23210, r23208, r23209, MPFR_RNDN); mpfr_set_d(r23211, z, MPFR_RNDN); mpfr_set_d(r23212, t, MPFR_RNDN); mpfr_mul(r23213, r23211, r23212, MPFR_RNDN); mpfr_add(r23214, r23210, r23213, MPFR_RNDN); mpfr_set_d(r23215, a, MPFR_RNDN); mpfr_set_d(r23216, b, MPFR_RNDN); mpfr_mul(r23217, r23215, r23216, MPFR_RNDN); mpfr_add(r23218, r23214, r23217, MPFR_RNDN); mpfr_set_d(r23219, c, MPFR_RNDN); mpfr_set_d(r23220, i, MPFR_RNDN); mpfr_mul(r23221, r23219, r23220, MPFR_RNDN); mpfr_add(r23222, r23218, r23221, MPFR_RNDN); return mpfr_get_d(r23222, MPFR_RNDN); }