(DO NOT USE) Smart-temporary Complex Double-precision General Dence Matrix Class More...

#include <_zgematrix.hpp>

Public Member Functions
	_zgematrix ()

	_zgematrix (const _zgematrix &)

	~_zgematrix ()

comple &	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...

comple *	array
	1D array to store matrix data More...

comple **	darray
	array of pointers of column head addresses More...

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

_zgematrix	t (const _zgematrix &)

_zgematrix	i (const _zgematrix &)

_zgematrix	conj (const _zgematrix &)

_zgematrix	conjt (const _zgematrix &)

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

comple	damax (const _zgematrix &)

const _zgematrix &	operator+ (const _zgematrix &)

_zgematrix	operator- (const _zgematrix &)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

_zgematrix	operator/ (const _zgematrix &, const comple &)

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

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

Detailed Description

(DO NOT USE) Smart-temporary Complex Double-precision General Dence Matrix Class

Definition at line 3 of file _zgematrix.hpp.

Constructor & Destructor Documentation

_zgematrix::_zgematrix ( )

inline

_zgematrix constructor without arguments

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

References array, darray, m, and n.

 {CPPL_VERBOSE_REPORT;
   //////// initialize ////////
   m =0;
   n =0;
   array =NULL;
   darray =NULL;
 }

_zgematrix::_zgematrix ( const _zgematrix & mat )

inline

_zgematrix copy constructor

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

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

 {CPPL_VERBOSE_REPORT;
   //////// initialize ////////
   m =mat.m;
   n =mat.n;
   array =mat.array;
   darray =mat.darray;
   
   mat.nullify();
 }

_zgematrix::~_zgematrix ( )

inline

zgematrix destructor

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

References array, and darray.

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

Member Function Documentation

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

inline

operator() for object

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

References darray, i, m, and n.

Referenced by write().

 {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
   
   //return array[i+m*j];
   return darray[j][i];
 }

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

inline

Definition at line 41 of file _zgematrix-io.hpp.

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

 {CPPL_VERBOSE_REPORT;
   std::ofstream ofs(filename, std::ios::trunc);
   ofs.setf(std::cout.flags());
   ofs.precision(std::cout.precision());
   ofs.width(std::cout.width());
   ofs.fill(std::cout.fill());
   
   ofs << "#zgematrix" << " " << m << " " << n << std::endl;
   for(CPPL_INT i=0; i<m; i++){
     for(CPPL_INT j=0; j<n; j++){
       ofs << operator()(i,j) << " ";
     }
     ofs << std::endl;
   }
   
   ofs.close();
   destroy();
 }

void _zgematrix::nullify ( ) const

inline

nullify all the matrix data

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

References array, darray, m, and n.

Referenced by _zgematrix(), zgematrix::shallow_copy(), and zgematrix::zgematrix().

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

void _zgematrix::destroy ( ) const

inline

destroy all the matrix data

Definition at line 13 of file _zgematrix-misc.hpp.

References array, and darray.

Referenced by conjt(), damax(), idamax(), operator*(), zgematrix::operator*=(), operator+(), zgematrix::operator+=(), operator-(), zgematrix::operator-=(), operator<<(), t(), and write().

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

Friends And Related Function Documentation

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

friend

Definition at line 63 of file zgematrix-io.hpp.

 {CPPL_VERBOSE_REPORT;
   for(CPPL_INT i=0; i<mat.m; i++){
     for(CPPL_INT j=0; j<mat.n; j++){
       s << " " << mat(i,j);
     }
     s << std::endl;
   }
   return s;
 }

_zgematrix t ( const _zgematrix & mat )

friend

return transposed zgematrix

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

 {CPPL_VERBOSE_REPORT;
   zgematrix newmat(mat.n,mat.m);
   
   for(CPPL_INT i=0; i<newmat.m; i++){
     for(CPPL_INT j=0; j<newmat.n; j++){
       newmat(i,j) =mat(j,i);
     }
   }
   
   mat.destroy();
   return _(newmat);
 }

_zgematrix i ( const _zgematrix & mat )

friend

return its inverse matrix

Definition at line 19 of file _zgematrix-calc.hpp.

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

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

_zgematrix conj ( const _zgematrix & mat )

friend

return its conjugate matrix

