#include #include #include #include #include char *name = "Linear.Matrix:det33 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 j) { float r17559 = x; float r17560 = y; float r17561 = z; float r17562 = r17560 * r17561; float r17563 = t; float r17564 = a; float r17565 = r17563 * r17564; float r17566 = r17562 - r17565; float r17567 = r17559 * r17566; float r17568 = b; float r17569 = c; float r17570 = r17569 * r17561; float r17571 = i; float r17572 = r17571 * r17564; float r17573 = r17570 - r17572; float r17574 = r17568 * r17573; float r17575 = r17567 - r17574; float r17576 = j; float r17577 = r17569 * r17563; float r17578 = r17571 * r17560; float r17579 = r17577 - r17578; float r17580 = r17576 * r17579; float r17581 = r17575 + r17580; return r17581; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i, double j) { double r17582 = x; double r17583 = y; double r17584 = z; double r17585 = r17583 * r17584; double r17586 = t; double r17587 = a; double r17588 = r17586 * r17587; double r17589 = r17585 - r17588; double r17590 = r17582 * r17589; double r17591 = b; double r17592 = c; double r17593 = r17592 * r17584; double r17594 = i; double r17595 = r17594 * r17587; double r17596 = r17593 - r17595; double r17597 = r17591 * r17596; double r17598 = r17590 - r17597; double r17599 = j; double r17600 = r17592 * r17586; double r17601 = r17594 * r17583; double r17602 = r17600 - r17601; double r17603 = r17599 * r17602; double r17604 = r17598 + r17603; return r17604; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i, float j) { float r17605 = t; float r17606 = -5.169060029703922e+185f; bool r17607 = r17605 <= r17606; float r17608 = c; float r17609 = r17608 * r17605; float r17610 = i; float r17611 = y; float r17612 = r17610 * r17611; float r17613 = r17609 - r17612; float r17614 = j; float r17615 = z; float r17616 = r17611 * r17615; float r17617 = a; float r17618 = r17605 * r17617; float r17619 = r17616 - r17618; float r17620 = x; float r17621 = r17619 * r17620; float r17622 = fma(r17613, r17614, r17621); float r17623 = b; float r17624 = r17615 * r17608; float r17625 = r17617 * r17610; float r17626 = r17624 - r17625; float r17627 = r17623 * r17626; float r17628 = r17622 - r17627; float r17629 = 4.693685328487849e+35f; bool r17630 = r17605 <= r17629; float r17631 = r17630 ? r17628 : r17628; float r17632 = r17607 ? r17628 : r17631; return r17632; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i, double j) { double r17633 = t; double r17634 = -5.169060029703922e+185; bool r17635 = r17633 <= r17634; double r17636 = c; double r17637 = r17636 * r17633; double r17638 = i; double r17639 = y; double r17640 = r17638 * r17639; double r17641 = r17637 - r17640; double r17642 = j; double r17643 = z; double r17644 = r17639 * r17643; double r17645 = a; double r17646 = r17633 * r17645; double r17647 = r17644 - r17646; double r17648 = x; double r17649 = r17647 * r17648; double r17650 = fma(r17641, r17642, r17649); double r17651 = b; double r17652 = r17643 * r17636; double r17653 = r17645 * r17638; double r17654 = r17652 - r17653; double r17655 = r17651 * r17654; double r17656 = r17650 - r17655; double r17657 = 4.693685328487849e+35; bool r17658 = r17633 <= r17657; double r17659 = r17658 ? r17656 : r17656; double r17660 = r17635 ? r17656 : r17659; return r17660; } 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 r17661, r17662, r17663, r17664, r17665, r17666, r17667, r17668, r17669, r17670, r17671, r17672, r17673, r17674, r17675, r17676, r17677, r17678, r17679, r17680, r17681, r17682, r17683; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r17661); mpfr_init(r17662); mpfr_init(r17663); mpfr_init(r17664); mpfr_init(r17665); mpfr_init(r17666); mpfr_init(r17667); mpfr_init(r17668); mpfr_init(r17669); mpfr_init(r17670); mpfr_init(r17671); mpfr_init(r17672); mpfr_init(r17673); mpfr_init(r17674); mpfr_init(r17675); mpfr_init(r17676); mpfr_init(r17677); mpfr_init(r17678); mpfr_init(r17679); mpfr_init(r17680); mpfr_init(r17681); mpfr_init(r17682); mpfr_init(r17683); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r17661, x, MPFR_RNDN); mpfr_set_d(r17662, y, MPFR_RNDN); mpfr_set_d(r17663, z, MPFR_RNDN); mpfr_mul(r17664, r17662, r17663, MPFR_RNDN); mpfr_set_d(r17665, t, MPFR_RNDN); mpfr_set_d(r17666, a, MPFR_RNDN); mpfr_mul(r17667, r17665, r17666, MPFR_RNDN); mpfr_sub(r17668, r17664, r17667, MPFR_RNDN); mpfr_mul(r17669, r17661, r17668, MPFR_RNDN); mpfr_set_d(r17670, b, MPFR_RNDN); mpfr_set_d(r17671, c, MPFR_RNDN); mpfr_mul(r17672, r17671, r17663, MPFR_RNDN); mpfr_set_d(r17673, i, MPFR_RNDN); mpfr_mul(r17674, r17673, r17666, MPFR_RNDN); mpfr_sub(r17675, r17672, r17674, MPFR_RNDN); mpfr_mul(r17676, r17670, r17675, MPFR_RNDN); mpfr_sub(r17677, r17669, r17676, MPFR_RNDN); mpfr_set_d(r17678, j, MPFR_RNDN); mpfr_mul(r17679, r17671, r17665, MPFR_RNDN); mpfr_mul(r17680, r17673, r17662, MPFR_RNDN); mpfr_sub(r17681, r17679, r17680, MPFR_RNDN); mpfr_mul(r17682, r17678, r17681, MPFR_RNDN); mpfr_add(r17683, r17677, r17682, MPFR_RNDN); return mpfr_get_d(r17683, MPFR_RNDN); } static mpfr_t r17684, r17685, r17686, r17687, r17688, r17689, r17690, r17691, r17692, r17693, r17694, r17695, r17696, r17697, r17698, r17699, r17700, r17701, r17702, r17703, r17704, r17705, r17706, r17707, r17708, r17709, r17710, r17711; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17684); mpfr_init_set_str(r17685, "-5.169060029703922e+185", 10, MPFR_RNDN); mpfr_init(r17686); mpfr_init(r17687); mpfr_init(r17688); mpfr_init(r17689); mpfr_init(r17690); mpfr_init(r17691); mpfr_init(r17692); mpfr_init(r17693); mpfr_init(r17694); mpfr_init(r17695); mpfr_init(r17696); mpfr_init(r17697); mpfr_init(r17698); mpfr_init(r17699); mpfr_init(r17700); mpfr_init(r17701); mpfr_init(r17702); mpfr_init(r17703); mpfr_init(r17704); mpfr_init(r17705); mpfr_init(r17706); mpfr_init(r17707); mpfr_init_set_str(r17708, "4.693685328487849e+35", 10, MPFR_RNDN); mpfr_init(r17709); mpfr_init(r17710); mpfr_init(r17711); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r17684, t, MPFR_RNDN); ; mpfr_set_si(r17686, mpfr_cmp(r17684, r17685) <= 0, MPFR_RNDN); mpfr_set_d(r17687, c, MPFR_RNDN); mpfr_mul(r17688, r17687, r17684, MPFR_RNDN); mpfr_set_d(r17689, i, MPFR_RNDN); mpfr_set_d(r17690, y, MPFR_RNDN); mpfr_mul(r17691, r17689, r17690, MPFR_RNDN); mpfr_sub(r17692, r17688, r17691, MPFR_RNDN); mpfr_set_d(r17693, j, MPFR_RNDN); mpfr_set_d(r17694, z, MPFR_RNDN); mpfr_mul(r17695, r17690, r17694, MPFR_RNDN); mpfr_set_d(r17696, a, MPFR_RNDN); mpfr_mul(r17697, r17684, r17696, MPFR_RNDN); mpfr_sub(r17698, r17695, r17697, MPFR_RNDN); mpfr_set_d(r17699, x, MPFR_RNDN); mpfr_mul(r17700, r17698, r17699, MPFR_RNDN); mpfr_fma(r17701, r17692, r17693, r17700, MPFR_RNDN); mpfr_set_d(r17702, b, MPFR_RNDN); mpfr_mul(r17703, r17694, r17687, MPFR_RNDN); mpfr_mul(r17704, r17696, r17689, MPFR_RNDN); mpfr_sub(r17705, r17703, r17704, MPFR_RNDN); mpfr_mul(r17706, r17702, r17705, MPFR_RNDN); mpfr_sub(r17707, r17701, r17706, MPFR_RNDN); ; mpfr_set_si(r17709, mpfr_cmp(r17684, r17708) <= 0, MPFR_RNDN); if (mpfr_get_si(r17709, MPFR_RNDN)) { mpfr_set(r17710, r17707, MPFR_RNDN); } else { mpfr_set(r17710, r17707, MPFR_RNDN); }; if (mpfr_get_si(r17686, MPFR_RNDN)) { mpfr_set(r17711, r17707, MPFR_RNDN); } else { mpfr_set(r17711, r17710, MPFR_RNDN); }; return mpfr_get_d(r17711, MPFR_RNDN); } static mpfr_t r17712, r17713, r17714, r17715, r17716, r17717, r17718, r17719, r17720, r17721, r17722, r17723, r17724, r17725, r17726, r17727, r17728, r17729, r17730, r17731, r17732, r17733, r17734, r17735, r17736, r17737, r17738, r17739; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17712); mpfr_init_set_str(r17713, "-5.169060029703922e+185", 10, MPFR_RNDN); mpfr_init(r17714); mpfr_init(r17715); mpfr_init(r17716); mpfr_init(r17717); mpfr_init(r17718); mpfr_init(r17719); mpfr_init(r17720); mpfr_init(r17721); mpfr_init(r17722); mpfr_init(r17723); mpfr_init(r17724); mpfr_init(r17725); mpfr_init(r17726); mpfr_init(r17727); mpfr_init(r17728); mpfr_init(r17729); mpfr_init(r17730); mpfr_init(r17731); mpfr_init(r17732); mpfr_init(r17733); mpfr_init(r17734); mpfr_init(r17735); mpfr_init_set_str(r17736, "4.693685328487849e+35", 10, MPFR_RNDN); mpfr_init(r17737); mpfr_init(r17738); mpfr_init(r17739); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i, double j) { mpfr_set_d(r17712, t, MPFR_RNDN); ; mpfr_set_si(r17714, mpfr_cmp(r17712, r17713) <= 0, MPFR_RNDN); mpfr_set_d(r17715, c, MPFR_RNDN); mpfr_mul(r17716, r17715, r17712, MPFR_RNDN); mpfr_set_d(r17717, i, MPFR_RNDN); mpfr_set_d(r17718, y, MPFR_RNDN); mpfr_mul(r17719, r17717, r17718, MPFR_RNDN); mpfr_sub(r17720, r17716, r17719, MPFR_RNDN); mpfr_set_d(r17721, j, MPFR_RNDN); mpfr_set_d(r17722, z, MPFR_RNDN); mpfr_mul(r17723, r17718, r17722, MPFR_RNDN); mpfr_set_d(r17724, a, MPFR_RNDN); mpfr_mul(r17725, r17712, r17724, MPFR_RNDN); mpfr_sub(r17726, r17723, r17725, MPFR_RNDN); mpfr_set_d(r17727, x, MPFR_RNDN); mpfr_mul(r17728, r17726, r17727, MPFR_RNDN); mpfr_fma(r17729, r17720, r17721, r17728, MPFR_RNDN); mpfr_set_d(r17730, b, MPFR_RNDN); mpfr_mul(r17731, r17722, r17715, MPFR_RNDN); mpfr_mul(r17732, r17724, r17717, MPFR_RNDN); mpfr_sub(r17733, r17731, r17732, MPFR_RNDN); mpfr_mul(r17734, r17730, r17733, MPFR_RNDN); mpfr_sub(r17735, r17729, r17734, MPFR_RNDN); ; mpfr_set_si(r17737, mpfr_cmp(r17712, r17736) <= 0, MPFR_RNDN); if (mpfr_get_si(r17737, MPFR_RNDN)) { mpfr_set(r17738, r17735, MPFR_RNDN); } else { mpfr_set(r17738, r17735, MPFR_RNDN); }; if (mpfr_get_si(r17714, MPFR_RNDN)) { mpfr_set(r17739, r17735, MPFR_RNDN); } else { mpfr_set(r17739, r17738, MPFR_RNDN); }; return mpfr_get_d(r17739, MPFR_RNDN); }