get_timestep_pseudo()

template<class GG >

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

Definition at line 1046 of file cntr_herm_pseudo_impl.hpp.

Referenced by get_timestep_matrix_pseudo(), and increment_timestep_matrix_pseudo().

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

Here is the caller graph for this function: