From 6571e840a295a680a883673dc0818e5709b3abdb Mon Sep 17 00:00:00 2001
From: vyudu <vincent.duyuan@gmail.com>
Date: Thu, 11 Jul 2024 16:16:27 -0400
Subject: [PATCH] fix docstring

---
 src/network_analysis.jl | 8 +-------
 1 file changed, 1 insertion(+), 7 deletions(-)

diff --git a/src/network_analysis.jl b/src/network_analysis.jl
index e9f5af1160..abab4f66ac 100644
--- a/src/network_analysis.jl
+++ b/src/network_analysis.jl
@@ -260,17 +260,12 @@ end
 Construct a directed simple graph where nodes correspond to reaction complexes and directed
 edges to reactions converting between two complexes.
 
-Notes:
-- Requires the `incidencemat` to already be cached in `rn` by a previous call to
-  `reactioncomplexes`.
-
 For example,
 ```julia
 sir = @reaction_network SIR begin
     β, S + I --> 2I
     ν, I --> R
 end
-complexes,incidencemat = reactioncomplexes(sir)
 incidencematgraph(sir)
 ```
 """
@@ -498,8 +493,7 @@ isweaklyreversible(rn, subnets)
 function isweaklyreversible(rn::ReactionSystem, subnets)
     nps = get_networkproperties(rn)
     isempty(nps.incidencemat) && reactioncomplexes(rn)
-    imat = nps.incidencemat
-    sparseig = issparse(imat)
+    sparseig = issparse(nps.incidencemat)
 
     for subnet in subnets
         subnps = get_networkproperties(subnet)