ParticleUpdatePredicates.h

Go to the documentation of this file.

#pragma once
 
#include "peano4/datamanagement/CellMarker.h"
#include "peano4/datamanagement/VertexMarker.h"
#include "toolbox/particles/MultiscaleTransitions.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 ParticleA, typename ParticleB = ParticleA>
    bool alwaysUpdateParticlePairs(
      const peano4::datamanagement::CellMarker& marker,
      const ParticleA&                           localParticle,
      const ParticleB&                           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 ParticleA, typename ParticleB = ParticleA>
    bool localParticleCanBeUpdatedInCellKernelFromAnyOtherParticle(
      const peano4::datamanagement::CellMarker& marker,
      const ParticleA&                          localParticle,
      const ParticleB&                          activeParticle
    ) {
      return (void*)(&localParticle) != (void*)(&activeParticle)
             and not localParticle.getCellHasUpdatedParticle()
             and marker.isContained(
               localParticle.getX(),
               toolbox::particles::internal::
                   relativeGrabOwnershipSpatialSortingTolerance(marker)
                 * tarch::la::NUMERICAL_ZERO_DIFFERENCE
             )
             and localParticle.getParallelState(
                 ) == ParticleA::ParallelState::Local;
    }
 
    template <typename ParticleA, typename ParticleB = ParticleA>
    bool updateFromAnyOtherParticle(
        const peano4::datamanagement::CellMarker& marker, const ParticleA& localParticle, const ParticleB& activeParticle
    ) {
        return (void*)(&localParticle) != (void*)(&activeParticle);
    }
 
    template <typename ParticleA, typename ParticleB = ParticleA>
    bool localParticleCanBeUpdatedInCellKernelFromAnyOtherParticleWithinIterationRange(
      const peano4::datamanagement::CellMarker& marker,
      const ParticleA&                          localParticle,
      const ParticleB&                          activeParticle
    ) {
      return (void*)(&localParticle) != (void*)(&activeParticle)
             and not localParticle.getCellHasUpdatedParticle()
             and marker.isContained(
               localParticle.getX(),
               toolbox::particles::internal::
                   relativeGrabOwnershipSpatialSortingTolerance(marker)
                 * tarch::la::NUMERICAL_ZERO_DIFFERENCE
             )
             and localParticle.getParallelState(
                 ) == ParticleA::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