Complex Double-precision General Sparse Matrix Class. More...

#include <zgsmatrix.hpp>

Public Member Functions
	zgsmatrix ()

	zgsmatrix (const zgsmatrix &)

	zgsmatrix (const _zgsmatrix &)

	zgsmatrix (const CPPL_INT &, const CPPL_INT &, const CPPL_INT=0)

	zgsmatrix (const char *)

	~zgsmatrix ()

_zgematrix	to_zgematrix () const

comple	operator() (const CPPL_INT &, const CPPL_INT &) const

comple &	operator() (const CPPL_INT &, const CPPL_INT &)

zgsmatrix &	put (const CPPL_INT &, const CPPL_INT &, const comple &)

zgsmatrix &	del (const CPPL_INT, const CPPL_INT)

zgsmatrix &	del (const CPPL_INT)

void	write (const char *) const

void	read (const char *)

void	clear ()

zgsmatrix &	zero ()

zgsmatrix &	identity ()

void	chsign ()

void	copy (const zgsmatrix &)

void	shallow_copy (const _zgsmatrix &)

zgsmatrix &	resize (const CPPL_INT &, const CPPL_INT &, const CPPL_INT=0, const CPPL_INT=0)

void	stretch (const CPPL_INT &, const CPPL_INT &)

void	expand (const CPPL_INT &)

bool	isListed (const CPPL_INT &, const CPPL_INT &)

CPPL_INT	number (const CPPL_INT &, const CPPL_INT &)

void	checkup ()

void	diet (const double=DBL_MIN)

_zrovector	row (const CPPL_INT &) const

_zcovector	col (const CPPL_INT &) const

zgsmatrix &	operator= (const zgsmatrix &)

zgsmatrix &	operator= (const _zgsmatrix &)

zgsmatrix &	operator+= (const zgsmatrix &)

zgsmatrix &	operator+= (const _zgsmatrix &)

zgsmatrix &	operator-= (const zgsmatrix &)

zgsmatrix &	operator-= (const _zgsmatrix &)

zgsmatrix &	operator*= (const zgsmatrix &)

zgsmatrix &	operator*= (const _zgsmatrix &)

zgsmatrix &	operator*= (const double &)

zgsmatrix &	operator*= (const comple &)

zgsmatrix &	operator/= (const double &)

zgsmatrix &	operator/= (const comple &)

Public Attributes
CPPL_INT	m
	matrix row size More...

CPPL_INT	n
	matrix column size More...

std::vector< zcomponent >	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 zgsmatrix &)

void	swap (zgsmatrix &, zgsmatrix &)

_zgsmatrix	_ (zgsmatrix &)

_zgsmatrix	t (const zgsmatrix &)

_zgsmatrix	conj (const zgsmatrix &)

_zgsmatrix	conjt (const zgsmatrix &)

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

comple	damax (const zgsmatrix &)

const zgsmatrix &	operator+ (const zgsmatrix &)

_zgsmatrix	operator- (const zgsmatrix &)

_zgematrix	operator+ (const zgsmatrix &, const zgematrix &)

_zgematrix	operator+ (const zgsmatrix &, const _zgematrix &)

_zgematrix	operator+ (const zgsmatrix &, const zhematrix &)

_zgematrix	operator+ (const zgsmatrix &, const _zhematrix &)

_zgematrix	operator+ (const zgsmatrix &, const zgbmatrix &)

_zgematrix	operator+ (const zgsmatrix &, const _zgbmatrix &)

_zgsmatrix	operator+ (const zgsmatrix &, const zgsmatrix &)

_zgsmatrix	operator+ (const zgsmatrix &, const _zgsmatrix &)

_zgsmatrix	operator+ (const zgsmatrix &, const zhsmatrix &)

_zgsmatrix	operator+ (const zgsmatrix &, const _zhsmatrix &)

_zgematrix	operator- (const zgsmatrix &, const zgematrix &)

_zgematrix	operator- (const zgsmatrix &, const _zgematrix &)

_zgematrix	operator- (const zgsmatrix &, const zhematrix &)

_zgematrix	operator- (const zgsmatrix &, const _zhematrix &)

_zgematrix	operator- (const zgsmatrix &, const zgbmatrix &)

_zgematrix	operator- (const zgsmatrix &, const _zgbmatrix &)

_zgsmatrix	operator- (const zgsmatrix &, const zgsmatrix &)

_zgsmatrix	operator- (const zgsmatrix &, const _zgsmatrix &)

_zgsmatrix	operator- (const zgsmatrix &, const zhsmatrix &)

_zgsmatrix	operator- (const zgsmatrix &, const _zhsmatrix &)

_zcovector	operator* (const zgsmatrix &, const zcovector &)

_zcovector	operator* (const zgsmatrix &, const _zcovector &)

_zgematrix	operator* (const zgsmatrix &, const zgematrix &)

_zgematrix	operator* (const zgsmatrix &, const _zgematrix &)

_zgematrix	operator* (const zgsmatrix &, const zhematrix &)

_zgematrix	operator* (const zgsmatrix &, const _zhematrix &)

_zgematrix	operator* (const zgsmatrix &, const zgbmatrix &)

_zgematrix	operator* (const zgsmatrix &, const _zgbmatrix &)

_zgsmatrix	operator* (const zgsmatrix &, const zgsmatrix &)

_zgsmatrix	operator* (const zgsmatrix &, const _zgsmatrix &)

_zgsmatrix	operator* (const zgsmatrix &, const zhsmatrix &)

_zgsmatrix	operator* (const zgsmatrix &, const _zhsmatrix &)

