CopyParticleAttribs.H

Go to the documentation of this file.

/* Copyright 2020 Andrew Myers
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
#ifndef WARPX_PARTICLES_PUSHER_COPYPARTICLEATTRIBS_H_
#define WARPX_PARTICLES_PUSHER_COPYPARTICLEATTRIBS_H_
 
#include "Particles/WarpXParticleContainer.H"
 
#include <AMReX_REAL.H>
 
#include <limits>
 
 
struct CopyParticleAttribs
{
    using TmpParticles = WarpXParticleContainer::TmpParticles;
 
    GetParticlePosition m_get_position;
 
    const amrex::ParticleReal* AMREX_RESTRICT uxp = nullptr;
    const amrex::ParticleReal* AMREX_RESTRICT uyp = nullptr;
    const amrex::ParticleReal* AMREX_RESTRICT uzp = nullptr;
 
    amrex::ParticleReal* AMREX_RESTRICT xpold = nullptr;
    amrex::ParticleReal* AMREX_RESTRICT ypold = nullptr;
    amrex::ParticleReal* AMREX_RESTRICT zpold = nullptr;
 
    amrex::ParticleReal* AMREX_RESTRICT uxpold = nullptr;
    amrex::ParticleReal* AMREX_RESTRICT uypold = nullptr;
    amrex::ParticleReal* AMREX_RESTRICT uzpold = nullptr;
 
    CopyParticleAttribs () = default;
 
    CopyParticleAttribs (const WarpXParIter& a_pti, TmpParticles& tmp_particle_data,
                         int a_offset = 0) noexcept
    {
        if (tmp_particle_data.size() == 0) return;
 
        auto& attribs = a_pti.GetAttribs();
 
        uxp = attribs[PIdx::ux].dataPtr() + a_offset;
        uyp = attribs[PIdx::uy].dataPtr() + a_offset;
        uzp = attribs[PIdx::uz].dataPtr() + a_offset;
 
        const auto lev = a_pti.GetLevel();
        const auto index = a_pti.GetPairIndex();
        xpold  = tmp_particle_data[lev].at(index)[TmpIdx::xold ].dataPtr() + a_offset;
        ypold  = tmp_particle_data[lev].at(index)[TmpIdx::yold ].dataPtr() + a_offset;
        zpold  = tmp_particle_data[lev].at(index)[TmpIdx::zold ].dataPtr() + a_offset;
        uxpold = tmp_particle_data[lev].at(index)[TmpIdx::uxold].dataPtr() + a_offset;
        uypold = tmp_particle_data[lev].at(index)[TmpIdx::uyold].dataPtr() + a_offset;
        uzpold = tmp_particle_data[lev].at(index)[TmpIdx::uzold].dataPtr() + a_offset;
 
        m_get_position = GetParticlePosition(a_pti, a_offset);
    }
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
    void operator() (const int i) const noexcept
    {
        AMREX_ASSERT(uxp != nullptr);
        AMREX_ASSERT(uyp != nullptr);
        AMREX_ASSERT(uzp != nullptr);
 
        AMREX_ASSERT(xpold != nullptr);
        AMREX_ASSERT(ypold != nullptr);
        AMREX_ASSERT(zpold != nullptr);
 
        AMREX_ASSERT(uxpold != nullptr);
        AMREX_ASSERT(uypold != nullptr);
        AMREX_ASSERT(uzpold != nullptr);
 
        amrex::ParticleReal x, y, z;
        m_get_position(i, x, y, z);
 
        xpold[i] = x;
        ypold[i] = y;
        zpold[i] = z;
 
        uxpold[i] = uxp[i];
        uypold[i] = uyp[i];
        uzpold[i] = uzp[i];
    }
};
 
#endif // WARPX_PARTICLES_PUSHER_COPYPARTICLEATTRIBS_H_