Definition at line 44 of file _zgematrix-calc.hpp.

 {CPPL_VERBOSE_REPORT;
   for(CPPL_INT i=0; i<mat.m; i++){
     for(CPPL_INT j=0; j<mat.n; j++){
       mat(i,j) =std::conj(mat(i,j));
     }
   }
   
   return mat;
 }

_zgematrix conjt ( const _zgematrix & mat )

friend

return its conjugate transposed matrix

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

 {CPPL_VERBOSE_REPORT;
   zgematrix newmat(mat.n,mat.m);
   
   for(CPPL_INT i=0; i<newmat.m; i++){
     for(CPPL_INT j=0; j<newmat.n; j++){
       newmat(i,j) =std::conj(mat(j,i));
     }
   }
   
   mat.destroy();
   return _(newmat);
 }

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

friend

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

Definition at line 78 of file _zgematrix-calc.hpp.

 {CPPL_VERBOSE_REPORT;
   CPPL_INT size =mat.m*mat.n;
   CPPL_INT inc =1;
   CPPL_INT index =izamax_(&size, mat.array, &inc) -1;
   i =index%mat.m;
   j =index/mat.m;
   
   mat.destroy();
 }

comple damax ( const _zgematrix & mat )

friend

return its largest absolute value

Definition at line 91 of file _zgematrix-calc.hpp.

 {CPPL_VERBOSE_REPORT;
   CPPL_INT size =mat.m*mat.n;
   CPPL_INT inc =1;
   comple val =mat.array[izamax_(&size, mat.array, &inc) -1];
   
   mat.destroy();
   return val;
 }

const _zgematrix& operator+ ( const _zgematrix & mat )

friend

+_zgematrix operator

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

 {CPPL_VERBOSE_REPORT;
   return mat;
 }

_zgematrix operator- ( const _zgematrix & mat )

friend

-_zgematrix operator

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

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

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

friend

_zgematrix+zgematrix operator

Definition at line 3 of file _zgematrix-zgematrix.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
   
   const CPPL_INT size =matA.m*matA.n;
   for(CPPL_INT i=0; i<size; i++){
     matA.array[i] += matB.array[i];
   }
   
   return matA;
 }

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

friend

_zgematrix+_zgematrix operator

Definition at line 3 of file _zgematrix-_zgematrix.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
   
   const CPPL_INT size =matA.m*matA.n;
   for(CPPL_INT i=0; i<size; i++){
     matA.array[i] += matB.array[i];
   }
 
   matB.destroy();
   return matA;
 }

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

friend

_zgematrix+zhematrix operator

Definition at line 3 of file _zgematrix-zhematrix.hpp.

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(matA.n!=matB.n || matA.m!=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
 
   for(CPPL_INT i=0; i<matB.n; i++){
     for(CPPL_INT j=0; j<matB.n; j++){
       matA(i,j) += matB(i,j);
     }
   }
   
   return matA;
 }

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

friend

_zgematrix+_zhematrix operator

Definition at line 3 of file _zgematrix-_zhematrix.hpp.

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(matA.n!=matB.n || matA.m!=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
   
   for(CPPL_INT i=0; i<matB.n; i++){
     for(CPPL_INT j=0; j<matB.n; j++){
       matA(i,j) += matB(i,j);
     }
   }
   
   matB.destroy();
   return matA;
 }

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

friend

_zgematrix+zgbmatrix operator

Definition at line 3 of file _zgematrix-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
   
   for(CPPL_INT i=0; i<matB.m; i++){
     const CPPL_INT jmax =std::min(matB.n,i+matB.ku+1);
     for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax; j++){
       matA(i,j) += matB(i,j);
     }
   }
   
   return matA;
 }

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

friend

_zgematrix+_zgbmatrix operator

Definition at line 3 of file _zgematrix-_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
   
   for(CPPL_INT i=0; i<matB.m; i++){
     const CPPL_INT jmax =std::min(matB.n,i+matB.ku+1);
     for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax; j++){
       matA(i,j)+=matB(i,j);
     }
   }
   
   matB.destroy();
   return matA;
 }

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

friend

_zgematrix+zgsmatrix operator