_zgsmatrix	operator* (const zgsmatrix &, const double &)

_zgsmatrix	operator* (const zgsmatrix &, const comple &)

_zgsmatrix	operator/ (const zgsmatrix &, const double &)

_zgsmatrix	operator/ (const zgsmatrix &, const comple &)

_zgsmatrix	operator* (const double &, const zgsmatrix &)

_zgsmatrix	operator* (const comple &, const zgsmatrix &)

Detailed Description

Complex Double-precision General Sparse Matrix Class.

Definition at line 3 of file zgsmatrix.hpp.

Constructor & Destructor Documentation

zgsmatrix::zgsmatrix ( )

inline

zgsmatrix constructor without arguments

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

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

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

zgsmatrix::zgsmatrix ( const zgsmatrix & mat )

inline

zgsmatrix copy constructor

Definition at line 19 of file zgsmatrix-constructor.hpp.

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

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

zgsmatrix::zgsmatrix ( const _zgsmatrix & mat )

inline

zgsmatrix constructor to cast _zgsmatrix

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

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

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

zgsmatrix::zgsmatrix	(	const CPPL_INT &	_m,
		const CPPL_INT &	_n,
		const CPPL_INT	_c = `0`
	)

inline

zgsmatrix constructor with size specification

Definition at line 52 of file zgsmatrix-constructor.hpp.

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

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if( _m<0 || _n<0 || _c<0 ){
     ERROR_REPORT;
     std::cerr << "Matrix sizes and the length of arrays must be positive integers. " << std::endl
               << "Your input was (" << _m << "," << _n << "," << _c << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   //////// initialize ////////
   m =_m;
   n =_n;
   data.resize(0);
   data.reserve(_c);
   rows.resize(m);
   cols.resize(n);
 }

zgsmatrix::zgsmatrix ( const char * filename )

inline

zgsmatrix constructor with filename

Definition at line 74 of file zgsmatrix-constructor.hpp.

References cols, data, read(), and rows.

 {CPPL_VERBOSE_REPORT;
   data.clear();
   rows.clear();
   cols.clear();
   
   //// read ////
   read(filename);
 }

zgsmatrix::~zgsmatrix ( )

inline

zgsmatrix destructor

Definition at line 90 of file zgsmatrix-constructor.hpp.

References cols, data, and rows.

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

Member Function Documentation

_zgematrix zgsmatrix::to_zgematrix ( ) const

inline

convert to _zgematrix

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

References _, data, zcomponent::i, zcomponent::j, m, n, zcomponent::v, and zero().

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

 {CPPL_VERBOSE_REPORT;
   zgematrix newmat( zgematrix(m,n).zero() );
   
   const size_t data_size =data.size();
   for(size_t c=0; c<data_size; c++){
     const zcomponent& z =data[c];
     newmat(z.i,z.j) =z.v;
   }
   
   return _(newmat);
 }

comple zgsmatrix::operator()	(	const CPPL_INT &	i,
		const CPPL_INT &	j
	)		const

inline

operator() for const object

Definition at line 3 of file zgsmatrix-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 << "x" << n << "." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   //// search (i,j) component ////
   const std::vector<CPPL_INT>::const_iterator rows_i_end =rows[i].end();
   for(std::vector<CPPL_INT>::const_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 comple(0.0,0.0);
 }

comple & zgsmatrix::operator()	(	const CPPL_INT &	i,
		const CPPL_INT &	j
	)

inline

operator() for const object

Definition at line 26 of file zgsmatrix-io.hpp.

References cols, 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 << "x" << 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 not found ////////
   rows[i].push_back(CPPL_INT(data.size()));
   cols[j].push_back(CPPL_INT(data.size()));
   data.push_back( zcomponent(i,j, comple(0.,0.)) );
   return data.back().v;
 }

zgsmatrix & zgsmatrix::put	(	const CPPL_INT &	i,
		const CPPL_INT &	j,
		const comple &	v
	)

inline

put value without isListed check

Definition at line 56 of file zgsmatrix-io.hpp.

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

Referenced by identity(), read(), _dgsmatrix::to_zgsmatrix(), _zhsmatrix::to_zgsmatrix(), dgsmatrix::to_zgsmatrix(), and zhsmatrix::to_zgsmatrix().

 {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 << "x" << n << "." << std::endl;
     exit(1);
   }
   
   if( isListed(i,j) ){
     ERROR_REPORT;
     std::cerr << "The required component is already listed." << std::endl
               << "Your input was (" << i << "," << j << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   //// push ////
   rows[i].push_back(CPPL_INT(data.size()));
   cols[j].push_back(CPPL_INT(data.size()));
   data.push_back(zcomponent(i,j,v));
   
   return *this;
 }

zgsmatrix & zgsmatrix::del	(	const CPPL_INT	i,
		const CPPL_INT	j
	)

inline

delete the entry of a component

Definition at line 88 of file zgsmatrix-io.hpp.

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

