incr_timestep() [3/4]

template<typename T >

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

Adds a herm_matrix with given weight to the herm_matrix at given time step.

Purpose

Performs the operation \(C \rightarrow C + \alpha g\), where \(C\) is the herm_matrix, \(g\) is a time slice taken from a herm_matrix 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.
g	[herm_matrix] The herm_matrix \(g\) which is added to the herm_matrix at the time step.
alpha	[complex<T>] The weight in front of \(g\).

Definition at line 2669 of file cntr_herm_matrix_impl.hpp.

                                                        {
    int m;
    cplx *x, *x0;
    assert(tstp >= -1 && tstp <= nt_ && "tstp >= -1 && tstp <= nt_");
    assert(g.nt_ >= tstp && g.ntau_ == ntau_ && g.size1_ == size1_
       && "g.nt_ >= tstp && g.ntau_ == ntau_ && g.size1_ == size1_");
    if (tstp == -1) {
        x0 = matptr(0);
        x = g.matptr(0);
        for (m = 0; m <= ntau_; m++) {
            element_incr<T, LARGESIZE>(size1_, x0 + m * element_size_, alpha,
                                       x + m * element_size_);
        }
    } else {
        x0 = retptr(tstp, 0);
        x = g.retptr(tstp, 0);
        for (m = 0; m <= tstp; m++) {
            element_incr<T, LARGESIZE>(size1_, x0 + m * element_size_, alpha,
                                       x + m * element_size_);
        }
        x0 = tvptr(tstp, 0);
        x = g.tvptr(tstp, 0);
        for (m = 0; m <= ntau_; m++) {
            element_incr<T, LARGESIZE>(size1_, x0 + m * element_size_, alpha,
                                       x + m * element_size_);
        }
        x0 = lesptr(0, tstp);
        x = g.lesptr(0, tstp);
        for (m = 0; m <= tstp; m++) {
            element_incr<T, LARGESIZE>(size1_, x0 + m * element_size_, alpha,
                                       x + m * element_size_);
        }
    }
}