Adding wrappers for device math functions and polar decomposition

cmake/Config.cmake

@@ -2,7 +2,7 @@
 set( PREPROCESSOR_DEFINES UMPIRE
                           CHAI
                           CUDA
-			  HIP
+                          HIP
                           TOTALVIEW_OUTPUT
                           CALIPER )
 

host-configs/LLNL/dane-clang-14-mpm.cmake

@@ -0,0 +1,5 @@
+set(CONFIG_NAME "dane-clang-14-mpm" CACHE PATH "")
+
+set(GEOS_TPL_DIR /usr/WS1/crook5/thirdPartyLibs/install-dane-clang-14-mpm-release CACHE PATH "")
+
+include(${CMAKE_CURRENT_LIST_DIR}/llnl-cpu-clang-14.cmake)

src/ArrayOfArrays.hpp

@@ -484,4 +484,21 @@ class ArrayOfArrays : protected ArrayOfArraysView< T, INDEX_TYPE, false, BUFFER_
   }
 };
 
+/**
+ * @brief True if the template type is an ArrayOfArrays.
+ */
+template< class >
+constexpr bool isArrayOfArrays = false;
+
+/**
+ * @tparam T The type contained in the ArrayOfArrays.
+ * @tparam INDEX_TYPE The integral type used as an index.
+ * @tparam BUFFER_TYPE The type used to manage the underlying allocation.
+ * @brief Specialization of isArrayOfArrays for the ArrayOfArrays class.
+ */
+template< typename T,
+          typename INDEX_TYPE,
+          template< typename > class BUFFER_TYPE >
+constexpr bool isArrayOfArrays< ArrayOfArrays< T, INDEX_TYPE, BUFFER_TYPE > > = true;
+
 } /* namespace LvArray */

src/ArrayOfArraysView.hpp

@@ -249,7 +249,7 @@ class ArrayOfArraysView
 
   /**
    * @brief Move assignment operator..
-   * @param src the SparsityPatternView to be moved from.
+   * @param src the ArrayOfArraysView to be moved from.
    * @return *this.
    */
   LVARRAY_HOST_DEVICE
@@ -1087,4 +1087,23 @@ class ArrayOfArraysView
   }
 };
 
+/**
+ * @brief True if the template type is an ArrayOfArraysView.
+ */
+template< class >
+constexpr bool isArrayOfArraysView = false;
+
+/**
+ * @tparam T The type contained in the ArrayOfArraysView.
+ * @tparam INDEX_TYPE The integral type used as an index.
+ * @tparam CONST_SIZES True iff the size of each array is constant.
+ * @tparam BUFFER_TYPE The type used to manager the underlying allocation.
+ * @brief Specialization of isArrayOfArraysView for the ArrayOfArraysView class.
+ */
+template< typename T,
+          typename INDEX_TYPE,
+          bool CONST_SIZES,
+          template< typename > class BUFFER_TYPE >
+constexpr bool isArrayOfArraysView< ArrayOfArraysView< T, INDEX_TYPE, CONST_SIZES, BUFFER_TYPE > > = true;
+
 } /* namespace LvArray */

src/ChaiBuffer.hpp

