ParticleUpdatePredicates.h

Go to the documentation of this file.

#pragma once
 
#include "toolbox/particles/MultiscaleTransitions.h"
#include "peano4/datamanagement/CellMarker.h"
#include "peano4/datamanagement/VertexMarker.h"
 
namespace swift2 {
  namespace kernels {
    template <typename Particle>
    using UpdateParticlePairWithinCellPredicate = std::function<bool(const peano4::datamanagement::CellMarker& marker, const Particle& localParticle, const Particle& activeParticle)>;
 
    template <typename Particle>
    using UpdateParticleAssignedToCellPredicate = std::function<bool(const peano4::datamanagement::CellMarker& marker, const bool& localParticleIsContainedInCell, const Particle& localParticle)>;
 
    template <typename Particle>
    using UpdateParticleAssignedToVertexPredicate = std::function<bool(const peano4::datamanagement::VertexMarker& marker, const Particle& localParticle)>;
 
 
    template <typename Particle>
    bool alwaysUpdateInVertexKernel(const peano4::datamanagement::VertexMarker&  marker, const Particle& localParticle) {
      return true;
    }
 
 
    template <typename Particle>
    bool alwaysUpdateInCellKernel(const peano4::datamanagement::CellMarker&  marker, const bool& localParticleIsContainedInCell, const Particle& localParticle) {
      return true;
    }
 
 
    template <typename Particle>
    bool alwaysUpdateParticlePairs(const peano4::datamanagement::CellMarker&  marker, const Particle& localParticle, const Particle& activeParticle) {
      return true;
    }
 
 
    template <typename Particle>
    bool particleIsLocal(const peano4::datamanagement::VertexMarker& marker, const Particle& localParticle) {
      return localParticle.getParallelState() == Particle::ParallelState::Local;
    }
 
 
    template <typename Particle>
    bool localParticleCanBeUpdatedAndMovedInVertexKernel(const peano4::datamanagement::VertexMarker& marker, const Particle& localParticle) {
      return (
        (localParticle.getMoveState() == Particle::MoveState::NotMovedYet)
        & particleIsLocal(marker, localParticle)
      );
    }
 
    template <typename Particle>
    bool localParticleCanBeUpdatedInCellKernelFromAnyOtherParticle(
      const peano4::datamanagement::CellMarker& marker, const Particle& localParticle, const Particle& activeParticle
    ) {
      return &localParticle != &activeParticle
         and not localParticle.getCellHasUpdatedParticle()
         and marker.isContained(
           localParticle.getX(),
           toolbox::particles::internal::relativeGrabOwnershipSpatialSortingTolerance(marker) * tarch::la::NUMERICAL_ZERO_DIFFERENCE
         )
         and localParticle.getParallelState() == Particle::ParallelState::Local;
    }
 
 
    template <typename Particle>
    bool localParticleCanBeUpdatedInCellKernelFromAnyOtherParticleWithinIterationRange(
      const peano4::datamanagement::CellMarker& marker, const Particle& localParticle, const Particle& activeParticle
    ) {
      return &localParticle != &activeParticle
         and not localParticle.getCellHasUpdatedParticle()
         and marker.isContained(
           localParticle.getX(),
           toolbox::particles::internal::relativeGrabOwnershipSpatialSortingTolerance(marker) * tarch::la::NUMERICAL_ZERO_DIFFERENCE
         )
         and localParticle.getParallelState() == Particle::ParallelState::Local;
    }
 
 
    template <typename Particle>
    bool localParticleCanBeUpdatedInCellKernel(
      const peano4::datamanagement::CellMarker& marker, const bool& localParticleIsContainedInCell, const Particle& localParticle
    ) {
      assertionEquals2(
        localParticleIsContainedInCell,
        marker.isContained(
          localParticle.getX(),
          toolbox::particles::internal::relativeGrabOwnershipSpatialSortingTolerance(marker) * tarch::la::NUMERICAL_ZERO_DIFFERENCE
        ),
        marker.toString(), localParticle.toString()
      );
 
      return not localParticle.getCellHasUpdatedParticle()
         and localParticleIsContainedInCell
         and localParticle.getParallelState() == Particle::ParallelState::Local;
    }
 
 
    template <typename Particle>
    bool localParticleCanBeUpdatedInVertexKernel(const peano4::datamanagement::VertexMarker& marker, const Particle& localParticle) {
      return particleIsLocal(marker, localParticle);
    }
  } // namespace kernels
} // namespace swift2