distributed_timestep_array() [3/3]

template<typename T >

cntr::distributed_timestep_array< T >::distributed_timestep_array	(	int	n,
		int	nt,
		int	ntau,
		int	size,
		int	sig,
		bool	mpi
	)

Initializes the distributed_timestep_array class.

Purpose

Initializes the distributed_timestep_array class including n blocks of the class herm_matrix_ti* mestep_view, which serves as an interface to herm_matrix_timestep without copying the data. At * the initialization the pointers don't point anywhere yet.

Parameters

n	Number of blocks
nt	The maximum allowed number of timestep (set by size of data_)
ntau	Number of points on Matsubara axis
size	Matrix rank of the herm_matrix_timestep function
sig	Set sig = -1 for fermions or sig = +1 for bosons
mpi	If 'true' use MPI, otherwise one task with tid_=0

Definition at line 87 of file cntr_distributed_timestep_array_impl.hpp.

References cntr::distributed_array< T >::n().

                                                                                                                              {
    assert(-1<=nt && 0<=ntau &&  sig*sig==1 && 1<=size);
    int size_tstp=(ntau+1+2*(nt+1))*size*size;
    int maxlen=size_tstp;
    data_=cntr::distributed_array<std::complex<T> >(n,maxlen,mpi);
    n_=data_.n();   
    tid_=data_.tid();
    ntasks_=data_.ntasks(); 
    tstp_=-2;
    nt_=nt;
    ntau_=ntau;
    size_=size;
    sig_=sig;
    // they all point nowhere, since tstp=-2
    G_=std::vector<cntr::herm_matrix_timestep_view<T> >(n_);
}

Here is the call graph for this function: