62   CMatrixType m2(cols, cols2);  in lu_non_invertible()  local
94   m2 = CMatrixType::Random(cols,cols2);  in lu_non_invertible()
95   m3 = m1*m2;  in lu_non_invertible()
96   m2 = CMatrixType::Random(cols,cols2);  in lu_non_invertible()
98   m2.block(0,0,m2.rows(),m2.cols()) = lu.solve(m3);  in lu_non_invertible()
99   VERIFY_IS_APPROX(m3, m1*m2);  in lu_non_invertible()
103   m2 = m1.transpose()*m3;  in lu_non_invertible()
105   lu.template _solve_impl_transposed<false>(m2, m3);  in lu_non_invertible()
106   VERIFY_IS_APPROX(m2, m1.transpose()*m3);  in lu_non_invertible()
108   m3 = lu.transpose().solve(m2);  in lu_non_invertible()
109   VERIFY_IS_APPROX(m2, m1.transpose()*m3);  in lu_non_invertible()
113   m2 = m1.adjoint()*m3;  in lu_non_invertible()
115   lu.template _solve_impl_transposed<true>(m2, m3);  in lu_non_invertible()
116   VERIFY_IS_APPROX(m2, m1.adjoint()*m3);  in lu_non_invertible()
118   m3 = lu.adjoint().solve(m2);  in lu_non_invertible()
119   VERIFY_IS_APPROX(m2, m1.adjoint()*m3);  in lu_non_invertible()
132   MatrixType m1(size, size), m2(size, size), m3(size, size);  in lu_invertible()  local
149   m2 = lu.solve(m3);  in lu_invertible()
150   VERIFY_IS_APPROX(m3, m1*m2);  in lu_invertible()
152   VERIFY_IS_APPROX(m2, m1_inverse*m3);  in lu_invertible()
161   lu.template _solve_impl_transposed<false>(m3, m2);  in lu_invertible()
162   VERIFY_IS_APPROX(m3, m1.transpose()*m2);  in lu_invertible()
164   m3 = lu.transpose().solve(m2);  in lu_invertible()
165   VERIFY_IS_APPROX(m2, m1.transpose()*m3);  in lu_invertible()
168   lu.template _solve_impl_transposed<true>(m3, m2);  in lu_invertible()
169   VERIFY_IS_APPROX(m3, m1.adjoint()*m2);  in lu_invertible()
171   m3 = lu.adjoint().solve(m2);  in lu_invertible()
172   VERIFY_IS_APPROX(m2, m1.adjoint()*m3);  in lu_invertible()
188   MatrixType m1(size, size), m2(size, size), m3(size, size);  in lu_partial_piv()  local
195   m2 = plu.solve(m3);  in lu_partial_piv()
196   VERIFY_IS_APPROX(m3, m1*m2);  in lu_partial_piv()
198   VERIFY_IS_APPROX(m2, m1_inverse*m3);  in lu_partial_piv()
206   plu.template _solve_impl_transposed<false>(m3, m2);  in lu_partial_piv()
207   VERIFY_IS_APPROX(m3, m1.transpose()*m2);  in lu_partial_piv()
209   m3 = plu.transpose().solve(m2);  in lu_partial_piv()
210   VERIFY_IS_APPROX(m2, m1.transpose()*m3);  in lu_partial_piv()
213   plu.template _solve_impl_transposed<true>(m3, m2);  in lu_partial_piv()
214   VERIFY_IS_APPROX(m3, m1.adjoint()*m2);  in lu_partial_piv()
216   m3 = plu.adjoint().solve(m2);  in lu_partial_piv()
217   VERIFY_IS_APPROX(m2, m1.adjoint()*m3);  in lu_partial_piv()