ENH: SFINAE is_vectorspace in adiosWritePrimitives

io/reduced-foam-extend/src/foam/db/IOstreams/A2streams/adiosWritePrimitives.C

@@ -28,7 +28,16 @@ License
 #include "adiosWriting.H"
 #include "adiosFileStream.H"
 
-#include "vector.H"
+#include "foamString.H"
+#include "labelList.H"
+
+Foam::List<Foam::label> Foam::create2DList( const Foam::label& val0,
+                                            const Foam::label& val1 ) {
+    Foam::List<Foam::label> ret{2};
+    ret[0] = val0;
+    ret[1] = val1;
+    return ret;
+};
 
 void Foam::adiosWritePrimitives( const Foam::string type,
                                  const Foam::string blockId,
@@ -50,18 +59,13 @@ void Foam::adiosWritePrimitives( const Foam::string type,
 
 void Foam::adiosWritePrimitives( const Foam::string type,
                                  const Foam::string blockId,
-                                 const Foam::label count,
-                                 const Foam::Vector<Foam::scalar>* buf ) {
-    Foam::labelList iShape{2};
-    iShape[0] = count;
-    iShape[1] = buf[0].size();
-    Foam::labelList iStart{2};
-    iStart[0] = iStart[1] = 0;
-    Foam::labelList iCount{ iShape };
-
+                                 const Foam::List<label> shape,
+                                 const Foam::List<label> start,
+                                 const Foam::List<label> count,
+                                 const Foam::scalar* buf ) {
     auto adiosStreamPtr = adiosWriting{}.createStream();
     adiosStreamPtr->open( type );
-    adiosStreamPtr->transfer( blockId, iShape, iStart, iCount, reinterpret_cast<const Foam::scalar*>( buf ) );
+    adiosStreamPtr->transfer( blockId, shape, start, count, buf );
 }
 
 // ************************************************************************* //

io/reduced-foam-extend/src/foam/db/IOstreams/A2streams/adiosWritePrimitives.H

@@ -26,25 +26,37 @@ License
 #ifndef adiosWritePrimitives_H
 #define adiosWritePrimitives_H
 
-#include "foamString.H"
 #include "label.H"
-#include "uLabel.H"
 #include "scalar.H"
+#include "direction.H"
 
-#include <typeinfo>
+#include <type_traits> // defines false_type, true_typ, enable_if, void_t
+#include <utility> // defines declval
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam
 {
 
-template<typename T> class Vector;
+// is_vectorspace
+// primary template
+template<typename, typename = std::void_t<>>
+struct is_vectorspace : std::false_type {};
 
+// partial specialization, maybe SFINAE'd away
 template<typename T>
-void adiosWritePrimitives( const string type,
-                           const string name,
-                           const label count,
-                           const T* buf );
+struct is_vectorspace< T,
+                       std::void_t<
+                           decltype( std::declval<T>().size() ),
+                           decltype( std::declval<T>().component( std::declval<direction>() ) ),
+                           decltype( std::declval<T>()[ std::declval<direction>() ] )
+                       >
+                     > : std::true_type {};
+
+class string;
+template<typename T> class List;
+
+List<label> create2DList( const label& val0, const label& val1 );
 
 void adiosWritePrimitives( const string type,
                            const string blockId,
@@ -58,14 +70,16 @@ void adiosWritePrimitives( const string type,
 
 void adiosWritePrimitives( const string type,
                            const string blockId,
-                           const label count,
-                           const Vector<scalar>* buf );
+                           const List<label> shape,
+                           const List<label> start,
+                           const List<label> count,
+                           const scalar* buf );
 
 } // End namespace Foam
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-template<typename T>
+template< typename T, std::enable_if_t< !Foam::is_vectorspace<T>::value, void >* = nullptr >
 void Foam::adiosWritePrimitives( const Foam::string type,
                                  const Foam::string name,
                                  const Foam::label count,
@@ -75,6 +89,22 @@ void Foam::adiosWritePrimitives( const Foam::string type,
     std::cout << name << "\n";
 }
 
+template< typename T, std::enable_if_t< Foam::is_vectorspace<T>::value, void >* = nullptr >
+void Foam::adiosWritePrimitives( const Foam::string type,
+                                 const Foam::string name,
+                                 const Foam::label count,
+                                 const T* buf ) {
+    Foam::List<Foam::label> shapeList = create2DList( count, buf[0].size() );
+    Foam::List<Foam::label> startList = create2DList( 0, 0 );
+    Foam::List<Foam::label> countList{ shapeList };
+    adiosWritePrimitives( type,
+                          name,
+                          shapeList,
+                          startList,
+                          countList,
+                          reinterpret_cast<const Foam::scalar*>( buf ) );
+}
+
 #endif
 
 // ************************************************************************* //