CPPLapack Code

//=============================================================================
/*! convert to _zgematrix */
inline _zgematrix zgbmatrix::to_zgematrix() const
{CPPL_VERBOSE_REPORT;
  zgematrix newmat( zgematrix(m,n).zero() );
  
  for(long i=0; i<m; i++){
    const long jmax =std::min(n,i+ku+1);
    for(long j=std::max(0l,i-kl); j<jmax; j++){
      newmat(i,j) =(*this)(i,j);
    }
  }
  
  return _(newmat);
}