Container for instances of ThreePointSusceptibility. More...

#include <ThreePointSusceptibilityContainer.hpp>

Inheritance diagram for Pomerol::ThreePointSusceptibilityContainer:

Collaboration diagram for Pomerol::ThreePointSusceptibilityContainer:

[legend]

Public Member Functions
template<typename... IndexTypes>
	ThreePointSusceptibilityContainer (Channel channel, 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, FreqVec2 const &freqs={}, MPI_Comm const &comm=MPI_COMM_WORLD)

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

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

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

bool	isInContainer (IndexCombination4 const &Indices) const

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

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

ElementWithPermFreq< ThreePointSusceptibility > &	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...

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< ThreePointSusceptibility >	createElement (IndexCombination4 const &Indices) const

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

Protected Attributes
Channel	channel
	Channel of \(\chi^{(3)}\). More...

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< ThreePointSusceptibility, ThreePointSusceptibilityContainer >
ParticleIndex	NumIndices
	Each of the four indices can change in the range [0; NumIndices[. More...

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

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

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

Friends
class	IndexContainer4< ThreePointSusceptibility, ThreePointSusceptibilityContainer >

Detailed Description

Container for instances of ThreePointSusceptibility.

This class stores multiple \((i, j, k, l)\)-elements of a 3-point susceptibility.

Definition at line 44 of file ThreePointSusceptibilityContainer.hpp.

Constructor & Destructor Documentation

◆ ThreePointSusceptibilityContainer()

template<typename... IndexTypes>

Pomerol::ThreePointSusceptibilityContainer::ThreePointSusceptibilityContainer	(	Channel	channel,
		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]	channel	Channel of the 3-point susceptibility.
[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 62 of file ThreePointSusceptibilityContainer.hpp.

Member Function Documentation

◆ computeAll()

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

Compute all prepared elements \(\chi^{(3)}_{ijkl}\).

Parameters

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

Precondition: prepareAll() has been called.

◆ createElement()

std::shared_ptr<ThreePointSusceptibility> Pomerol::ThreePointSusceptibilityContainer::createElement ( IndexCombination4 const & Indices ) const

protected

Create a single element \(\chi^{(3)}_{ijkl}\).

Parameters

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

◆ prepareAll()

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

Prepare a set of elements \(\chi^{(3)}_{ijkl}\).

Parameters

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