_dgsmatrix Class Reference

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

#include <_dgsmatrix.hpp>

Public Member Functions
	_dgsmatrix ()
	_dgsmatrix (const _dgsmatrix &)
	~_dgsmatrix ()
_zgsmatrix	to_zgsmatrix () 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...
std::vector< dcomponent >	data
	matrix data More...
std::vector< std::vector < CPPL_INT > >	rows
	array of vector to store the entry information of component for each row More...
std::vector< std::vector < CPPL_INT > >	cols
	array of vector to store the entry information of component for each column More...

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

Detailed Description

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

Definition at line 3 of file _dgsmatrix.hpp.

Constructor & Destructor Documentation

_dgsmatrix::_dgsmatrix ( )

inline

_dgsmatrix constructor without arguments

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

References cols, data, m, n, and rows.

{CPPL_VERBOSE_REPORT;
  m =0;
  n =0;
  data.clear();
  rows.clear();
  cols.clear();
}

_dgsmatrix::_dgsmatrix ( const _dgsmatrix &  mat )

inline

_dgsmatrix copy constructor

Definition at line 14 of file _dgsmatrix-constructor.hpp.

References cols, data, m, n, nullify(), and rows.

{CPPL_VERBOSE_REPORT;
  m =mat.m;
  n =mat.n;
  data.swap(mat.data);
  rows.swap(mat.rows);
  cols.swap(mat.cols);
  
  mat.nullify();
}

_dgsmatrix::~_dgsmatrix ( )

inline

_dgsmatrix destructor

Definition at line 31 of file _dgsmatrix-constructor.hpp.

References cols, data, and rows.

{CPPL_VERBOSE_REPORT;
  data.clear();
  rows.clear();
  cols.clear();
}

Member Function Documentation

_zgsmatrix _dgsmatrix::to_zgsmatrix ( ) const

inline

cast to _zgsmatrix

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

References _(), data, destroy(), m, n, and zgsmatrix::put().

{CPPL_VERBOSE_REPORT;
  zgsmatrix newmat(m,n,CPPL_INT(data.size()));
  
  const std::vector<dcomponent>::const_iterator data_end =data.end();
  for(std::vector<dcomponent>::const_iterator it=data.begin(); it!=data_end; it++){
    newmat.put(it->i, it->j, comple(it->v,0.));
  }
  
  destroy();
  return _(newmat);
}

_dgematrix _dgsmatrix::to_dgematrix ( ) const

inline

convert to _dgematrix

Definition at line 22 of file _dgsmatrix-cast.hpp.

References _(), data, destroy(), m, n, and dgematrix::zero().

Referenced by operator+(), and operator-().

{CPPL_VERBOSE_REPORT;
  dgematrix newmat(m,n);
  newmat.zero();
  
  const std::vector<dcomponent>::const_iterator data_end =data.end();
  for(std::vector<dcomponent>::const_iterator it=data.begin(); it!=data_end; it++){
    newmat(it->i,it->j) = it->v;
  }
  
  destroy();
  return _(newmat);
}

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

inline

operator() for const object

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

References data, i(), m, n, and rows.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if( i<0 || j<0 || m<=i || n<=j ){
    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 << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  //////// search (i,j) component ////////
  const std::vector<CPPL_INT>::iterator rows_i_end =rows[i].end();
  for(std::vector<CPPL_INT>::iterator p=rows[i].begin(); p!=rows_i_end; p++){
    if(data[*p].j==j){ return data[*p].v; }
  }
  
  //////// (i,j) component was not found ////////
  return 0.0;
}

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

inline

Definition at line 53 of file _dgsmatrix-io.hpp.

References data, destroy(), m, and n.

{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 << "#dgsmatrix " << m << " " << n << " " << data.size() << std::endl;
  
  const std::vector<dcomponent>::const_iterator data_end =data.end();
  for(std::vector<dcomponent>::const_iterator it=data.begin(); it!=data_end; it++){
    ofs << it->i << " " << it->j << " " << it->v << std::endl;
  }
  
  ofs.close();
  destroy();
}

void _dgsmatrix::nullify ( ) const

inline

nullify all the matrix data

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

References cols, data, m, n, and rows.

