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

#include <_zhematrix.hpp>

Public Member Functions
	_zhematrix ()

	_zhematrix (const _zhematrix &)

	~_zhematrix ()

_zgematrix	to_zgematrix () const

zhecomplex	operator() (const CPPL_INT &, const CPPL_INT &) const

void	write (const char *) const

void	nullify () const

void	destroy () const

void	complete () const

Public Attributes
CPPL_INT const &	m
	matrix row and 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 zhematrix &)

_zhematrix	t (const _zhematrix &)

_zgematrix	i (const _zhematrix &)

const _zhematrix &	operator+ (const _zhematrix &)

_zhematrix	operator- (const _zhematrix &)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Detailed Description

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

Definition at line 3 of file _zhematrix.hpp.

Constructor & Destructor Documentation

_zhematrix::_zhematrix ( )

inline

_zhematrix constructor without arguments

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

References array, darray, and n.

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

_zhematrix::_zhematrix ( const _zhematrix & mat )

inline

_zhematrix copy constructor

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

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

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

_zhematrix::~_zhematrix ( )

inline

zhematrix destructor

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

References array, and darray.

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

Member Function Documentation

_zgematrix _zhematrix::to_zgematrix ( ) const

inline

convert to _zgematrix

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

References _(), destroy(), i, and n.

Referenced by operator*(), operator+(), and operator/().

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

zhecomplex _zhematrix::operator()	(	const CPPL_INT &	i,
		const CPPL_INT &	j
	)		const

inline

operator() for object

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

References darray, and n.

Referenced by write().

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

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

inline

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

References destroy(), i, imag(), 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 << "#zhematrix " << n << std::endl;
   for(CPPL_INT i=0; i<n; i++){
     for(CPPL_INT j=0; j<=i; j++ ){
       ofs << operator()(i,j) << " ";
     }
     ofs << std::endl;
     
 #ifdef  CPPL_DEBUG
     if(std::fabs(std::imag(operator()(i,i))) > DBL_MIN){
       WARNING_REPORT;
       std::cerr << "The " << i << "th diagonal component of the zhematrix is not a real number." << std::endl;
     }
 #endif//CPPL_DEBUG
   }
   
   ofs.close();
   destroy();
 }

void _zhematrix::nullify ( ) const

inline

nullify all the matrix data

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

References array, darray, and n.

Referenced by _zhematrix(), zhematrix::shallow_copy(), and zhematrix::zhematrix().

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

void _zhematrix::destroy ( ) const

inline

destroy all the matrix data

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

References array, and darray.

Referenced by operator*(), zgematrix::operator*=(), operator+(), zhematrix::operator+=(), zgematrix::operator+=(), operator-(), zhematrix::operator-=(), zgematrix::operator-=(), operator<<(), to_zgematrix(), and write().

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

void _zhematrix::complete ( ) const

inline

complete the upper-right components

Definition at line 22 of file _zhematrix-misc.hpp.

References conj(), darray, i, and n.

Referenced by operator*().

 {CPPL_VERBOSE_REPORT;
   for(CPPL_INT i=0; i<n; i++){
     for(CPPL_INT j=0; j<i; j++){
       darray[i][j] =std::conj(darray[j][i]);
     }
   }
 }

Friends And Related Function Documentation

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

friend

Definition at line 65 of file zhematrix-io.hpp.

 {CPPL_VERBOSE_REPORT;
   for(CPPL_INT i=0; i<mat.n; i++){
     for(CPPL_INT j=0; j<mat.n; j++){ 
       if(i>j){ s << " " << mat(i,j) << " "; }
       else if(i==j){ s << " " << std::real(mat(i,i)) << " "; }
       else{ s << "{" << std::conj(mat(j,i)) << "} "; }
     }
     s << std::endl;
     
 #ifdef  CPPL_DEBUG
     if(std::fabs(std::imag(mat(i,i))) > DBL_MIN){
       WARNING_REPORT;
       std::cerr << "The " << i << "th diagonal component of the zhematrix is not a real number." << std::endl;
     }
 #endif//CPPL_DEBUG
   }
   
   return s;
 }

_zhematrix t ( const _zhematrix & mat )

friend

return transposed _zhematrix

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

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

_zgematrix i ( const _zhematrix & mat )

friend

return its inverse matrix

Definition at line 15 of file _zhematrix-calc.hpp.

Referenced by complete(), to_zgematrix(), and write().

 {CPPL_VERBOSE_REPORT;
   zhematrix mat_cp(mat);
   zgematrix mat_inv(mat_cp.n,mat_cp.n);
   mat_inv.identity();
   mat_cp.zhesv(mat_inv);
   
   return _(mat_inv);
 }

const _zhematrix& operator+ ( const _zhematrix & mat )

friend

+_zhematrix operator

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

 {CPPL_VERBOSE_REPORT;
   return mat;
 }

_zhematrix operator- ( const _zhematrix & mat )

friend

-_zhematrix operator

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

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

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

friend

_zhematrix+zgematrix operator

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

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

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

friend

_zgematrix+zhematrix operator

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

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

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

friend

_zhematrix+zhematrix operator

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

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(matA.n!=matB.n){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a summation." << std::endl
               << "Your input was (" << matA.n << "x" << matA.n << ") + (" << matB.n << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   for(CPPL_INT j=0; j<matA.n; j++){
     for(CPPL_INT i=j; i<matA.n; i++){
       matA.darray[j][i] += matB.darray[j][i];
     }
   }
   
   return matA;
 }

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

friend

