#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 r10360 = x; float r10361 = 18.0; float r10362 = r10360 * r10361; float r10363 = y; float r10364 = r10362 * r10363; float r10365 = z; float r10366 = r10364 * r10365; float r10367 = t; float r10368 = r10366 * r10367; float r10369 = a; float r10370 = 4.0; float r10371 = r10369 * r10370; float r10372 = r10371 * r10367; float r10373 = r10368 - r10372; float r10374 = b; float r10375 = c; float r10376 = r10374 * r10375; float r10377 = r10373 + r10376; float r10378 = r10360 * r10370; float r10379 = i; float r10380 = r10378 * r10379; float r10381 = r10377 - r10380; float r10382 = j; float r10383 = 27.0; float r10384 = r10382 * r10383; float r10385 = k; float r10386 = r10384 * r10385; float r10387 = r10381 - r10386; return r10387; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r10388 = x; double r10389 = 18.0; double r10390 = r10388 * r10389; double r10391 = y; double r10392 = r10390 * r10391; double r10393 = z; double r10394 = r10392 * r10393; double r10395 = t; double r10396 = r10394 * r10395; double r10397 = a; double r10398 = 4.0; double r10399 = r10397 * r10398; double r10400 = r10399 * r10395; double r10401 = r10396 - r10400; double r10402 = b; double r10403 = c; double r10404 = r10402 * r10403; double r10405 = r10401 + r10404; double r10406 = r10388 * r10398; double r10407 = i; double r10408 = r10406 * r10407; double r10409 = r10405 - r10408; double r10410 = j; double r10411 = 27.0; double r10412 = r10410 * r10411; double r10413 = k; double r10414 = r10412 * r10413; double r10415 = r10409 - r10414; return r10415; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j, float k) { float r10416 = t; float r10417 = x; float r10418 = 18.0; float r10419 = r10417 * r10418; float r10420 = y; float r10421 = r10419 * r10420; float r10422 = z; float r10423 = r10421 * r10422; float r10424 = r10416 * r10423; float r10425 = a; float r10426 = -r10425; float r10427 = 4.0; float r10428 = r10427 * r10416; float r10429 = b; float r10430 = c; float r10431 = r10429 * r10430; float r10432 = fma(r10426, r10428, r10431); float r10433 = r10424 + r10432; float r10434 = r10427 * r10417; float r10435 = i; float r10436 = r10434 * r10435; float r10437 = r10433 - r10436; float r10438 = 27.0; float r10439 = j; float r10440 = r10438 * r10439; float r10441 = k; float r10442 = r10440 * r10441; float r10443 = r10437 - r10442; float r10444 = -4.970360939493975e+300; bool r10445 = r10443 <= r10444; float r10446 = r10416 * r10417; float r10447 = r10446 * r10422; float r10448 = r10420 * r10418; float r10449 = r10427 * r10425; float r10450 = -r10416; float r10451 = fma(r10449, r10450, r10431); float r10452 = fma(r10447, r10448, r10451); float r10453 = cbrt(r10452); float r10454 = r10416 * r10422; float r10455 = r10448 * r10454; float r10456 = fma(r10455, r10417, r10432); float r10457 = cbrt(r10456); float r10458 = r10457 * r10457; float r10459 = r10453 * r10458; float r10460 = r10459 - r10436; float r10461 = r10460 - r10442; float r10462 = 5.074392852067346e+299; bool r10463 = r10443 <= r10462; float r10464 = r10419 * r10454; float r10465 = fma(r10420, r10464, r10432); float r10466 = r10435 * r10427; float r10467 = r10417 * r10466; float r10468 = r10467 + r10442; float r10469 = r10465 - r10468; float r10470 = r10463 ? r10443 : r10469; float r10471 = r10445 ? r10461 : r10470; return r10471; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) { double r10472 = t; double r10473 = x; double r10474 = 18.0; double r10475 = r10473 * r10474; double r10476 = y; double r10477 = r10475 * r10476; double r10478 = z; double r10479 = r10477 * r10478; double r10480 = r10472 * r10479; double r10481 = a; double r10482 = -r10481; double r10483 = 4.0; double r10484 = r10483 * r10472; double r10485 = b; double r10486 = c; double r10487 = r10485 * r10486; double r10488 = fma(r10482, r10484, r10487); double r10489 = r10480 + r10488; double r10490 = r10483 * r10473; double r10491 = i; double r10492 = r10490 * r10491; double r10493 = r10489 - r10492; double r10494 = 27.0; double r10495 = j; double r10496 = r10494 * r10495; double r10497 = k; double r10498 = r10496 * r10497; double r10499 = r10493 - r10498; double r10500 = -4.970360939493975e+300; bool r10501 = r10499 <= r10500; double r10502 = r10472 * r10473; double r10503 = r10502 * r10478; double r10504 = r10476 * r10474; double r10505 = r10483 * r10481; double r10506 = -r10472; double r10507 = fma(r10505, r10506, r10487); double r10508 = fma(r10503, r10504, r10507); double r10509 = cbrt(r10508); double r10510 = r10472 * r10478; double r10511 = r10504 * r10510; double r10512 = fma(r10511, r10473, r10488); double r10513 = cbrt(r10512); double r10514 = r10513 * r10513; double r10515 = r10509 * r10514; double r10516 = r10515 - r10492; double r10517 = r10516 - r10498; double r10518 = 5.074392852067346e+299; bool r10519 = r10499 <= r10518; double r10520 = r10475 * r10510; double r10521 = fma(r10476, r10520, r10488); double r10522 = r10491 * r10483; double r10523 = r10473 * r10522; double r10524 = r10523 + r10498; double r10525 = r10521 - r10524; double r10526 = r10519 ? r10499 : r10525; double r10527 = r10501 ? r10517 : r10526; return r10527; } 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 r10528, r10529, r10530, r10531, r10532, r10533, r10534, r10535, r10536, r10537, r10538, r10539, r10540, r10541, r10542, r10543, r10544, r10545, r10546, r10547, r10548, r10549, r10550, r10551, r10552, r10553, r10554, r10555; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10528); mpfr_init_set_str(r10529, "18.0", 10, MPFR_RNDN); mpfr_init(r10530); mpfr_init(r10531); mpfr_init(r10532); mpfr_init(r10533); mpfr_init(r10534); mpfr_init(r10535); mpfr_init(r10536); mpfr_init(r10537); mpfr_init_set_str(r10538, "4.0", 10, MPFR_RNDN); mpfr_init(r10539); mpfr_init(r10540); mpfr_init(r10541); mpfr_init(r10542); mpfr_init(r10543); mpfr_init(r10544); mpfr_init(r10545); mpfr_init(r10546); mpfr_init(r10547); mpfr_init(r10548); mpfr_init(r10549); mpfr_init(r10550); mpfr_init_set_str(r10551, "27.0", 10, MPFR_RNDN); mpfr_init(r10552); mpfr_init(r10553); mpfr_init(r10554); mpfr_init(r10555); } 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(r10528, x, MPFR_RNDN); ; mpfr_mul(r10530, r10528, r10529, MPFR_RNDN); mpfr_set_d(r10531, y, MPFR_RNDN); mpfr_mul(r10532, r10530, r10531, MPFR_RNDN); mpfr_set_d(r10533, z, MPFR_RNDN); mpfr_mul(r10534, r10532, r10533, MPFR_RNDN); mpfr_set_d(r10535, t, MPFR_RNDN); mpfr_mul(r10536, r10534, r10535, MPFR_RNDN); mpfr_set_d(r10537, a, MPFR_RNDN); ; mpfr_mul(r10539, r10537, r10538, MPFR_RNDN); mpfr_mul(r10540, r10539, r10535, MPFR_RNDN); mpfr_sub(r10541, r10536, r10540, MPFR_RNDN); mpfr_set_d(r10542, b, MPFR_RNDN); mpfr_set_d(r10543, c, MPFR_RNDN); mpfr_mul(r10544, r10542, r10543, MPFR_RNDN); mpfr_add(r10545, r10541, r10544, MPFR_RNDN); mpfr_mul(r10546, r10528, r10538, MPFR_RNDN); mpfr_set_d(r10547, i, MPFR_RNDN); mpfr_mul(r10548, r10546, r10547, MPFR_RNDN); mpfr_sub(r10549, r10545, r10548, MPFR_RNDN); mpfr_set_d(r10550, j, MPFR_RNDN); ; mpfr_mul(r10552, r10550, r10551, MPFR_RNDN); mpfr_set_d(r10553, k, MPFR_RNDN); mpfr_mul(r10554, r10552, r10553, MPFR_RNDN); mpfr_sub(r10555, r10549, r10554, MPFR_RNDN); return mpfr_get_d(r10555, MPFR_RNDN); } static mpfr_t r10556, r10557, r10558, r10559, r10560, r10561, r10562, r10563, r10564, r10565, r10566, r10567, r10568, r10569, r10570, r10571, r10572, r10573, r10574, r10575, r10576, r10577, r10578, r10579, r10580, r10581, r10582, r10583, r10584, r10585, r10586, r10587, r10588, r10589, r10590, r10591, r10592, r10593, r10594, r10595, r10596, r10597, r10598, r10599, r10600, r10601, r10602, r10603, r10604, r10605, r10606, r10607, r10608, r10609, r10610, r10611; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r10556); mpfr_init(r10557); mpfr_init_set_str(r10558, "18.0", 10, MPFR_RNDN); mpfr_init(r10559); mpfr_init(r10560); mpfr_init(r10561); mpfr_init(r10562); mpfr_init(r10563); mpfr_init(r10564); mpfr_init(r10565); mpfr_init(r10566); mpfr_init_set_str(r10567, "4.0", 10, MPFR_RNDN); mpfr_init(r10568); mpfr_init(r10569); mpfr_init(r10570); mpfr_init(r10571); mpfr_init(r10572); mpfr_init(r10573); mpfr_init(r10574); mpfr_init(r10575); mpfr_init(r10576); mpfr_init(r10577); mpfr_init_set_str(r10578, "27.0", 10, MPFR_RNDN); mpfr_init(r10579); mpfr_init(r10580); mpfr_init(r10581); mpfr_init(r10582); mpfr_init(r10583); mpfr_init_set_str(r10584, "-4.970360939493975e+300", 10, MPFR_RNDN); mpfr_init(r10585); mpfr_init(r10586); mpfr_init(r10587); mpfr_init(r10588); mpfr_init(r10589); mpfr_init(r10590); mpfr_init(r10591); mpfr_init(r10592); mpfr_init(r10593); mpfr_init(r10594); mpfr_init(r10595); mpfr_init(r10596); mpfr_init(r10597); mpfr_init(r10598); mpfr_init(r10599); mpfr_init(r10600); mpfr_init(r10601); mpfr_init_set_str(r10602, "5.074392852067346e+299", 10, MPFR_RNDN); mpfr_init(r10603); mpfr_init(r10604); mpfr_init(r10605); mpfr_init(r10606); mpfr_init(r10607); mpfr_init(r10608); mpfr_init(r10609); mpfr_init(r10610); mpfr_init(r10611); } 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(r10556, t, MPFR_RNDN); mpfr_set_d(r10557, x, MPFR_RNDN); ; mpfr_mul(r10559, r10557, r10558, MPFR_RNDN); mpfr_set_d(r10560, y, MPFR_RNDN); mpfr_mul(r10561, r10559, r10560, MPFR_RNDN); mpfr_set_d(r10562, z, MPFR_RNDN); mpfr_mul(r10563, r10561, r10562, MPFR_RNDN); mpfr_mul(r10564, r10556, r10563, MPFR_RNDN); mpfr_set_d(r10565, a, MPFR_RNDN); mpfr_neg(r10566, r10565, MPFR_RNDN); ; mpfr_mul(r10568, r10567, r10556, MPFR_RNDN); mpfr_set_d(r10569, b, MPFR_RNDN); mpfr_set_d(r10570, c, MPFR_RNDN); mpfr_mul(r10571, r10569, r10570, MPFR_RNDN); mpfr_fma(r10572, r10566, r10568, r10571, MPFR_RNDN); mpfr_add(r10573, r10564, r10572, MPFR_RNDN); mpfr_mul(r10574, r10567, r10557, MPFR_RNDN); mpfr_set_d(r10575, i, MPFR_RNDN); mpfr_mul(r10576, r10574, r10575, MPFR_RNDN); mpfr_sub(r10577, r10573, r10576, MPFR_RNDN); ; mpfr_set_d(r10579, j, MPFR_RNDN); mpfr_mul(r10580, r10578, r10579, MPFR_RNDN); mpfr_set_d(r10581, k, MPFR_RNDN); mpfr_mul(r10582, r10580, r10581, MPFR_RNDN); mpfr_sub(r10583, r10577, r10582, MPFR_RNDN); ; mpfr_set_si(r10585, mpfr_cmp(r10583, r10584) <= 0, MPFR_RNDN); mpfr_mul(r10586, r10556, r10557, MPFR_RNDN); mpfr_mul(r10587, r10586, r10562, MPFR_RNDN); mpfr_mul(r10588, r10560, r10558, MPFR_RNDN); mpfr_mul(r10589, r10567, r10565, MPFR_RNDN); mpfr_neg(r10590, r10556, MPFR_RNDN); mpfr_fma(r10591, r10589, r10590, r10571, MPFR_RNDN); mpfr_fma(r10592, r10587, r10588, r10591, MPFR_RNDN); mpfr_cbrt(r10593, r10592, MPFR_RNDN); mpfr_mul(r10594, r10556, r10562, MPFR_RNDN); mpfr_mul(r10595, r10588, r10594, MPFR_RNDN); mpfr_fma(r10596, r10595, r10557, r10572, MPFR_RNDN); mpfr_cbrt(r10597, r10596, MPFR_RNDN); mpfr_mul(r10598, r10597, r10597, MPFR_RNDN); mpfr_mul(r10599, r10593, r10598, MPFR_RNDN); mpfr_sub(r10600, r10599, r10576, MPFR_RNDN); mpfr_sub(r10601, r10600, r10582, MPFR_RNDN); ; mpfr_set_si(r10603, mpfr_cmp(r10583, r10602) <= 0, MPFR_RNDN); mpfr_mul(r10604, r10559, r10594, MPFR_RNDN); mpfr_fma(r10605, r10560, r10604, r10572, MPFR_RNDN); mpfr_mul(r10606, r10575, r10567, MPFR_RNDN); mpfr_mul(r10607, r10557, r10606, MPFR_RNDN); mpfr_add(r10608, r10607, r10582, MPFR_RNDN); mpfr_sub(r10609, r10605, r10608, MPFR_RNDN); if (mpfr_get_si(r10603, MPFR_RNDN)) { mpfr_set(r10610, r10583, MPFR_RNDN); } else { mpfr_set(r10610, r10609, MPFR_RNDN); }; if (mpfr_get_si(r10585, MPFR_RNDN)) { mpfr_set(r10611, r10601, MPFR_RNDN); } else { mpfr_set(r10611, r10610, MPFR_RNDN); }; return mpfr_get_d(r10611, MPFR_RNDN); } static mpfr_t r10612, r10613, r10614, r10615, r10616, r10617, r10618, r10619, r10620, r10621, r10622, r10623, r10624, r10625, r10626, r10627, r10628, r10629, r10630, r10631, r10632, r10633, r10634, r10635, r10636, r10637, r10638, r10639, r10640, r10641, r10642, r10643, r10644, r10645, r10646, r10647, r10648, r10649, r10650, r10651, r10652, r10653, r10654, r10655, r10656, r10657, r10658, r10659, r10660, r10661, r10662, r10663, r10664, r10665, r10666, r10667; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r10612); mpfr_init(r10613); mpfr_init_set_str(r10614, "18.0", 10, MPFR_RNDN); mpfr_init(r10615); mpfr_init(r10616); mpfr_init(r10617); mpfr_init(r10618); mpfr_init(r10619); mpfr_init(r10620); mpfr_init(r10621); mpfr_init(r10622); mpfr_init_set_str(r10623, "4.0", 10, MPFR_RNDN); mpfr_init(r10624); mpfr_init(r10625); mpfr_init(r10626); mpfr_init(r10627); mpfr_init(r10628); mpfr_init(r10629); mpfr_init(r10630); mpfr_init(r10631); mpfr_init(r10632); mpfr_init(r10633); mpfr_init_set_str(r10634, "27.0", 10, MPFR_RNDN); mpfr_init(r10635); mpfr_init(r10636); mpfr_init(r10637); mpfr_init(r10638); mpfr_init(r10639); mpfr_init_set_str(r10640, "-4.970360939493975e+300", 10, MPFR_RNDN); mpfr_init(r10641); mpfr_init(r10642); mpfr_init(r10643); mpfr_init(r10644); mpfr_init(r10645); mpfr_init(r10646); mpfr_init(r10647); mpfr_init(r10648); mpfr_init(r10649); mpfr_init(r10650); mpfr_init(r10651); mpfr_init(r10652); mpfr_init(r10653); mpfr_init(r10654); mpfr_init(r10655); mpfr_init(r10656); mpfr_init(r10657); mpfr_init_set_str(r10658, "5.074392852067346e+299", 10, MPFR_RNDN); mpfr_init(r10659); mpfr_init(r10660); mpfr_init(r10661); mpfr_init(r10662); mpfr_init(r10663); mpfr_init(r10664); mpfr_init(r10665); mpfr_init(r10666); mpfr_init(r10667); } 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(r10612, t, MPFR_RNDN); mpfr_set_d(r10613, x, MPFR_RNDN); ; mpfr_mul(r10615, r10613, r10614, MPFR_RNDN); mpfr_set_d(r10616, y, MPFR_RNDN); mpfr_mul(r10617, r10615, r10616, MPFR_RNDN); mpfr_set_d(r10618, z, MPFR_RNDN); mpfr_mul(r10619, r10617, r10618, MPFR_RNDN); mpfr_mul(r10620, r10612, r10619, MPFR_RNDN); mpfr_set_d(r10621, a, MPFR_RNDN); mpfr_neg(r10622, r10621, MPFR_RNDN); ; mpfr_mul(r10624, r10623, r10612, MPFR_RNDN); mpfr_set_d(r10625, b, MPFR_RNDN); mpfr_set_d(r10626, c, MPFR_RNDN); mpfr_mul(r10627, r10625, r10626, MPFR_RNDN); mpfr_fma(r10628, r10622, r10624, r10627, MPFR_RNDN); mpfr_add(r10629, r10620, r10628, MPFR_RNDN); mpfr_mul(r10630, r10623, r10613, MPFR_RNDN); mpfr_set_d(r10631, i, MPFR_RNDN); mpfr_mul(r10632, r10630, r10631, MPFR_RNDN); mpfr_sub(r10633, r10629, r10632, MPFR_RNDN); ; mpfr_set_d(r10635, j, MPFR_RNDN); mpfr_mul(r10636, r10634, r10635, MPFR_RNDN); mpfr_set_d(r10637, k, MPFR_RNDN); mpfr_mul(r10638, r10636, r10637, MPFR_RNDN); mpfr_sub(r10639, r10633, r10638, MPFR_RNDN); ; mpfr_set_si(r10641, mpfr_cmp(r10639, r10640) <= 0, MPFR_RNDN); mpfr_mul(r10642, r10612, r10613, MPFR_RNDN); mpfr_mul(r10643, r10642, r10618, MPFR_RNDN); mpfr_mul(r10644, r10616, r10614, MPFR_RNDN); mpfr_mul(r10645, r10623, r10621, MPFR_RNDN); mpfr_neg(r10646, r10612, MPFR_RNDN); mpfr_fma(r10647, r10645, r10646, r10627, MPFR_RNDN); mpfr_fma(r10648, r10643, r10644, r10647, MPFR_RNDN); mpfr_cbrt(r10649, r10648, MPFR_RNDN); mpfr_mul(r10650, r10612, r10618, MPFR_RNDN); mpfr_mul(r10651, r10644, r10650, MPFR_RNDN); mpfr_fma(r10652, r10651, r10613, r10628, MPFR_RNDN); mpfr_cbrt(r10653, r10652, MPFR_RNDN); mpfr_mul(r10654, r10653, r10653, MPFR_RNDN); mpfr_mul(r10655, r10649, r10654, MPFR_RNDN); mpfr_sub(r10656, r10655, r10632, MPFR_RNDN); mpfr_sub(r10657, r10656, r10638, MPFR_RNDN); ; mpfr_set_si(r10659, mpfr_cmp(r10639, r10658) <= 0, MPFR_RNDN); mpfr_mul(r10660, r10615, r10650, MPFR_RNDN); mpfr_fma(r10661, r10616, r10660, r10628, MPFR_RNDN); mpfr_mul(r10662, r10631, r10623, MPFR_RNDN); mpfr_mul(r10663, r10613, r10662, MPFR_RNDN); mpfr_add(r10664, r10663, r10638, MPFR_RNDN); mpfr_sub(r10665, r10661, r10664, MPFR_RNDN); if (mpfr_get_si(r10659, MPFR_RNDN)) { mpfr_set(r10666, r10639, MPFR_RNDN); } else { mpfr_set(r10666, r10665, MPFR_RNDN); }; if (mpfr_get_si(r10641, MPFR_RNDN)) { mpfr_set(r10667, r10657, MPFR_RNDN); } else { mpfr_set(r10667, r10666, MPFR_RNDN); }; return mpfr_get_d(r10667, MPFR_RNDN); }