_dgbmatrix Class Reference

(DO NOT USE) Smart-temporary Real Double-precision General Band Matrix Class More...

#include <_dgbmatrix.hpp>

Public Member Functions
	_dgbmatrix ()
	_dgbmatrix (const _dgbmatrix &)
	~_dgbmatrix ()
_zgbmatrix	to_zgbmatrix () const
_dgematrix	to_dgematrix () const
double &	operator() (const CPPL_INT &, const CPPL_INT &) const
void	write (const char *) const
void	nullify () const
void	destroy () const

Public Attributes
CPPL_INT	m
	matrix row size More...
CPPL_INT	n
	matrix column size More...
CPPL_INT	kl
	lower band width More...
CPPL_INT	ku
	upper band width More...
double *	array
	1D array to store matrix data More...
double **	darray
	array of pointers of column head addresses More...

Friends
std::ostream &	operator<< (std::ostream &, const _dgbmatrix &)
_dgbmatrix	t (const _dgbmatrix &)
_dgematrix	i (const _dgbmatrix &)
const _dgbmatrix &	operator+ (const _dgbmatrix &)
_dgbmatrix	operator- (const _dgbmatrix &)
_dgematrix	operator+ (const _dgbmatrix &, const dgematrix &)
_dgematrix	operator+ (const _dgbmatrix &, const _dgematrix &)
_dgematrix	operator+ (const _dgbmatrix &, const dsymatrix &)
_dgematrix	operator+ (const _dgbmatrix &, const _dsymatrix &)
_dgbmatrix	operator+ (const _dgbmatrix &, const dgbmatrix &)
_dgbmatrix	operator+ (const _dgbmatrix &, const _dgbmatrix &)
_dgematrix	operator+ (const _dgbmatrix &, const dgsmatrix &)
_dgematrix	operator+ (const _dgbmatrix &, const _dgsmatrix &)
_dgematrix	operator+ (const _dgbmatrix &, const dssmatrix &)
_dgematrix	operator+ (const _dgbmatrix &, const _dssmatrix &)
_dgematrix	operator- (const _dgbmatrix &, const dgematrix &)
_dgematrix	operator- (const _dgbmatrix &, const _dgematrix &)
_dgematrix	operator- (const _dgbmatrix &, const dsymatrix &)
_dgematrix	operator- (const _dgbmatrix &, const _dsymatrix &)
_dgbmatrix	operator- (const _dgbmatrix &, const dgbmatrix &)
_dgbmatrix	operator- (const _dgbmatrix &, const _dgbmatrix &)
_dgematrix	operator- (const _dgbmatrix &, const dgsmatrix &)
_dgematrix	operator- (const _dgbmatrix &, const _dgsmatrix &)
_dgematrix	operator- (const _dgbmatrix &, const dssmatrix &)
_dgematrix	operator- (const _dgbmatrix &, const _dssmatrix &)
_dcovector	operator* (const _dgbmatrix &, const dcovector &)
_dcovector	operator* (const _dgbmatrix &, const _dcovector &)
_dgematrix	operator* (const _dgbmatrix &, const dgematrix &)
_dgematrix	operator* (const _dgbmatrix &, const _dgematrix &)
_dgematrix	operator* (const _dgbmatrix &, const dsymatrix &)
_dgematrix	operator* (const _dgbmatrix &, const _dsymatrix &)
_dgbmatrix	operator* (const _dgbmatrix &, const dgbmatrix &)
_dgbmatrix	operator* (const _dgbmatrix &, const _dgbmatrix &)
_dgematrix	operator* (const _dgbmatrix &, const dgsmatrix &)
_dgematrix	operator* (const _dgbmatrix &, const _dgsmatrix &)
_dgematrix	operator* (const _dgbmatrix &, const dssmatrix &)
_dgematrix	operator* (const _dgbmatrix &, const _dssmatrix &)
_dgbmatrix	operator* (const _dgbmatrix &, const double &)
_dgbmatrix	operator/ (const _dgbmatrix &, const double &)
_dgbmatrix	operator* (const double &, const _dgbmatrix &)

Detailed Description

(DO NOT USE) Smart-temporary Real Double-precision General Band Matrix Class

Definition at line 3 of file _dgbmatrix.hpp.

Constructor & Destructor Documentation

_dgbmatrix::_dgbmatrix ( )

inline

_dgbmatrix constructor

Definition at line 3 of file _dgbmatrix-constructor.hpp.

References array, darray, kl, ku, m, and n.