Referenced by diet(), and stretch().

 {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 << "x" << 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){//exists
       //// save position ////
       CPPL_INT c =*p;
       CPPL_INT C =CPPL_INT(data.size()-1);
       
       //// data translation ////
       data[c]=data.back();
       data.pop_back();
       
       //// remove from List ////
       rows[i].erase(p);
       const std::vector<CPPL_INT>::iterator cols_j_end =cols[j].end();
       for(std::vector<CPPL_INT>::iterator q=cols[j].begin(); q!=cols_j_end; q++){
         if(*q==c){ cols[j].erase(q); break; }
       }
       
       //// modify the entry of translated component ////
       CPPL_INT I(data[c].i), J(data[c].j);
       const std::vector<CPPL_INT>::iterator rows_I_end =rows[I].end();
       for(std::vector<CPPL_INT>::iterator q=rows[I].begin(); q!=rows_I_end; q++){
         if(*q==C){ *q=c; break; }
       }
       const std::vector<CPPL_INT>::iterator cols_J_end =cols[J].end();
       for(std::vector<CPPL_INT>::iterator q=cols[J].begin(); q!=cols_J_end; q++){
         if(*q==C){ *q=c; break; }
       }
       return *this;
     }
   }
   
 #ifdef  CPPL_DEBUG
   std::cerr << "# [NOTE] zgsmatrix::del(CPPL_INT&, CPPL_INT&): The required component was not listed. Your input was (" << i << "," << j << ")." << std::endl;
 #endif//CPPL_DEBUG
   
   return *this;
 }

zgsmatrix & zgsmatrix::del ( const CPPL_INT c )

inline

delete the entry of an element

Definition at line 141 of file zgsmatrix-io.hpp.

References cols, data, i(), and rows.

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if( c<0 || c>=CPPL_INT(data.size()) ){
     ERROR_REPORT;
     std::cerr << "The required element is out of the matrix volume." << std::endl
               << "Your input was (" << c << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   if( c==CPPL_INT(data.size()-1) ){//if c is the last element
     CPPL_INT i(data[c].i), j(data[c].j);
     const std::vector<CPPL_INT>::iterator rows_i_end =rows[i].end();
     for(std::vector<CPPL_INT>::iterator q=rows[i].begin(); q!=rows_i_end; q++){
       if(*q==c){ rows[i].erase(q); break; }
     }
     const std::vector<CPPL_INT>::iterator cols_j_end =cols[j].end();
     for(std::vector<CPPL_INT>::iterator q=cols[j].begin(); q!=cols_j_end; q++){
       if(*q==c){ cols[j].erase(q); break; }
     }
     data.pop_back();
   }
 
   else{//if c is NOT the last element
     //// data translation ////
     CPPL_INT C =CPPL_INT(data.size()-1);
     CPPL_INT i(data[c].i), j(data[c].j), I(data.back().i), J(data.back().j);
     data[c]=data.back();
     //// remove entry of component ////
     const std::vector<CPPL_INT>::iterator rows_i_end =rows[i].end();
     for(std::vector<CPPL_INT>::iterator q=rows[i].begin(); q!=rows_i_end; q++){
       if(*q==c){ rows[i].erase(q); break; }
     }
     const std::vector<CPPL_INT>::iterator cols_j_end =cols[j].end();
     for(std::vector<CPPL_INT>::iterator q=cols[j].begin(); q!=cols_j_end; q++){
       if(*q==c){ cols[j].erase(q); break; }
     }
     //// modify the entry of translated component ////
     const std::vector<CPPL_INT>::iterator rows_I_end =rows[I].end();
     for(std::vector<CPPL_INT>::iterator q=rows[I].begin(); q!=rows_I_end; q++){
       if(*q==C){ *q=c; break; }
     }
     const std::vector<CPPL_INT>::iterator cols_J_end =cols[J].end();
     for(std::vector<CPPL_INT>::iterator q=cols[J].begin(); q!=cols_J_end; q++){
       if(*q==C){ *q=c; break; }
     }
     //// pop_back ////
     data.pop_back();
   }
   
   return *this;
 }

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

inline

Definition at line 223 of file zgsmatrix-io.hpp.

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

void zgsmatrix::read ( const char * filename )

inline

Definition at line 242 of file zgsmatrix-io.hpp.

References data, i(), m, n, put(), and resize().

Referenced by zgsmatrix().

 {CPPL_VERBOSE_REPORT;
   std::ifstream s( filename );
   if(!s){
     ERROR_REPORT;
     std::cerr << "The file \"" << filename << "\" can not be opened." << std::endl;
     exit(1);
   }
 
   std::string id;
   s >> id;
   if( id != "zgsmatrix" && id != "#zgsmatrix" ){
     ERROR_REPORT;
     std::cerr << "The type name of the file \"" << filename << "\" is not zgsmatrix." << std::endl
               << "Its type name was " << id << " ." << std::endl;
     exit(1);
   }
   
   //////// m, n, size ////////
   size_t size;
   s >> m >> n >> size;
   resize(m, n);
   data.reserve(size);
   
   //////// i, j, v ////////
   CPPL_INT i, j;
   comple v;
   for(size_t k=0; k<size; k++){
     s >> i >> j >> v;
     put(i,j, v);
   }
   
   //// check ////
   s >> i;
   if(!s.eof()){
     ERROR_REPORT;
     std::cerr << "There is something is wrong with the file \"" << filename << " ." << std::endl
               << "Most likely, there are too many data components over the context." << std::endl;
     exit(1);
   }
   s.close();
 }

void zgsmatrix::clear ( )

inline

clear all the matrix data and set the sizes 0

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

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

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

zgsmatrix & zgsmatrix::zero ( )

inline

change the matrix into a zero matrix

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

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

Referenced by col(), identity(), row(), to_zgematrix(), _zhsmatrix::to_zgsmatrix(), and zhsmatrix::to_zgsmatrix().

 {CPPL_VERBOSE_REPORT;
   data.resize(0);
   for(CPPL_INT i=0; i<m; i++){ rows[i].resize(0); }
   for(CPPL_INT j=0; j<n; j++){ cols[j].resize(0); }
   return *this;
 }

zgsmatrix & zgsmatrix::identity ( )

inline

