Sur-invariants back but commented

rate_independence.pl

@@ -435,6 +435,7 @@ add_merge_reaction(M, Reaction1, Reaction2):-
 
 % Dynamic predicate for storing structural information of the graph
 :- dynamic(p_invariant/1).
+%% :- dynamic(sur_invariant/1).
 :- dynamic(output/1).
 :- dynamic(siphon/1).
 
@@ -445,6 +446,7 @@ 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).'),
   retractall(p_invariant(_)),
+  %% retractall(sur_invariant(_)),
   retractall(output(_)),
   retractall(siphon(_)),
   with_output_to(atom(_), invariants:find_invar_aux(4, is_place, is_transition, '#=<')),
@@ -452,10 +454,12 @@ test_rate_independence_invariants :-
     retract(invariants:base_mol(B)),
     (
       normalize(B, BB),
+      %% assertz(sur_invariant(BB)),
+      %% debug(rate, "sur_invariant: ~w", [BB]),
       (
-        B = [Out]
+        BB = [Out]
       ->
-        debug(rate, "output: ~w", [B]),
+        debug(rate, "output: ~w", [BB]),
         assertz(output(Out))
       ;
         true