get_timestep() [1/2]

template<typename T >

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

Definition at line 2156 of file cntr_herm_matrix_impl.hpp.

Referenced by cntr::Reduce_timestep().

                                                                           {
    int len = (2 * (tstp + 1) + ntau_ + 1) * element_size_;
    cplx *x;
    assert(tstp >= -1 && tstp <= nt_ && "tstp >= -1 && tstp <= nt_");
    if (timestep.total_size_ < len)
        timestep.resize(tstp, ntau_, size1_);
    x = timestep.data_;
    timestep.tstp_ = tstp;
    timestep.ntau_ = ntau_;
    timestep.size1_ = size1_;
    if (tstp == -1) {
        memcpy(x, mat_, sizeof(cplx) * (ntau_ + 1) * element_size_);
    } else {
        memcpy(x, retptr(tstp, 0), sizeof(cplx) * (tstp + 1) * element_size_);
        memcpy(x + (tstp + 1) * element_size_, tvptr(tstp, 0),
               sizeof(cplx) * (ntau_ + 1) * element_size_);
        memcpy(x + (tstp + 1 + ntau_ + 1) * element_size_, lesptr(0, tstp),
               sizeof(cplx) * (tstp + 1) * element_size_);
    }
}

Here is the caller graph for this function: