CD_ParticleOps.H

Go to the documentation of this file.

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
#ifndef CD_PARTICLEOPS_H
#define CD_PARTICLEOPS_H
 
// Chombo includes
#include <RefCountedPtr.H>
#include <BaseIF.H>
#include <RealVect.H>
 
// Our includes
#include <CD_EBAMRData.H>
#include <CD_ParticleMap.H>
#include <CD_ParticleContainer.H>
#include <CD_NamespaceHeader.H>
 
class ParticleOps
{
public:
  static inline IntVect
  getParticleCellIndex(const RealVect& a_particlePosition, const RealVect& a_probLo, const Real& a_dx) noexcept;
 
  static inline IntVect
  getParticleCellIndex(const RealVect& a_particlePosition, const RealVect& a_probLo, const RealVect& a_dx) noexcept;
 
  template <typename P, typename Traits>
  static inline void
  getPhysicalParticlesPerCell(EBAMRCellData& a_ppc, const ParticleContainer<P, Traits>& a_src) noexcept;
 
  template <typename P, typename Traits>
  static inline void
  getComputationalParticlesPerCell(EBAMRCellData& a_ppc, const ParticleContainer<P, Traits>& a_src) noexcept;
 
  static inline IntVect
  getParticleGridCell(const RealVect& a_particlePosition, const RealVect& a_probLo, const RealVect& a_dx) noexcept;
 
  static inline bool
  domainIntersection(const RealVect& a_oldPos,
                     const RealVect& a_newPos,
                     const RealVect& a_probLo,
                     const RealVect& a_probHi,
                     Real&           a_s);
 
  static inline bool
  ebIntersectionBisect(const RefCountedPtr<BaseIF>& a_impFunc,
                       const RealVect&              a_oldPos,
                       const RealVect&              a_newPos,
                       const Real&                  a_bisectStep,
                       Real&                        a_s);
 
  static inline bool
  ebIntersectionRaycast(const RefCountedPtr<BaseIF>& a_impFunc,
                        const RealVect&              a_oldPos,
                        const RealVect&              a_newPos,
                        const Real&                  a_tolerance,
                        Real&                        a_s);
 
  template <typename P, typename Traits>
  static inline void
  copyDestructive(ParticleContainer<P, Traits>& a_dst, ParticleContainer<P, Traits>& a_src) noexcept;
 
  template <typename P, typename Traits>
  static inline Real
  sum(const ParticleContainer<P, Traits>& a_particles) noexcept;
 
  template <typename P, typename Traits>
  static inline void
  setData(ParticleContainer<P, Traits>&                                    a_particles,
          const std::function<void(ParticleSoA<P, Traits>&, std::size_t)>& a_functor) noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#include <CD_ParticleOpsImplem.H>
 
#endif