fill_from_lapack_upper() [2/2]

void fill_from_lapack_upper	(	const std::vector< double > &	lapack_numbers,
		int	numrows,
		int	numcols,
		DMatRR &	upper )

Definition at line 74 of file lapack.cpp.

{
  (void) numcols;
  // At this point, upper should be a zero matrix of size (min(numrows, numcols) x numcols)
  assert(upper.numRows() == std::min(numrows, numcols));
  assert(upper.numColumns() == numcols);
  assert(MatrixOps::isZero(upper));
 
  const double* U = lapack_numbers.data();
  for (size_t c = 0; c < upper.numColumns(); c++, U += numrows)
    for (size_t r = 0; r <= c; ++r)
      {
        if (r >= upper.numRows()) break;
        upper.entry(r, c) = U[r];
      }
}

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

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