change the matrix into an identity matrix

Definition at line 24 of file zgsmatrix-misc.hpp.

References i(), m, n, put(), and zero().

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(m!=n){
     ERROR_REPORT;
     std::cerr << "Only square matrix can be a identity matrix." << std::endl
               << "The matrix size was " << m << "x" << n << "." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zero();
   for(CPPL_INT i=0; i<m; i++){
     put(i,i, comple(1.,0.));
   }
   return *this;
 }

void zgsmatrix::chsign ( )

inline

change sign(+/-) of the matrix

Definition at line 44 of file zgsmatrix-misc.hpp.

References data.

Referenced by operator-().

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

void zgsmatrix::copy ( const zgsmatrix & mat )

inline

make a deep copy of the matrix

Definition at line 54 of file zgsmatrix-misc.hpp.

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

Referenced by operator=(), and zgsmatrix().

 {CPPL_VERBOSE_REPORT;
   m =mat.m;
   n =mat.n;
   data =mat.data;
   rows =mat.rows;
   cols =mat.cols;
 }

void zgsmatrix::shallow_copy ( const _zgsmatrix & mat )

inline

make a shallow copy of the matrix
This function is not designed to be used in project codes.

Definition at line 66 of file zgsmatrix-misc.hpp.

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

Referenced by operator=().

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

zgsmatrix & zgsmatrix::resize	(	const CPPL_INT &	_m,
		const CPPL_INT &	_n,
		const CPPL_INT	_c = `0`,
		const CPPL_INT	_l = `0`
	)

inline

resize the matrix

Definition at line 83 of file zgsmatrix-misc.hpp.

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

Referenced by read().

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if( _m<0 || _n<0 || _c<0 ){
     ERROR_REPORT;
     std::cerr << "Matrix sizes and the length of arrays must be positive integers. " << std::endl
               << "Your input was (" << _m << "," << _n << "," << _c << "," << _l << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   m =_m;
   n =_n;
   data.resize(0);
   data.reserve(_c);
   rows.resize(m);
   for(CPPL_INT i=0; i<m; i++){
     rows[i].resize(0);
     rows[i].reserve(_l);
   }
   cols.resize(n);
   for(CPPL_INT i=0; i<n; i++){
     cols[i].resize(0);
     cols[i].reserve(_l);
   }
   
   return *this;
 }

void zgsmatrix::stretch	(	const CPPL_INT &	dm,
		const CPPL_INT &	dn
	)

inline

stretch the matrix size

Definition at line 114 of file zgsmatrix-misc.hpp.

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

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if( m+dm<0 || n+dn<0 ){
     ERROR_REPORT;
     std::cerr << "The new matrix size must be larger than zero. " << std::endl
               << "Your input was (" << dm << ", " << dn << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   //////// zero ////////
   if(dm==0 && dn==0){ return; }
   
   //////// non-zero ////////
   m +=dm;
   n +=dn;
   
   //// for rows ////
   if(dm<0){
     //// delete components over the new size ////
     const std::vector<zcomponent>::reverse_iterator data_rend =data.rend();
     for(std::vector<zcomponent>::reverse_iterator it=data.rbegin(); it!=data_rend; it++){
       if( it->i>=m ){ del( CPPL_INT(data_rend-it-1) ); }
     }
     //// shrink rows ////
     for(CPPL_INT i=0; i<-dm; i++){
       rows.pop_back();
     }
   }
   else{//dm>0
     //// expand rows ////
     for(CPPL_INT i=0; i<dm; i++){
       rows.push_back( std::vector<CPPL_INT>(0) );
     }
   }
 
   //// for cols ////
   if(dn<0){
     //// delete components over the new size ////
     const std::vector<zcomponent>::reverse_iterator data_rend =data.rend();
     for(std::vector<zcomponent>::reverse_iterator it=data.rbegin(); it!=data_rend; it++){
       if( it->j>=n ){ del( CPPL_INT(data_rend-it-1) ); }
     }
     for(CPPL_INT j=0; j<-dn; j++){
       cols.pop_back();
     }
   }
   else{//dn>0
     //// expand cols ////
     for(CPPL_INT j=0; j<dn; j++){
       cols.push_back( std::vector<CPPL_INT>(0) );
     }
   }
 }

void zgsmatrix::expand ( const CPPL_INT & )

inline

bool zgsmatrix::isListed	(	const CPPL_INT &	i,
		const CPPL_INT &	j
	)

inline

check if the component is listed

Definition at line 172 of file zgsmatrix-misc.hpp.

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

Referenced by put().

 {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 << "x" << n << "." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   const std::vector<CPPL_INT>::const_iterator rows_i_end =rows[i].end();
   for(std::vector<CPPL_INT>::const_iterator p=rows[i].begin(); p!=rows_i_end; p++){
     if( data[*p].j==j ){ return 1; }
   }
   
   return 0;
 }

CPPL_INT zgsmatrix::number	(	const CPPL_INT &	i,
		const CPPL_INT &	j
	)

inline

return the element number of the component

Definition at line 193 of file zgsmatrix-misc.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 << "x" << n << "." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   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 *p; }
   }
   
   return -1;
 }

void zgsmatrix::checkup ( )

inline

health checkup

Definition at line 224 of file zgsmatrix-misc.hpp.

