Method inv

src/complex/matrix.rs:1841–1857

(&self)

    }

    pub fn inv(&self) -> ComplexMatrix {
        let (p, q, l, u) = self.extract();
        let mut m = complex_inv_u(u) * complex_inv_l(l);
        // Q = Q1 Q2 Q3 ..
        for (idx1, idx2) in q.into_iter().enumerate().rev() {
            unsafe {
                m.swap(idx1, idx2, Shape::Row);
            }
        }
        // P = Pn-1 .. P3 P2 P1
        for (idx1, idx2) in p.into_iter().enumerate().rev() {
            unsafe {
                m.swap(idx1, idx2, Shape::Col);
            }
        }
        m
    }
}

/// MATLAB like eye - Identity matrix