Add the stats option in numerical simulation to control the display of simulation time

numerical_simulation.pl

@@ -118,7 +118,7 @@ filter(only_extrema).
 :- initial(option(stochastic_conversion: 100)).
 :- initial(option(stochastic_thresholding: 1000)).
 :- initial(option(filter: no_filter)).
-
+:- initial(option(stats: no)).
 
 
 :- devdoc('\\section{Commands}').
@@ -126,6 +126,7 @@ filter(only_extrema).
 numerical_simulation :-
   biocham_command,
   option(time, time, Time, 'time horizon of the numerical integration'),
+  option(stats, yesno, Stats, 'display computation time'),
   option(method, method, Method, 'method for the numerical solver'),
   option(
     error_epsilon_absolute, number, _ErrorEpsilonAbsolute,
@@ -223,7 +224,13 @@ numerical_simulation :-
     retractall(last_method_rsbk)
   ),
   statistics(runtime,[_,T]),
-  format("Simulation time: ~pms~n",[T]).
+  (
+    Stats = yes
+  ->
+    format("Simulation time: ~pms~n",[T])
+  ;
+    true
+  ).
 
 
 continue(Time) :-