Container for instances of TwoParticleGF. More...

#include <TwoParticleGFContainer.hpp>

Inheritance diagram for Pomerol::TwoParticleGFContainer:

Collaboration diagram for Pomerol::TwoParticleGFContainer:

[legend]

Public Member Functions
template<typename... IndexTypes>
	TwoParticleGFContainer (IndexClassification< IndexTypes... > const &IndexInfo, StatesClassification const &S, Hamiltonian const &H, DensityMatrix const &DM, FieldOperatorContainer const &Ops)

void	prepareAll (std::set< IndexCombination4 > const &Indices={})

std::map< IndexCombination4, std::vector< ComplexType > >	computeAll (bool clearTerms=false, FreqVec3 const &freqs={}, MPI_Comm const &comm=MPI_COMM_WORLD, bool split=true)

Public Member Functions inherited from Pomerol::IndexContainer4< TwoParticleGF, TwoParticleGFContainer >
	IndexContainer4 (TwoParticleGFContainer const &Source, IndexClassification< IndexTypes... > const &IndexInfo)

void	fill (std::set< IndexCombination4 > Indices=std::set< IndexCombination4 >())

ElementWithPermFreq< TwoParticleGF > &	create (IndexCombination4 const &Indices)

bool	isInContainer (IndexCombination4 const &Indices) const

bool	isInContainer (ParticleIndex Index1, ParticleIndex Index2, ParticleIndex Index3, ParticleIndex Index4) const

ElementWithPermFreq< TwoParticleGF > &	operator() (IndexCombination4 const &Indices)

ElementWithPermFreq< TwoParticleGF > &	operator() (ParticleIndex Index1, ParticleIndex Index2, ParticleIndex Index3, ParticleIndex Index4)

Public Member Functions inherited from Pomerol::Thermal
	Thermal (RealType beta)

Data Fields
RealType	ReduceResonanceTolerance = 1e-8
	A difference in energies with magnitude below this value is treated as zero. More...

RealType	CoefficientTolerance = 1e-16
	Minimal magnitude of the coefficient of a term for it to be taken into account. More...

RealType	MultiTermCoefficientTolerance = 1e-5

Data Fields inherited from Pomerol::Thermal
const RealType	beta
	Inverse temperature \(\beta\). More...

const ComplexType	MatsubaraSpacing
	Spacing between (imaginary) Matsubara frequencies, \(i\pi/\beta\). More...

Protected Member Functions
std::shared_ptr< TwoParticleGF >	createElement (IndexCombination4 const &Indices) const

Protected Member Functions inherited from Pomerol::IndexContainer4< TwoParticleGF, TwoParticleGFContainer >
std::set< IndexCombination4 >	enumerateIndices () const

Protected Attributes
StatesClassification const &	S
	Information about invariant subspaces of the Hamiltonian. More...

Hamiltonian const &	H
	The Hamiltonian. More...

DensityMatrix const &	DM
	Many-body density matrix \(\hat\rho\). More...

FieldOperatorContainer const &	Operators
	A set of creation/annihilation operators \(c^\dagger\)/ \(c\). More...

Protected Attributes inherited from Pomerol::IndexContainer4< TwoParticleGF, TwoParticleGFContainer >
ParticleIndex	NumIndices
	Each of the four indices can change in the range [0; NumIndices[. More...

TwoParticleGFContainer const &	Source
	Stored elements are created by calling Source.createElement(Indices). More...

std::map< IndexCombination4, ElementWithPermFreq< TwoParticleGF > >	ElementsMap
	Sparse storage for the decorated elements. More...

std::map< IndexCombination4, std::shared_ptr< TwoParticleGF > >	NonTrivialElements
	Sparse storage for the plain (non-decorated) elements. More...

Friends
class	IndexContainer4< TwoParticleGF, TwoParticleGFContainer >

Detailed Description

Container for instances of TwoParticleGF.

This class stores multiple \((i,j,k,l)\)-elements of a fermionic two-particle Matsubara Green's function

\[ \chi_{ijkl}(\omega_{n_1},\omega_{n_2};\omega_{n_3},\omega_{n_1}+\omega_{n_2}-\omega_{n_3}) = \int_0^\beta Tr[\mathcal{T}_\tau \hat\rho c_i(\tau_1)c_j(\tau_2)c^\dagger_k(\tau_3)c^\dagger_l(0)] e^{i\omega_{n_1}\tau_1+i\omega_{n_2}\tau_2-i\omega_{n_3}\tau_3} d\tau_1 d\tau_2 d\tau_3. \]

Definition at line 51 of file TwoParticleGFContainer.hpp.

Constructor & Destructor Documentation

◆ TwoParticleGFContainer()

template<typename... IndexTypes>

Pomerol::TwoParticleGFContainer::TwoParticleGFContainer	(	IndexClassification< IndexTypes... > const &	IndexInfo,
		StatesClassification const &	S,
		Hamiltonian const &	H,
		DensityMatrix const &	DM,
		FieldOperatorContainer const &	Ops
	)

inline

Constructor.

Template Parameters

IndexTypes Types of indices carried by the creation and annihilation operators.

Parameters

[in]	IndexInfo	Map for fermionic operator index tuples.
[in]	S	Information about invariant subspaces of the Hamiltonian.
[in]	H	The Hamiltonian.
[in]	DM	Many-body density matrix \(\hat\rho\).
[in]	Ops	A set of creation/annihilation operators \(c^\dagger\)/ \(c\).

Definition at line 69 of file TwoParticleGFContainer.hpp.

Member Function Documentation

◆ computeAll()

std::map<IndexCombination4, std::vector<ComplexType> > Pomerol::TwoParticleGFContainer::computeAll	(	bool	clearTerms = `false`,
		FreqVec3 const &	freqs = `{}`,
		MPI_Comm const &	comm = `MPI_COMM_WORLD`,
		bool	split = `true`
	)

Compute all prepared elements \(\chi_{ijkl}\).

Parameters

[in]	clearTerms	If true, computed TwoParticleGFPart's of all elements will be destroyed immediately after filling the precomputed value cache.
[in]	freqs	List of frequency triplets \((\omega_{n_1},\omega_{n_2},\omega_{n_3})\) for value pre-computation.
[in]	comm	MPI communicator used to parallelize the computation.
[in]	split	Enable MPI parallelization.

Precondition: prepareAll() has been called.

◆ createElement()

std::shared_ptr<TwoParticleGF> Pomerol::TwoParticleGFContainer::createElement ( IndexCombination4 const & Indices ) const

protected

Create a single element \(\chi_{ijkl}\).

Parameters

[in] Indices Index combination \((i,j,k,l)\).

◆ prepareAll()

void Pomerol::TwoParticleGFContainer::prepareAll ( std::set< IndexCombination4 > const & Indices = {} )

Prepare a set of elements \(\chi_{ijkl}\).

Parameters

[in] Indices Set of index combinations of the elements \(\chi_{ijkl}\) to be prepared. An empty set results in creation of elements for all possible index combinations \((i,j,k,l)\).