Added call to a routine filling the fields required by some nonlinear solvers.

At some point this task might be delegated to the preprocessor.

(cherry picked from commit ca230d0902eadeb89ab705cf66790ffa2f135192)

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);