herm_matrix() [2/5]

template<typename T >

cntr::herm_matrix< T >::herm_matrix	(	int	nt,
		int	ntau,
		int	size1 = 1,
		int	sig = -1
	)

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

Purpose

Initializes the herm_matrix class for a 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	Matrix rank of the contour function
sig	Set sig = -1 for fermions or sig = +1 for bosons.

Definition at line 62 of file cntr_herm_matrix_impl.hpp.

                                                                {
    assert(size1 >= 0 && nt >= -1 && sig * sig == 1 && ntau >= 0);
    nt_ = nt;
    ntau_ = ntau;
    sig_ = sig;
    size1_ = size1;
    size2_ = size1;
    element_size_ = size1 * size1;
    if (size1 > 0) {
        mat_ = new 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_];
    } else {
        les_ = 0;
        tv_ = 0;
        ret_ = 0;
    }
}