107       return m_lu;  in matrixLU()
286                                       : NumTraits<Scalar>::epsilon() * m_lu.diagonalSize();  in threshold()
302         result += (abs(m_lu.coeff(i,i)) > premultiplied_threshold);  in rank()
353       return isInjective() && (m_lu.rows() == m_lu.cols());  in isInvertible()
366 …  eigen_assert(m_lu.rows() == m_lu.cols() && "You can't take the inverse of a non-square matrix!");  in inverse()
368                (*this, MatrixType::Identity(m_lu.rows(), m_lu.cols()));  in inverse()
373     inline Index rows() const { return m_lu.rows(); }  in rows()
374     inline Index cols() const { return m_lu.cols(); }  in cols()
377     MatrixType m_lu;
395   : m_lu(rows, cols),  in FullPivLU()
407   : m_lu(matrix.rows(), matrix.cols()),  in FullPivLU()
425   m_lu = matrix;  in compute()
447     biggest_in_corner = m_lu.bottomRightCorner(rows-k, cols-k)  in compute()
474       m_lu.row(k).swap(m_lu.row(row_of_biggest_in_corner));  in compute()
478       m_lu.col(k).swap(m_lu.col(col_of_biggest_in_corner));  in compute()
486       m_lu.col(k).tail(rows-k-1) /= m_lu.coeff(k,k);  in compute()
488 …m_lu.block(k+1,k+1,rows-k-1,cols-k-1).noalias() -= m_lu.col(k).tail(rows-k-1) * m_lu.row(k).tail(c…  in compute()
510 …eigen_assert(m_lu.rows() == m_lu.cols() && "You can't take the determinant of a non-square matrix!…  in determinant()
511   return Scalar(m_det_pq) * Scalar(m_lu.diagonal().prod());  in determinant()
521   const Index smalldim = (std::min)(m_lu.rows(), m_lu.cols());  in reconstructedMatrix()
523   MatrixType res(m_lu.rows(),m_lu.cols());  in reconstructedMatrix()
525   res = m_lu.leftCols(smalldim)  in reconstructedMatrix()
527       * m_lu.topRows(smalldim)  in reconstructedMatrix()