References data, m, and n.

 {CPPL_VERBOSE_REPORT;
   //////////////// data ////////////////
   //////// check i,j ////////
   const std::vector<zcomponent>::const_iterator data_end =data.end();
   for(std::vector<zcomponent>::const_iterator it=data.begin(); it!=data_end; it++){
     if( it->i>=m || it->j>=n ){
       ERROR_REPORT;
       std::cerr << "A component, (" << it->i << ", " << it->j << "), is out of matrix size." << std::endl;
       exit(1);
     }
   }
   //////// check double listing ////////
   
   //////////////// rows ////////////////
   //////////////// cols ////////////////
   
   //////////////// NOTE ////////////////
   std::cerr << "# [NOTE] zgsmatrix::checkup(): This sparse matrix is fine." << std::endl;
 }

void zgsmatrix::diet ( const double eps = DBL_MIN )

inline

erase components less than DBL_MIN

Definition at line 214 of file zgsmatrix-misc.hpp.

References data, and del().

 {CPPL_VERBOSE_REPORT;
   const std::vector<zcomponent>::reverse_iterator data_rend =data.rend();
   for(std::vector<zcomponent>::reverse_iterator it=data.rbegin(); it!=data_rend; it++){
     if( fabs(it->v.real())<eps && fabs(it->v.imag())<eps ){ del( CPPL_INT(data_rend-it-1) ); }
   }
 }

_zrovector zgsmatrix::row ( const CPPL_INT & _m ) const

inline

get row of the matrix

Definition at line 251 of file zgsmatrix-misc.hpp.

References _, data, m, n, rows, and zero().

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if( _m<0 || _m>m ){
     ERROR_REPORT;
     std::cerr << "Input row number must be between 0 and " << m << "." << std::endl
               << "Your input was " << _m << "." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zrovector vec( zrovector(n).zero() );
   
   const std::vector<CPPL_INT>::const_iterator rows__m_end =rows[_m].end();
   for(std::vector<CPPL_INT>::const_iterator p=rows[_m].begin(); p!=rows__m_end; p++){
     vec(data[*p].j) =data[*p].v;
   }
   
   return _(vec);
 }

_zcovector zgsmatrix::col ( const CPPL_INT & _n ) const

inline

get column of the matrix

Definition at line 274 of file zgsmatrix-misc.hpp.

References _, cols, data, i(), m, n, and zero().

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if( _n<0 || _n>n ){
     ERROR_REPORT;
     std::cerr << "Input row number must be between 0 and " << n << "." << std::endl
               << "Your input was " << _n << "." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zcovector vec( zcovector(m).zero() );
   
   const std::vector<CPPL_INT>::const_iterator cols__n_end =cols[_n].end();
   for(std::vector<CPPL_INT>::const_iterator p=cols[_n].begin(); p!=cols__n_end; p++){
     vec(data[*p].i) =data[*p].v;
   }
   
   return _(vec);
 }

zgsmatrix & zgsmatrix::operator= ( const zgsmatrix & mat )

inline

zgsmatrix=zgsmatrix operator

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

References copy(), and data.

 {CPPL_VERBOSE_REPORT;
   if(&data!=&mat.data){ // if it is NOT self substitution
     copy(mat);
   }
   return *this;
 }

zgsmatrix & zgsmatrix::operator= ( const _zgsmatrix & mat )

inline

zgsmatrix=_zgsmatrix operator

Definition at line 3 of file zgsmatrix-_zgsmatrix.hpp.

References shallow_copy().

 {CPPL_VERBOSE_REPORT;
   shallow_copy(mat);
   return *this;
 }

zgsmatrix & zgsmatrix::operator+= ( const zgsmatrix & mat )

inline

zgsmatrix+=zgsmatrix operator

Definition at line 17 of file zgsmatrix-zgsmatrix.hpp.

References data, m, and n.

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(n!=mat.n || m!=mat.m){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a summation." << std::endl
               << "Your input was (" << m << "x" << n << ") += (" << mat.m << "x" << mat.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   const std::vector<zcomponent>::const_iterator mat_data_end =mat.data.end();
   for(std::vector<zcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
     (*this)(it->i,it->j) +=it->v;
   }
   
   return *this;
 }

zgsmatrix & zgsmatrix::operator+= ( const _zgsmatrix & mat )

inline

zgsmatrix+=_zgsmatrix operator

Definition at line 15 of file zgsmatrix-_zgsmatrix.hpp.

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

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(n!=mat.n || m!=mat.m){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a summation." << std::endl
               << "Your input was (" << m << "x" << n << ") += (" << mat.m << "x" << mat.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   const std::vector<zcomponent>::const_iterator mat_data_end =mat.data.end();
   for(std::vector<zcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
     (*this)(it->i,it->j) += it->v;
   }
   
   mat.destroy();
   return *this;
 }

zgsmatrix & zgsmatrix::operator-= ( const zgsmatrix & mat )

inline

zgsmatrix-=zgsmatrix operator

Definition at line 38 of file zgsmatrix-zgsmatrix.hpp.

References data, m, and n.

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(n!=mat.n || m!=mat.m){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a sutraction." << std::endl
               << "Your input was (" << m << "x" << n << ") -= (" << mat.m << "x" << mat.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   const std::vector<zcomponent>::const_iterator mat_data_end =mat.data.end();
   for(std::vector<zcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
     (*this)(it->i,it->j) -=it->v;
   }
   
   return *this;
 }

zgsmatrix & zgsmatrix::operator-= ( const _zgsmatrix & mat )

inline

zgsmatrix-=_zgsmatrix operator

Definition at line 37 of file zgsmatrix-_zgsmatrix.hpp.

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

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(n!=mat.n || m!=mat.m){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a sutraction." << std::endl
               << "Your input was (" << m << "x" << n << ") -= (" << mat.m << "x" << mat.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   const std::vector<zcomponent>::const_iterator mat_data_end =mat.data.end();
   for(std::vector<zcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
     (*this)(it->i,it->j) -= it->v;
   }
   
   mat.destroy();
   return *this;
 }