{CPPL_VERBOSE_REPORT;
  m =0;
  n =0;
  kl =0;
  ku =0;
  array =NULL;
  darray =NULL;
}

_dgbmatrix::_dgbmatrix ( const _dgbmatrix &  mat )

inline

_dgbmatrix copy constructor

Definition at line 15 of file _dgbmatrix-constructor.hpp.

References array, darray, kl, ku, m, n, and nullify().

{CPPL_VERBOSE_REPORT;
  m =mat.m;
  n =mat.n;
  kl =mat.kl;
  ku =mat.ku;
  array =mat.array;
  darray =mat.darray;
  
  mat.nullify();
}

_dgbmatrix::~_dgbmatrix ( )

inline

_dgbmatrix destructor

Definition at line 33 of file _dgbmatrix-constructor.hpp.

References array, and darray.

{CPPL_VERBOSE_REPORT;
  delete[] array;
  delete[] darray;
}

Member Function Documentation

_zgbmatrix _dgbmatrix::to_zgbmatrix ( ) const

inline

cast to _zgbmatrix

Definition at line 3 of file _dgbmatrix-cast.hpp.

References _(), zgbmatrix::array, array, destroy(), i, kl, ku, m, and n.

{CPPL_VERBOSE_REPORT;
  zgbmatrix newmat(m,n,kl,ku);
  
  for(CPPL_INT i=0; i<(kl+ku+1)*n; i++){
    newmat.array[i] =comple(array[i],0.0);
  }
  
  destroy();
  return _(newmat);
}

_dgematrix _dgbmatrix::to_dgematrix ( ) const

inline

convert to _dgematrix

Definition at line 21 of file _dgbmatrix-cast.hpp.

References _(), destroy(), i, kl, ku, m, and n.

{CPPL_VERBOSE_REPORT;
  dgematrix newmat( dgematrix(m,n).zero() );
  
  for(CPPL_INT i=0; i<m; i++){
    const CPPL_INT jmax =std::min(n,i+ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-kl); j<jmax; j++){
      newmat(i,j) =(*this)(i,j);
    }
  }
  
  destroy();
  return _(newmat);
}

double & _dgbmatrix::operator() ( const CPPL_INT &  i, const CPPL_INT &  j  ) const

inline

operator() for const object

Definition at line 3 of file _dgbmatrix-io.hpp.

References darray, i, kl, ku, m, and n.

Referenced by write().

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if( i<0 || j<0 || m<=i || n<=j || i-j>kl || j-i>ku ){
    ERROR_REPORT;
    std::cerr << "The required component is out of the matrix size." << std::endl
              << "Your input is (" << i << "," << j << "), whereas the matrix size is (" << m << "," << n << " , " << kl << "," << ku << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  //return array[ku+i+(kl+ku)*j];
  return darray[j][ku-j+i];
}

void _dgbmatrix::write ( const char *  filename ) const

inline

Definition at line 42 of file _dgbmatrix-io.hpp.

References destroy(), i, kl, ku, m, n, and operator()().

{CPPL_VERBOSE_REPORT;
  std::ofstream ofs(filename, std::ios::trunc);
  ofs.setf(std::cout.flags());
  ofs.precision(std::cout.precision());
  ofs.width(std::cout.width());
  ofs.fill(std::cout.fill());
  
  ofs << "#dgbmatrix" << " " << m << " " << n << " " << kl << " " << ku << std::endl;
  for(CPPL_INT i=0; i<m; i++){
    const CPPL_INT jmax =std::min(n,i+ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-kl); j<jmax; j++){
      ofs << operator()(i,j) << " ";
    }
    ofs << std::endl;
  }
  
  ofs.close();
  destroy();
}

void _dgbmatrix::nullify ( ) const

inline

nullify all the matrix data

Definition at line 3 of file _dgbmatrix-misc.hpp.

References array, darray, kl, ku, m, and n.

Referenced by _dgbmatrix(), dgbmatrix::dgbmatrix(), and dgbmatrix::shallow_copy().

{CPPL_VERBOSE_REPORT;
  m=0;
  n=0;
  kl=0;
  ku=0;
  array=NULL;
  darray=NULL;
}

void _dgbmatrix::destroy ( ) const

inline

destroy all the matrix data

Definition at line 15 of file _dgbmatrix-misc.hpp.

References array, and darray.

