set_timestep() [2/2]

template<typename T >

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

Sets all components at time step tstp to the components of a given herm_matrix_timestep.

Purpose

Sets all components of the herm_matrix at time step tstp to the components of given herm_matrix_timestep. If tstp = -1, only the Matsubara component will be copied.

Parameters

tstp	[int] The time step at which the components are set.
timestep	[herm_matrix_timestep] The herm_matrix_timestep from which the time step is copied.

Definition at line 2137 of file cntr_herm_matrix_impl.hpp.

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