herm_matrix() [3/5]

template<typename T >

cntr::herm_matrix< T >::herm_matrix	(	int	nt,
		int	ntau,
		int	size1,
		int	size2,
		int	sig
	)

Initializes the herm_matrix class for a general matrix two-time contour function.

Purpose

Initializes the herm_matrix class for a general (square or non-square) matrix two-time contour function. If nt = 0, memory will be allocated for the Matsubara component only.

Parameters

nt	Number of time steps
ntau	Number of points on Matsubara axis
size1	Number of matrix rows
size2	Number of matrix columns
sig	Set sig = -1 for fermions or sig = +1 for bosons.

Definition at line 111 of file cntr_herm_matrix_impl.hpp.

                                                                                            {
   assert(size1>=0 && size2>=0 && nt>=-1 && sig*sig==1 && ntau>=0);
   nt_=nt;
   ntau_=ntau;
   sig_=sig;
   size1_=size1;
   size2_=size2;
   element_size_=size1*size2;
   if(size1>0){
       mat_ = new cplx [(ntau_+1)*element_size_];
       memset(mat_, 0, sizeof(cplx)*(ntau_+1)*element_size_);
   }else{
       mat_=0;
   }
   if(nt>=0 && size1>0){
       les_ = new cplx [((nt_+1)*(nt_+2))/2*element_size_];
       ret_ = new cplx [((nt_+1)*(nt_+2))/2*element_size_];
       tv_ = new cplx [(nt_+1)*(ntau_+1)*element_size_];
       memset(les_, 0, sizeof(cplx)*((nt_+1)*(nt_+2))/2*element_size_);
       memset(ret_, 0, sizeof(cplx)*((nt_+1)*(nt_+2))/2*element_size_);
       memset(tv_, 0, sizeof(cplx)*(nt_+1)*(ntau_+1)*element_size_);
   }else{
       les_=0;
       tv_=0;
       ret_=0;
   }
}