Referenced by operator*(), dgbmatrix::operator*=(), dgematrix::operator*=(), operator+(), dgbmatrix::operator+=(), dgematrix::operator+=(), operator-(), dgbmatrix::operator-=(), dgematrix::operator-=(), operator<<(), t(), to_dgematrix(), to_zgbmatrix(), and write().

{CPPL_VERBOSE_REPORT;
  delete [] array;
  delete [] darray;
  array=NULL;
  darray=NULL;
}

Friends And Related Function Documentation

std::ostream& operator<< ( std::ostream &  s, const _dgbmatrix &  mat  )

friend

Definition at line 23 of file _dgbmatrix-io.hpp.

{CPPL_VERBOSE_REPORT;
  for(CPPL_INT i=0; i<mat.m; i++){
    for(CPPL_INT j=0; j<mat.n; j++){
      if( i-j>mat.kl || j-i>mat.ku ){ s << " x"; }
      else{ s << " " << mat(i,j); }
    }
    s << std::endl;
  }
  
  mat.destroy();
  return s;
}

_dgbmatrix t ( const _dgbmatrix &  mat )

friend

return its transposed dgbmatrix

Definition at line 3 of file _dgbmatrix-calc.hpp.

{CPPL_VERBOSE_REPORT;
  dgbmatrix newmat(mat.n, mat.m, mat.ku, mat.kl);
  
  for(CPPL_INT i=0; i<newmat.m; i++){
    const CPPL_INT jmax =std::min(newmat.n,i+newmat.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-newmat.kl); j<jmax; j++){
      newmat(i,j) =mat(j,i);
    }
  }
  
  mat.destroy();
  return _(newmat);
}

_dgematrix i ( const _dgbmatrix &  mat )

friend

return its inverse matrix

Definition at line 20 of file _dgbmatrix-calc.hpp.

Referenced by operator()(), to_dgematrix(), to_zgbmatrix(), and write().

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(mat.m!=mat.n){
    ERROR_REPORT;
    std::cerr << "This matrix is not square and has no inverse matrix." << std::endl
              << "Your input was (" << mat.m << "x" << mat.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgbmatrix mat_cp(mat);
  dgematrix mat_inv(mat_cp.m,mat_cp.n);
  mat_inv.identity();
  mat_cp.dgbsv(mat_inv);

  return _(mat_inv);
}

const _dgbmatrix& operator+ ( const _dgbmatrix &  mat )

friend

+_dgbmatrix operator

Definition at line 3 of file _dgbmatrix-unary.hpp.

{CPPL_VERBOSE_REPORT;
  return mat;
}

_dgbmatrix operator- ( const _dgbmatrix &  mat )

friend

-_dgbmatrix operator

Definition at line 10 of file _dgbmatrix-unary.hpp.

{CPPL_VERBOSE_REPORT;
  for(CPPL_INT i=0; i<(mat.kl+mat.ku+1)*mat.n; i++){
    mat.array[i]=-mat.array[i];
  }
  
  return mat;
}

_dgematrix operator+ ( const _dgbmatrix &  matA, const dgematrix &  matB  )

friend

_dgbmatrix+dgematrix operator

Definition at line 3 of file _dgbmatrix-dgematrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG

  dgematrix newmat(matB);
  for(CPPL_INT i=0; i<matA.m; i++){
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      newmat(i,j) += matA(i,j);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

_dgematrix operator+ ( const _dgbmatrix &  matA, const _dgematrix &  matB  )

friend

_dgbmatrix+_dgematrix operator

Definition at line 3 of file _dgbmatrix-_dgematrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  for(CPPL_INT i=0; i<matA.m; i++){
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      matB(i,j) += matA(i,j);
    }
  }
  
  matA.destroy();
  return matB;
}

_dgematrix operator+ ( const _dgbmatrix &  matA, const dsymatrix &  matB  )

friend

_dgbmatrix+dsymatrix operator

