setUpperLower()

void DMatLinAlg< M2::ARingZZpFlint >::setUpperLower ( const Mat & LU, Mat & lower, Mat & upper )

private

Definition at line 108 of file dmat-lu.hpp.

{
  size_t min = std::min(LU.numRows(), LU.numColumns());
  lower.resize(LU.numRows(), min);
  upper.resize(min, LU.numColumns());
 
  // At this point, lower and upper should be zero matrices.
  assert(MatrixOps::isZero(lower));
  assert(MatrixOps::isZero(upper));
 
  for (size_t c = 0; c < LU.numColumns(); c++)
    {
      if (c < min) ring().set_from_long(lower.entry(c, c), 1);
      for (size_t r = 0; r < LU.numRows(); r++)
        {
          if (r <= c)
            ring().set(upper.entry(r, c), LU.entry(r, c));
          else if (c < lower.numRows())
            {
              ring().set(lower.entry(r, c), LU.entry(r, c));
            }
        }
    }
}