fixes

Author: vyudu

src/network_analysis.jl

@@ -221,12 +221,12 @@ function fluxmat(rn::ReactionSystem, pmap::Dict = Dict(); sparse=false)
     nc = length(rcmap)
     nr = length(rates)
     mtype = eltype(rates) <: Symbolics.BasicSymbolic ? Num : eltype(rates)
-    fluxmat = if sparse
+    mat = if sparse
         fluxmat(SparseMatrixCSC{mtype, Int}, rcmap, rates)
     else
         fluxmat(Matrix{mtype}, rcmap, rates)
     end
-    mtype == Num ? Matrix{Any}(fluxmat) : fluxmat
+    mtype == Num ? Matrix{Any}(mat) : mat
 end
 
 function fluxmat(::Type{SparseMatrixCSC{T, Int}}, rcmap, rates) where T
@@ -250,8 +250,7 @@ function fluxmat(::Type{Matrix{T}}, rcmap, rates) where T
     nc = length(rcmap)
     K = zeros(T, nr, nc)
     for (i, (complex, rxs)) in enumerate(rcmap)
-        for (rx, dir) in rxs
-            dir == -1 && (K[rx, i] = rates[rx])
+        for (rx, dir) in rxs dir == -1 && (K[rx, i] = rates[rx])
         end
     end
     K

test/network_analysis/network_properties.jl

@@ -379,7 +379,7 @@ let
               ]
     @test isequal(Φ, truevec)
 
-    K = fluxmat(MAPK)
+    K = Catalyst.fluxmat(MAPK)
     # Construct matrix from incidence matrix
     mat = zeros(Num, 30, 26)
     D = incidencemat(MAPK)
@@ -396,7 +396,7 @@ let
     @test isequal(K[3, 2], MAPK.k₃)
     @test all(==(0), vcat(K[3,1], K[3,3:end]))
     @test count(k -> !isequal(k, 0), K) == length(reactions(MAPK))
-    K = fluxmat(MAPK; sparse = true)
+    K = Catalyst.fluxmat(MAPK; sparse = true)
     @test Catalyst.issparse(K)
 
     A_k = Catalyst.laplacianmat(MAPK)