Definition at line 3 of file _dgbmatrix-dsymatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.m!=matB.n || matA.n!=matB.n){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.n << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat(matB.to_dgematrix());
  for(CPPL_INT i=0; i<matA.m; i++){
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      newmat(i,j) += matA(i,j);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

_dgematrix operator+ ( const _dgbmatrix &  matA, const _dsymatrix &  matB  )

friend

_dgbmatrix+_dsymatrix operator

Definition at line 3 of file _dgbmatrix-_dsymatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.m!=matB.n || matA.n!=matB.n){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.n << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG

  dgematrix newmat(matB.n, matB.n);
  for(CPPL_INT i=0; i<matA.m; i++){
    for(CPPL_INT j=0; j<matB.n; j++){
      newmat(i,j) = matB(i,j);
    }
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      newmat(i,j) += matA(i,j);
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgbmatrix operator+ ( const _dgbmatrix &  matA, const dgbmatrix &  matB  )

friend

_dgbmatrix+dgbmatrix operator

Definition at line 3 of file _dgbmatrix-dgbmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  if(matA.kl>=matB.kl && matA.ku>=matB.ku){
    for(CPPL_INT i=0; i<matB.m; i++){
      const CPPL_INT jmax =std::min(matB.n,i+matB.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax; j++){
        matA(i,j) += matB(i,j);
      }
    }
    
    return matA;
  }
  
  else{
    dgbmatrix newmat(matA.m,matA.n,std::max(matA.kl,matB.kl),std::max(matA.ku,matB.ku));
    newmat.zero();
    
    for(CPPL_INT i=0; i<matA.m; i++){
      const CPPL_INT jmax1 =std::min(matA.n,i+matA.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax1; j++){
        newmat(i,j) += matA(i,j);
      }
      const CPPL_INT jmax2 =std::min(matB.n,i+matB.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax2; j++){
        newmat(i,j) += matB(i,j);
      }
    }
    
    matA.destroy();
    return _(newmat);
  }
}

_dgbmatrix operator+ ( const _dgbmatrix &  matA, const _dgbmatrix &  matB  )

friend

_dgbmatrix+_dgbmatrix operator

Definition at line 3 of file _dgbmatrix-_dgbmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  if(matA.kl>matB.kl && matA.ku>matB.ku){
    for(CPPL_INT i=0; i<matB.m; i++){
      const CPPL_INT jmax =std::min(matB.n,i+matB.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax; j++){
        matA(i,j) += matB(i,j);
      }
    }
    
    matB.destroy();
    return matA;
  }

  else if(matB.kl>matA.kl && matB.ku>matA.ku){
    for(CPPL_INT i=0; i<matA.m; i++){
      const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
        matB(i,j) += matA(i,j);
      }
    }
    
    matA.destroy();
    return matB;
  }
  
  else{
    dgbmatrix newmat(matA.m,matA.n,std::max(matA.kl,matB.kl),std::max(matA.ku,matB.ku));
    newmat.zero();
    
    for(CPPL_INT i=0; i<matA.m; i++){
      const CPPL_INT jmax1 =std::min(matA.n,i+matA.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax1; j++){
        newmat(i,j) += matA(i,j);
      }
      const CPPL_INT jmax2 =std::min(matB.n,i+matB.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax2; j++){
        newmat(i,j) += matB(i,j);
      }
    }
    
    matA.destroy();
    matB.destroy();
    return _(newmat);
  }
}

_dgematrix operator+ ( const _dgbmatrix &  matA, const dgsmatrix &  matB  )

friend

_dgbmatrix+dgsmatrix operator

