fill_from_lapack_upper() [1/2]

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

Definition at line 100 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 += 2*numrows)
    for (size_t r = 0; r <= c; ++r)
      {
        if (r >= upper.numRows()) break;
        upper.ring().set_from_doubles(upper.entry(r, c), U[2*r], U[2*r+1]);
      }
}

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