Lines Matching refs:lu
52 LU<MatrixType> lu(m1); in lu_non_invertible() local
53 typename LU<MatrixType>::KernelResultType m1kernel = lu.kernel(); in lu_non_invertible()
54 typename LU<MatrixType>::ImageResultType m1image = lu.image(); in lu_non_invertible()
56 VERIFY(rank == lu.rank()); in lu_non_invertible()
57 VERIFY(cols - lu.rank() == lu.dimensionOfKernel()); in lu_non_invertible()
58 VERIFY(!lu.isInjective()); in lu_non_invertible()
59 VERIFY(!lu.isInvertible()); in lu_non_invertible()
60 VERIFY(lu.isSurjective() == (lu.rank() == rows)); in lu_non_invertible()
62 VERIFY(m1image.lu().rank() == rank); in lu_non_invertible()
65 VERIFY(sidebyside.lu().rank() == rank); in lu_non_invertible()
69 lu.solve(m3, &m2); in lu_non_invertible()
95 LU<MatrixType> lu(m1); in lu_invertible() local
96 VERIFY(0 == lu.dimensionOfKernel()); in lu_invertible()
97 VERIFY(size == lu.rank()); in lu_invertible()
98 VERIFY(lu.isInjective()); in lu_invertible()
99 VERIFY(lu.isSurjective()); in lu_invertible()
100 VERIFY(lu.isInvertible()); in lu_invertible()
101 VERIFY(lu.image().lu().isInvertible()); in lu_invertible()
103 lu.solve(m3, &m2); in lu_invertible()
105 VERIFY_IS_APPROX(m2, lu.inverse()*m3); in lu_invertible()
107 VERIFY(lu.solve(m3, &m2)); in lu_invertible()