A functor to handle lennard-jones interactions between two particles (molecules) This functor uses the SIMD abstraction library Google Highway to provide architecture independent vectorization. More...

#include <LJFunctorHWY.h>

Inheritance diagram for mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >:

[legend]

Collaboration diagram for mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >:

[legend]

Public Member Functions
	LJFunctorHWY ()=delete
	Deleted default constructor.

	LJFunctorHWY (double cutoff, std::optional< std::reference_wrapper< ParticlePropertiesLibrary< double, size_t > > > particlePropertiesLibrary=std::nullopt)
	Constructor for Functor with mixing enabled/disabled.

std::string	getName () final
	Returns name of functor.

bool	isRelevantForTuning () final
	Specifies whether the functor should be considered for the auto-tuning process.

bool	allowsNewton3 () final
	Specifies whether the functor is capable of Newton3-like functors.

bool	allowsNonNewton3 () final
	Specifies whether the functor is capable of non-Newton3-like functors.

bool	isVecPatternAllowed (const VectorizationPattern vecPattern) final
	Specifies whether the functor is capable of using the specified Vectorization Pattern in the SoA functor.

void	AoSFunctor (Particle_T &i, Particle_T &j, bool newton3) final
	PairwiseFunctor for arrays of structures (AoS).

void	SoAFunctorSingle (autopas::SoAView< SoAArraysType > soa, const bool newton3) final
	PairwiseFunctor for structure of arrays (SoA)

void	SoAFunctorPair (autopas::SoAView< SoAArraysType > soa1, autopas::SoAView< SoAArraysType > soa2, bool newton3) final
	PairwiseFunctor for structure of arrays (SoA)

void	SoAFunctorVerlet (autopas::SoAView< SoAArraysType > soa, const size_t indexFirst, const std::vector< size_t, autopas::AlignedAllocator< size_t > > &neighborList, bool newton3) final
	PairwiseFunctor for structure of arrays (SoA) for neighbor lists.

void	initTraversal () final
	Reset the global values.

void	endTraversal (const bool newton3) final
	Accumulates global values, e.g.

double	getPotentialEnergy () const
	Get the potential Energy.

double	getVirial () const
	Get the virial.

void	setParticleProperties (const double epsilon24, const double sigmaSquare)
	Sets the particle properties constants for this functor.

void	setVecPattern (const VectorizationPattern vecPattern) final
	Setter for the vectorization pattern to be used.

Public Member Functions inherited from autopas::PairwiseFunctor< Particle_T, CRTP_T >
	PairwiseFunctor (double cutoff)
	Constructor.

virtual void	AoSFunctor (Particle_T &i, Particle_T &j, bool newton3)
	PairwiseFunctor for arrays of structures (AoS).

virtual void	SoAFunctorSingle (SoAView< SoAArraysType > soa, bool newton3)
	PairwiseFunctor for structure of arrays (SoA)

virtual void	SoAFunctorVerlet (SoAView< SoAArraysType > soa, const size_t indexFirst, const std::vector< size_t, AlignedAllocator< size_t > > &neighborList, bool newton3)
	PairwiseFunctor for structure of arrays (SoA) for neighbor lists.

virtual void	SoAFunctorPair (SoAView< SoAArraysType > soa1, SoAView< SoAArraysType > soa2, bool newton3)
	PairwiseFunctor for structure of arrays (SoA)

bool	isVecPatternAllowed (const VectorizationPatternOption::Value vecPattern) override
	Specifies whether the functor is capable of using the specified Vectorization Pattern in the SoA functor.

Public Member Functions inherited from autopas::Functor< Particle_T, CRTP_T >
	Functor (double cutoff)
	Constructor.

virtual void	initTraversal ()
	This function is called at the start of each traversal.

virtual void	endTraversal (bool newton3)
	This function is called at the end of each traversal.

template<class ParticleCell >
void	SoALoader (ParticleCell &cell, SoA< SoAArraysType > &soa, size_t offset, bool skipSoAResize)
	Copies the AoS data of the given cell in the given soa.

