green_equilibrium_mat()

template<typename T , class dos_function >

void cntr::green_equilibrium_mat	(	herm_matrix< T > &	G,
		dos_function &	dos,
		double	beta,
		int	limit = 100,
		int	nn = 20,
		double	mu = 0.0
	)

Definition at line 438 of file cntr_equilibrium_impl.hpp.

References cntr::dos_wrapper< dos >::beta_, fourier::adft_func::dft(), cntr::dos_wrapper< dos >::mu_, cntr::herm_matrix< T >::ntau(), fourier::adft_func::sample(), cntr::herm_matrix< T >::sig(), cntr::herm_matrix< T >::size1(), cntr::dos_wrapper< dos >::tau_, and cntr::dos_wrapper< dos >::x_.

Referenced by green_equilibrium(), and green_equilibrium_mat_bethe().

{
  typedef std::complex<double> cplx;
  int ntau=G.ntau(),m,size1=G.size1();
  int sign=G.sig();
  double dtau;
  cplx res,err;
  fourier::adft_func adft;
  dos_wrapper<dos_function> dos1(dos,sign,mu);
  dos1.beta_=beta;
  dtau=beta/ntau;
  dos1.mu_=mu;
  dos1.x_=mat;
  for(m=0;m<=ntau;m++){
    dos1.tau_=m*dtau;
    adft.sample(0.0,dos.lo_,dos.hi_,dos1,nn,limit);
    adft.dft(0.0,res,err);
    element_set<T,LARGESIZE>(size1,G.matptr(m),(std::complex<T>)(res));
  }
}

Here is the call graph for this function:

Here is the caller graph for this function: