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,Data.size());
  for(std::vector<dcomponent>::const_iterator it=Data.begin(); it!=Data.end(); it++){
    newmat.fput(it->i, it->j, std::complex<double>(it->v,0.0));
  }
  
  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(M,N);
  newmat.zero();
  for(std::vector<dcomponent>::const_iterator it=Data.begin(); it!=Data.end(); it++){
    newmat(it->i, it->j) =it->v;
    newmat(it->j, it->i) =it->v;
  }
  
  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(N);
  newmat.zero();
  for(std::vector<dcomponent>::const_iterator it=Data.begin(); it!=Data.end(); it++){
    newmat(it->i, it->j) =it->v;
  }
  
  return _(newmat);
}

//============================================================================
/*! convert to _dgsmatrix */
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,Data.size()*2) );
  newmat.zero();
  for(std::vector<dcomponent>::const_iterator it=Data.begin(); it!=Data.end(); it++){
    newmat.fput(it->i, it->j, it->v);
    if(it->i!=it->j){
      newmat.fput(it->j, it->i, it->v);
    }
  }
  
  return _(newmat);
}