Referenced by _dgsmatrix(), dgsmatrix::dgsmatrix(), and dgsmatrix::shallow_copy().

{CPPL_VERBOSE_REPORT;
  m=0;
  n=0;
  data.clear();
  rows.clear();
  cols.clear();
}

void _dgsmatrix::destroy ( ) const

inline

destroy all the matrix data

Definition at line 14 of file _dgsmatrix-misc.hpp.

References cols, data, and rows.

Referenced by damax(), idamax(), operator*(), dgsmatrix::operator*=(), operator+(), dgsmatrix::operator+=(), operator-(), dgsmatrix::operator-=(), operator<<(), to_dgematrix(), to_zgsmatrix(), and write().

{CPPL_VERBOSE_REPORT;
  data.clear();
  rows.clear();
  cols.clear();
}

Friends And Related Function Documentation

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

friend

Definition at line 29 of file _dgsmatrix-io.hpp.

{CPPL_VERBOSE_REPORT;
  for(CPPL_INT i=0; i<mat.m; i++){
    for(CPPL_INT j=0; j<mat.n; j++){
      std::vector<CPPL_INT>::iterator q;
      const std::vector<CPPL_INT>::iterator mat_rows_i_end =mat.rows[i].end();
      for(q=mat.rows[i].begin(); q!=mat_rows_i_end; q++){
        if(mat.data[*q].j==j){ break; }
      }
      if(q!=mat_rows_i_end){ s << " " << mat.data[*q].v; }
      else{ s << " x"; }
    }
    s << std::endl;
  }
  
  mat.destroy();
  return s;
}

_dgsmatrix t ( const dgsmatrix &  mat )

friend

return transposed dgsmatrix

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

{CPPL_VERBOSE_REPORT;
  dgsmatrix newmat =mat;
  
  std::swap(newmat.m,newmat.n);
  std::swap(newmat.rows,newmat.cols);
  const std::vector<dcomponent>::iterator newmat_data_end =newmat.data.end();
  for(std::vector<dcomponent>::iterator it=newmat.data.begin(); it!=newmat_data_end; it++){
    std::swap(it->i,it->j);
  }
  
  return _(newmat);
}

void idamax ( CPPL_INT &  i, CPPL_INT &  j, const dgsmatrix &  mat  )

friend

search the index of element having the largest absolute value in 0-based numbering system

Definition at line 23 of file dgsmatrix-calc.hpp.

{CPPL_VERBOSE_REPORT;
  //////// exception ////////
  if(mat.data.size()==0){
    WARNING_REPORT;
    std::cerr << "The dgsmatrix is a zero matrix." << std::endl;
    return;
  }
  
  //////// find ////////
  std::vector<dcomponent>::const_iterator itx(mat.data.begin());
  double vmax =0.;
  
  const std::vector<dcomponent>::const_iterator mat_data_end =mat.data.end();
  for(std::vector<dcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
    if( vmax < fabs(it->v) ){
      vmax =fabs(it->v);
      itx =it;
    }
  }
  
  i =itx->i;
  j =itx->j;
}

double damax ( const dgsmatrix &  mat )

friend

return its largest absolute value

Definition at line 50 of file dgsmatrix-calc.hpp.

{CPPL_VERBOSE_REPORT;
  //////// exception ////////
  if(mat.data.size()==0){
    return 0.;
  }
  
  //////// find ////////
  std::vector<dcomponent>::const_iterator itx(mat.data.begin());
  double vmax =0.;
  
  const std::vector<dcomponent>::const_iterator mat_data_end =mat.data.end();
  for(std::vector<dcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
    if( vmax < fabs(it->v) ){
      vmax =fabs(it->v);
      itx =it;
    }
  }

  return itx->v;
}

const _dgsmatrix& operator+ ( const _dgsmatrix &  mat )

friend

+_dgsmatrix operator

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

{CPPL_VERBOSE_REPORT;
  return mat;
}

_dgsmatrix operator- ( const _dgsmatrix &  mat )

friend

-_dgsmatrix operator

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

{CPPL_VERBOSE_REPORT;
  const std::vector<dcomponent>::iterator mat_data_end =mat.data.end();
  for(std::vector<dcomponent>::iterator it=mat.data.begin(); it!=mat_data_end; it++){
    it->v = -it->v;
  }
  
  return mat;
}

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

