#include #include #include #include #include char *name = "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r17381 = x; float r17382 = 18.0f; float r17383 = r17381 * r17382; float r17384 = y; float r17385 = r17383 * r17384; float r17386 = z; float r17387 = r17385 * r17386; float r17388 = t; float r17389 = r17387 * r17388; float r17390 = a; float r17391 = 4.0f; float r17392 = r17390 * r17391; float r17393 = r17392 * r17388; float r17394 = r17389 - r17393; float r17395 = b; float r17396 = c; float r17397 = r17395 * r17396; float r17398 = r17394 + r17397; float r17399 = r17381 * r17391; float r17400 = i; float r17401 = r17399 * r17400; float r17402 = r17398 - r17401; float r17403 = j; float r17404 = 27.0f; float r17405 = r17403 * r17404; float r17406 = k; float r17407 = r17405 * r17406; float r17408 = r17402 - r17407; return r17408; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r17409 = x; double r17410 = 18.0; double r17411 = r17409 * r17410; double r17412 = y; double r17413 = r17411 * r17412; double r17414 = z; double r17415 = r17413 * r17414; double r17416 = t; double r17417 = r17415 * r17416; double r17418 = a; double r17419 = 4.0; double r17420 = r17418 * r17419; double r17421 = r17420 * r17416; double r17422 = r17417 - r17421; double r17423 = b; double r17424 = c; double r17425 = r17423 * r17424; double r17426 = r17422 + r17425; double r17427 = r17409 * r17419; double r17428 = i; double r17429 = r17427 * r17428; double r17430 = r17426 - r17429; double r17431 = j; double r17432 = 27.0; double r17433 = r17431 * r17432; double r17434 = k; double r17435 = r17433 * r17434; double r17436 = r17430 - r17435; return r17436; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r17437 = 18.0f; float r17438 = y; float r17439 = t; float r17440 = x; float r17441 = z; float r17442 = r17440 * r17441; float r17443 = r17439 * r17442; float r17444 = r17438 * r17443; float r17445 = 1.0f; float r17446 = pow(r17444, r17445); float r17447 = r17437 * r17446; float r17448 = 4.0f; float r17449 = a; float r17450 = r17439 * r17449; float r17451 = i; float r17452 = r17451 * r17440; float r17453 = r17450 + r17452; float r17454 = r17448 * r17453; float r17455 = r17447 - r17454; float r17456 = c; float r17457 = b; float r17458 = r17456 * r17457; float r17459 = j; float r17460 = 27.0f; float r17461 = k; float r17462 = r17460 * r17461; float r17463 = r17459 * r17462; float r17464 = r17458 - r17463; float r17465 = r17455 + r17464; return r17465; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r17466 = 18.0; double r17467 = y; double r17468 = t; double r17469 = x; double r17470 = z; double r17471 = r17469 * r17470; double r17472 = r17468 * r17471; double r17473 = r17467 * r17472; double r17474 = 1.0; double r17475 = pow(r17473, r17474); double r17476 = r17466 * r17475; double r17477 = 4.0; double r17478 = a; double r17479 = r17468 * r17478; double r17480 = i; double r17481 = r17480 * r17469; double r17482 = r17479 + r17481; double r17483 = r17477 * r17482; double r17484 = r17476 - r17483; double r17485 = c; double r17486 = b; double r17487 = r17485 * r17486; double r17488 = j; double r17489 = 27.0; double r17490 = k; double r17491 = r17489 * r17490; double r17492 = r17488 * r17491; double r17493 = r17487 - r17492; double r17494 = r17484 + r17493; return r17494; } 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 r17495, r17496, r17497, r17498, r17499, r17500, r17501, r17502, r17503, r17504, r17505, r17506, r17507, r17508, r17509, r17510, r17511, r17512, r17513, r17514, r17515, r17516, r17517, r17518, r17519, r17520, r17521, r17522; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17495); mpfr_init_set_str(r17496, "18.0", 10, MPFR_RNDN); mpfr_init(r17497); mpfr_init(r17498); mpfr_init(r17499); mpfr_init(r17500); mpfr_init(r17501); mpfr_init(r17502); mpfr_init(r17503); mpfr_init(r17504); mpfr_init_set_str(r17505, "4.0", 10, MPFR_RNDN); mpfr_init(r17506); mpfr_init(r17507); mpfr_init(r17508); mpfr_init(r17509); mpfr_init(r17510); mpfr_init(r17511); mpfr_init(r17512); mpfr_init(r17513); mpfr_init(r17514); mpfr_init(r17515); mpfr_init(r17516); mpfr_init(r17517); mpfr_init_set_str(r17518, "27.0", 10, MPFR_RNDN); mpfr_init(r17519); mpfr_init(r17520); mpfr_init(r17521); mpfr_init(r17522); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { mpfr_set_d(r17495, x, MPFR_RNDN); ; mpfr_mul(r17497, r17495, r17496, MPFR_RNDN); mpfr_set_d(r17498, y, MPFR_RNDN); mpfr_mul(r17499, r17497, r17498, MPFR_RNDN); mpfr_set_d(r17500, z, MPFR_RNDN); mpfr_mul(r17501, r17499, r17500, MPFR_RNDN); mpfr_set_d(r17502, t, MPFR_RNDN); mpfr_mul(r17503, r17501, r17502, MPFR_RNDN); mpfr_set_d(r17504, a, MPFR_RNDN); ; mpfr_mul(r17506, r17504, r17505, MPFR_RNDN); mpfr_mul(r17507, r17506, r17502, MPFR_RNDN); mpfr_sub(r17508, r17503, r17507, MPFR_RNDN); mpfr_set_d(r17509, b, MPFR_RNDN); mpfr_set_d(r17510, c, MPFR_RNDN); mpfr_mul(r17511, r17509, r17510, MPFR_RNDN); mpfr_add(r17512, r17508, r17511, MPFR_RNDN); mpfr_mul(r17513, r17495, r17505, MPFR_RNDN); mpfr_set_d(r17514, i, MPFR_RNDN); mpfr_mul(r17515, r17513, r17514, MPFR_RNDN); mpfr_sub(r17516, r17512, r17515, MPFR_RNDN); mpfr_set_d(r17517, j, MPFR_RNDN); ; mpfr_mul(r17519, r17517, r17518, MPFR_RNDN); mpfr_set_d(r17520, k, MPFR_RNDN); mpfr_mul(r17521, r17519, r17520, MPFR_RNDN); mpfr_sub(r17522, r17516, r17521, MPFR_RNDN); return mpfr_get_d(r17522, MPFR_RNDN); } static mpfr_t r17523, r17524, r17525, r17526, r17527, r17528, r17529, r17530, r17531, r17532, r17533, r17534, r17535, r17536, r17537, r17538, r17539, r17540, r17541, r17542, r17543, r17544, r17545, r17546, r17547, r17548, r17549, r17550, r17551; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17523, "18.0", 10, MPFR_RNDN); mpfr_init(r17524); mpfr_init(r17525); mpfr_init(r17526); mpfr_init(r17527); mpfr_init(r17528); mpfr_init(r17529); mpfr_init(r17530); mpfr_init_set_str(r17531, "1", 10, MPFR_RNDN); mpfr_init(r17532); mpfr_init(r17533); mpfr_init_set_str(r17534, "4.0", 10, MPFR_RNDN); mpfr_init(r17535); mpfr_init(r17536); mpfr_init(r17537); mpfr_init(r17538); mpfr_init(r17539); mpfr_init(r17540); mpfr_init(r17541); mpfr_init(r17542); mpfr_init(r17543); mpfr_init(r17544); mpfr_init(r17545); mpfr_init_set_str(r17546, "27.0", 10, MPFR_RNDN); mpfr_init(r17547); mpfr_init(r17548); mpfr_init(r17549); mpfr_init(r17550); mpfr_init(r17551); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { ; mpfr_set_d(r17524, y, MPFR_RNDN); mpfr_set_d(r17525, t, MPFR_RNDN); mpfr_set_d(r17526, x, MPFR_RNDN); mpfr_set_d(r17527, z, MPFR_RNDN); mpfr_mul(r17528, r17526, r17527, MPFR_RNDN); mpfr_mul(r17529, r17525, r17528, MPFR_RNDN); mpfr_mul(r17530, r17524, r17529, MPFR_RNDN); ; mpfr_pow(r17532, r17530, r17531, MPFR_RNDN); mpfr_mul(r17533, r17523, r17532, MPFR_RNDN); ; mpfr_set_d(r17535, a, MPFR_RNDN); mpfr_mul(r17536, r17525, r17535, MPFR_RNDN); mpfr_set_d(r17537, i, MPFR_RNDN); mpfr_mul(r17538, r17537, r17526, MPFR_RNDN); mpfr_add(r17539, r17536, r17538, MPFR_RNDN); mpfr_mul(r17540, r17534, r17539, MPFR_RNDN); mpfr_sub(r17541, r17533, r17540, MPFR_RNDN); mpfr_set_d(r17542, c, MPFR_RNDN); mpfr_set_d(r17543, b, MPFR_RNDN); mpfr_mul(r17544, r17542, r17543, MPFR_RNDN); mpfr_set_d(r17545, j, MPFR_RNDN); ; mpfr_set_d(r17547, k, MPFR_RNDN); mpfr_mul(r17548, r17546, r17547, MPFR_RNDN); mpfr_mul(r17549, r17545, r17548, MPFR_RNDN); mpfr_sub(r17550, r17544, r17549, MPFR_RNDN); mpfr_add(r17551, r17541, r17550, MPFR_RNDN); return mpfr_get_d(r17551, MPFR_RNDN); } static mpfr_t r17552, r17553, r17554, r17555, r17556, r17557, r17558, r17559, r17560, r17561, r17562, r17563, r17564, r17565, r17566, r17567, r17568, r17569, r17570, r17571, r17572, r17573, r17574, r17575, r17576, r17577, r17578, r17579, r17580; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17552, "18.0", 10, MPFR_RNDN); mpfr_init(r17553); mpfr_init(r17554); mpfr_init(r17555); mpfr_init(r17556); mpfr_init(r17557); mpfr_init(r17558); mpfr_init(r17559); mpfr_init_set_str(r17560, "1", 10, MPFR_RNDN); mpfr_init(r17561); mpfr_init(r17562); mpfr_init_set_str(r17563, "4.0", 10, MPFR_RNDN); mpfr_init(r17564); mpfr_init(r17565); mpfr_init(r17566); mpfr_init(r17567); mpfr_init(r17568); mpfr_init(r17569); mpfr_init(r17570); mpfr_init(r17571); mpfr_init(r17572); mpfr_init(r17573); mpfr_init(r17574); mpfr_init_set_str(r17575, "27.0", 10, MPFR_RNDN); mpfr_init(r17576); mpfr_init(r17577); mpfr_init(r17578); mpfr_init(r17579); mpfr_init(r17580); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { ; mpfr_set_d(r17553, y, MPFR_RNDN); mpfr_set_d(r17554, t, MPFR_RNDN); mpfr_set_d(r17555, x, MPFR_RNDN); mpfr_set_d(r17556, z, MPFR_RNDN); mpfr_mul(r17557, r17555, r17556, MPFR_RNDN); mpfr_mul(r17558, r17554, r17557, MPFR_RNDN); mpfr_mul(r17559, r17553, r17558, MPFR_RNDN); ; mpfr_pow(r17561, r17559, r17560, MPFR_RNDN); mpfr_mul(r17562, r17552, r17561, MPFR_RNDN); ; mpfr_set_d(r17564, a, MPFR_RNDN); mpfr_mul(r17565, r17554, r17564, MPFR_RNDN); mpfr_set_d(r17566, i, MPFR_RNDN); mpfr_mul(r17567, r17566, r17555, MPFR_RNDN); mpfr_add(r17568, r17565, r17567, MPFR_RNDN); mpfr_mul(r17569, r17563, r17568, MPFR_RNDN); mpfr_sub(r17570, r17562, r17569, MPFR_RNDN); mpfr_set_d(r17571, c, MPFR_RNDN); mpfr_set_d(r17572, b, MPFR_RNDN); mpfr_mul(r17573, r17571, r17572, MPFR_RNDN); mpfr_set_d(r17574, j, MPFR_RNDN); ; mpfr_set_d(r17576, k, MPFR_RNDN); mpfr_mul(r17577, r17575, r17576, MPFR_RNDN); mpfr_mul(r17578, r17574, r17577, MPFR_RNDN); mpfr_sub(r17579, r17573, r17578, MPFR_RNDN); mpfr_add(r17580, r17570, r17579, MPFR_RNDN); return mpfr_get_d(r17580, MPFR_RNDN); }