@@ -340,7 +340,7 @@ class ChaiBuffer
 
     if( size > 0 )
     {
-      LVARRAY_ERROR_IF_NE_MSG( space, MemorySpace::host, "Calling reallocate with a non-zero current size is not yet supporeted for the GPU." );
+      LVARRAY_ERROR_IF_NE_MSG( space, MemorySpace::host, "Calling reallocate with a non-zero current size is not yet supported for the GPU." );
       std::ptrdiff_t const overlapAmount = std::min( newCapacity, size );
       arrayManipulation::uninitializedMove( newPointer, overlapAmount, m_pointer );
       arrayManipulation::destroy( m_pointer, size );

src/fixedSizeSquareMatrixOps.hpp

@@ -38,6 +38,42 @@ LVARRAY_HOST_DEVICE CONSTEXPR_WITHOUT_BOUNDS_CHECK inline
 auto determinant( MATRIX const & matrix )
 { return internal::SquareMatrixOps< M >::determinant( matrix ); }
 
+/**
+ * @brief Compute the cofactor of the source matrix @p srcMatrix and store the result in @p dstMatrix.
+ * @tparam M The size of the matrices @p dstMatrix and @p srcMatrix.
+ * @tparam DST_MATRIX The type of @p dstMatrix.
+ * @tparam SRC_MATRIX The type of @p srcMatrix.
+ * @param dstMatrix The M x M matrix to write the cofactor to.
+ * @param srcMatrix The M x M matrix to take the cofactor of.
+ * @note @p srcMatrix can contain integers but @p dstMatrix must contain floating point values.
+ */
+template< std::ptrdiff_t M, typename DST_MATRIX, typename SRC_MATRIX >
+LVARRAY_HOST_DEVICE CONSTEXPR_WITHOUT_BOUNDS_CHECK inline
+auto cofactor( DST_MATRIX && LVARRAY_RESTRICT_REF dstMatrix,
+               SRC_MATRIX const & LVARRAY_RESTRICT_REF srcMatrix )
+{
+  static_assert( std::is_floating_point< std::decay_t< decltype( dstMatrix[ 0 ][ 0 ] ) > >::value,
+                 "The destination matrix must be contain floating point values." );
+  internal::SquareMatrixOps< M >::cofactor( std::forward< DST_MATRIX >( dstMatrix ), srcMatrix );
+}
+
+/**
+ * @brief Compute the cofactor of the matrix @p matrix overwritting it.
+ * @tparam M The size of the matrix @p matrix.
+ * @tparam MATRIX The type of @p matrix.
+ * @param matrix The M x M matrix to take the cofactor of and overwrite.
+ * @return The determinant.
+ * @note @p matrix must contain floating point values.
+ */
+template< std::ptrdiff_t M, typename MATRIX >
+LVARRAY_HOST_DEVICE CONSTEXPR_WITHOUT_BOUNDS_CHECK inline
+auto cofactor( MATRIX && matrix )
+{
+  static_assert( std::is_floating_point< std::decay_t< decltype( matrix[ 0 ][ 0 ] ) > >::value,
+                 "The matrix must be contain floating point values." );
+  internal::SquareMatrixOps< M >::cofactor( std::forward< MATRIX >( matrix ) );
+}
+
 /**
  * @brief Invert the source matrix @p srcMatrix and store the result in @p dstMatrix.
  * @tparam M The size of the matrices @p dstMatrix and @p srcMatrix.
@@ -354,6 +390,24 @@ void symmetricToDense( DST_MATRIX && dstMatrix, SRC_SYM_MATRIX const & srcSymMat
                                                            srcSymMatrix );
 }
 
+/**
+ * @brief Determine the polar decomposition of the matrix @p srcMatrix
+ * @tparam M The size of @p R and @p srcMatrix.
+ * @tparam DST_MATRIX The type of @p R.
+ * @tparam MATRIX The type of @p srcMatrix.
+ * @param R The resultant rotation matrix.
+ * @param matrix The matrix to be decomposed.
+ * @details The polar decomposition returns a rotation matrix such that @p R . U = V . @p R = @p srcMatrix.
+ *   This is done using Higham's iterative algorithm.
+ */
+template< std::ptrdiff_t M, typename DST_MATRIX, typename MATRIX >
+LVARRAY_HOST_DEVICE constexpr inline
+void polarDecomposition( DST_MATRIX && R, MATRIX const & srcMatrix )
+{
+  return internal::SquareMatrixOps< M >::polarDecomposition( std::forward< DST_MATRIX >( R ),
+                                                             srcMatrix );
+}
+
 ///@}
 
 } // namespace tensorOps