fill_lower_and_upper() [2/2]

void fill_lower_and_upper	(	const std::vector< double > &	lapack_numbers,
		DMatRR &	lower,
		DMatRR &	upper )

Definition at line 126 of file lapack.cpp.

{
  int nrows = static_cast<int>(lower.numRows());
  int ncols = static_cast<int>(upper.numColumns());
  int min = static_cast<int>(lower.numColumns());
  assert(min == static_cast<int>(upper.numRows()));
         
  // At this point, lower and upper should be zero matrices.
  assert(MatrixOps::isZero(lower));
  assert(MatrixOps::isZero(upper));
 
  const double* U = lapack_numbers.data();
  for (size_t c = 0; c < upper.numColumns(); c++)
    {
      for (size_t r = 0; r <= c and r <= min; r++)
        upper.entry(r, c) = *U++;
      lower.entry(c, c) = 1;
      for (size_t r = c+1 ; r <= lower.numColumns(); ++r)
        lower.entry(r,c) = *U++;
    }
}

References DMat< ACoeffRing >::entry(), MatrixOps::isZero(), DMat< ACoeffRing >::numColumns(), and DMat< ACoeffRing >::numRows().

Referenced by Lapack::LU(), and Lapack::LU().