friend

_dgsmatrix+dgematrix operator

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

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.m!=matB.m || 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.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat(matB);
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    newmat(it->i,it->j) += it->v;
  }
  
  matA.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix+_dgematrix operator

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

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.m!=matB.m || 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.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    matB(it->i,it->j) += it->v;
  }
  
  matA.destroy();
  return matB;
}

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

friend

_dgsmatrix+dsymatrix operator

Definition at line 3 of file _dgsmatrix-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() );
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    newmat(it->i,it->j) += it->v;
  }
  
  matA.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix+_dsymatrix operator

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

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.m!=matB.m || 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.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat( matB.to_dgematrix() );
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    newmat(it->i,it->j) += it->v;
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix+dgbmatrix operator

Definition at line 3 of file _dgsmatrix-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
  
  dgematrix newmat( matA.to_dgematrix() );
  
  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++){
      newmat(i,j)+=matB(i,j);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix+_dgbmatrix operator

Definition at line 3 of file _dgsmatrix-_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
  
  dgematrix newmat( matA.to_dgematrix() );
  
  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++){
      newmat(i,j)+=matB(i,j);
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgsmatrix operator+ ( const _dgsmatrix &  matA, const dgsmatrix &  matB  )

friend

_dgsmatrix+dgsmatrix operator

Definition at line 3 of file _dgsmatrix-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
  
  dgsmatrix newmat(matA);
  
  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++){
    newmat(it->i,it->j) += it->v;
  }
  
  return _(newmat);
}

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

friend

_dgsmatrix+_dgsmatrix operator

Definition at line 3 of file _dgsmatrix-_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
  
  dgsmatrix newmat(matA);
  
  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++){
    newmat(it->i,it->j) += it->v;
  }
  
  matB.destroy();
  return _(newmat);
}

_dgsmatrix operator+ ( const _dgsmatrix &  , const dssmatrix &    )

friend

_dgsmatrix operator+ ( const _dgsmatrix &  , const _dssmatrix &    )

friend

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

friend

_dgsmatrix-dgematrix operator

Definition at line 27 of file _dgsmatrix-dgematrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.m!=matB.m || matA.n!=matB.n){
    ERROR_REPORT;
    std::cerr << "These two matrises can not make a subtraction." << 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);
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    newmat(it->i,it->j) += it->v;
  }
  
  matA.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix-_dgematrix operator

Definition at line 25 of file _dgsmatrix-_dgematrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.m!=matB.m || matA.n!=matB.n){
    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
  
  //// change sign ////
  for(int i=0; i<matB.m*matB.n; i++){
    matB.array[i] = -matB.array[i];
  }
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    matB(it->i,it->j) += it->v;
  }
  
  matA.destroy();
  return matB;
}

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

friend

_dgsmatrix-dsymatrix operator

Definition at line 27 of file _dgsmatrix-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 subtraction." << 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() );
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    newmat(it->i,it->j) += it->v;
  }
  
  matA.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix-_dsymatrix operator

Definition at line 28 of file _dgsmatrix-_dsymatrix.hpp.

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if(matA.m!=matB.m || matA.n!=matB.n){
    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

  //// shallow copy to dgematrix  ////
  dgematrix newmat( (-matB).to_dgematrix() );
  
  //// add ////
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    newmat(it->i,it->j) += it->v;
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix-dgbmatrix operator

Definition at line 29 of file _dgsmatrix-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
  
  dgematrix newmat( matA.to_dgematrix() );
  
  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++){
      newmat(i,j)-=matB(i,j);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix-_dgbmatrix operator

Definition at line 30 of file _dgsmatrix-_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
  
  dgematrix newmat( matA.to_dgematrix() );
  
  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++){
      newmat(i,j)-=matB(i,j);
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgsmatrix operator- ( const _dgsmatrix &  matA, const dgsmatrix &  matB  )

friend

_dgsmatrix-dgsmatrix operator

Definition at line 26 of file _dgsmatrix-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 subtraction." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") - (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgsmatrix newmat(matA);
  
  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++){
    newmat(it->i,it->j) -= it->v;
  }
  
  return _(newmat);
}

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

friend

