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

#include <_dssmatrix.hpp>

Public Member Functions
	_dssmatrix ()

	_dssmatrix (const _dssmatrix &)

	~_dssmatrix ()

_zhsmatrix	to_zhsmatrix () const

_dgematrix	to_dgematrix () const

_dsymatrix	to_dsymatrix () const

_dgsmatrix	to_dgsmatrix () 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 const &	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 > >	line
	vector of vector to store the entry information of component for each row and column More...

Friends
std::ostream &	operator<< (std::ostream &, const _dssmatrix &)

_dssmatrix	t (const dssmatrix &)

void	idamax (CPPL_INT &, CPPL_INT &, const dssmatrix &)

double	damax (const dssmatrix &)

const _dssmatrix &	operator+ (const _dssmatrix &)

_dssmatrix	operator- (const _dssmatrix &)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

_dssmatrix	operator/ (const _dssmatrix &, const double &)

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

Detailed Description

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

Definition at line 3 of file _dssmatrix.hpp.

Constructor & Destructor Documentation

_dssmatrix::_dssmatrix ( )

inline

_dssmatrix constructor without arguments

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

References data, line, and n.

   :m(n)
 {CPPL_VERBOSE_REPORT;
   n =0;
   data.clear();
   line.clear();
 }

_dssmatrix::_dssmatrix ( const _dssmatrix & mat )

inline

_dssmatrix copy constructor

Definition at line 13 of file _dssmatrix-constructor.hpp.

References data, line, n, and nullify().

   :m(n)
 {CPPL_VERBOSE_REPORT;
   n =mat.n;
   data.swap(mat.data);
   line.swap(mat.line);
   
   mat.nullify();
 }

_dssmatrix::~_dssmatrix ( )

inline

_dssmatrix destructor

Definition at line 29 of file _dssmatrix-constructor.hpp.

References data, and line.

 {CPPL_VERBOSE_REPORT;
   data.clear();
   line.clear();
 }

Member Function Documentation

_zhsmatrix _dssmatrix::to_zhsmatrix ( ) const

inline

cast to _zhsmatrix

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

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

 {CPPL_VERBOSE_REPORT;
   zhsmatrix newmat(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.0));
   }
   
   destroy();
   return _(newmat);
 }

_dgematrix _dssmatrix::to_dgematrix ( ) const

inline

convert to _dgematrix

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

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

 {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;
     newmat(it->j, it->i) =it->v;
   }
   
   destroy();
   return _(newmat);
 }

_dsymatrix _dssmatrix::to_dsymatrix ( ) const

inline

convert to _dsymatrix

Definition at line 39 of file _dssmatrix-cast.hpp.

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

 {CPPL_VERBOSE_REPORT;
   dsymatrix newmat(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);
 }

_dgsmatrix _dssmatrix::to_dgsmatrix ( ) const

inline

convert to _dgsmatrix

Definition at line 55 of file _dssmatrix-cast.hpp.

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

 {CPPL_VERBOSE_REPORT;
   dgsmatrix newmat(m,n,CPPL_INT(data.size()*2));
   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.put(it->i, it->j, it->v);
     if(it->i!=it->j){
       newmat.put(it->j, it->i, it->v);
     }
   }
   
   destroy();
   return _(newmat);
 }

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

inline

operator() for const object

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

References data, line, and n.

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

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

inline

Definition at line 66 of file _dssmatrix-io.hpp.

References data, destroy(), 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 << "#dssmatrix " << 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 _dssmatrix::nullify ( ) const

inline

nullify all the matrix data

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

References data, line, and n.

Referenced by _dssmatrix(), dssmatrix::dssmatrix(), and dssmatrix::shallow_copy().

 {CPPL_VERBOSE_REPORT;
   n=0;
   data.clear();
   line.clear();
 }

void _dssmatrix::destroy ( ) const

inline

destroy all the matrix data

Definition at line 12 of file _dssmatrix-misc.hpp.

References data, and line.

Referenced by damax(), idamax(), operator*(), operator+(), dssmatrix::operator+=(), operator-(), dssmatrix::operator-=(), operator<<(), to_dgematrix(), to_dgsmatrix(), to_dsymatrix(), to_zhsmatrix(), and write().

 {CPPL_VERBOSE_REPORT;
   data.clear();
   line.clear();
 }

Friends And Related Function Documentation

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

friend

Definition at line 31 of file _dssmatrix-io.hpp.

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

_dssmatrix t ( const dssmatrix & mat )

friend

return transposed dssmatrix

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

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   WARNING_REPORT;
   std::cerr << "This function call has no effect since the matrix is symmetric." << std::endl;
 #endif//CPPL_DEBUG
 
   dssmatrix newmat(mat);
   return _(newmat);
 }

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

friend

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

Definition at line 21 of file dssmatrix-calc.hpp.

 {CPPL_VERBOSE_REPORT;
   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 dssmatrix & mat )

friend

return its largest absolute value

Definition at line 40 of file dssmatrix-calc.hpp.

 {CPPL_VERBOSE_REPORT;
   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 _dssmatrix& operator+ ( const _dssmatrix & mat )

friend

+_dssmatrix operator

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

 {CPPL_VERBOSE_REPORT;
   return mat;
 }

_dssmatrix operator- ( const _dssmatrix & mat )

friend

-_dssmatrix operator

Definition at line 10 of file _dssmatrix-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 _dssmatrix &	,
		const dgematrix &
	)

friend

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

_dssmatrix operator+	(	const _dssmatrix &	matA,
		const dssmatrix &	matB
	)

friend

_dssmatrix+dssmatrix operator

Definition at line 3 of file _dssmatrix-dssmatrix.hpp.

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

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

friend

_dssmatrix+_dssmatrix operator

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

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

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

_dssmatrix operator-	(	const _dssmatrix &	matA,
		const dssmatrix &	matB
	)

friend

_dssmatrix-dssmatrix operator

Definition at line 26 of file _dssmatrix-dssmatrix.hpp.

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(matA.n!=matB.n){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a subtraction." << std::endl
               << "Your input was (" << matA.n << "x" << matA.n << ") - (" << matB.n << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
 
   dssmatrix 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);
 }

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

friend

_dssmatrix-_dssmatrix operator

Definition at line 27 of file _dssmatrix-_dssmatrix.hpp.

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(matA.n!=matB.n){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a subtraction." << std::endl
               << "Your input was (" << matA.n << "x" << matA.n << ") - (" << matB.n << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   dssmatrix 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);
 }

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

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

friend

_dssmatrix*double operator

Definition at line 3 of file _dssmatrix-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;
 }

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

friend

_dssmatrix/double operator

Definition at line 15 of file _dssmatrix-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;
 }

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

friend

double*_dssmatrix operator

Definition at line 3 of file double-_dssmatrix.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;
 }