ADMconstraints: update admconstraints.cxx

EinsteinToolkit · Nov 18, 2024 · e0141d5 · e0141d5
1 parent 49cf5bf
commit e0141d5
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/ADMconstraints/src/admconstraints.cxx b/ADMconstraints/src/admconstraints.cxx
@@ -6,6 +6,11 @@ namespace ADMconstraints {
 using namespace Arith;
 using namespace Loop;
 
+template <typename T>
+CCTK_DEVICE CCTK_HOST CCTK_ATTRIBUTE_ALWAYS_INLINE inline T Power(T x, int y) {
+  return (y == 2) ? Arith::pow2(x) : Arith::pown(x, y);
+}
+
 extern "C" void ADMconstraints(CCTK_ARGUMENTS) {
   DECLARE_CCTK_ARGUMENTSX_ADMconstraints;
   DECLARE_CCTK_PARAMETERS;