_dgsmatrix-_dgsmatrix operator

Definition at line 27 of file _dgsmatrix-_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 subtraction." << std::endl
              << "Your input was (" << matA.m << "x" << matA.n << ") - (" << matB.m << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgsmatrix newmat(matA);
  
  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++){
    newmat(it->i,it->j) -= it->v;
  }
  
  matB.destroy();
  return _(newmat);
}

_dgsmatrix operator- ( const _dgsmatrix &  , const dssmatrix &    )

friend

_dgsmatrix operator- ( const _dgsmatrix &  , const _dssmatrix &    )

friend

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

friend

_dgsmatrix*dcovector operator

Definition at line 3 of file _dgsmatrix-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);
  newvec.zero();
  
  const std::vector<dcomponent>::const_iterator mat_data_end =mat.data.end();
  for(std::vector<dcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
    newvec(it->i) += it->v*vec(it->j);
  }
  
  mat.destroy();
  return _(newvec);
}

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

friend

_dgsmatrix*_dcovector operator

Definition at line 3 of file _dgsmatrix-_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);
  newvec.zero();
  
  const std::vector<dcomponent>::const_iterator mat_data_end =mat.data.end();
  for(std::vector<dcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
    newvec(it->i) += it->v * vec(it->j);
  }
  
  mat.destroy();
  vec.destroy();
  return _(newvec);
}

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

friend

_dgsmatrix*dgematrix operator

Definition at line 51 of file _dgsmatrix-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.n << "x" << matB.n << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  dgematrix newmat(matA.m, matB.n);
  newmat.zero();
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    for(CPPL_INT j=0; j<matB.n; j++){
      newmat(it->i,j) += it->v*matB(it->j,j);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix*_dgematrix operator

Definition at line 52 of file _dgsmatrix-_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();
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    for(CPPL_INT j=0; j<matB.n; j++){
      newmat(it->i,j) += it->v *matB(it->j,j);
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix*dsymatrix operator

Definition at line 51 of file _dgsmatrix-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();
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    for(CPPL_INT i=0; i<matB.n; i++){
      newmat(it->i,i) += it->v*matB(it->j,i);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix*_dsymatrix operator

Definition at line 55 of file _dgsmatrix-_dsymatrix.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 matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    for(CPPL_INT i=0; i<matB.n; i++){
      newmat(it->i,i) += it->v*matB(it->j,i);
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix*dgbmatrix operator

Definition at line 55 of file _dgsmatrix-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
  
  dgematrix newmat( matA.m, matB.n );
  newmat.zero();
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    const CPPL_INT jmax =std::min(matB.n,it->j+matB.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),it->j-matB.kl); j<jmax; j++){
      newmat(it->i,j) += it->v*matB(it->j,j);
    }
  }
  
  matA.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix*_dgbmatrix operator

Definition at line 57 of file _dgsmatrix-_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
  
  dgematrix newmat( matA.m, matB.n );
  newmat.zero();
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    const CPPL_INT jmax =std::min(matB.n,it->j+matB.ku+1);
    for(CPPL_INT j=std::max(CPPL_INT(0),it->j-matB.kl); j<jmax; j++){
      newmat(it->i,j) += it->v*matB(it->j,j);
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgsmatrix operator* ( const _dgsmatrix &  matA, const dgsmatrix &  matB  )

friend

_dgsmatrix*dgsmatrix operator

Definition at line 49 of file _dgsmatrix-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
  
  dgsmatrix newmat( matA.m, matB.n );
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    CPPL_INT k =it->j;
    const std::vector<CPPL_INT>::const_iterator matB_rows_k_end =matB.rows[k].end();
    for(std::vector<CPPL_INT>::const_iterator p=matB.rows[k].begin(); p!=matB_rows_k_end; p++){
      newmat(it->i,matB.data[*p].j) += it->v*matB.data[*p].v;
    }
  }
  
  matA.destroy();
  return _(newmat);
}

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

friend

_dgsmatrix*_dgsmatrix operator