zgsmatrix & zgsmatrix::operator*= ( const zgsmatrix & mat )

inline

zgsmatrix*=zgsmatrix operator

Definition at line 59 of file zgsmatrix-zgsmatrix.hpp.

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

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(n!=mat.m){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a product." << std::endl
               << "Your input was (" << m << "x" << n << ") *= (" << mat.m << "x" << mat.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zgsmatrix newmat(m, mat.n);
   
   const std::vector<zcomponent>::const_iterator data_end =data.end();
   for(std::vector<zcomponent>::const_iterator it=data.begin(); it!=data_end; it++){
     CPPL_INT k =it->j;
     const std::vector<CPPL_INT>::const_iterator mat_rows_k_end =mat.rows[k].end();
     for(std::vector<CPPL_INT>::const_iterator p=mat.rows[k].begin(); p!=mat_rows_k_end; p++){
       newmat(it->i,mat.data[*p].j) +=it->v*mat.data[*p].v;
     }
   }
   
   swap(*this,newmat);
   return *this;
 }

zgsmatrix & zgsmatrix::operator*= ( const _zgsmatrix & mat )

inline

zgsmatrix*=_zgsmatrix operator

Definition at line 59 of file zgsmatrix-_zgsmatrix.hpp.

References data, _zgsmatrix::data, _zgsmatrix::destroy(), _zgsmatrix::m, m, n, _zgsmatrix::n, _zgsmatrix::rows, and swap.

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(n!=mat.m){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a product." << std::endl
               << "Your input was (" << m << "x" << n << ") *= (" << mat.m << "x" << mat.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zgsmatrix newmat( m, mat.n, 0 );
   
   const std::vector<zcomponent>::const_iterator data_end =data.end();
   for(std::vector<zcomponent>::const_iterator it=data.begin(); it!=data_end; it++){
     CPPL_INT k =it->j;
     const std::vector<CPPL_INT>::iterator mat_rows_k_end =mat.rows[k].end();
     for(std::vector<CPPL_INT>::iterator p=mat.rows[k].begin(); p!=mat_rows_k_end; p++){
       newmat(it->i,mat.data[*p].j) += it->v*mat.data[*p].v;
     }
   }
   
   swap(*this,newmat);
   mat.destroy();
   return *this;
 }

zgsmatrix & zgsmatrix::operator*= ( const double & d )

inline

zgsmatrix*=double operator

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

References data.

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

zgsmatrix & zgsmatrix::operator*= ( const comple & d )

inline

zgsmatrix*=comple operator

Definition at line 3 of file zgsmatrix-complex.hpp.

References data.

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

zgsmatrix & zgsmatrix::operator/= ( const double & d )

inline

zgsmatrix/=double operator

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

References data.

 {CPPL_VERBOSE_REPORT;
   const std::vector<zcomponent>::iterator data_end =data.end();
   for(std::vector<zcomponent>::iterator it=data.begin(); it!=data_end; it++){
     it->v /=d;
   }
   return *this;
 }

zgsmatrix & zgsmatrix::operator/= ( const comple & d )

inline

zgsmatrix/=comple operator

Definition at line 15 of file zgsmatrix-complex.hpp.

References data.

 {CPPL_VERBOSE_REPORT;
   const std::vector<zcomponent>::iterator data_end =data.end();
   for(std::vector<zcomponent>::iterator it=data.begin(); it!=data_end; it++){
     it->v /=d;
   }
   
   return *this;
 }

Friends And Related Function Documentation

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

friend

