10 #ifndef IMPACTX_TO_FIXED_S_H
11 #define IMPACTX_TO_FIXED_S_H
52 amrex::ParticleReal & x,
53 amrex::ParticleReal & y,
54 amrex::ParticleReal & z,
55 amrex::ParticleReal & px,
56 amrex::ParticleReal & py,
57 amrex::ParticleReal & pz)
const
59 using namespace amrex::literals;
62 amrex::ParticleReal
const argd = 1.0_prt +
pow(
m_pzd, 2);
64 amrex::ParticleReal
const ptdf = argd > 0.0_prt ? -
sqrt(argd) : -1.0_prt;
73 amrex::ParticleReal
const arg = 1.0_prt +
pow(
m_pzd + pz, 2) +
pow(px, 2) +
pow(py, 2);
75 amrex::ParticleReal
const ptf =
arg > 0.0_prt ? -
sqrt(
arg) : -1.0_prt;
80 y = y - py * z / (
m_pzd + pz);
83 t = ptf * z / (
m_pzd + pz);
#define AMREX_ASSERT_WITH_MESSAGE(EX, MSG)
#define AMREX_FORCE_INLINE
#define AMREX_GPU_HOST_DEVICE
T_ParticleType ParticleType
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE T arg(const GpuComplex< T > &a_z) noexcept
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE GpuComplex< T > pow(const GpuComplex< T > &a_z, const T &a_y) noexcept
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE GpuComplex< T > sqrt(const GpuComplex< T > &a_z) noexcept
@ t
fixed t as the independent variable