CPPLapack Code

//============================================================================
/*! cast to _zgbmatrix */
inline _zgbmatrix _dgbmatrix::to_zgbmatrix() const
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] _dgbmatrix::to_zgbmatrix() const" << std::endl;
#endif//CPPL_VERBOSE
  
  zgbmatrix newmat(M,N,KL,KU);
  for(long i=0; i<(KL+KU+1)*N; i++){
    newmat.array[i] =std::complex<double>(Array[i],0.0);
  }
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] _dgbmatrix::to_zgbmatrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
  
  destroy();
  return _(newmat);
}

///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////

//============================================================================
/*! convert to _dgematrix */
inline _dgematrix _dgbmatrix::to_dgematrix() const
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] _dgbmatrix::to_dgematrix() const" << std::endl;
#endif//CPPL_VERBOSE

  dgematrix newmat( dgematrix(M,N).zero() );
  for(long i=0; i<M; i++){
    for(long j=std::max(long(0),i-KL); j<std::min(N,i+KU+1); j++){
      newmat(i,j) =(*this)(i,j);
    }
  }

#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] _dgbmatrix::to_dgematrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_DEBUG
  destroy();
  return _(newmat);
}