Definition at line 200 of file zgsmatrix-io.hpp.

 {CPPL_VERBOSE_REPORT;
   for(CPPL_INT i=0; i<mat.m; i++){
     const std::vector<CPPL_INT>::const_iterator mat_rows_i_end =mat.rows[i].end();
     for(CPPL_INT j=0; j<mat.n; j++){
       std::vector<CPPL_INT>::const_iterator q;
       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;
   }
   
   return s;
 }

void swap	(	zgsmatrix &	A,
		zgsmatrix &	B
	)

friend

swap two matrices

Definition at line 301 of file zgsmatrix-misc.hpp.

Referenced by operator*=().

 {CPPL_VERBOSE_REPORT;
   std::swap(A.n,B.n);
   std::swap(A.m,B.m);
   std::swap(A.data,B.data);
   std::swap(A.rows,B.rows);
   std::swap(A.cols,B.cols);
 }

_zgsmatrix _ ( zgsmatrix & mat )

friend

convert user object to smart-temporary object

Definition at line 312 of file zgsmatrix-misc.hpp.

Referenced by col(), row(), and to_zgematrix().

 {CPPL_VERBOSE_REPORT;
   _zgsmatrix newmat;
   
   //////// shallow copy ////////
   newmat.n =mat.n;
   newmat.m =mat.m;
   std::swap(newmat.data,mat.data);
   std::swap(newmat.rows,mat.rows);
   std::swap(newmat.cols,mat.cols);
 
   //////// nullify ////////
   mat.m =0;
   mat.n =0;
   
   return newmat;
 }

_zgsmatrix t ( const zgsmatrix & mat )

friend

return transposed zgsmatrix

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

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

_zgsmatrix conj ( const zgsmatrix & mat )

friend

return its conjugate matrix

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

 {CPPL_VERBOSE_REPORT;
   zgsmatrix newmat(mat);
   
   const std::vector<zcomponent>::iterator newmat_data_end =newmat.data.end();
   for(std::vector<zcomponent>::iterator it=newmat.data.begin(); it!=newmat_data_end; it++){
     it->v =std::conj(it->v);
   }
   
   return _(newmat);
 }

_zgsmatrix conjt ( const zgsmatrix & mat )

friend

return its conjugate transposed matrix

Definition at line 37 of file zgsmatrix-calc.hpp.

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

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

friend

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

Definition at line 57 of file zgsmatrix-calc.hpp.

 {CPPL_VERBOSE_REPORT;
   std::vector<zcomponent>::const_iterator itx(mat.data.begin());
   double vmax =0.;
   
   const std::vector<zcomponent>::const_iterator mat_data_end =mat.data.end();
   for(std::vector<zcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
     if( vmax < norm(it->v) ){
       vmax =norm(it->v);
       itx =it;
     }
   }
   
   i=itx->i;
   j=itx->j;
 }

comple damax ( const zgsmatrix & mat )

friend

return its largest absolute value

Definition at line 76 of file zgsmatrix-calc.hpp.

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

const zgsmatrix& operator+ ( const zgsmatrix & mat )

friend

+zgsmatrix operator

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

 {CPPL_VERBOSE_REPORT;
   return mat;
 }

_zgsmatrix operator- ( const zgsmatrix & mat )

friend

-zgsmatrix operator

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

 {CPPL_VERBOSE_REPORT;
   zgsmatrix newmat(mat);
   
   const std::vector<zcomponent>::iterator newmat_data_end =newmat.data.end();
   for(std::vector<zcomponent>::iterator it=newmat.data.begin(); it!=newmat_data_end; it++){
     it->v =-it->v;
   }
   
   return _(newmat);
 }

_zgematrix operator+	(	const zgsmatrix &	matA,
		const zgematrix &	matB
	)

friend

zgsmatrix+zgematrix operator

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

_zgematrix operator+	(	const zgsmatrix &	matA,
		const _zgematrix &	matB
	)

friend

zgsmatrix+_zgematrix operator

Definition at line 3 of file zgsmatrix-_zgematrix.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<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
     matB(it->i,it->j) += it->v;
   }
   
   return matB;
 }

_zgematrix operator+	(	const zgsmatrix &	matA,
		const zhematrix &	matB
	)

friend

zgsmatrix+zhematrix operator

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

_zgematrix operator+	(	const zgsmatrix &	matA,
		const _zhematrix &	matB
	)

friend

zgsmatrix+_zhematrix operator

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

_zgematrix operator+	(	const zgsmatrix &	matA,
		const zgbmatrix &	matB
	)

friend

zgsmatrix+zgbmatrix operator

Definition at line 3 of file zgsmatrix-zgbmatrix.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
   
   zgematrix newmat( matA.to_zgematrix() );
   
   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);
     }
   }
   
   return _(newmat);
 }

_zgematrix operator+	(	const zgsmatrix &	matA,
		const _zgbmatrix &	matB
	)

friend

zgsmatrix+_zgbmatrix operator

Definition at line 3 of file zgsmatrix-_zgbmatrix.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
   
   zgematrix newmat( matA.to_zgematrix() );
   
   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);
     }
   }
   
   matB.destroy();
   return _(newmat);
 }

_zgsmatrix operator+	(	const zgsmatrix &	matA,
		const zgsmatrix &	matB
	)

friend

zgsmatrix+zgsmatrix operator

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

_zgsmatrix operator+	(	const zgsmatrix &	matA,
		const _zgsmatrix &	matB
	)

friend

zgsmatrix+_zgsmatrix operator

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

_zgsmatrix operator+	(	const zgsmatrix &	,
		const zhsmatrix &
	)

friend

_zgsmatrix operator+	(	const zgsmatrix &	,
		const _zhsmatrix &
	)

friend

_zgematrix operator-	(	const zgsmatrix &	matA,
		const zgematrix &	matB
	)

friend

zgsmatrix-zgematrix operator

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

_zgematrix operator-	(	const zgsmatrix &	matA,
		const _zgematrix &	matB
	)

friend

zgsmatrix-_zgematrix operator

Definition at line 24 of file zgsmatrix-_zgematrix.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
   
   //// change sign ////
   for(CPPL_INT i=0; i<matB.m*matB.n; i++){
     matB.array[i]=-matB.array[i];
   }
   
   //// add ////
   const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
     matB(it->i,it->j) += it->v;
   }
   
   return matB;
 }

_zgematrix operator-	(	const zgsmatrix &	matA,
		const zhematrix &	matB
	)

friend

zgsmatrix-zhematrix operator

Definition at line 26 of file zgsmatrix-zhematrix.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
 
   //// shallow copy to zgematrix ////
   zgematrix newmat( (-matB).to_zgematrix() );
   
   //// add ////
   const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
     newmat(it->i,it->j) += it->v;
   }
   
   return _(newmat);
 }

_zgematrix operator-	(	const zgsmatrix &	matA,
		const _zhematrix &	matB
	)

friend

zgsmatrix-_zhematrix operator

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

_zgematrix operator-	(	const zgsmatrix &	matA,
		const zgbmatrix &	matB
	)

friend

zgsmatrix-zgbmatrix operator

Definition at line 28 of file zgsmatrix-zgbmatrix.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
   
   zgematrix newmat( matA.to_zgematrix() );
   
   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);
     }
   }
   
   return _(newmat);
 }

_zgematrix operator-	(	const zgsmatrix &	matA,
		const _zgbmatrix &	matB
	)

friend

zgsmatrix-_zgbmatrix operator

