ERF::Interpolation Namespace Reference

Functions
template<typename ParticleType >
AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	interpolateFields (const ParticleType &p, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &plo, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxi, const amrex::Array4< const amrex::Real > *field_arrays, amrex::ParticleReal *result_values, int num_fields, int is_periodic_z=1, const amrex::Array4< amrex::Real > *zheight=nullptr)
	Helper function to interpolate multiple field values at a particle position. More...

Function Documentation

interpolateFields()

template<typename ParticleType >

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void ERF::Interpolation::interpolateFields	(	const ParticleType &	p,
		const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &	plo,
		const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &	dxi,
		const amrex::Array4< const amrex::Real > *	field_arrays,
		amrex::ParticleReal *	result_values,
		int	num_fields,
		int	is_periodic_z = 1,
		const amrex::Array4< amrex::Real > *	zheight = nullptr
	)

Helper function to interpolate multiple field values at a particle position.

This function interpolates values from multiple field arrays at a particle's position using the cloud-in-cell (CIC) interpolation method. It supports both periodic and non-periodic z-boundaries.

Template Parameters

ParticleType

The type of particle

Parameters

p	The particle
plo	The domain lower bounds
dxi	The inverse of cell spacing
field_arrays	Array of field arrays to interpolate from
result_values	Output array to store interpolated values
num_fields	Number of fields to interpolate
is_periodic_z	Whether the z-boundary is periodic
zheight	Height field array (only used for non-periodic z-boundary)

  {
    // Iterate over each field and interpolate its value at the particle position
    for (int i = 0; i < num_fields; ++i) {
        if (is_periodic_z != 0) {
            cic_interpolate(p, plo, dxi, field_arrays[i], &result_values[i], 1);
        } else {
            cic_interpolate_mapped_z(p, plo, dxi, field_arrays[i], *zheight, &result_values[i], 1);
        }
    }
}