Interpolate.H

Go to the documentation of this file.

/* Copyright 2019-2022 Axel Huebl, Remi Lehe
 *
 * This file is part of WarpX.
 *
 * License: BSD-3-Clause-LBNL
 */
#ifndef ABLASTR_INTERPOLATE_H
#define ABLASTR_INTERPOLATE_H
 
#include <AMReX_Array4.H>
#include <AMReX_IntVect.H>
#include <AMReX_MFInterp_C.H>
 
#include <array>
#include <optional>
 
namespace ablastr::fields::details {

    struct PoissonInterpCPtoFP
    {
        PoissonInterpCPtoFP(
                amrex::Array4<amrex::Real> const phi_fp_arr,
                amrex::Array4<amrex::Real const> const phi_cp_arr,
                amrex::IntVect const refratio)
        : m_phi_fp_arr(phi_fp_arr), m_phi_cp_arr(phi_cp_arr), m_refratio(refratio)
        {}
 
        AMREX_GPU_DEVICE AMREX_FORCE_INLINE
        void
        operator() (int i, int j, int k) const noexcept
        {
            amrex::mf_nodebilin_interp(i, j, k, 0, m_phi_fp_arr, 0, m_phi_cp_arr,
                                       0, m_refratio);
        }
 
        amrex::Array4<amrex::Real> m_phi_fp_arr;
        amrex::Array4<amrex::Real const> m_phi_cp_arr;
        amrex::IntVect m_refratio;
    };
} // namespace ablastr::fields::details
 
#endif // ABLASTR_INTERPOLATE_H