_zhematrix+_zhematrix operator

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

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(matA.n!=matB.n){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a summation." << std::endl
               << "Your input was (" << matA.n << "x" << matA.n << ") + (" << matB.n << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   for(CPPL_INT j=0; j<matA.n; j++){
     for(CPPL_INT i=j; i<matA.n; i++){
       matA.darray[j][i] += matB.darray[j][i];
     }
   }
   
   matB.destroy();
   return matA;
 }

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

friend

_zhematrix+zgbmatrix operator

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

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

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

friend

_zhematrix+_zgbmatrix operator

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

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

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

friend

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

friend

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

friend

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

friend

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

friend

_zhematrix-zgematrix operator

Definition at line 28 of file _zhematrix-zgematrix.hpp.

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

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

friend

_zhematrix-zgematrix operator

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

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

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

friend

_zhematrix-zhematrix operator

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

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(matA.n!=matB.n){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a subtraction." << std::endl
               << "Your input was (" << matA.n << "x" << matA.n << ") - (" << matB.n << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   for(CPPL_INT j=0; j<matA.n; j++){
     for(CPPL_INT i=j; i<matA.n; i++){
       matA.darray[j][i] -= matB.darray[j][i];
     }
   }
   
   return matA;
 }

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

friend

_zhematrix-_zhematrix operator

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

 {CPPL_VERBOSE_REPORT;
 #ifdef  CPPL_DEBUG
   if(matA.n!=matB.n){
     ERROR_REPORT;
     std::cerr << "These two matrises can not make a subtraction." << std::endl
               << "Your input was (" << matA.n << "x" << matA.n << ") - (" << matB.n << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
 
   for(CPPL_INT j=0; j<matA.n; j++){
     for(CPPL_INT i=j; i<matA.n; i++){
       matA.darray[j][i] -= matB.darray[j][i];
     }
   }
   
   matB.destroy();
   return matA;
 }

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

friend

_zhematrix-zgbmatrix operator

Definition at line 32 of file _zhematrix-zgbmatrix.hpp.

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

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

friend

_zhematrix-_zgbmatrix operator

Definition at line 33 of file _zhematrix-_zgbmatrix.hpp.

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

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

friend

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

friend

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

friend

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

friend

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

friend

_zhematrix*zcovector operator

Definition at line 3 of file _zhematrix-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.n << "x" << mat.n << ") * (" << vec.l << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zcovector newvec(mat.n);
   char uplo ='l';
   comple alpha =comple(1.,0.);
   CPPL_INT inc =1;
   comple beta =comple(0.,0.);
   
   zhemv_( &uplo, &mat.n, &alpha, mat.array, &mat.n, vec.array, &inc, &beta, newvec.array, &inc );
   
   mat.destroy();
   return _(newvec);
 }

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

friend

_zhematrix*_zcovector operator

Definition at line 3 of file _zhematrix-_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.n << "x" << mat.n << ") * (" << vec.l << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zcovector newvec(mat.n);
   char uplo ='l';
   comple alpha =comple(1.,0.);
   CPPL_INT inc =1;
   comple beta =comple(0.,0.);
 
   zhemv_( &uplo, &mat.n, &alpha, mat.array, &mat.n, vec.array, &inc, &beta, newvec.array, &inc );
   
   mat.destroy();
   vec.destroy();
   return _(newvec);
 }

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

friend

_zhematrix*zgematrix operator

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

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

friend

_zgematrix*zgematrix operator

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

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

friend

_zhematrix*zhematrix operator

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

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

friend

_zhematrix*_zhematrix operator

Definition at line 49 of file _zhematrix-_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.n << "x" << matA.n << ") * (" << matB.n << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   matB.complete();
   
   zgematrix newmat( matA.n, matB.n );
   char side ='l';
   char uplo ='l';
   comple alpha =comple(1.,0.);
   comple beta =comple(0.,0.);
   
   zhemm_( &side, &uplo, &matA.n, &matB.n, &alpha, matA.array, &matA.n, matB.array, &matB.n, &beta, newmat.array, &newmat.m );
   
   matA.destroy();
   matB.destroy();
   return _(newmat);
 }

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

friend

_zgematrix*zgbmatrix operator

Definition at line 61 of file _zhematrix-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.n << "x" << matA.n << ") * (" << matB.m << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zgematrix newmat( matA.n, 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 _zhematrix &	matA,
		const _zgbmatrix &	matB
	)

friend

_zhematrix*_zgbmatrix operator

Definition at line 63 of file _zhematrix-_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.n << "x" << matA.n << ") * (" << matB.m << "x" << matB.n << ")." << std::endl;
     exit(1);
   }
 #endif//CPPL_DEBUG
   
   zgematrix newmat( matA.n, 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 _zhematrix &	,
		const zgsmatrix &
	)

friend

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

friend

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

friend

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

friend

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

friend

_zhematrix*double operator

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

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

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

friend

_zhematrix*comple operator

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

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

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

friend

_zhematrix/double operator

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

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

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

friend

zhematrix/comple operator

Definition at line 16 of file _zhematrix-complex.hpp.

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

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

friend

double*_zhematrix operator

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

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

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

friend

comple*_zhematrix operator

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

 {CPPL_VERBOSE_REPORT;
   zgematrix newmat =mat.to_zgematrix();
   
   CPPL_INT size =mat.n*mat.n;
   CPPL_INT inc =1;
   zscal_(&size, &d, newmat.array, &inc);
   
   mat.destroy();
   return _(newmat);
 }