Fix requirements for simd type

Author: niermann999

common/include/algebra/concepts.hpp

@@ -19,10 +19,19 @@ namespace algebra::concepts {
 template <typename T>
 concept arithmetic = std::is_arithmetic_v<T>;
 
+/// Arithmetic types
+template <typename T>
+concept enumeration = std::is_enum_v<T>;
+
 // Value concept: Single entry
 template <typename T>
 concept value = algebra::concepts::arithmetic<std::decay_t<T>>;
 
+// Element concept: Single entry, not necessarily for vector/matrix operations
+template <typename T>
+concept element = algebra::concepts::value<T> ||
+                  algebra::concepts::enumeration<std::decay_t<T>>;
+
 /// Scalar concept: Elements of vectors/matrices (can be simd vectors)
 template <typename T>
 concept scalar = !algebra::traits::is_matrix<T> &&

common/include/algebra/type_traits.hpp

@@ -33,13 +33,13 @@ using index_t = typename index<M>::type;
 template <class M>
 struct value {};
 
-template <class M>
-using value_t = typename value<M>::type;
-
 template <typename T>
 requires(std::is_arithmetic_v<T>&& std::is_scalar_v<T>) struct value<T> {
   using type = T;
 };
+
+template <class M>
+using value_t = typename value<M>::type;
 /// @}
 
 /// Scalar type that is used with the matrix (can be multiple values in SoA)

frontend/array_cmath/include/algebra/array_cmath.hpp

@@ -130,7 +130,7 @@ struct array {
   /// Define scalar type
   using value_type = V;
 
-  template <concepts::value T>
+  template <concepts::element T>
   using simd = T;
 
   using boolean = bool;

frontend/eigen_eigen/include/algebra/eigen_eigen.hpp

@@ -94,7 +94,7 @@ struct eigen {
   /// Define scalar type
   using value_type = V;
 
-  template <concepts::value T>
+  template <concepts::element T>
   using simd = T;
 
   using boolean = bool;

frontend/eigen_generic/include/algebra/eigen_generic.hpp

@@ -106,7 +106,7 @@ struct eigen_generic {
   /// Define scalar type
   using value_type = V;
 
-  template <concepts::value T>
+  template <concepts::element T>
   using simd = T;
 
   using boolean = bool;

frontend/fastor_fastor/include/algebra/fastor_fastor.hpp

@@ -101,7 +101,7 @@ struct fastor {
   /// Define scalar type
   using value_type = V;
 
-  template <concepts::value T>
+  template <concepts::element T>
   using simd = T;
 
   using boolean = bool;

frontend/smatrix_generic/include/algebra/smatrix_generic.hpp

@@ -97,7 +97,7 @@ struct smatrix_generic {
   /// Define scalar type
   using value_type = V;
 
-  template <concepts::value T>
+  template <concepts::element T>
   using simd = T;
 
   using boolean = bool;

frontend/smatrix_smatrix/include/algebra/smatrix_smatrix.hpp

@@ -92,7 +92,7 @@ struct smatrix {
   /// Define scalar type
   using value_type = V;
 
-  template <concepts::value T>
+  template <concepts::element T>
   using simd = T;
 
   using boolean = bool;

frontend/vc_aos/include/algebra/vc_aos.hpp

@@ -119,7 +119,7 @@ struct vc_aos {
   /// Define scalar type
   using value_type = V;
 
-  template <concepts::value T>
+  template <concepts::element T>
   using simd = T;
 
   using boolean = bool;

frontend/vc_aos_generic/include/algebra/vc_aos_generic.hpp

@@ -105,7 +105,7 @@ struct vc_aos_generic {
   /// Define scalar type
   using value_type = V;
 
-  template <concepts::value T>
+  template <concepts::element T>
   using simd = T;
 
   using boolean = bool;