Merge branch 'development' into derive_vector_multifab

Author: WeiqunZhang

.github/workflows/cuda.yml

@@ -20,6 +20,7 @@ jobs:
     - uses: actions/checkout@v5
     - name: Dependencies
       run: |
+        .github/workflows/dependencies/ubuntu_free_disk_space.sh
         .github/workflows/dependencies/dependencies_nvcc.sh 12.6
         .github/workflows/dependencies/dependencies_ccache.sh
     - name: Set Up Cache
@@ -82,6 +83,7 @@ jobs:
     - uses: actions/checkout@v5
     - name: Dependencies
       run: |
+        .github/workflows/dependencies/ubuntu_free_disk_space.sh
         .github/workflows/dependencies/dependencies_llvm_cuda11_clang15.sh
         .github/workflows/dependencies/dependencies_ccache.sh
     - name: Set Up Cache
@@ -193,6 +195,7 @@ jobs:
     - uses: actions/checkout@v5
     - name: Dependencies
       run: |
+        .github/workflows/dependencies/ubuntu_free_disk_space.sh
         .github/workflows/dependencies/dependencies_nvcc_2404.sh 13.0
         .github/workflows/dependencies/dependencies_ccache.sh
     - name: Set Up Cache

Src/Base/AMReX_Math.H

@@ -325,6 +325,34 @@ T comp_ellint_2 (T k)
     return Kcomp*sum_val;
 }
 
+//! Return inverse square root of x
+AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
+double rsqrt (double x)
+{
+    double r;
+#if defined(AMREX_USE_SYCL)
+    AMREX_IF_ON_DEVICE(( r = sycl::rsqrt(x); ))
+#else
+    AMREX_IF_ON_DEVICE(( r = ::rsqrt(x); ))
+#endif
+    AMREX_IF_ON_HOST(( r = 1. / std::sqrt(x); ))
+    return r;
+}
+
+//! Return inverse square root of x
+AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
+float rsqrt (float x)
+{
+    float r;
+#if defined(AMREX_USE_SYCL)
+    AMREX_IF_ON_DEVICE(( r = sycl::rsqrt(x); ))
+#else
+    AMREX_IF_ON_DEVICE(( r = ::rsqrtf(x); ))
+#endif
+    AMREX_IF_ON_HOST(( r = 1.F / std::sqrt(x); ))
+    return r;
+}
+
 /***************************************************************************************************
  * Copyright (c) 2017 - 2024 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
  * SPDX-License-Identifier: BSD-3-Clause

Src/Particle/AMReX_ParticleContainer.H

@@ -1370,14 +1370,14 @@ public:
     /** Create an empty particle container
      *
      * This creates a new AMReX particle container type with same compile-time
-     * and run-time attributes. But, it can change its allocator. This is
-     * helpful when creating temporary particle buffers for filter operations
+     * and run-time attributes. But, it can change its allocator (default: same allocator).
+     * This is helpful when creating temporary particle buffers for filter operations
      * and device-to-host copies.
      *
-     * @tparam Allocator AMReX allocator, e.g., amrex::PinnedArenaAllocator
+     * @tparam Allocator AMReX allocator, e.g., amrex::PinnedArenaAllocator, default: same allocator as the creating PC
      * @return an empty particle container
      * */
-    template <template<class> class NewAllocator=amrex::DefaultAllocator>
+    template <template<class> class NewAllocator=Allocator>
     ContainerLike<NewAllocator>
     make_alike () const
     {