CPPLapack Code

//============================================================================
/*! cast to _zgsmatrix */
inline _zgsmatrix dgsmatrix::to_zgsmatrix() const
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dgsmatrix::to_zgsmatrix() const" << std::endl;
#endif//CPPL_VERBOSE
  
  zgsmatrix newmat(M,N,CAP);
  for(long c=0; c<VOL; c++){
    newmat.fput(Indx[c], Jndx[c], std::complex<double>(Array[c],0.0));
  }
  
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dgsmatrix::to_zgsmatrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_DEBUG
  
  return _(newmat);
}

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

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

  dgematrix newmat( dgematrix(M,N).zero() );
  for(long c=0; c<VOL; c++){
    newmat(Indx[c],Jndx[c]) =Array[c];
  }

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

  return _(newmat);
}