Added some edge cases in reactionsystems and networkanalysis

Signed-off-by: sivasathyaseeelan <[email protected]>

Author: sivasathyaseeelan

test/network_analysis/network_properties.jl

@@ -180,6 +180,17 @@ let
     rates = Dict(zip(parameters(rn4), k))
     @test Catalyst.iscomplexbalanced(rn4, rates) == true
 end
+
+let
+    rn = @reaction_network begin
+        (k1, k2), C1 <--> C2
+        (k3, k4), C2 <--> C3
+        (k5, k6), C3 <--> C1
+    end
+
+    incorrect_params = Dict(:k1 => 0.5)
+    @test_throws ErrorException Catalyst.iscomplexbalanced(rn, incorrect_params)
+end
 
 ### Tests Reversibility ###
 
@@ -708,6 +719,13 @@ let
     Catalyst.ratematrix(rn, rates_vec) == rate_mat
     Catalyst.ratematrix(rn, rates_tup) == rate_mat
     Catalyst.ratematrix(rn, rates_dict) == rate_mat
+
+    # Tests that throws error in rate matrix.
+    incorrect_param_dict = Dict(:k1 => 1.0)
+
+    @test_throws ErrorException Catalyst.ratematrix(rn, 123)
+    @test_throws ErrorException Catalyst.ratematrix(rn, incorrect_param_dict)
+
     @test_throws Exception Catalyst.iscomplexbalanced(rn, rates_invalid)
 end
 
@@ -739,6 +757,17 @@ let
     @test Catalyst.robustspecies(EnvZ_OmpR) == [6]
 end
 
+let
+    # Define a reaction network with bi-directional reactions
+    non_deficient_network = @reaction_network begin
+        (k1, k2), A <--> B
+        (k3, k4), B <--> C 
+    end
+
+    # Test: Check that the error is raised for networks with deficiency != 1
+    @test_throws ErrorException Catalyst.robustspecies(non_deficient_network)
+end
+
 
 ### Complex and detailed balance tests
 

test/reactionsystem_core/reactionsystem.jl

@@ -107,6 +107,9 @@ end
 let
     @named rs2 = ReactionSystem(rxs, t)
     @test Catalyst.isequivalent(rs, rs2)
+
+    # Test with a type mismatch
+    @test Catalyst.isequivalent(rs, "Not a ReactionSystem") == false
 end
 
 # Defaults test.