Definition at line 3 of file _zgematrix-zgsmatrix.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 size_t matB_data_size =matB.data.size();
   for(size_t c=0; c<matB_data_size; c++){
     const zcomponent& z =matB.data[c];
     matA(z.i,z.j) += z.v;
   }
   
   return matA;
 }

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

friend

_zgematrix+_zgsmatrix operator

Definition at line 3 of file _zgematrix-_zgsmatrix.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 size_t matB_data_size =matB.data.size();
   for(size_t c=0; c<matB_data_size; c++){
     const zcomponent& z =matB.data[c];
     matA(z.i,z.j) += z.v;
   }
   
   matB.destroy();
   return matA;
 }

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

friend

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

friend

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

friend

_zgematrix-zgematrix operator

Definition at line 24 of file _zgematrix-zgematrix.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
 
   const CPPL_INT size =matA.m*matA.n;
   for(CPPL_INT i=0; i<size; i++){
     matA.array[i] -= matB.array[i];
   }
   
   return matA;
 }

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

friend

_zgematrix-_zgematrix operator

Definition at line 25 of file _zgematrix-_zgematrix.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
   
   const CPPL_INT size =matA.m*matA.n;
   for(CPPL_INT i=0; i<size; i++){
     matA.array[i]-=matB.array[i];
   }
   
   matB.destroy();
   return matA;
 }

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

friend

_zgematrix-zhematrix operator

Definition at line 25 of file _zgematrix-zhematrix.hpp.

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(matA.n!=matB.n || matA.m!=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
   
   for(CPPL_INT i=0; i<matB.n; i++){
     for(CPPL_INT j=0; j<matB.n; j++){
       matA(i,j) -= matB(i,j);
     }
   }
   
   return matA;
 }

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

friend

_zgematrix-_zhematrix operator

Definition at line 26 of file _zgematrix-_zhematrix.hpp.

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(matA.n!=matB.n || matA.m!=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
   
   for(CPPL_INT i=0; i<matB.n; i++){
     for(CPPL_INT j=0; j<matB.n; j++){
       matA(i,j) -= matB(i,j);
     }
   }
   
   matB.destroy();
   return matA;
 }

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

friend

_zgematrix-zgbmatrix operator

Definition at line 26 of file _zgematrix-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
   
   for(CPPL_INT i=0; i<matB.m; i++){
     const CPPL_INT jmax =std::min(matB.n,i+matB.ku+1);
     for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax; j++){
       matA(i,j) -= matB(i,j);
     }
   }
   
   return matA;
 }

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

friend

_zgematrix-_zgbmatrix operator

Definition at line 27 of file _zgematrix-_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
   
   for(CPPL_INT i=0; i<matB.m; i++){
     const CPPL_INT jmax =std::min(matB.n,i+matB.ku+1);
     for(CPPL_INT j=std::max(CPPL_INT(0),i-matB.kl); j<jmax; j++){
       matA(i,j)-=matB(i,j);
     }
   }
   
   matB.destroy();
   return matA;
 }

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

friend

_zgematrix-zgsmatrix operator

Definition at line 25 of file _zgematrix-zgsmatrix.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
   
   const size_t matB_data_size =matB.data.size();
   for(size_t c=0; c<matB_data_size; c++){
     const zcomponent& z =matB.data[c];
     matA(z.i,z.j) -= z.v;
   }
   
   return matA;
 }

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

friend

_zgematrix-_zgsmatrix operator

Definition at line 26 of file _zgematrix-_zgsmatrix.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
   
   const size_t matB_data_size =matB.data.size();
   for(size_t c=0; c<matB_data_size; c++){
     const zcomponent& z =matB.data[c];
     matA(z.i,z.j) -= z.v;
   }
   
   matB.destroy();
   return matA;
 }

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

friend

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

friend

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

friend

_zgematrix*zcovector operator

Definition at line 3 of file _zgematrix-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);
   char trans ='n';
   comple alpha =comple(1.,0.);
   CPPL_INT inc =1;
   comple beta =comple(0.,0.);
   
   zgemv_( &trans, &mat.m, &mat.n, &alpha, mat.array, &mat.m, vec.array, &inc, &beta, newvec.array, &inc );
   
   mat.destroy();
   return _(newvec);
 }

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

friend

_zgematrix*_zcovector operator