template<typename ParticleCell >
void	SoAExtractor (ParticleCell &cell, SoA< SoAArraysType > &soa, size_t offset)
	Copies the data stored in the soa back into the cell.

virtual bool	allowsNewton3 ()=0
	Specifies whether the functor is capable of Newton3-like functors.

virtual bool	allowsNonNewton3 ()=0
	Specifies whether the functor is capable of non-Newton3-like functors.

virtual bool	isVecPatternAllowed (const VectorizationPatternOption::Value vecPattern)=0
	Specifies whether the functor is capable of using the specified Vectorization Pattern in the SoA functor.

virtual bool	isRelevantForTuning ()=0
	Specifies whether the functor should be considered for the auto-tuning process.

virtual std::string	getName ()=0
	Returns name of functor.

double	getCutoff () const
	Getter for the functor's cutoff.

virtual void	setVecPattern (const VectorizationPatternOption::Value vecPattern)
	Setter for the vectorization pattern to be used.

virtual size_t	getNumFLOPs () const
	Get the number of FLOPs.

virtual double	getHitRate () const
	Get the hit rate.

Static Public Member Functions
static constexpr auto	getNeededAttr ()
	Get attributes needed for computation.

static constexpr auto	getNeededAttr (std::false_type)
	Get attributes needed for computation without N3 optimization.

static constexpr auto	getComputedAttr ()
	Get attributes computed by this functor.

static constexpr bool	getMixing ()

Static Public Member Functions inherited from autopas::Functor< Particle_T, CRTP_T >
static constexpr std::array< typename Particle_T::AttributeNames, 0 >	getNeededAttr ()
	Get attributes needed for computation.

static constexpr std::array< typename Particle_T::AttributeNames, 0 >	getNeededAttr (std::false_type)
	Get attributes needed for computation without N3 optimization.

static constexpr std::array< typename Particle_T::AttributeNames, 0 >	getComputedAttr ()
	Get attributes computed by this functor.

Additional Inherited Members
Public Types inherited from autopas::PairwiseFunctor< Particle_T, CRTP_T >
using	SoAArraysType = typename Particle_T::SoAArraysType
	Structure of the SoAs defined by the particle.

Public Types inherited from autopas::Functor< Particle_T, CRTP_T >
using	SoAArraysType = typename Particle_T::SoAArraysType
	Structure of the SoAs defined by the particle.

using	Functor_T = CRTP_T
	Make the Implementation type template publicly available.

Detailed Description

template<class Particle_T, bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>
class mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >

A functor to handle lennard-jones interactions between two particles (molecules) This functor uses the SIMD abstraction library Google Highway to provide architecture independent vectorization.

Template Parameters

Particle_T	The type of particle
applyShift	Switch for the lj potential to be truncated shifted
useMixing	Switch for the functor to be used with multiple particle types. If set to false, _epsilon and _sigma need to be set and the constructor with PPL can be omitted.
useNewton3	Switch for the functor to support newton3 on, off or both. See FunctorN3Modes for possible values.
calculateGlobals	Defines whether the global values are to be calculated (energy, virial).
relevantForTuning	Whether the auto-tuner should consider this functor.
countFLOPs	counts FLOPs and hitrate. Not implemented for this functor. Please use the AutoVec functor.

Constructor & Destructor Documentation

◆ LJFunctorHWY()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::LJFunctorHWY	(	double	cutoff,
		std::optional< std::reference_wrapper< ParticlePropertiesLibrary< double, size_t > > >	particlePropertiesLibrary = `std::nullopt`
	)

inlineexplicit

Constructor for Functor with mixing enabled/disabled.

When using this functor without particlePropertiesLibrary it is necessary to call setParticleProperties() to set internal constants.

Parameters

cutoff
particlePropertiesLibrary

Member Function Documentation

◆ allowsNewton3()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

bool mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::allowsNewton3 ( )

