get_timestep() [2/2]

template<typename T >

void cntr::herm_matrix< T >::get_timestep	(	int	tstp,
		herm_matrix< T > &	timestep
	)		const

Definition at line 2179 of file cntr_herm_matrix_impl.hpp.

                                                                    {

    assert(tstp >= -1 && tstp <= nt_ && tstp <= g1.nt() && "tstp >= -1 && tstp <= nt_ && tstp <= g1.nt()");
    assert(g1.size1() == size1_ && "g1.size1() == size1_");
    assert(g1.ntau() == ntau_ && "g1.ntau() == ntau_");
    if (tstp == -1) {
        memcpy(g1.mat_, mat_, sizeof(cplx) * (ntau_ + 1) * element_size_);
    } else {
        memcpy(g1.retptr(tstp, 0), retptr(tstp, 0),
               sizeof(cplx) * (tstp + 1) * element_size_);
        memcpy(g1.tvptr(tstp, 0), tvptr(tstp, 0),
               sizeof(cplx) * (ntau_ + 1) * element_size_);
        memcpy(g1.lesptr(0, tstp), lesptr(0, tstp),
               sizeof(cplx) * (tstp + 1) * element_size_);
    }
}