Definition at line 3 of file _zgematrix-_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);
   char trans ='n';
   comple alpha =comple(1.,0.);
   CPPL_INT inc =1;
   comple beta =comple(0.,0.);
   
   zgemv_( &trans, &mat.m, &mat.n, &alpha, mat.array, &mat.m, vec.array, &inc, &beta, newvec.array, &inc );
   
   mat.destroy();
   vec.destroy();
   return _(newvec);
 }

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

friend

_zgematrix*zgematrix operator

Definition at line 45 of file _zgematrix-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.m << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zgematrix newmat( matA.m, matB.n );
   char transa ='n';
   char transb ='n';
   comple alpha =comple(1.,0.);
   comple beta =comple(0.,0.);
   
   zgemm_( &transa, &transb, &matA.m, &matB.n, &matA.n, &alpha, matA.array, &matA.m, matB.array, &matB.m, &beta, newmat.array, &matA.m );
   
   matA.destroy();
   return _(newmat);
 }

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

friend

_zgematrix*_zgematrix operator

Definition at line 47 of file _zgematrix-_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.m << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zgematrix newmat( matA.m, matB.n );
   char transa ='n';
   char transb ='n';
   comple alpha =comple(1.,0.);
   comple beta =comple(0.,0.);
 
   zgemm_( &transa, &transb, &matA.m, &matB.n, &matA.n, &alpha, matA.array, &matA.m, matB.array, &matB.m, &beta, newmat.array, &matA.m );
   
   matA.destroy();
   matB.destroy();
   return _(newmat);
 }

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

friend

_zgematrix*zhematrix operator

Definition at line 47 of file _zgematrix-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 );
   char side ='R';
   char uplo ='l';
   comple alpha =comple(1.,0.);
   comple beta =comple(0.,0.);
   
   zhemm_( &side, &uplo, &newmat.m, &newmat.n, &alpha, matB.array, &matB.n, matA.array, &matA.m, &beta, newmat.array, &newmat.m );
   
   matA.destroy();
   return _(newmat);
 }

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

friend

_zgematrix*_zhematrix operator

Definition at line 49 of file _zgematrix-_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 );
   char side ='R';
   char uplo ='l';
   comple alpha =comple(1.,0.);
   comple beta =comple(0.,0.);
 
   zhemm_( &side, &uplo, &newmat.m, &newmat.n, &alpha, matB.array, &matB.n, matA.array, &matA.m, &beta, newmat.array, &newmat.m );
   
   matA.destroy();
   matB.destroy();
   return _(newmat);
 }

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

friend

_zgematrix*zgbmatrix operator

Definition at line 49 of file _zgematrix-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();
   
   for(CPPL_INT i=0; i<newmat.m; i++){
     for(CPPL_INT j=0; j<newmat.n; j++){
       const CPPL_INT kmax =std::min(matB.m,j+matB.kl+1);
       for(CPPL_INT k=std::max(CPPL_INT(0),j-matB.ku); k<kmax; k++){
         newmat(i,j) += matA(i,k)*matB(k,j);
       }
     }
   }
   
   matA.destroy();
   return _(newmat);
 }

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

friend

_zgematrix*_zgbmatrix operator

Definition at line 51 of file _zgematrix-_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();
   
   for(CPPL_INT i=0; i<newmat.m; i++){
     for(CPPL_INT j=0; j<newmat.n; j++){
       const CPPL_INT kmax =std::min(matB.m,j+matB.kl+1);
       for(CPPL_INT k=std::max(CPPL_INT(0),j-matB.ku); k<kmax; k++){
         newmat(i,j)+=matA(i,k)*matB(k,j);
       }
     }
   }
   
   matA.destroy();
   matB.destroy();
   return _(newmat);
 }

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

friend

_zgematrix*zgsmatrix operator

Definition at line 47 of file _zgematrix-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.n << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zgematrix newmat(matA.m, matB.n);
   newmat.zero();
   
   const size_t matB_data_size =matB.data.size();
   for(size_t c=0; c<matB_data_size; c++){
     const zcomponent& z =matB.data[c];
     for(CPPL_INT i=0; i<matA.m; i++){
       newmat(i,z.j) += matA(i,z.i)*z.v;
     }
   }
   
   matA.destroy();
   return _(newmat);
 }

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

friend

_zgematrix*_zgsmatrix operator

