incr_timestep() [1/4]

template<typename T >

void cntr::herm_matrix< T >::incr_timestep	(	int	tstp,
		herm_matrix_timestep< T > &	timestep,
		cplx	alpha
	)

Adds a herm_matrix_timestep with given weight to the herm_matrix.

Purpose

Performs the operation \(C \rightarrow C + \alpha A\), where \(C\) is the herm_matrix, \(A\) is a time slice described by a herm_matrix_timestep and \(\alpha\) is a complex weight. The operation is performed at given time step tstp.

Parameters

tstp	[int] The time step where the herm_matrix is increased.
timestep	[herm_matrix_timestep] The herm_matrix_timestep which is added to the herm_matrix.
alpha	[complex<T>] The weight in front of timestep.

Definition at line 2567 of file cntr_herm_matrix_impl.hpp.

Referenced by cntr::vie2_timestep_sin(), and cntr::vie2_timestep_sin_omp().

                                                        {
    int i, len;
    cplx *x, *x0;
    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) {
        len = (ntau_ + 1) * element_size_;
        x0 = matptr(0);
        x = timestep.data_;
        HERM_MATRIX_INCR_TSTP
    } else {
        len = (tstp + 1) * element_size_;
        x0 = retptr(tstp, 0);
        x = timestep.data_;
        HERM_MATRIX_INCR_TSTP
        len = (ntau_ + 1) * element_size_;
        x0 = tvptr(tstp, 0);
        x = timestep.data_ + (tstp + 1) * element_size_;
        HERM_MATRIX_INCR_TSTP
        len = (tstp + 1) * element_size_;
        x0 = lesptr(0, tstp);
        x = timestep.data_ + (tstp + 1 + ntau_ + 1) * element_size_;
        HERM_MATRIX_INCR_TSTP
    }
}

Here is the caller graph for this function: