diff --git a/src/ModFile.cc b/src/ModFile.cc
index a32dc93b..48b4489d 100644
--- a/src/ModFile.cc
+++ b/src/ModFile.cc
@@ -957,6 +957,9 @@ ModFile::writeMOutput(const string &basename, bool clear_all, bool clear_global,
 
   symbol_table.writeOutput(mOutputFile);
 
+  // Fill the fields related to solve_algo={12,14} if possible.
+  mOutputFile << "M_ = setup_solvers(M_);" << endl;
+
   var_model_table.writeOutput(basename, mOutputFile);
   trend_component_model_table.writeOutput(basename, mOutputFile);