pub fn solve_owned( operator: DMatrix<Complex>, rhs: DMatrix<Complex>, ) -> Result<SolveResidual>
Solve while reusing the owned operator buffer for its LU factors.