Kokkos Portability Layer

Adds portability via Kokkos:

• Adds Kokkos as a dependency
• Updates Cabana::AoSoA template parameters to use syntax and structure identical to Kokkos::View
• Adds a Cabana::InnerArraySize template parameter which allows users to manually set the inner array size or will be set to default option based on the Kokkos::ExecutionSpace
• Adds a Cabana::parallel_for interface which will dispatch to Kokkos::parallel_for with a variety of Kokkos::RangePolicy and Kokkos::MDRangePolicy options based on the execution space and the struct/array sizes.
• Adds a portable testing infrastructure which currently support Serial, OpenMP, and CudaUVM executions and allocations.
• Adds Cabana::PerformanceTraits which can be used to define reasonably performant settings for different execution spaces. These can always be overridden by the user.

Future work that will not be performed in this MR includes:

• Utilizing the Kokkos memory tracking infrastructure instead of std::shared_ptr
• Add a deep copy capability to Cabana::AoSoA (see #6 (closed))
• Add CudaSpace testing support after addition of deep copy.

core/src/Cabana_Parallel.hpp

@@ -50,7 +50,7 @@ inline void parallel_for( const ExecutionPolicy& exec_policy,
     // each thread loops over the inner arrays.
     std::size_t array_size = begin.a();
     auto functor_wrapper =
-        KOKKOS_LAMBDA( std::size_t s )
+        KOKKOS_LAMBDA( const std::size_t s )
         {
             std::size_t i_begin = (s == s_begin) ? begin.i() : 0;
             std::size_t i_end = ((s == s_end - 1) && (end.i() != 0))
@@ -99,7 +99,7 @@ inline void parallel_for( const ExecutionPolicy& exec_policy,
         // Create a wrapper for the functor. Each struct is given a thread and
         // each thread loops over the inner arrays.
         auto functor_wrapper =
-            KOKKOS_LAMBDA( std::size_t i )
+            KOKKOS_LAMBDA( const std::size_t i )
             {
                 Index idx( array_size, s, i );
                 functor( idx );
@@ -140,7 +140,7 @@ inline void parallel_for( const ExecutionPolicy& exec_policy,
 
     // Create a wrapper for the functor.
     auto functor_wrapper =
-        KOKKOS_LAMBDA( std::size_t s, std::size_t i )
+        KOKKOS_LAMBDA( const std::size_t s, const std::size_t i )
         {
             Index idx( array_size, s, i );
             if ( idx >= begin && idx < end ) functor( idx );