incr_timestep() [2/2]

template<typename T >

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

Increase the value of the herm_matrix_timestep by \(\alpha g(t)\)

Purpose

Increase the value of herm_matrix_timestep by a value of \(\alpha g(t) \), where \( g(t)\) is a herm_matrix and \(\alpha \) is a complex number. If \(t>-1\) then ret,les,tv components are set, otherwise mat. Works for scalar or square-matrix contour objects.

Parameters

tstp	time step tstp, identical to the time step of herm_matrix_timestep e(dummy argument in release mode)
g	The herm_matrix which is added
alpha	Constant multiplication factor

Definition at line 2019 of file cntr_herm_matrix_timestep_impl.hpp.

References cntr::herm_matrix< T >::nt(), cntr::herm_matrix< T >::ntau(), and cntr::herm_matrix< T >::size1().

                                                                                {
   assert(tstp == tstp_);
   assert(tstp_ <= g.nt() && ntau_ == g.ntau() && size1_ == g.size1());
   int i, len;
   cplx *x, *x0;
   if (tstp_ == -1) {
      len = (ntau_ + 1) * element_size_;
      x = g.matptr(0);
      x0 = data_;
      HERM_MATRIX_INCR_TSTP
   } else {
      len = (tstp_ + 1) * element_size_;
      x = g.retptr(tstp_, 0);
      x0 = data_;
      HERM_MATRIX_INCR_TSTP
      len = (ntau_ + 1) * element_size_;
      x = g.tvptr(tstp_, 0);
      x0 = data_ + (tstp_ + 1) * element_size_;
      HERM_MATRIX_INCR_TSTP
      len = (tstp_ + 1) * element_size_;
      x = g.lesptr(0, tstp_);
      x0 = data_ + (tstp_ + 1 + ntau_ + 1) * element_size_;
      HERM_MATRIX_INCR_TSTP
   }
}

Here is the call graph for this function: