Show time for rate_independence

rate_independence.pl

@@ -24,6 +24,7 @@
 test_rate_independence :-
   biocham_command,
   doc('Tests graphical sufficient conditions for rate independence of the current model for all output species (assuming well-formed kinetics).'),
+  statistics(walltime, _),
   reaction_graph,
   get_current_graph(Id),
   (
@@ -37,7 +38,10 @@ test_rate_independence :-
     nb_getval(rate_dependence_reason, Failed),
     Failed =.. [Stage, Place],
     format("not ~w-free (on ~w).~n", [Stage, Place])
-  ).
+  ),
+  statistics(walltime, [_, Time]),
+  TTime is Time / 1000,
+  format("Time: ~p s~n", [TTime]).
 
 
 				% to generalize to influence models too
@@ -445,6 +449,7 @@ add_merge_reaction(M, Reaction1, Reaction2):-
 test_rate_independence_invariants :-
   biocham_command,
   doc('Tests invariant sufficient conditions for rate independence of the current model for all output species (assuming well-formed kinetics).'),
+  statistics(walltime, _),
   retractall(p_invariant(_)),
   %% retractall(sur_invariant(_)),
   retractall(output(_)),
@@ -489,8 +494,11 @@ test_rate_independence_invariants :-
       true
     )
   ;
-    write('Undecided: no outputs!')
-  ).
+    writeln('Undecided: no outputs!')
+  ),
+  statistics(walltime, [_, Time]),
+  TTime is Time / 1000,
+  format("Time: ~p s~n", [TTime]).
 
 
 %! test_rate_independence_invariants(+Output, -Id) is det.
@@ -503,9 +511,9 @@ test_rate_independence_invariants(Output, _Id) :-
     p_invariant(P),
     member(Output, P)
   ->
-    format("The output ~w is rate independent.", [Output])
+    format("The output ~w is rate independent.~n", [Output])
   ;
-    format("Undecided, ~w not persistent or not covered by any P-invariant.", [Output]),
+    format("Undecided, ~w not persistent or not covered by any P-invariant.~n", [Output]),
     fail
   ).