set_matrixelement() [2/6]

template<typename T >

void cntr::herm_matrix< T >::set_matrixelement	(	int	tstp,
		int	i1,
		int	i2,
		herm_matrix< T > &	g
	)

Sets given matrix elements of all components at time step tstp to the components of a given herm_matrix of scaler type.

Purpose

Sets the (sub-) matrix elements (i1,i2) of all components at time step tstp to the components at time step tstp of a given herm_matrix.

Parameters

tstp	The time step at which the components are set.
i1	First index for submatrix.
i2	Second index for submatrix.
g	The herm_matrix from which the time step is copied.

Definition at line 2264 of file cntr_herm_matrix_impl.hpp.

                                                          {
    int i, sij = i1 * size2_ + i2;
    assert(tstp <= g.nt() && tstp <= nt_ && "tstp <= g.nt() && tstp <= nt_");
    assert(0 <= i1 && i1 < size1_ && 0 <= i2 && i2 < size2_ && g.size1() == 1
           && "0 <= i1 && i1 < size1_ && 0 <= i2 && i2 < size2_ && g.size1() == 1");
    if (tstp == -1) {
        for (i = 0; i <= ntau_; i++)
            matptr(i)[sij] = *(g.matptr(i));
    } else {
        for (i = 0; i <= tstp; i++)
            retptr(tstp, i)[sij] = *(g.retptr(tstp, i));
        for (i = 0; i <= ntau_; i++)
            tvptr(tstp, i)[sij] = *(g.tvptr(tstp, i));
        for (i = 0; i <= tstp; i++)
            lesptr(i, tstp)[sij] = *(g.lesptr(i, tstp));
    }
}