inlinefinalvirtual

Specifies whether the functor is capable of Newton3-like functors.

If the functor provides an interface to soa or aos functions that utilize Newton's third law of motion (actio = reactio) to reduce the computational complexity this function should return true. If this is not the case this function should return false.

Returns: true if and only if this functor provides an interface to Newton3-like functions.

Implements autopas::Functor< Particle_T, CRTP_T >.

◆ allowsNonNewton3()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

bool mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::allowsNonNewton3 ( )

inlinefinalvirtual

Specifies whether the functor is capable of non-Newton3-like functors.

If the functor provides an interface to soa or aos functions that do not utilize Newton's third law of motion (actio = reactio) this function should return true. If this is not the case this function should return false.

Returns: true if and only if this functor provides an interface to functions that do not utilize Newton3.

Implements autopas::Functor< Particle_T, CRTP_T >.

◆ AoSFunctor()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

void mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::AoSFunctor	(	Particle_T &	i,
		Particle_T &	j,
		bool	newton3
	)

inlinefinalvirtual

PairwiseFunctor for arrays of structures (AoS).

This functor should calculate the forces or any other pair-wise interaction between two particles. This should include a cutoff check if needed!

Parameters

i	Particle i
j	Particle j
newton3	defines whether or whether not to use newton 3

Reimplemented from autopas::PairwiseFunctor< Particle_T, CRTP_T >.

◆ endTraversal()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

void mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::endTraversal ( const bool newton3 )

inlinefinalvirtual

Accumulates global values, e.g.

upot and virial.

Parameters

newton3

Reimplemented from autopas::Functor< Particle_T, CRTP_T >.

◆ getComputedAttr()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

static constexpr auto mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::getComputedAttr ( )

inlinestaticconstexpr

Get attributes computed by this functor.

Returns: Attributes computed by this functor.

Todo:: C++20: make this function virtual

◆ getMixing()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

static constexpr bool mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::getMixing ( )

inlinestaticconstexpr

Returns: useMixing

◆ getName()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

std::string mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::getName ( )

inlinefinalvirtual

Returns name of functor.

Intended for use with the iteration logger, to differentiate between calls to computeInteractions using different functors in the logs.

Returns: name of functor.

Implements autopas::Functor< Particle_T, CRTP_T >.

◆ getNeededAttr() [1/2]

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

static constexpr auto mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::getNeededAttr ( )

inlinestaticconstexpr

Get attributes needed for computation.

Returns: Attributes needed for computation.

Todo:: C++20: make this function virtual

◆ getNeededAttr() [2/2]

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

static constexpr auto mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::getNeededAttr ( std::false_type )

inlinestaticconstexpr

Get attributes needed for computation without N3 optimization.

Returns: Attributes needed for computation.

Todo:: C++20: make this function virtual

◆ getPotentialEnergy()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

double mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::getPotentialEnergy ( ) const

inline

Get the potential Energy.

Returns: the potential Energy

◆ getVirial()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

double mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::getVirial ( ) const

inline

Get the virial.

Returns: the virial

◆ initTraversal()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

void mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::initTraversal ( )

inlinefinalvirtual

Reset the global values.

Will set the global values to zero to prepare for the next iteration.

Reimplemented from autopas::Functor< Particle_T, CRTP_T >.

◆ isRelevantForTuning()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

bool mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::isRelevantForTuning ( )

inlinefinalvirtual

Specifies whether the functor should be considered for the auto-tuning process.

Returns: true if and only if this functor is relevant for auto-tuning.

Implements autopas::Functor< Particle_T, CRTP_T >.

◆ isVecPatternAllowed()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

bool mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::isVecPatternAllowed ( const VectorizationPattern vecPattern )

inlinefinal

Specifies whether the functor is capable of using the specified Vectorization Pattern in the SoA functor.

Thi functor can handle p1xVec, p2xVecDiv2, pVecDiv2x2, pVecx1

Parameters

vecPattern

