set_matrixelement() [4/6]

template<typename T >

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

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

Purpose

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

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.
j1	First index for submatrix of g.
j2	Second index for submatrix of g.

Definition at line 2365 of file cntr_herm_matrix_impl.hpp.

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