fill_lower_and_upper() [1/2]

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

Definition at line 154 of file lapack.cpp.

{
  auto& ring = lower.ring();
  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++)
        {
          double re = *U++;
          double im = *U++;
          ring.set_from_doubles(upper.entry(r, c), re, im);
          // upper.entry(r, c).re = *U++;
          // upper.entry(r, c).im = *U++;
        }
      ring.set_from_long(lower.entry(c, c), 1);
      for (size_t r = c+1 ; r <= lower.numColumns(); ++r)
        {
          double re = *U++;
          double im = *U++;
          ring.set_from_doubles(lower.entry(r, c), re, im);
          // lower.entry(r,c).re = *U++;
          // lower.entry(r,c).im = *U++;
        }
    }
}

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