Returns: whether the functor is capable of using the specified Vectorization Pattern

◆ setParticleProperties()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

void mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::setParticleProperties	(	const double	epsilon24,
		const double	sigmaSquare
	)

inline

Sets the particle properties constants for this functor.

This is only necessary if no particlePropertiesLibrary is used.

Parameters

epsilon24
sigmaSquare

◆ setVecPattern()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

void mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::setVecPattern ( const VectorizationPattern vecPattern )

inlinefinal

Setter for the vectorization pattern to be used.

Parameters

vecPattern

◆ SoAFunctorPair()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

void mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::SoAFunctorPair	(	autopas::SoAView< SoAArraysType >	soa1,
		autopas::SoAView< SoAArraysType >	soa2,
		bool	newton3
	)

inlinefinalvirtual

PairwiseFunctor for structure of arrays (SoA)

This functor should calculate the forces or any other pair-wise interaction between all particles of soa1 and soa2. This should include a cutoff check if needed!

Parameters

soa1	First structure of arrays.
soa2	Second structure of arrays.
newton3	defines whether or whether not to use newton 3

Reimplemented from autopas::PairwiseFunctor< Particle_T, CRTP_T >.

◆ SoAFunctorSingle()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

void mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::SoAFunctorSingle	(	autopas::SoAView< SoAArraysType >	soa,
		const bool	newton3
	)

inlinefinalvirtual

PairwiseFunctor for structure of arrays (SoA)

This functor should calculate the forces or any other pair-wise interaction between all particles in an soa. This should include a cutoff check if needed!

Parameters

soa	Structure of arrays
newton3	defines whether or whether not to use newton 3 This functor will always do a newton3 like traversal of the soa.

Reimplemented from autopas::PairwiseFunctor< Particle_T, CRTP_T >.

◆ SoAFunctorVerlet()

template<class Particle_T , bool applyShift = false, bool useMixing = false, autopas::FunctorN3Modes useNewton3 = autopas::FunctorN3Modes::Both, bool calculateGlobals = false, bool countFLOPs = false, bool relevantForTuning = true>

void mdLib::LJFunctorHWY< Particle_T, applyShift, useMixing, useNewton3, calculateGlobals, countFLOPs, relevantForTuning >::SoAFunctorVerlet	(	autopas::SoAView< SoAArraysType >	soa,
		const size_t	indexFirst,
		const std::vector< size_t, autopas::AlignedAllocator< size_t > > &	neighborList,
		bool	newton3
	)

inlinefinalvirtual

PairwiseFunctor for structure of arrays (SoA) for neighbor lists.

This functor should calculate the forces or any other pair-wise interaction between the particle in the SoA with index indexFirst and all particles with indices in the neighborList. This should include a cutoff check if needed!

Parameters

soa	Structure of arrays
indexFirst	The index of the first particle for each interaction
neighborList	The list of neighbors
newton3	defines whether or whether not to use newton 3

Note: If you want to parallelize this by openmp, please ensure that there are no dependencies, i.e. introduce colors and specify iFrom and iTo accordingly.

Reimplemented from autopas::PairwiseFunctor< Particle_T, CRTP_T >.

The documentation for this class was generated from the following file:

/github/workspace/applicationLibrary/molecularDynamics/molecularDynamicsLibrary/LJFunctorHWY.h

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ LJFunctorHWY()

Member Function Documentation

◆ allowsNewton3()

◆ allowsNonNewton3()

◆ AoSFunctor()

◆ endTraversal()

◆ getComputedAttr()

◆ getMixing()

◆ getName()

◆ getNeededAttr() [1/2]

◆ getNeededAttr() [2/2]

◆ getPotentialEnergy()

◆ getVirial()

◆ initTraversal()

◆ isRelevantForTuning()

◆ isVecPatternAllowed()

◆ setParticleProperties()

◆ setVecPattern()

◆ SoAFunctorPair()

◆ SoAFunctorSingle()

◆ SoAFunctorVerlet()