Definition at line 29 of file zgsmatrix-_zgbmatrix.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
   
   zgematrix newmat( matA.to_zgematrix() );
   
   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);
     }
   }
   
   matB.destroy();
   return _(newmat);
 }

_zgsmatrix operator-	(	const zgsmatrix &	matA,
		const zgsmatrix &	matB
	)

friend

zgsmatrix-zgsmatrix operator

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

_zgsmatrix operator-	(	const zgsmatrix &	matA,
		const _zgsmatrix &	matB
	)

friend

zgsmatrix-_zgsmatrix operator

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

_zgsmatrix operator-	(	const zgsmatrix &	,
		const zhsmatrix &
	)

friend

_zgsmatrix operator-	(	const zgsmatrix &	,
		const _zhsmatrix &
	)

friend

_zcovector operator*	(	const zgsmatrix &	mat,
		const zcovector &	vec
	)

friend

zgsmatrix*zcovector operator

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

_zcovector operator*	(	const zgsmatrix &	mat,
		const _zcovector &	vec
	)

friend

zgsmatrix*_zcovector operator

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

_zgematrix operator*	(	const zgsmatrix &	matA,
		const zgematrix &	matB
	)

friend

zgsmatrix*zgematrix operator

Definition at line 49 of file zgsmatrix-zgematrix.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
   
   zgematrix newmat(matA.m, matB.n);
   newmat.zero();
   
   const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::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);
     }
   }
   
   return _(newmat);
 }

_zgematrix operator*	(	const zgsmatrix &	matA,
		const _zgematrix &	matB
	)

friend

zgsmatrix*_zgematrix operator

Definition at line 51 of file zgsmatrix-_zgematrix.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
   
   zgematrix newmat(matA.m, matB.n);
   newmat.zero();
   
   const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::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);
     }
   }
   
   matB.destroy();
   return _(newmat);
 }

_zgematrix operator*	(	const zgsmatrix &	matA,
		const zhematrix &	matB
	)

friend

zgsmatrix*zhematrix operator

Definition at line 51 of file zgsmatrix-zhematrix.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
   
   zgematrix newmat(matA.m, matB.n);
   newmat.zero();
   
   const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::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);
     }
   }
   
   return _(newmat);
 }

_zgematrix operator*	(	const zgsmatrix &	matA,
		const _zhematrix &	matB
	)

friend

zgsmatrix*_zhematrix operator

Definition at line 51 of file zgsmatrix-_zhematrix.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
   
   zgematrix newmat(matA.m, matB.n);
   newmat.zero();
   
   const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::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);
     }
   }
   
   matB.destroy();
   return _(newmat);
 }

_zgematrix operator*	(	const zgsmatrix &	matA,
		const zgbmatrix &	matB
	)

friend

zgsmatrix*zgbmatrix operator

Definition at line 53 of file zgsmatrix-zgbmatrix.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
   
   zgematrix newmat( matA.m, matB.n );
   newmat.zero();
   
   const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::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);
     }
   }
   
   return _(newmat);
 }

_zgematrix operator*	(	const zgsmatrix &	matA,
		const _zgbmatrix &	matB
	)

friend

zgsmatrix*_zgbmatrix operator

Definition at line 55 of file zgsmatrix-_zgbmatrix.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
   
   zgematrix newmat( matA.m, matB.n );
   newmat.zero();
   
   const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::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);
     }
   }
   
   matB.destroy();
   return _(newmat);
 }

_zgsmatrix operator*	(	const zgsmatrix &	matA,
		const zgsmatrix &	matB
	)

friend

zgsmatrix*zgsmatrix operator

Definition at line 137 of file zgsmatrix-zgsmatrix.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
   
   zgsmatrix newmat(matA.m, matB.n);
   
   const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::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;
     }
   }
   
   return _(newmat);
 }

_zgsmatrix operator*	(	const zgsmatrix &	matA,
		const _zgsmatrix &	matB
	)

friend

zgsmatrix*_zgsmatrix operator

Definition at line 139 of file zgsmatrix-_zgsmatrix.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
 
   zgsmatrix newmat(matA.m, matB.n);
   
   const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
   for(std::vector<zcomponent>::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;
     }
   }
   
   matB.destroy();
   return _(newmat);
 }

_zgsmatrix operator*	(	const zgsmatrix &	,
		const zhsmatrix &
	)

friend

_zgsmatrix operator*	(	const zgsmatrix &	,
		const _zhsmatrix &
	)

friend

_zgsmatrix operator*	(	const zgsmatrix &	mat,
		const double &	d
	)

friend

zgsmatrix*double operator

Definition at line 30 of file zgsmatrix-double.hpp.

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

_zgsmatrix operator*	(	const zgsmatrix &	mat,
		const comple &	d
	)

friend

zgsmatrix*comple operator

Definition at line 31 of file zgsmatrix-complex.hpp.

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

_zgsmatrix operator/	(	const zgsmatrix &	mat,
		const double &	d
	)

friend

zgsmatrix/double operator

Definition at line 44 of file zgsmatrix-double.hpp.

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

_zgsmatrix operator/	(	const zgsmatrix &	mat,
		const comple &	d
	)

friend

zgsmatrix/comple operator

Definition at line 45 of file zgsmatrix-complex.hpp.

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

_zgsmatrix operator*	(	const double &	d,
		const zgsmatrix &	mat
	)

friend

double*zgsmatrix operator

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

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

_zgsmatrix operator*	(	const comple &	d,
		const zgsmatrix &	mat
	)

friend

complex*zgsmatrix operator

Definition at line 3 of file complex-zgsmatrix.hpp.

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