Also write the results as binary data, to alloc EOC measurements

src/gradient-flow.cc

@@ -253,6 +253,11 @@ int main (int argc, char *argv[]) try
 #endif
   vtkWriter.write("gradientflow_result_0");
 
+  // Write the corresponding coefficient vector: verbatim in binary, to be completely lossless
+  std::ofstream outFile("gradientflow_result_0.data", std::ios_base::binary);
+  GenericVector::writeBinary(outFile, xEmbedded);
+  outFile.close();
+
   ///////////////////////////////////////////////////////
   //   Time loop
   ///////////////////////////////////////////////////////
@@ -291,6 +296,12 @@ int main (int argc, char *argv[]) try
     vtkWriter.addVertexData(vtkFunction(xFunction), VTK::FieldInfo("orientation", VTK::FieldInfo::Type::scalar, xEmbedded[0].size()));
 #endif
     vtkWriter.write("gradientflow_result_" + std::to_string(i+1));
+
+    // Write the corresponding coefficient vector: verbatim in binary, to be completely lossless
+    std::ofstream outFile("gradientflow_result_" + std::to_string(i+1) + ".data", std::ios_base::binary);
+    GenericVector::writeBinary(outFile, xEmbedded);
+    outFile.close();
+
   }
 
   return 0;