herm_matrix() [4/5]

template<typename T >

cntr::herm_matrix< T >::herm_matrix

(

const herm_matrix< T > &

g

)

Initializes the herm_matrix class with the same layout as a given herm_matrix.

Purpose

Initializes the herm_matrix class with the same number of time steps nt, number of points on the imaginary branch ntau, matrix rank size1 and bosonic/fermionic symmetry sig. Works for scalar or square-matrix contour objects only.

Parameters

g	The herm_matrix according to which the class should be initialized

Definition at line 157 of file cntr_herm_matrix_impl.hpp.

                                                {
    nt_ = g.nt_;
    ntau_ = g.ntau_;
    sig_ = g.sig_;
    size1_ = g.size1_;
    size2_ = g.size1_;
    element_size_ = size1_ * size1_;
    if (size1_ > 0) {
        mat_ = new cplx[(ntau_ + 1) * element_size_];
        memcpy(mat_, g.mat_, 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_];
        memcpy(les_, g.les_,
               sizeof(cplx) * ((nt_ + 1) * (nt_ + 2)) / 2 * element_size_);
        memcpy(ret_, g.ret_,
               sizeof(cplx) * ((nt_ + 1) * (nt_ + 2)) / 2 * element_size_);
        memcpy(tv_, g.tv_,
               sizeof(cplx) * (nt_ + 1) * (ntau_ + 1) * element_size_);
    } else {
        les_ = 0;
        ret_ = 0;
        tv_ = 0;
    }
}