CPPLapack Code

//============================================================================
/*! zgbmatrix constructor */
inline zgbmatrix::zgbmatrix()
  : m(M), n(N), kl(KL), ku(KU), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgbmatrix::zgbmatrix()" << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  M =0;
  N =0;
  KL =0;
  KU =0;
  Array =NULL;
  Darray =NULL;
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgbmatrix::zgbmatrix() "
            << "A new 0x0:0,0 matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

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

//============================================================================
/*! zgbmatrix copy constructor */
inline zgbmatrix::zgbmatrix(const zgbmatrix& mat)
  : m(M), n(N), kl(KL), ku(KU), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgbmatrix::zgbmatrix(const zgbmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  M =mat.M;
  N =mat.N;
  KL =mat.KL;
  KU =mat.KU;
  Array =new std::complex<double>[(KL+KU+1)*N];
  Darray =new std::complex<double>*[N];
  for(int i=0; i<N; i++){ Darray[i] =&Array[i*(KL+KU+1)]; }

  //////// copy ////////
  zcopy_((KL+KU+1)*N, mat.Array, 1, Array, 1);
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgbmatrix::zgbmatrix(const zgbmatrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

//============================================================================
/*! zgbmatrix constructor to cast _zgbmatrix */
inline zgbmatrix::zgbmatrix(const _zgbmatrix& mat)
  : m(M), n(N), kl(KL), ku(KU), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgbmatrix::zgbmatrix(const _zgbmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  M =mat.M;
  N =mat.N;
  KL =mat.KL;
  KU =mat.KU;
  Array =mat.Array;
  Darray =mat.Darray;
  
  mat.nullify();
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgbmatrix::zgbmatrix(const zgbmatrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

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

//============================================================================
/*! zgbmatrix constructor with size specification */
inline zgbmatrix::zgbmatrix(const long& _m, const long& _n,
                            const long& _kl, const long& _ku)
  : m(M), n(N), kl(KL), ku(KU), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgbmatrix::zgbmatrix(const long&, const long&, const long&, const long&)" << std::endl;
#endif//CPPL_VERBOSE
  
#ifdef  CPPL_DEBUG
  if( _m<0 || _n<0 || _kl<0 || _ku<0 || _m<_kl || _n<_ku ){
    std::cerr << "[ERROR] zgbmatrix::zgbmatrix"
              << "(const long, const long, const long const long)" << std::endl
              << "It is impossible to make a matrix you ordered. " << std::endl
              << "Your input was (" << _m << "," << _n << ","<< _ku << "," << _kl << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  //////// initialize ////////
  M =_m;
  N =_n;
  KL =_kl;
  KU =_ku;
  Array =new std::complex<double>[(KL+KU+1)*N];
  Darray =new std::complex<double>*[N];
  for(int i=0; i<N; i++){ Darray[i] =&Array[i*(KL+KU+1)]; }

#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgbmatrix::zgbmatrix(const long&, const long&, const long&,const long&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

//============================================================================
/*! zgbmatrix constructor with filename */
inline zgbmatrix::zgbmatrix(const char* filename)
  : m(M), n(N), kl(KL), ku(KU), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgbmatrix::zgbmatrix(const char*)" << std::endl;
#endif//CPPL_VERBOSE
  
  Array =NULL;
  Darray =NULL;
  
  //// read ////
  read(filename);
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgbmatrix::zgbmatrix(const char*) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

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

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