increment_timestep_matrix_pseudo()

template<typename T , class Matrix >

void cntr::increment_timestep_matrix_pseudo	(	int	tstp,
		matrix< T > &	G,
		std::complex< T >	alpha,
		std::vector< Matrix > &	data
	)

Definition at line 1094 of file cntr_herm_pseudo_impl.hpp.

References get_timestep_pseudo(), and set_timestep_pseudo().

                                                                                                        {
  int i,l,m,sg,dim,ntau=G.ntau(),nmax=(tstp<0 ? ntau+1 : ntau+2*tstp+3);
  std::complex<T> *data1;
  assert((int)data.size()>=nmax);
  dim=G.size1();
  sg=dim*dim;
  data1=new std::complex<T> [nmax*sg];
  get_timestep_pseudo(tstp,G,data1);
  if(alpha.real()==1 && alpha.imag()==0) for(i=0;i<nmax;i++) for(l=0;l<dim;l++) for(m=0;m<dim;m++) data1[i*sg+l*dim+m] += data[i](l,m);
  else for(i=0;i<nmax;i++) for(l=0;l<dim;l++) for(m=0;m<dim;m++) data1[i*sg+l*dim+m] += alpha*data[i](l,m);
  set_timestep_pseudo(tstp,G,data1);
  delete [] data1;
}

Here is the call graph for this function: