resize_nt()

template<typename T >

void cntr::herm_matrix< T >::resize_nt

(

int

nt

)

Resizes herm_matrix object with respect to the number of time points nt. Internal routine; see top-level interface resize.

Purpose

Resizes herm_matrix class with respect to number of time steps nt. If nt >= 0 real-time components are resized, otherwise deallocated (i.e. only the Matsubara component is kept in memory).

Parameters

nt	New number of time steps.

Definition at line 361 of file cntr_herm_matrix_impl.hpp.

                                     {
    int nt1 = (nt_ > nt ? nt : nt_);
    cplx *ret, *les, *tv;
    assert(nt >= -1);
    nt_ = nt;
    if (size1_ == 0)
        return;
    if (nt_ >= 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_];
        if (nt1 >= 0) {
            memcpy(les, les_, sizeof(cplx) * ((nt1 + 1) * (nt1 + 2)) / 2 *
                                  element_size_);
            memcpy(ret, ret_, sizeof(cplx) * ((nt1 + 1) * (nt1 + 2)) / 2 *
                                  element_size_);
            memcpy(tv, tv_,
                   sizeof(cplx) * (nt1 + 1) * (ntau_ + 1) * element_size_);
        }
    } else {
        les = 0;
        ret = 0;
        tv = 0;
    }
    delete[] les_;
    delete[] ret_;
    delete[] tv_;
    les_ = les;
    ret_ = ret;
    tv_ = tv;
}