Definition at line 51 of file _dgsmatrix-_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
  
  dgsmatrix newmat(matA.m, matB.n);
  
  const std::vector<dcomponent>::const_iterator matA_data_end =matA.data.end();
  for(std::vector<dcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
    CPPL_INT k =it->j;
    const std::vector<CPPL_INT>::iterator matB_rows_k_end =matB.rows[k].end();
    for(std::vector<CPPL_INT>::iterator p=matB.rows[k].begin(); p!=matB_rows_k_end; p++){
      newmat(it->i,matB.data[*p].j) += it->v*matB.data[*p].v;
    }
  }
  
  matA.destroy();
  matB.destroy();
  return _(newmat);
}

_dgsmatrix operator* ( const _dgsmatrix &  , const dssmatrix &    )

friend

_dgsmatrix operator* ( const _dgsmatrix &  , const _dssmatrix &    )

friend

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

friend

_dgsmatrix*double operator

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

{CPPL_VERBOSE_REPORT;
  const std::vector<dcomponent>::iterator mat_data_end =mat.data.end();
  for(std::vector<dcomponent>::iterator it=mat.data.begin(); it!=mat_data_end; it++){
    it->v *= d;
  }
  
  return mat;
}

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

friend

_dgsmatrix/double operator

Definition at line 15 of file _dgsmatrix-double.hpp.

{CPPL_VERBOSE_REPORT;
  const std::vector<dcomponent>::iterator mat_data_end =mat.data.end();
  for(std::vector<dcomponent>::iterator it=mat.data.begin(); it!=mat_data_end; it++){
    it->v /= d;
  }
  
  return mat;
}

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

friend

double*_dgsmatrix operator

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

{CPPL_VERBOSE_REPORT;
  const std::vector<dcomponent>::iterator mat_data_end =mat.data.end();
  for(std::vector<dcomponent>::iterator it=mat.data.begin(); it!=mat_data_end; it++){
    it->v *= d;
  }
  
  return mat;
}

Member Data Documentation

CPPL_INT _dgsmatrix::m

mutable

matrix row size

Definition at line 9 of file _dgsmatrix.hpp.

Referenced by _(), _dgsmatrix(), dgsmatrix::dgsmatrix(), nullify(), operator()(), operator*(), dgsmatrix::operator*=(), operator+(), dgsmatrix::operator+=(), operator-(), dgsmatrix::operator-=(), operator<<(), dgsmatrix::shallow_copy(), to_dgematrix(), to_zgsmatrix(), and write().

CPPL_INT _dgsmatrix::n

mutable

matrix column size

Definition at line 10 of file _dgsmatrix.hpp.

Referenced by _(), _dgsmatrix(), dgsmatrix::dgsmatrix(), nullify(), operator()(), operator*(), dgsmatrix::operator*=(), operator+(), dgsmatrix::operator+=(), operator-(), dgsmatrix::operator-=(), operator<<(), dgsmatrix::shallow_copy(), to_dgematrix(), to_zgsmatrix(), and write().

std::vector<dcomponent> _dgsmatrix::data

mutable

matrix data

Definition at line 11 of file _dgsmatrix.hpp.

Referenced by _(), _dgsmatrix(), damax(), destroy(), dgsmatrix::dgsmatrix(), idamax(), nullify(), operator()(), operator*(), dgsmatrix::operator*=(), operator+(), dgsmatrix::operator+=(), operator-(), dgsmatrix::operator-=(), operator/(), operator<<(), dgsmatrix::shallow_copy(), t(), to_dgematrix(), to_zgsmatrix(), write(), and ~_dgsmatrix().

std::vector< std::vector<CPPL_INT> > _dgsmatrix::rows

mutable

array of vector to store the entry information of component for each row

Definition at line 12 of file _dgsmatrix.hpp.

Referenced by _(), _dgsmatrix(), destroy(), dgsmatrix::dgsmatrix(), nullify(), operator()(), operator*(), dgsmatrix::operator*=(), operator<<(), dgsmatrix::shallow_copy(), t(), and ~_dgsmatrix().

std::vector< std::vector<CPPL_INT> > _dgsmatrix::cols

mutable

array of vector to store the entry information of component for each column

Definition at line 13 of file _dgsmatrix.hpp.

Referenced by _(), _dgsmatrix(), destroy(), dgsmatrix::dgsmatrix(), nullify(), dgsmatrix::shallow_copy(), t(), and ~_dgsmatrix().

---

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

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