cell_batch_particle_cache.hpp

Go to the documentation of this file.

#pragma once
 
#include <meltpooldg/particles/obstacle_data_structure.hpp>
#include <meltpooldg/particles/particle_accessor.hpp>
#include <meltpooldg/utilities/matrix_free_util.hpp>
 
#include <boost/container/small_vector.hpp>
 
#include <vector>
 
namespace MeltPoolDG
{
  template <int dim, typename number, typename ObstacleType>
  class MatrixFreeCellBatchParticleCache
  {
  public:
    MatrixFreeCellBatchParticleCache(const MatrixFreeContext<dim, number> mf_context)
      : mf_context(mf_context)
    {}
 
    void
    update(CellListParticleHandler<dim, number, ObstacleType> &particle_handler,
           const typename CellListParticleHandler<dim, number, ObstacleType>::NotifyEvent event);
 
    inline DEAL_II_ALWAYS_INLINE const std::vector<DEMParticleAccessor<dim, number>> &
    obstacles_in_cell_batch(const unsigned int cell_batch_id) const;
 
    inline DEAL_II_ALWAYS_INLINE std::vector<DEMParticleAccessor<dim, number>>                              &
    obstacles_in_cell_batch(const unsigned int cell_batch_id);
 
  private:
    const MatrixFreeContext<dim, number> mf_context;
 
    std::vector<std::vector<DEMParticleAccessor<dim, number>>> cell_batch_to_particle_cache;
  };
 
 
  template <int dim, typename number, typename ObstacleType>
  const std::vector<DEMParticleAccessor<dim, number>> &
  MatrixFreeCellBatchParticleCache<dim, number, ObstacleType>::obstacles_in_cell_batch(
    const unsigned int cell_batch_id) const
  {
    AssertIndexRange(cell_batch_id, cell_batch_to_particle_cache.size());
    return cell_batch_to_particle_cache[cell_batch_id];
  }
 
  template <int dim, typename number, typename ObstacleType>
  std::vector<DEMParticleAccessor<dim, number>> &
  MatrixFreeCellBatchParticleCache<dim, number, ObstacleType>::obstacles_in_cell_batch(
    const unsigned int cell_batch_id)
  {
    AssertIndexRange(cell_batch_id, cell_batch_to_particle_cache.size());
    return cell_batch_to_particle_cache[cell_batch_id];
  }
} // namespace MeltPoolDG