(const vec_t *)intmp.data(), (const vec_t *)D0, \
(const vec_t *)D1, (const vec_t *)D2, \
(vec_t *)out_d0.data(), (vec_t *)out_d1.data(), \
(vec_t *)out_d2.data(), Deriv0, Deriv1, Deriv2)
static NEK_FORCE_INLINE void PhysDerivTensor3DKernel(const unsigned int nq0, const unsigned int nq1, const unsigned int nq2, const simd_type *in, const simd_type *D0, const simd_type *D1, const simd_type *D2, simd_type *out_d0, simd_type *out_d1, simd_type *out_d2, bool Deriv0=true, bool Deriv1=true, bool Deriv2=true)