Utilities.hxx

Go to the documentation of this file.

 
#ifndef UTILITIES_HXX
#define UTILITIES_HXX
 
#include "AMReX_REAL.H"
#include "AMReX_Scan.H"
#include "cctk_Types.h"
#include <cctk.h>
 
AMREX_GPU_DEVICE AMREX_GPU_HOST
    AMREX_FORCE_INLINE CCTK_ATTRIBUTE_ALWAYS_INLINE amrex::GpuArray<CCTK_REAL, 3>
    SMatVecMul(amrex::GpuArray<CCTK_REAL, 6> A,
               amrex::GpuArray<CCTK_REAL, 3> V) {
 
  amrex::GpuArray<CCTK_REAL, 3> result = {
      A[0] * V[0] + A[1] * V[1] + A[2] * V[2],
      A[1] * V[0] + A[3] * V[1] + A[4] * V[2],
      A[2] * V[0] + A[4] * V[1] + A[5] * V[2]};
 
  return result;
} // SMatVecMul
 
AMREX_GPU_DEVICE
AMREX_GPU_HOST AMREX_FORCE_INLINE CCTK_ATTRIBUTE_ALWAYS_INLINE CCTK_REAL
VecVecMul(amrex::GpuArray<CCTK_REAL, 3> U, amrex::GpuArray<CCTK_REAL, 3> V) {
 
  return U[0] * V[0] + U[1] * V[1] + U[2] * V[2];
} // SMatVecMul
 
#endif // !UTILITIES_HXX