Logo Search packages:      
Sourcecode: feel++ version File versions  Download package

template<typename _MatrixType>
Index FullPivLU< _MatrixType >::rank (  ) const [inline]
the rank of the matrix of which *this is the LU decomposition.
This method has to determine which pivots should be considered nonzero. For that, it uses the threshold value that you can control by calling setThreshold(const RealScalar&).

Definition at line 306 of file FullPivLU.h.

References FullPivLU< _MatrixType >::threshold().

Referenced by FullPivLU< _MatrixType >::dimensionOfKernel(), FullPivLU< _MatrixType >::isInjective(), and FullPivLU< _MatrixType >::isSurjective().

      ei_assert(m_isInitialized && "LU is not initialized.");
      RealScalar premultiplied_threshold = ei_abs(m_maxpivot) * threshold();
      Index result = 0;
      for(Index i = 0; i < m_nonzero_pivots; ++i)
        result += (ei_abs(m_lu.coeff(i,i)) > premultiplied_threshold);
      return result;

Here is the call graph for this function:

Here is the caller graph for this function:

Generated by  Doxygen 1.6.0   Back to index