findPivot() [3/3]

template<class RingType>

size_t DMatLUinPlace< RingType >::findPivot ( size_t row, size_t col )

private

Definition at line 171 of file dmat-lu-inplace.hpp.

{
  // Look at elements A[row,col], A[row+1,col], ..., A[nrows-1, col]
  // Return the index r s.y. abs(A[r,col]) is maximum over all of these
 
  for (size_t i = row; i < mLU.numRows(); i++)
    {
      if (!ring().is_zero(mLU.entry(i, col))) return i;
    }
  return static_cast<size_t>(-1);
}

References mLU, and ring().

Referenced by computeLU().