Definition at line 49 of file _zgematrix-_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.n << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zgematrix newmat(matA.m, matB.n);
   newmat.zero();
   
   const size_t matB_data_size =matB.data.size();
   for(size_t c=0; c<matB_data_size; c++){
     const zcomponent& z =matB.data[c];
     for(CPPL_INT i=0; i<matA.m; i++){
       newmat(i,z.j) += matA(i,z.i)*z.v;
     }
   }
   
   matA.destroy();
   matB.destroy();
   return _(newmat);
 }

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

friend

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

friend

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

friend

_zgematrix*double operator

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

 {CPPL_VERBOSE_REPORT;
   CPPL_INT size =mat.m*mat.n;
   CPPL_INT inc =1;
   zdscal_(&size, &d, mat.array, &inc);
   return mat;
 }

_zgematrix operator*	(	const _zgematrix &	mat,
		const comple &	d
	)

friend

_zgematrix*comple operator

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

 {CPPL_VERBOSE_REPORT;
   CPPL_INT size =mat.m*mat.n;
   CPPL_INT inc =1;
   zscal_(&size, &d, mat.array, &inc);
   return mat;
 }

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

friend

_zgematrix/double operator

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

 {CPPL_VERBOSE_REPORT;
   CPPL_INT size =mat.m*mat.n;
   double dinv =1./d;
   CPPL_INT inc =1;
   zdscal_(&size, &dinv, mat.array, &inc);
   return mat;
 }

_zgematrix operator/	(	const _zgematrix &	mat,
		const comple &	d
	)

friend

_zgematrix/comple operator

Definition at line 13 of file _zgematrix-complex.hpp.

 {CPPL_VERBOSE_REPORT;
   CPPL_INT size =mat.m*mat.n;
   comple dinv =1./d;
   CPPL_INT inc =1;
   zscal_(&size, &dinv, mat.array, &inc);
   return mat;
 }

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

friend

double*_zgematrix operator

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

 {CPPL_VERBOSE_REPORT;
   CPPL_INT size =mat.m*mat.n;
   CPPL_INT inc =1;
   zdscal_(&size, &d, mat.array, &inc);
   return mat;
 }

_zgematrix operator*	(	const comple &	d,
		const _zgematrix &	mat
	)

friend

comple*_zgematrix operator

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

 {CPPL_VERBOSE_REPORT;
   CPPL_INT size =mat.m*mat.n;
   CPPL_INT inc =1;
   zscal_(&size, &d, mat.array, &inc);
   return mat;
 }

Member Data Documentation

CPPL_INT _zgematrix::m

mutable

matrix row size

Definition at line 9 of file _zgematrix.hpp.

Referenced by _(), _zgematrix(), conj(), conjt(), damax(), i(), idamax(), nullify(), operator()(), operator*(), zgematrix::operator*=(), operator+(), zgematrix::operator+=(), operator-(), zgematrix::operator-=(), operator/(), operator<<(), zgematrix::shallow_copy(), t(), write(), and zgematrix::zgematrix().

CPPL_INT _zgematrix::n

mutable

matrix column size

Definition at line 10 of file _zgematrix.hpp.

Referenced by _(), _zgematrix(), conj(), conjt(), damax(), i(), idamax(), nullify(), operator()(), operator*(), zgematrix::operator*=(), operator+(), zgematrix::operator+=(), operator-(), zgematrix::operator-=(), operator/(), operator<<(), zgematrix::shallow_copy(), t(), write(), and zgematrix::zgematrix().

comple* _zgematrix::array

mutable

1D array to store matrix data

Definition at line 11 of file _zgematrix.hpp.

Referenced by _(), _zgematrix(), damax(), destroy(), idamax(), nullify(), operator*(), zgematrix::operator*=(), operator+(), zgematrix::operator+=(), operator-(), zgematrix::operator-=(), operator/(), zgematrix::shallow_copy(), zgematrix::zgematrix(), and ~_zgematrix().

comple** _zgematrix::darray

mutable

array of pointers of column head addresses

Definition at line 12 of file _zgematrix.hpp.

Referenced by _(), _zgematrix(), destroy(), nullify(), operator()(), zgematrix::shallow_copy(), zgematrix::zgematrix(), and ~_zgematrix().

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

Public Member Functions

Public Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation