#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 r20429 = x; float r20430 = y; float r20431 = r20429 * r20430; float r20432 = z; float r20433 = t; float r20434 = r20432 * r20433; float r20435 = r20431 + r20434; float r20436 = a; float r20437 = b; float r20438 = r20436 * r20437; float r20439 = r20435 + r20438; float r20440 = c; float r20441 = i; float r20442 = r20440 * r20441; float r20443 = r20439 + r20442; return r20443; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r20444 = x; double r20445 = y; double r20446 = r20444 * r20445; double r20447 = z; double r20448 = t; double r20449 = r20447 * r20448; double r20450 = r20446 + r20449; double r20451 = a; double r20452 = b; double r20453 = r20451 * r20452; double r20454 = r20450 + r20453; double r20455 = c; double r20456 = i; double r20457 = r20455 * r20456; double r20458 = r20454 + r20457; return r20458; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r20459 = x; float r20460 = y; float r20461 = r20459 * r20460; float r20462 = z; float r20463 = t; float r20464 = r20462 * r20463; float r20465 = r20461 + r20464; float r20466 = a; float r20467 = b; float r20468 = r20466 * r20467; float r20469 = r20465 + r20468; float r20470 = c; float r20471 = i; float r20472 = r20470 * r20471; float r20473 = r20469 + r20472; return r20473; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r20474 = x; double r20475 = y; double r20476 = r20474 * r20475; double r20477 = z; double r20478 = t; double r20479 = r20477 * r20478; double r20480 = r20476 + r20479; double r20481 = a; double r20482 = b; double r20483 = r20481 * r20482; double r20484 = r20480 + r20483; double r20485 = c; double r20486 = i; double r20487 = r20485 * r20486; double r20488 = r20484 + r20487; return r20488; } 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 r20489, r20490, r20491, r20492, r20493, r20494, r20495, r20496, r20497, r20498, r20499, r20500, r20501, r20502, r20503; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r20489); mpfr_init(r20490); mpfr_init(r20491); mpfr_init(r20492); mpfr_init(r20493); mpfr_init(r20494); mpfr_init(r20495); mpfr_init(r20496); mpfr_init(r20497); mpfr_init(r20498); mpfr_init(r20499); mpfr_init(r20500); mpfr_init(r20501); mpfr_init(r20502); mpfr_init(r20503); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r20489, x, MPFR_RNDN); mpfr_set_d(r20490, y, MPFR_RNDN); mpfr_mul(r20491, r20489, r20490, MPFR_RNDN); mpfr_set_d(r20492, z, MPFR_RNDN); mpfr_set_d(r20493, t, MPFR_RNDN); mpfr_mul(r20494, r20492, r20493, MPFR_RNDN); mpfr_add(r20495, r20491, r20494, MPFR_RNDN); mpfr_set_d(r20496, a, MPFR_RNDN); mpfr_set_d(r20497, b, MPFR_RNDN); mpfr_mul(r20498, r20496, r20497, MPFR_RNDN); mpfr_add(r20499, r20495, r20498, MPFR_RNDN); mpfr_set_d(r20500, c, MPFR_RNDN); mpfr_set_d(r20501, i, MPFR_RNDN); mpfr_mul(r20502, r20500, r20501, MPFR_RNDN); mpfr_add(r20503, r20499, r20502, MPFR_RNDN); return mpfr_get_d(r20503, MPFR_RNDN); } static mpfr_t r20504, r20505, r20506, r20507, r20508, r20509, r20510, r20511, r20512, r20513, r20514, r20515, r20516, r20517, r20518; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r20504); mpfr_init(r20505); mpfr_init(r20506); mpfr_init(r20507); mpfr_init(r20508); mpfr_init(r20509); mpfr_init(r20510); mpfr_init(r20511); mpfr_init(r20512); mpfr_init(r20513); mpfr_init(r20514); mpfr_init(r20515); mpfr_init(r20516); mpfr_init(r20517); mpfr_init(r20518); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r20504, x, MPFR_RNDN); mpfr_set_d(r20505, y, MPFR_RNDN); mpfr_mul(r20506, r20504, r20505, MPFR_RNDN); mpfr_set_d(r20507, z, MPFR_RNDN); mpfr_set_d(r20508, t, MPFR_RNDN); mpfr_mul(r20509, r20507, r20508, MPFR_RNDN); mpfr_add(r20510, r20506, r20509, MPFR_RNDN); mpfr_set_d(r20511, a, MPFR_RNDN); mpfr_set_d(r20512, b, MPFR_RNDN); mpfr_mul(r20513, r20511, r20512, MPFR_RNDN); mpfr_add(r20514, r20510, r20513, MPFR_RNDN); mpfr_set_d(r20515, c, MPFR_RNDN); mpfr_set_d(r20516, i, MPFR_RNDN); mpfr_mul(r20517, r20515, r20516, MPFR_RNDN); mpfr_add(r20518, r20514, r20517, MPFR_RNDN); return mpfr_get_d(r20518, MPFR_RNDN); } static mpfr_t r20519, r20520, r20521, r20522, r20523, r20524, r20525, r20526, r20527, r20528, r20529, r20530, r20531, r20532, r20533; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r20519); mpfr_init(r20520); mpfr_init(r20521); mpfr_init(r20522); mpfr_init(r20523); mpfr_init(r20524); mpfr_init(r20525); mpfr_init(r20526); mpfr_init(r20527); mpfr_init(r20528); mpfr_init(r20529); mpfr_init(r20530); mpfr_init(r20531); mpfr_init(r20532); mpfr_init(r20533); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r20519, x, MPFR_RNDN); mpfr_set_d(r20520, y, MPFR_RNDN); mpfr_mul(r20521, r20519, r20520, MPFR_RNDN); mpfr_set_d(r20522, z, MPFR_RNDN); mpfr_set_d(r20523, t, MPFR_RNDN); mpfr_mul(r20524, r20522, r20523, MPFR_RNDN); mpfr_add(r20525, r20521, r20524, MPFR_RNDN); mpfr_set_d(r20526, a, MPFR_RNDN); mpfr_set_d(r20527, b, MPFR_RNDN); mpfr_mul(r20528, r20526, r20527, MPFR_RNDN); mpfr_add(r20529, r20525, r20528, MPFR_RNDN); mpfr_set_d(r20530, c, MPFR_RNDN); mpfr_set_d(r20531, i, MPFR_RNDN); mpfr_mul(r20532, r20530, r20531, MPFR_RNDN); mpfr_add(r20533, r20529, r20532, MPFR_RNDN); return mpfr_get_d(r20533, MPFR_RNDN); }