CPPLapack Code

//============================================================================
/*! dssmatrix constructor without arguments */
inline dssmatrix::dssmatrix()
  : n(N), cap(CAP), vol(VOL), array(Array), indx(Indx), jndx(Jndx), row(Row)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dssmatrix::dssmatrix()"
            << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  N =0;
  CAP =VOL =0;
  Array =NULL;
  Indx =Jndx =NULL;
  Row =NULL;
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dssmatrix::dssmatrix() "
            << "A new 0x0 matrix at " << Array
            << " has been made." << std::endl;
#endif//CPPL_DEBUG
}

//============================================================================
/*! dssmatrix copy constructor */
inline dssmatrix::dssmatrix(const dssmatrix& mat)
  : n(N), cap(CAP), vol(VOL), array(Array), indx(Indx), jndx(Jndx), row(Row)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dssmatrix::dssmatrix(const dssmatrix&)"
            << std::endl;
#endif//CPPL_VERBOSE

  Array =NULL;
  Indx =Jndx =NULL;
  Row =NULL;
  copy(mat);
  
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dssmatrix::dssmatrix(const dssmatrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_DEBUG
}

//============================================================================
/*! dssmatrix constructor to cast _dssmatrix */
inline dssmatrix::dssmatrix(const _dssmatrix& mat)
  : n(N), cap(CAP), vol(VOL), array(Array), indx(Indx), jndx(Jndx), row(Row)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dssmatrix::dssmatrix(const _dssmatrix&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  N =mat.N;
  CAP =mat.CAP;
  VOL =mat.VOL;
  Array =mat.Array;
  Indx =mat.Indx;
  Jndx =mat.Jndx;
  Row =mat.Row;
  
  mat.nullify();
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dssmatrix::dssmatrix(const _dssmatrix&) "
            << "A new matrix pointing at " << Array << " has been made."
            << std::endl;
#endif//CPPL_DEBUG
}

//============================================================================
/*! dssmatrix constructor with size specification */
inline dssmatrix::dssmatrix(const long& _n, const long& _c)
  : n(N), cap(CAP), vol(VOL), array(Array), indx(Indx), jndx(Jndx), row(Row)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dssmatrix::dssmatrix(const long&, const long&, const long&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
#ifdef  CPPL_DEBUG
  if( _n<0 || _c<0 ){
    std::cerr << "[ERROR] dssmatrix::dssmatrix"
              << "(const long&, const long&)" << std::endl
              << "Matrix sizes and the length of arrays "
              << "must be positive integers. " << std::endl
              << "Your input was (" << _n << "," << _c << ")."
              << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  //////// initialize ////////
  N =_n;
  CAP =_c;
  VOL =0;
  Array =new double[CAP];
  Indx =new long[CAP];
  Jndx =new long[CAP];
  Row =new std::vector< std::pair<long,long> >[N];
  
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dssmatrix(const long&, const long&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_DEBUG
}

//============================================================================
/*! dssmatrix constructor with filename */
inline dssmatrix::dssmatrix(const char* filename)
  : n(N), cap(CAP), vol(VOL), array(Array), indx(Indx), jndx(Jndx), row(Row)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dssmatrix::dssmatrix(const char*)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  Array =NULL;
  Indx =Jndx =NULL;
  Row =NULL;
  
  //// read ////
  read(filename);
  
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dssmatrix::dssmatrix(const char*) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_DEBUG
}

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

//============================================================================
/*! dssmatrix constructor from dgematrix */
inline dssmatrix::dssmatrix(const dsymatrix& mat)
  : n(N), cap(CAP), vol(VOL), array(Array), indx(Indx), jndx(Jndx), row(Row)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dssmatrix::dssmatrix(const dgematrix&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  N =mat.N;
  CAP =CPPL_SS_SECTOR;
  VOL =0;
  Array =new double[CAP];
  Indx =new long[CAP];
  Jndx =new long[CAP];
  Row =new std::vector< std::pair<long,long> >[N];
  
  //////// copy ////////
  for(long i=0; i<N; i++){
    for(long j=0; j<=i; j++){
      if(mat(i,j)!=0.0){ put(i,j,mat(i,j)); }
    }
  }
  
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dssmatrix::dssmatrix(const dssmatrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_DEBUG
}

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

//============================================================================
/*! dssmatrix destructor */
inline dssmatrix::~dssmatrix()
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dssmatrix::~dssmatrix()"
            << std::endl;
#endif//CPPL_VERBOSE
  
#ifdef  CPPL_DEBUG
  std::cerr << "# [NOTE] dssmatrix::~dssmatrix() "
            << "A matrix at " << Array << " is going to be deleted."
            << std::endl;
#endif//CPPL_DEBUG
  
  //////// delete Array ////////
  delete [] Array;
  delete [] Indx;
  delete [] Jndx;
  delete [] Row;
}