vie2_mat()

template<typename T >

void cntr::vie2_mat	(	herm_matrix< T > &	G,
		herm_matrix< T > &	F,
		herm_matrix< T > &	Fcc,
		herm_matrix< T > &	Q,
		T	beta,
		const int	SolveOrder,
		const int	method
	)

VIE solver \((1+F)*G=Q\) for a Green's function \(G\) on the Matsubara axis

Purpose

One solves the linear equation \((1+F)*G=Q\) for a hermitian matrix \(G(t, t^\prime)\) on a Matsubara axis, for given: input kernel \(F(t, t^\prime)\), its hermitian conjugate \(F^\ddagger(t, t^\prime)\), and the source term \(Q(t, t^\prime)\). There are 3 possible methods for solution: Fourier, steep, and fixpoint. Fixpoint method is choosen by default.

Parameters

&G	[herm_matrix] solution
&F	[herm_matrix] green's function on left-hand side
&Fcc	[herm_matrix] Complex conjugate of F
&Q	[herm_matrix] green's function on right-hand side
beta	[double] inverse temperature
SolveOrder	[int] integrator order
method	[const] Solution method on the Matsubara axis with 0: Fourier, 1: steep, 2: fixpoint

Definition at line 1283 of file cntr_vie2_impl.hpp.

                                                                       {

  const int fourier_order=3;
  int maxiter;
  T tol = 1.0e-12;

  switch(method) {
  case CNTR_MAT_FOURIER:
    vie2_mat_fourier(G, F, Fcc, Q, beta, fourier_order);
    break;
  case CNTR_MAT_CG:
    maxiter = 40;
    vie2_mat_steep(G, F, Fcc, Q, beta, integration::I<T>(SolveOrder), maxiter, tol);
    break;
  default:
    maxiter = 6;
    vie2_mat_fixpoint(G, F, Fcc, Q, beta, integration::I<T>(SolveOrder), maxiter);
    break;
  }

}