I don't remember

models.pl

@@ -334,7 +334,8 @@ check_model :-
       format("The model is NOT wellformed.~n", [])
     ),
     (
-      maplist(test_all_reactions(catalyst, test_strictness), List_Molecule)
+      maplist(test_all_reactions(catalyst, test_strictness), List_Molecule),
+      maplist(test_all_reactions(reactant, test_strictness), List_Molecule)
     ->
       format("The model is strict.~n", [])
     ;

models.plt

@@ -63,7 +63,8 @@ test('test_wf_other', []) :-
 test('test_strictness', []):-
   with_output_to(atom(_), (
     test_strictness(a, a*b+a^2.5/c),
-    \+(test_strictness(mapk, j0*mk/(1+(mapk/j1)^j2)/(j1+mk)))
+    \+(test_strictness(mapk, j0*mk/(1+(mapk/j1)^j2)/(j1+mk))),
+    \+(test_strictness(a, k1*a - k2*b))
   )).
 
 :- end_tests(models).