CPPLapack Code

Status: Beta
Brought to you by: yonishi
[r44]: / trunk / include / zgematrix- / zgematrix-constructor.hpp Maximize Restore History

404 lines (341 with data), 13.0 kB

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

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

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

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

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

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

//============================================================================
/*! zgematrix constructor to cast zgbmatrix */
inline zgematrix::zgematrix(const zgbmatrix& mat)
  : m(M), n(N), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgematrix::zgematrix(const zgbmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  M =mat.M;
  N =mat.N;
  Array =new std::complex<double>[M*N];
  Darray =new std::complex<double>*[N];
  for(int i=0; i<N; i++){ Darray[i] =&Array[i*M]; }
  
  //////// copy ////////
  zero();
  for(long i=0; i<mat.M; i++){
    for(long j=max(0,i-mat.KL); j<min(N,i+mat.KU+1); j++){
      operator()(i,j) =mat(i,j);
    }
  }
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgematrix::zgematrix(const zgbmatrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

//============================================================================
/*! zgematrix constructor to cast _zgbmatrix */
inline zgematrix::zgematrix(const _zgbmatrix& mat)
  : m(M), n(N), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgematrix::zgematrix(const _zgbmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  M =mat.M;
  N =mat.N;
  Array =new std::complex<double>[M*N];
  Darray =new std::complex<double>*[N];
  for(int i=0; i<N; i++){ Darray[i] =&Array[i*M]; }
  
  //////// copy ////////
  zero();
  for(long i=0; i<mat.M; i++){
    for(long j=max(0,i-mat.KL); j<min(N,i+mat.KU+1); j++){
      operator()(i,j) =mat(i,j);
    }
  }
  
  mat.destroy();
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgematrix::zgematrix(const zgbmatrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

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

//============================================================================
/*! zgematrix constructor to cast zhematrix */
inline zgematrix::zgematrix(const zhematrix& mat)
  : m(M), n(N), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgematrix::zgematrix(const zhematrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  M =mat.N;
  N =mat.N;
  Array =new std::complex<double>[N*N];
  Darray =new std::complex<double>*[N];
  for(int i=0; i<N; i++){ Darray[i] =&Array[i*M]; }
  
  //////// copy ////////
  zero();
  for(long i=0; i<mat.N; i++){
    for(long j=0; j<mat.N; j++ ){
      operator()(i,j) =mat(i,j);
    }
  }
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgematrix::zgematrix(const zhematrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

//============================================================================
/*! zgematrix constructor to cast _zhematrix */
inline zgematrix::zgematrix(const _zhematrix& mat)
  : m(M), n(N), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgematrix::zgematrix(const _zhematrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  M =mat.N;
  N =mat.N;
  Array =mat.Array;
  Darray =mat.Darray;
  //Darray =new std::complex<double>*[N];
  //for(int i=0; i<N; i++){ Darray[i] =&Array[i*M]; }
  
  //////// self copy to upper right ////////
  for(long i=0; i<mat.N; i++){ for(long j=0; j<i; j++){
    operator()(j,i) =std::conj(operator()(i,j));
    }}
  
  mat.nullify();
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgematrix::zgematrix(const _zhematrix&) "
            << "A matrix at " << Array << " has been transformed into " << "zgematrix." << std::endl;
#endif//CPPL_VERBOSE
}

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

//============================================================================
/*! zgematrix constructor to cast zgsmatrix */
inline zgematrix::zgematrix(const zgsmatrix& mat)
  : m(M), n(N), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgematrix::zgematrix(const zgsmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  M =mat.M;
  N =mat.N;
  Array =new std::complex<double>[M*N];
  Darray =new std::complex<double>*[N];
  for(int i=0; i<N; i++){ Darray[i] =&Array[i*M]; }
  
  //////// copy ////////
  zero();
  for(long c=0; c<mat.VOL; c++) {
    operator()(mat.Indx[c],mat.Jndx[c]) =mat.Array[c];
  }
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgematrix::zgematrix(const zgsmatrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

//============================================================================
/*! zgematrix constructor to cast _zgsmatrix */
inline zgematrix::zgematrix(const _zgsmatrix& mat)
  : m(M), n(N), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgematrix::zgematrix(const _zhematrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  M =mat.M;
  N =mat.N;
  Array =new std::complex<double>[M*N];
  Darray =new std::complex<double>*[N];
  for(int i=0; i<N; i++){ Darray[i] =&Array[i*M]; }
  
  //////// copy ////////
  zero();
  for(long c=0; c<mat.VOL; c++){
    operator()(mat.Indx[c],mat.Jndx[c]) =mat.Array[c];
  }
  
  mat.destroy();
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgematrix::zgematrix(const _zgsmatrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

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

//============================================================================
/*! zgematrix constructor to cast zhsmatrix */
inline zgematrix::zgematrix(const zhsmatrix& mat)
  : m(M), n(N), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgematrix::zgematrix(const zhsmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  M =mat.N;
  N =mat.N;
  Array =new std::complex<double>[M*N];
  Darray =new std::complex<double>*[N];
  for(int i=0; i<N; i++){ Darray[i] =&Array[i*M]; }
  
  //////// copy ////////
  zero();
  for(long c=0; c<mat.VOL; c++) {
    operator()(mat.Indx[c],mat.Jndx[c]) =mat.Array[c];
    operator()(mat.Jndx[c],mat.Indx[c]) =std::conj(mat.Array[c]);
  }
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgematrix::zgematrix(const zhsmatrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

//============================================================================
/*! zgematrix constructor to cast _zhsmatrix */
inline zgematrix::zgematrix(const _zhsmatrix& mat)
  : m(M), n(N), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgematrix::zgematrix(const _zhematrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  //////// initialize ////////
  M =mat.N;
  N =mat.N;
  Array =new std::complex<double>[M*N];
  Darray =new std::complex<double>*[N];
  for(int i=0; i<N; i++){ Darray[i] =&Array[i*M]; }
  
  //////// copy ////////
  zero();
  for(long c=0; c<mat.VOL; c++){
    operator()(mat.Indx[c],mat.Jndx[c]) =mat.Array[c];
    operator()(mat.Jndx[c],mat.Indx[c]) =std::conj(mat.Array[c]);
  }
  
  mat.destroy();
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgematrix::zgematrix(const _zhsmatrix&) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

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

//============================================================================
/*! zgematrix constructor with size specification */
inline zgematrix::zgematrix(const long& _m, const long& _n)
  : m(M), n(N), array(Array), darray(Darray)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] zgematrix::zgematrix(const long&, const long&)" << std::endl;
#endif//CPPL_VERBOSE
  
#ifdef  CPPL_DEBUG
  if( _m<0 || _n<0 ){
    std::cerr << "[ERROR] zgematrix::zgematrix(const long, const long)" << std::endl
              << "Matrix sizes must be positive integers. " << std::endl
              << "Your input was (" << _m << "," << _n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  //////// initialize ////////
  M =_m;
  N =_n;
  Array =new std::complex<double>[M*N];
  Darray =new std::complex<double>*[N];
  for(int i=0; i<N; i++){ Darray[i] =&Array[i*M]; }
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] zgematrix(long,long) "
            << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}

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

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

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