SmoothingLength.h

Go to the documentation of this file.

// This file is part of the SWIFT2 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/VertexMarker.h"
#include "swift2/kernels/ParticleUpdatePredicates.h"
#include "tarch/Assertions.h"
 
 
namespace swift2 {
  namespace kernels {
    namespace legacy {
      template <typename Particle>
      void hydro_update_smoothing_length_and_rerun_if_required(
        Particle* localParticle
      );
 
      template <typename Particle>
      bool updateSmoothingLengthAndRerunIfRequiredUpdateParticlePredicate(
        const peano4::datamanagement::VertexMarker& marker,
        const Particle&                             localParticle
      ) {
        return ::swift2::kernels::localParticleCanBeUpdatedInVertexKernel(marker, localParticle);
      }
 
      template <typename Particle>
      void updateSmoothingLengthAndRerunIfRequired(
        const peano4::datamanagement::VertexMarker& marker,
        Particle&                                   localParticle
      ) {
        assertion( updateSmoothingLengthAndRerunIfRequiredUpdateParticlePredicate(marker,localParticle) );
        hydro_update_smoothing_length_and_rerun_if_required(&localParticle);
        assertion2(
            localParticle.getSearchRadius() <= tarch::la::min(marker.h()),
            marker.toString(),
            localParticle.toString()
        );
        assertion2(
            localParticle.getSmoothingLength() <= localParticle.getSearchRadius(
            ),
            marker.toString(),
            localParticle.toString()
        );
      }
 
      template <typename Particle>
      bool resetSmoothingLengthIterationCounterUpdateParticlePredicate(
        const peano4::datamanagement::VertexMarker&  marker,
        const Particle&                              particle
      ) {
        return ::swift2::kernels::localParticleCanBeUpdatedInVertexKernel(marker, particle);
      }
 
      template <typename Particle>
      void resetSmoothingLengthIterationCounter(
        const peano4::datamanagement::VertexMarker& marker,
        Particle&                                   particle
      ) {
        assertion(resetSmoothingLengthIterationCounterUpdateParticlePredicate(marker, particle));
        particle.setSmoothingLengthConverged(false);
        particle.setSmoothingLengthIterCount(0);
      }
    } // namespace legacy
  }   // namespace kernels
} // namespace swift2
 
 
#include "SmoothingLength.cpph"