BoxLibParticle.h

Go to the documentation of this file.

#ifndef BOXLIB_PARTICLE_H
#define BOXLIB_PARTICLE_H

#include "AmrParticle/AmrParticleBase.h"

#include <AMReX_REAL.H>
#include <AMReX_IntVect.H>
#include <AMReX_Vector.H>
#include <AMReX_Utility.H>
#include <AMReX_Geometry.H>
#include <AMReX_RealBox.H>

template<class PLayout>
class BoxLibParticle : public virtual AmrParticleBase<PLayout>
{
public:
    typedef typename AmrParticleBase<PLayout>::ParticlePos_t             ParticlePos_t;
    typedef typename AmrParticleBase<PLayout>::ParticleIndex_t           ParticleIndex_t;
    typedef typename AmrParticleBase<PLayout>::SingleParticlePos_t       SingleParticlePos_t;
    typedef typename AmrParticleBase<PLayout>::AmrField_t                AmrField_t;
    // Array<std::unique_ptr<MultiFab> >
    typedef typename AmrParticleBase<PLayout>::AmrScalarFieldContainer_t AmrScalarFieldContainer_t;
    typedef typename AmrParticleBase<PLayout>::AmrVectorFieldContainer_t AmrVectorFieldContainer_t;
    typedef typename AmrParticleBase<PLayout>::AmrVectorField_t          AmrVectorField_t;
    typedef typename AmrParticleBase<PLayout>::ParticleLevelCounter_t    ParticleLevelCounter_t;
    
    typedef typename PLayout::AmrProcMap_t  AmrProcMap_t;
    typedef typename PLayout::AmrGrid_t     AmrGrid_t;
    typedef typename PLayout::AmrGeometry_t AmrGeometry_t;
    typedef typename PLayout::AmrIntVect_t  AmrIntVect_t;
    typedef typename PLayout::AmrBox_t      AmrBox_t;
    typedef typename PLayout::AmrReal_t     AmrReal_t;
    
    typedef amrex::FArrayBox                FArrayBox_t;
    
public:
    BoxLibParticle();
    
    BoxLibParticle(PLayout *layout);
    
    template <class FT, unsigned Dim, class PT>
    void scatter(ParticleAttrib<FT>& attrib, AmrScalarFieldContainer_t& f,
                 ParticleAttrib<Vektor<PT, Dim> >& pp,
                 int lbase, int lfine,
                 const ParticleAttrib<int>& pbin, int bin = -1);
    
    
    template <class FT, unsigned Dim, class PT>
    void scatter(ParticleAttrib<FT>& attrib, AmrField_t& f,
                 ParticleAttrib<Vektor<PT, Dim> >& pp,
                 const ParticleAttrib<int>& pbin, int bin = -1,
                 int level = 0);
    
    template <class FT, unsigned Dim, class PT>
    void gather(ParticleAttrib<FT>& attrib, AmrVectorFieldContainer_t& f,
                ParticleAttrib<Vektor<PT, Dim> >& pp,
                int lbase, int lfine);
    
private:
    /*
     * AMReX functions adjusted to work with Ippl
     */
    
    template <class AType>
    void AssignDensityFort(ParticleAttrib<AType> &pa,
                           AmrScalarFieldContainer_t& mf_to_be_filled, 
                           int lev_min, int ncomp, int finest_level,
                           const ParticleAttrib<int>& pbin, int bin = -1) const;
    
    template <class AType>
    void InterpolateFort(ParticleAttrib<AType> &pa,
                         AmrVectorFieldContainer_t& mesh_data, 
                         int lev_min, int lev_max);
    
    template <class AType>
    void InterpolateSingleLevelFort(ParticleAttrib<AType> &pa, AmrVectorField_t& mesh_data, int lev);
    
    
    template <class AType>
    void InterpolateMultiLevelFort(ParticleAttrib<AType> &pa,
                                   AmrVectorFieldContainer_t& mesh_data,
                                   int lev);
    
    
    template <class AType>
    void AssignCellDensitySingleLevelFort(ParticleAttrib<AType> &pa, AmrField_t& mf, int level,
                                          const ParticleAttrib<int>& pbin, int bin = -1,
                                          int ncomp=1, int particle_lvl_offset = 0) const;
    
private:
    IpplTimings::TimerRef AssignDensityTimer_m;
};


#include "BoxLibParticle.hpp"

#endif