set_timestep_pseudo()

template<class GG >

void cntr::set_timestep_pseudo	(	int	tstp,
		GG &	G,
		std::complex< double > *	timestep
	)

Definition at line 1011 of file cntr_herm_pseudo_impl.hpp.

Referenced by increment_timestep_matrix_pseudo(), and set_timestep_matrix_pseudo().

                                                                     {
  int sg=G.element_size(),ntau=G.ntau(),m,l=0;
  std::complex<double> *x1;
  x1=timestep;
  if(tstp>=0){
    for(m=0;m<=tstp;m++){
      // ret(tstp,m) = gtr(tstp,m)
      // = -G(m-,tstp-)^*
      G.element_set(G.retptr(tstp,tstp-m),x1,G);
      G.element_smul(G.retptr(tstp,tstp-m),-1);
      G.element_conj(G.retptr(tstp,tstp-m));
      x1+=sg;l++;
    }
    for(m=0;m<=ntau;m++){
      // x1 points to Gvt(m,tstp) = -bosefermi *Gtv(tstp,ntau-m)^*
      G.element_set(G.tvptr(tstp,ntau-m),x1,G);
      G.element_conj(G.tvptr(tstp,ntau-m));
      G.element_smul(G.tvptr(tstp,ntau-m),-G.sig());
      x1+=sg;l++;
    }
    for(m=0;m<=tstp;m++){
      // x1 points to G(m+,tstp+)=les(m,tstp)
      G.element_set(G.lesptr(m,tstp),x1,G);
      x1+=sg;l++;
    }
  }else{
    for(m=0;m<=ntau;m++){
      // x1 points to Gmat(m,0) = i*Gmat(m)
      G.element_set(G.matptr(m),x1,G);
      G.element_smul(G.matptr(m),std::complex<double>(0,-1));
      x1+=sg;l++;
    }
  }
}

Here is the caller graph for this function: