resize_discard()

template<typename T >

void cntr::herm_matrix< T >::resize_discard	(	int	nt,
		int	ntau,
		int	size1
	)

Discards the herm_matrix object and resize with respect to the number of time points nt, points on the Matsubara branch ntau and the square-matrix size size1.

Purpose

Discard herm_matrix class and then resize with respect to number of time steps nt, points on the Matsubara branch ntau and the maxtix size size1. If nt = 0, real-time components are deallocated (i.e. only the Matsubara component is kept in memory). Internal routine; see top-level interface resize.

Parameters

nt	New number of time steps.
ntau	New number of points on Matsubara axis
size1	New matrix rank of the contour function

Definition at line 313 of file cntr_herm_matrix_impl.hpp.

                                                               {
    assert(ntau >= 0 && nt >= -1 && size1 >= 0);
    delete[] les_;
    delete[] ret_;
    delete[] tv_;
    delete[] mat_;
    nt_ = nt;
    ntau_ = ntau;
    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;
    }
}