Definition at line 3 of file _dgbmatrix-dgsmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( matB.to_dgematrix() );
  
  for(CPPL_INT i=0; i<matA.m; i++){
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      newmat(i,j)+=matA(i,j);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

_dgematrix operator+ ( const _dgbmatrix &  matA, const _dgsmatrix &  matB  )

friend

_dgbmatrix+_dgsmatrix operator

Definition at line 3 of file _dgbmatrix-_dgsmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( matB.to_dgematrix() );
  
  for(CPPL_INT i=0; i<matA.m; i++){
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      newmat(i,j)+=matA(i,j);
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgematrix operator+ ( const _dgbmatrix &  , const dssmatrix &    )

friend

_dgematrix operator+ ( const _dgbmatrix &  , const _dssmatrix &    )

friend

_dgematrix operator- ( const _dgbmatrix &  matA, const dgematrix &  matB  )

friend

_dgbmatrix-dgematrix operator

Definition at line 28 of file _dgbmatrix-dgematrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat(-matB);
  for(CPPL_INT i=0; i<matA.m; i++){
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      newmat(i,j) += matA(i,j);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

_dgematrix operator- ( const _dgbmatrix &  matA, const _dgematrix &  matB  )

friend

_dgbmatrix-_dgematrix operator

Definition at line 27 of file _dgbmatrix-_dgematrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  for(CPPL_INT i=0; i<matA.m; i++){
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      matB(i,j) = matA(i,j)-matB(i,j);
    }
  }
  
  matA.destroy();
  return matB;
}

_dgematrix operator- ( const _dgbmatrix &  matA, const dsymatrix &  matB  )

friend

_dgbmatrix-dsymatrix operator

Definition at line 28 of file _dgbmatrix-dsymatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.m!=matB.n || matA.n!=matB.n){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.n << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( (-matB).to_dgematrix() );
  for(CPPL_INT i=0; i<matA.m; i++){
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      newmat(i,j) += matA(i,j);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

_dgematrix operator- ( const _dgbmatrix &  matA, const _dsymatrix &  matB  )

friend

_dgbmatrix-_dsymatrix operator

Definition at line 32 of file _dgbmatrix-_dsymatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.m!=matB.n || matA.n!=matB.n){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.n << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat(matB.n, matB.n);
  for(CPPL_INT i=0; i<newmat.m; i++){
    for(CPPL_INT j=0; j<newmat.n; j++){
      newmat(i,j) = -matB(i,j);
    }
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      newmat(i,j) += matA(i,j);
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgbmatrix operator- ( const _dgbmatrix &  matA, const dgbmatrix &  matB  )

friend

_dgbmatrix-dgbmatrix operator

Definition at line 47 of file _dgbmatrix-dgbmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a subtraction." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") - (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  if(matA.kl>=matB.kl && matA.ku>=matB.ku){
    for(CPPL_INT i=0; i<matB.m; i++){
      const CPPL_INT jmax =std::min(matB.n,i+matB.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax; j++){
        matA(i,j) -= matB(i,j);
      }
    }
    
    return matA;
  }
  
  else{
    dgbmatrix newmat(matA.m,matA.n,std::max(matA.kl,matB.kl),std::max(matA.ku,matB.ku));
    newmat.zero();
    
    for(CPPL_INT i=0; i<matA.m; i++){
      const CPPL_INT jmax1 =std::min(matA.n,i+matA.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax1; j++){
        newmat(i,j) += matA(i,j);
      }
      const CPPL_INT jmax2 =std::min(matB.n,i+matB.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax2; j++){
        newmat(i,j) -= matB(i,j);
      }
    }
    
    return _(newmat);
  }
}

_dgbmatrix operator- ( const _dgbmatrix &  matA, const _dgbmatrix &  matB  )

friend

_dgbmatrix-_dgbmatrix operator

Definition at line 61 of file _dgbmatrix-_dgbmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a subtraction." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") - (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  if(matA.kl>matB.kl && matA.ku>matB.ku){
    for(CPPL_INT i=0; i<matB.m; i++){
      const CPPL_INT jmax =std::min(matB.n,i+matB.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax; j++){
        matA(i,j) -= matB(i,j);
      }
    }
    
    matB.destroy();
    return matA;
  }
  
  else{
    dgbmatrix newmat(matA.m,matA.n,std::max(matA.kl,matB.kl),std::max(matA.ku,matB.ku));
    newmat.zero();
    
    for(CPPL_INT i=0; i<matA.m; i++){
      const CPPL_INT jmax1 =std::min(matA.n,i+matA.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax1; j++){
        newmat(i,j) += matA(i,j);
      }
      const CPPL_INT jmax2 =std::min(matB.n,i+matB.ku+1);
      for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax2; j++){
        newmat(i,j) -= matB(i,j);
      }
    }
    
    matA.destroy();
    matB.destroy();
    return _(newmat);
  }
}

_dgematrix operator- ( const _dgbmatrix &  matA, const dgsmatrix &  matB  )

friend

_dgbmatrix-dgsmatrix operator

