diff --git a/src/cosserat-continuum.cc b/src/cosserat-continuum.cc
index 6e1336705369c7c437b4a42b299015a9e9e5273e..f9880cfd097a151751115a1467e22a205ca4ebb9 100644
--- a/src/cosserat-continuum.cc
+++ b/src/cosserat-continuum.cc
@@ -438,7 +438,7 @@ int main (int argc, char *argv[]) try
       xEmbedded[i] = x[i].globalCoordinates();
 
     std::ofstream outFile("cosserat-continuum-result-" + std::to_string(numLevels) + ".data", std::ios_base::binary);
-    GenericVector::writeBinary(outFile, xEmbedded);
+    MatrixVector::Generic::writeBinary(outFile, xEmbedded);
     outFile.close();
 
     // finally: compute the average deformation of the Neumann boundary
diff --git a/src/gradient-flow.cc b/src/gradient-flow.cc
index aaef1da78d436bfba5d2730fe7654bce4933211b..6d3a79cc82075a4b8a5ab8673b0f137478198558 100644
--- a/src/gradient-flow.cc
+++ b/src/gradient-flow.cc
@@ -259,7 +259,7 @@ int main (int argc, char *argv[]) try
 
   // 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);
+  MatrixVector::Generic::writeBinary(outFile, xEmbedded);
   outFile.close();
 
   ///////////////////////////////////////////////////////
@@ -299,7 +299,7 @@ int main (int argc, char *argv[]) try
 
     // 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);
+    MatrixVector::Generic::writeBinary(outFile, xEmbedded);
     outFile.close();
 
   }