particles.h

Go to the documentation of this file.

// This file is part of the Peano project. For conditions of distribution and
// use, please see the copyright notice at www.peano-framework.org
#pragma once
 
 
#include "peano4/datamanagement/CellMarker.h"
#include "peano4/datamanagement/VertexEnumerator.h"
#include "tarch/la/Vector.h"
 
 
namespace toolbox {
  namespace particles {
    constexpr double SpatialDuplicateTolerance = 1e-6;
 
    template <typename Iterator>
    Iterator particleIsDuplicate(
      const typename Iterator::value_type& particle,
      Iterator                             begin,
      const Iterator&                      end
    );
 
    template <typename Particle>
    bool particleWillStayWithinComputationalDomain(
      const Particle&                              particle,
      const tarch::la::Vector<Dimensions, double>& domainOffset,
      const tarch::la::Vector<Dimensions, double>& domainWidth
    );
 
 
    template <typename Particle, typename ParticleSet>
    void insertParticlesIntoCell(
      const peano4::datamanagement::CellMarker&              marker,
      const std::vector<Particle*>&                          newParticles,
      peano4::datamanagement::VertexEnumerator<ParticleSet>& fineGridVertices,
      int                                                    spacetreeId
    );
 
 
    template <typename Particle, typename ParticleSet>
    void insertParticleIntoCell(
      const peano4::datamanagement::CellMarker&              marker,
      Particle*                                              newParticles,
      peano4::datamanagement::VertexEnumerator<ParticleSet>& fineGridVertices,
      int                                                    spacetreeId
    );
 
    void updateContainedInAdjacentCellProperty(
      const peano4::datamanagement::CellMarker& marker,
      const std::bitset<Dimensions>&            assignedVertex,
      auto&                                     particle
    );
 
    void updateContainedInAdjacentCellProperty(
      const peano4::datamanagement::VertexMarker& marker,
      auto&                                       particle
    );
 
    void updateContainedInAdjacentCellProperty(
      const tarch::la::Vector<Dimensions, double>& x,
      const tarch::la::Vector<Dimensions, double>& h,
      auto&                                        particle
    );
 
    bool particleAssignedToVertexWillBeLocal(
      const tarch::la::Vector<Dimensions, double>& x,
      const peano4::datamanagement::VertexMarker&  marker
    );
 
    tarch::la::Vector<Dimensions, double> mirrorParticleAlongPeriodicDomains(
      const tarch::la::Vector<Dimensions, double>& x,
      const peano4::datamanagement::VertexMarker&  marker,
      const tarch::la::Vector<Dimensions, double>  domainOffset,
      const tarch::la::Vector<Dimensions, double>  domainSize,
      const std::bitset<Dimensions>                periodicBC
    );
 
    tarch::la::Vector<Dimensions, double> applyPeriodicBoundaryConditions(
      const tarch::la::Vector<Dimensions, double>& x,
      const tarch::la::Vector<Dimensions, double>  domainOffset,
      const tarch::la::Vector<Dimensions, double>  domainSize,
      const std::bitset<Dimensions>                periodicBC
    );
  } // namespace particles
} // namespace toolbox
 
 
#include "particles.cpph"