Functor that returns forces acting on a particle surface, returns data in output for every particle in a row(described are return values for the first particle). More...

#include <latticeMomentumExchangeForce.h>

Inheritance diagram for olb::SuperLatticeParticleForce< T, DESCRIPTOR, PARTICLETYPE, BLOCKFUNCTOR >:

Collaboration diagram for olb::SuperLatticeParticleForce< T, DESCRIPTOR, PARTICLETYPE, BLOCKFUNCTOR >:

Public Member Functions
	SuperLatticeParticleForce (SuperLattice< T, DESCRIPTOR > &sLattice, const SuperGeometry< T, DESCRIPTOR::d > &superGeometry, particles::ParticleSystem< T, PARTICLETYPE > &particleSystem, const UnitConverter< T, DESCRIPTOR > &converter, Vector< bool, DESCRIPTOR::d > periodic=Vector< bool, DESCRIPTOR::d >(false), std::size_t iP0=0, const std::unordered_set< int > &ignoredMaterials=std::unordered_set< int >{}, const F f=[](auto &, const auto &, const auto &, const auto &){})

bool	operator() (T output[], const int input[]) override

Detailed Description

template<typename T, typename DESCRIPTOR, typename PARTICLETYPE, typename BLOCKFUNCTOR = BlockLatticeMomentumExchangeForce<T,DESCRIPTOR,PARTICLETYPE>>
class olb::SuperLatticeParticleForce< T, DESCRIPTOR, PARTICLETYPE, BLOCKFUNCTOR >

Functor that returns forces acting on a particle surface, returns data in output for every particle in a row(described are return values for the first particle).

Returns: output[0]-output[2] translational force - physical units; output[3]-output[5] torque - physical units; output[7] number of voxels

Definition at line 88 of file latticeMomentumExchangeForce.h.

Constructor & Destructor Documentation

◆ SuperLatticeParticleForce()

template<typename T , typename DESCRIPTOR , typename PARTICLETYPE , typename BLOCKFUNCTOR >

olb::SuperLatticeParticleForce< T, DESCRIPTOR, PARTICLETYPE, BLOCKFUNCTOR >::SuperLatticeParticleForce	(	SuperLattice< T, DESCRIPTOR > &	sLattice,
		const SuperGeometry< T, DESCRIPTOR::d > &	superGeometry,
		particles::ParticleSystem< T, PARTICLETYPE > &	particleSystem,
		const UnitConverter< T, DESCRIPTOR > &	converter,
		Vector< bool, DESCRIPTOR::d >	periodic = Vector<bool,DESCRIPTOR::d> (false),
		std::size_t	iP0 = 0,
		const std::unordered_set< int > &	ignoredMaterials = std::unordered_set<int>{},
		const F	f = [](auto&, const auto&, const auto&, const auto&){} )

Definition at line 35 of file latticeMomentumExchangeForce.hh.

  : SuperLatticePhysF<T,DESCRIPTOR>(sLattice,converter,
                                    (DESCRIPTOR::d+utilities::dimensions::convert<DESCRIPTOR::d>::rotation+1)*(particleSystem.size()-iP0))
{
  constexpr unsigned D = DESCRIPTOR::d;
 
  this->getName() = "physParticleForce";
  int maxC = this->_sLattice.getLoadBalancer().size();
  this->_blockF.reserve(maxC);
  const PhysR<T,D> min = particles::communication::getCuboidMin<T,D>(
      superGeometry.getCuboidGeometry());
  const PhysR<T,D> max = particles::communication::getCuboidMax<T,D>(
      superGeometry.getCuboidGeometry(), min);
 
  for (int iC = 0; iC < maxC; ++iC) {
    this->_blockF.emplace_back( new BLOCKFUNCTOR(
                                  this->_sLattice.getBlock(iC),
                                  superGeometry.getBlockGeometry(iC),
                                  particleSystem, converter, min, max,
                                  periodic, iP0, ignoredMaterials, f));
  }
}

References olb::SuperGeometry< T, D >::getBlockGeometry(), and olb::SuperStructure< T, D >::getCuboidGeometry().

Here is the call graph for this function:

Member Function Documentation

◆ operator()()

template<typename T , typename DESCRIPTOR , typename PARTICLETYPE , typename BLOCKFUNCTOR >

bool olb::SuperLatticeParticleForce< T, DESCRIPTOR, PARTICLETYPE, BLOCKFUNCTOR >::operator()	(	T	output[],
		const int	input[] )

override

Definition at line 65 of file latticeMomentumExchangeForce.hh.

{
  for (int iS=0; iS<this->getTargetDim(); ++iS) {
    output[iS] = 0.;
  }
  for (int iC = 0; iC < this->_sLattice.getLoadBalancer().size(); ++iC) {
    this->getBlockF(iC)(output, &input[1]);
  }
 
#ifdef PARALLEL_MODE_MPI
  for (int iS = 0; iS < this->getTargetDim(); ++iS) {
    singleton::mpi().reduceAndBcast(output[iS], MPI_SUM);
  }
#endif
  return true;
 
}

References olb::singleton::mpi(), and olb::singleton::MpiManager::reduceAndBcast().