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