CPPLapack Code

//============================================================================
/*! cast to _zhsmatrix */
inline _zhsmatrix dssmatrix::to_zhsmatrix() const
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dssmatrix::to_zhsmatrix() const" << std::endl;
#endif//CPPL_VERBOSE
  
  zhsmatrix newmat(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] dssmatrix::to_zhsmatrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_DEBUG
  
  return _(newmat);
}

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

//============================================================================
/*! convert to _dgematrix */
inline _dgematrix dssmatrix::to_dgematrix() const
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dssmatrix::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];
    newmat(Jndx[c],Indx[c]) =Array[c];
  }
  
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dssmatrix::to_dgematrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_DEBUG

  return _(newmat);
}

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

  dsymatrix newmat( dsymatrix(N).zero() );
  for(long c=0; c<VOL; c++){
    newmat(Indx[c],Jndx[c]) =Array[c];
  }
  
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dssmatrix::to_dsymatrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_DEBUG

  return _(newmat);
}

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

  dgsmatrix newmat( dgsmatrix(M,N,VOL*2).zero() );
  for(long c=0; c<VOL; c++){
    newmat.fput(Indx[c],Jndx[c], Array[c]);
    newmat.fput(Jndx[c],Indx[c], Array[c]);
  }
  
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dssmatrix::to_dgsmatrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_DEBUG

  return _(newmat);
}