(&self)
| 1821 | } |
| 1822 | |
| 1823 | pub fn det(&self) -> C64 { |
| 1824 | // sgn of perms |
| 1825 | let mut sgn_p = 1f64; |
| 1826 | let mut sgn_q = 1f64; |
| 1827 | for (i, &j) in self.p.iter().enumerate() { |
| 1828 | if i != j { |
| 1829 | sgn_p *= -1f64; |
| 1830 | } |
| 1831 | } |
| 1832 | for (i, &j) in self.q.iter().enumerate() { |
| 1833 | if i != j { |
| 1834 | sgn_q *= -1f64; |
| 1835 | } |
| 1836 | } |
| 1837 | |
| 1838 | self.u.diag().reduce(Complex::one(), |x, y| x * y) * sgn_p * sgn_q |
| 1839 | } |
| 1840 | |
| 1841 | pub fn inv(&self) -> ComplexMatrix { |
| 1842 | let (p, q, l, u) = self.extract(); |