Hadrons: random vector utility module I/O

Hadrons/Modules/MUtilities/RandomVectors.hpp

@@ -31,6 +31,7 @@ See the full license in the file "LICENSE" in the top level distribution directo
 #include <Hadrons/Global.hpp>
 #include <Hadrons/Module.hpp>
 #include <Hadrons/ModuleFactory.hpp>
+#include <Hadrons/A2AVectors.hpp>
 
 BEGIN_HADRONS_NAMESPACE
 
@@ -44,7 +45,9 @@ class RandomVectorsPar: Serializable
 public:
     GRID_SERIALIZABLE_CLASS_MEMBERS(RandomVectorsPar,
                                     unsigned int, size,
-                                    unsigned int, Ls);
+                                    unsigned int, Ls,
+                                    std::string,  output,
+                                    bool,         multiFile);
 };
 
 template <typename Field>
@@ -119,6 +122,11 @@ void TRandomVectors<Field>::execute(void)
     {
         random(rng4d(), vec[i]);
     }
+    // I/O if necessary
+    if (!par().output.empty())
+    {
+        A2AVectorsIo::write(par().output, vec, par().multiFile, vm().getTrajectory());
+    }
 }
 
 END_MODULE_NAMESPACE