Definition at line 29 of file _dgbmatrix-dgsmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( (-matB).to_dgematrix() );
  
  for(CPPL_INT i=0; i<matA.m; i++){
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      newmat(i,j)-=matA(i,j);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

_dgematrix operator- ( const _dgbmatrix &  matA, const _dgsmatrix &  matB  )

friend

_dgbmatrix-_dgsmatrix operator

Definition at line 30 of file _dgbmatrix-_dgsmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n || matA.m!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a summation." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( (-matB).to_dgematrix() );
  
  for(CPPL_INT i=0; i<matA.m; i++){
    const CPPL_INT jmax =std::min(matA.n,i+matA.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-matA.kl); j<jmax; j++){
      newmat(i,j)-=matA(i,j);
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgematrix operator- ( const _dgbmatrix &  , const dssmatrix &    )

friend

_dgematrix operator- ( const _dgbmatrix &  , const _dssmatrix &    )

friend

_dcovector operator* ( const _dgbmatrix &  mat, const dcovector &  vec  )

friend

_dgbmatrix*dcovector operator

Definition at line 3 of file _dgbmatrix-dcovector.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(mat.n!=vec.l){
    ERROR_REPORT;
    std::cerr << "These matrix and vector can not make a product." << std::endl
              << "Your input was (" << mat.m << "x" << mat.n << ") * (" << vec.l << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dcovector newvec(mat.m);
  char trans ='n';
  double alpha =1.;
  CPPL_INT lda =mat.kl+mat.ku+1;
  CPPL_INT inc =1;
  double beta =0.;
  
  dgbmv_( &trans, &mat.m, &mat.n, &mat.kl, &mat.ku, &alpha, mat.array, &lda, vec.array, &inc, &beta, newvec.array, &inc );
  
  mat.destroy();
  return _(newvec);
}

_dcovector operator* ( const _dgbmatrix &  mat, const _dcovector &  vec  )

friend

_dgbmatrix*_dcovector operator

Definition at line 3 of file _dgbmatrix-_dcovector.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(mat.n!=vec.l){
    ERROR_REPORT;
    std::cerr << "These matrix and vector can not make a product." << std::endl
              << "Your input was (" << mat.m << "x" << mat.n << ") * (" << vec.l << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dcovector newvec(mat.m);
  char trans ='n';
  double alpha =1.;
  CPPL_INT lda =mat.kl+mat.ku+1;
  CPPL_INT inc =1;
  double beta =0.;
  
  dgbmv_( &trans, &mat.m, &mat.n, &mat.kl, &mat.ku, &alpha, mat.array, &lda, vec.array, &inc, &beta, newvec.array, &inc );
  
  mat.destroy();
  vec.destroy();
  return _(newvec);
}

_dgematrix operator* ( const _dgbmatrix &  matA, const dgematrix &  matB  )

friend

_dgbmatrix*dgematrix operator

Definition at line 53 of file _dgbmatrix-dgematrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a product." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( matA.m, matB.n );
  newmat.zero();
  
  for(CPPL_INT i=0; i<newmat.m; i++){
    for(CPPL_INT j=0; j<newmat.n; j++){
      const CPPL_INT kmax =std::min(matA.n,i+matA.ku+1);
      for(CPPL_INT k=std::max(CPPL_INT(0),i-matA.kl); k<kmax; k++){
        newmat(i,j) += matA(i,k)*matB(k,j);
      }
    }
  }
  
  matA.destroy();
  return _(newmat);
}

_dgematrix operator* ( const _dgbmatrix &  matA, const _dgematrix &  matB  )

friend

_dgbmatrix*_dgematrix operator

Definition at line 51 of file _dgbmatrix-_dgematrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a product." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( matA.m, matB.n );
  newmat.zero();
  
  for(CPPL_INT i=0; i<newmat.m; i++){
    for(CPPL_INT j=0; j<newmat.n; j++){
      const CPPL_INT kmax =std::min(matA.n,i+matA.ku+1);
      for(CPPL_INT k=std::max(CPPL_INT(0),i-matA.kl); k<kmax; k++){
        newmat(i,j) += matA(i,k)*matB(k,j);
      }
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgematrix operator* ( const _dgbmatrix &  matA, const dsymatrix &  matB  )

friend

_dgbmatrix*dgematrix operator

Definition at line 53 of file _dgbmatrix-dsymatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a product." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.n << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( matA.m, matB.n );
  newmat.zero();
  
  for(CPPL_INT i=0; i<newmat.m; i++){
    for(CPPL_INT j=0; j<newmat.n; j++){
      const CPPL_INT kmax =std::min(matA.n,i+matA.ku+1);
      for(CPPL_INT k=std::max(CPPL_INT(0),i-matA.kl); k<kmax; k++){
        newmat(i,j) += matA(i,k)*matB(k,j);
      }
    }
  }
  
  matA.destroy();
  return _(newmat);
}

_dgematrix operator* ( const _dgbmatrix &  matA, const _dsymatrix &  matB  )

friend

_dgbmatrix*_dsymatrix operator

Definition at line 61 of file _dgbmatrix-_dsymatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.n){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a product." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.n << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( matA.m, matB.n );
  newmat.zero();
  
  for(CPPL_INT i=0; i<newmat.m; i++){
    for(CPPL_INT j=0; j<newmat.n; j++){
      const CPPL_INT kmax =std::min(matA.n,i+matA.ku+1);
      for(CPPL_INT k=std::max(CPPL_INT(0),i-matA.kl); k<kmax; k++){
        newmat(i,j) += matA(i,k)*matB(k,j);
      }
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgbmatrix operator* ( const _dgbmatrix &  matA, const dgbmatrix &  matB  )

friend

_dgbmatrix*dgbmatrix operator

Definition at line 90 of file _dgbmatrix-dgbmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a product." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgbmatrix newmat( matA.m, matB.n, std::min(matA.kl+matB.kl,matA.m-1), std::min(matA.ku+matB.ku,matB.n-1) );
  newmat.zero();
  
  for(CPPL_INT i=0; i<newmat.m; i++){
    const CPPL_INT jmax =std::min(newmat.n,i+newmat.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-newmat.kl); j<jmax; j++){
      const CPPL_INT kmax =std::min( std::min(matA.n,i+matA.ku+1), std::min(matB.m,j+matB.kl+1) );
      for(CPPL_INT k=std::max( std::max(CPPL_INT(0),i-matA.kl), std::max(CPPL_INT(0),j-matB.ku) ); k<kmax; k++){
        newmat(i,j) += matA(i,k)*matB(k,j);
      }
    }
  }
  
  matA.destroy();
  return _(newmat);
}

_dgbmatrix operator* ( const _dgbmatrix &  matA, const _dgbmatrix &  matB  )

friend

_dgbmatrix*_dgbmatrix operator

Definition at line 107 of file _dgbmatrix-_dgbmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a product." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgbmatrix newmat( matA.m, matB.n, std::min(matA.kl+matB.kl,matA.m-1), std::min(matA.ku+matB.ku,matB.n-1) );
  newmat.zero();
  
  for(CPPL_INT i=0; i<newmat.m; i++){
    const CPPL_INT jmax =std::min(newmat.n,i+newmat.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),i-newmat.kl); j<jmax; j++){
      const CPPL_INT kmax =std::min( std::min(matA.n,i+matA.ku+1), std::min(matB.m,j+matB.kl+1) );
      for(CPPL_INT k=std::max( std::max(CPPL_INT(0),i-matA.kl), std::max(CPPL_INT(0),j-matB.ku) ); k<kmax; k++){
        newmat(i,j) += matA(i,k)*matB(k,j);
      }
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgematrix operator* ( const _dgbmatrix &  matA, const dgsmatrix &  matB  )

friend

_dgbmatrix*dgsmatrix operator

Definition at line 55 of file _dgbmatrix-dgsmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a product." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( matA.m, matB.n );
  newmat.zero();
  
  const std::vector<dcomponent>::const_iterator matB_data_end =matB.data.end();
  for(std::vector<dcomponent>::const_iterator it=matB.data.begin(); it!=matB_data_end; it++){
    const CPPL_INT imax =std::min(matA.m,it->i+matA.kl);
    for(CPPL_INT i=std::max(CPPL_INT(0),it->i-(matA.ku+1)); i<imax; i++){
      newmat(i,it->j) += matA(i,it->i)*it->v;
    }
  }
  
  matA.destroy();
  return _(newmat);
}

_dgematrix operator* ( const _dgbmatrix &  matA, const _dgsmatrix &  matB  )

friend

_dgbmatrix*_dgsmatrix operator

Definition at line 57 of file _dgbmatrix-_dgsmatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.n!=matB.m){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a product." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( matA.m, matB.n );
  newmat.zero();
  
  const std::vector<dcomponent>::const_iterator matB_data_end =matB.data.end();
  for(std::vector<dcomponent>::const_iterator it=matB.data.begin(); it!=matB_data_end; it++){
    const CPPL_INT imax =std::min(matA.m,it->i+matA.kl);
    for(CPPL_INT i=std::max(CPPL_INT(0),it->i-(matA.ku+1)); i<imax; i++){
      newmat(i,it->j) += matA(i,it->i)*it->v;
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgematrix operator* ( const _dgbmatrix &  , const dssmatrix &    )

friend

_dgematrix operator* ( const _dgbmatrix &  , const _dssmatrix &    )

friend

_dgbmatrix operator* ( const _dgbmatrix &  mat, const double &  d  )

friend

_dgbmatrix*double operator

Definition at line 3 of file _dgbmatrix-double.hpp.

{CPPL_VERBOSE_REPORT;
  CPPL_INT size =(mat.kl+mat.ku+1)*mat.n;
  CPPL_INT inc =1;
  dscal_(&size, &d, mat.array, &inc);
  return mat;
}

_dgbmatrix operator/ ( const _dgbmatrix &  mat, const double &  d  )

friend

_dgbmatrix/double operator

Definition at line 13 of file _dgbmatrix-double.hpp.

{CPPL_VERBOSE_REPORT;
  CPPL_INT size =(mat.kl+mat.ku+1)*mat.n;
  double dinv =1./d;
  CPPL_INT inc =1;
  dscal_(&size, &dinv, mat.array, &inc);
  return mat;
}

_dgbmatrix operator* ( const double &  d, const _dgbmatrix &  mat  )

friend

double*_dgbmatrix operator

Definition at line 3 of file double-_dgbmatrix.hpp.

{CPPL_VERBOSE_REPORT;
  CPPL_INT size =(mat.kl+mat.ku+1)*mat.n;
  CPPL_INT inc =1;
  dscal_(&size, &d, mat.array, &inc);
  return mat;
}

Member Data Documentation

CPPL_INT _dgbmatrix::m

mutable

matrix row size

Definition at line 9 of file _dgbmatrix.hpp.

Referenced by _(), _dgbmatrix(), dgbmatrix::dgbmatrix(), i(), nullify(), operator()(), operator*(), dgbmatrix::operator*=(), dgematrix::operator*=(), operator+(), dgbmatrix::operator+=(), dgematrix::operator+=(), operator-(), dgbmatrix::operator-=(), dgematrix::operator-=(), operator<<(), dgbmatrix::shallow_copy(), t(), to_dgematrix(), to_zgbmatrix(), and write().

CPPL_INT _dgbmatrix::n

mutable

matrix column size

Definition at line 10 of file _dgbmatrix.hpp.

Referenced by _(), _dgbmatrix(), dgbmatrix::dgbmatrix(), i(), nullify(), operator()(), operator*(), dgbmatrix::operator*=(), dgematrix::operator*=(), operator+(), dgbmatrix::operator+=(), dgematrix::operator+=(), operator-(), dgbmatrix::operator-=(), dgematrix::operator-=(), operator/(), operator<<(), dgbmatrix::shallow_copy(), t(), to_dgematrix(), to_zgbmatrix(), and write().

CPPL_INT _dgbmatrix::kl

mutable

lower band width

Definition at line 11 of file _dgbmatrix.hpp.

Referenced by _(), _dgbmatrix(), dgbmatrix::dgbmatrix(), nullify(), operator()(), operator*(), dgbmatrix::operator*=(), dgematrix::operator*=(), operator+(), dgbmatrix::operator+=(), dgematrix::operator+=(), operator-(), dgbmatrix::operator-=(), dgematrix::operator-=(), operator/(), operator<<(), dgbmatrix::shallow_copy(), t(), to_dgematrix(), to_zgbmatrix(), and write().

CPPL_INT _dgbmatrix::ku

mutable

upper band width

Definition at line 12 of file _dgbmatrix.hpp.

Referenced by _(), _dgbmatrix(), dgbmatrix::dgbmatrix(), nullify(), operator()(), operator*(), dgbmatrix::operator*=(), dgematrix::operator*=(), operator+(), dgbmatrix::operator+=(), dgematrix::operator+=(), operator-(), dgbmatrix::operator-=(), dgematrix::operator-=(), operator/(), operator<<(), dgbmatrix::shallow_copy(), t(), to_dgematrix(), to_zgbmatrix(), and write().

double* _dgbmatrix::array

mutable

1D array to store matrix data

Definition at line 13 of file _dgbmatrix.hpp.

Referenced by _(), _dgbmatrix(), destroy(), dgbmatrix::dgbmatrix(), nullify(), operator*(), operator-(), operator/(), dgbmatrix::shallow_copy(), to_zgbmatrix(), and ~_dgbmatrix().

double** _dgbmatrix::darray

mutable

array of pointers of column head addresses

Definition at line 14 of file _dgbmatrix.hpp.

Referenced by _(), _dgbmatrix(), destroy(), dgbmatrix::dgbmatrix(), nullify(), operator()(), dgbmatrix::shallow_copy(), and ~_dgbmatrix().

---

The documentation for this class was generated from the following files:

• _dgbmatrix.hpp
• _dgbmatrix-cast.hpp
• _dgbmatrix-constructor.hpp
• _dgbmatrix-io.hpp
• _dgbmatrix-misc.hpp