Only allocate IncrementalCycleTracker once in tearing

```julia
julia> @time sysRed, = structural_simplify(sysEx, (inputs=(), outputs=()));
  5.951818 seconds (31.59 M allocations: 2.248 GiB, 10.44% gc time)
```

Author: YingboMa

src/structural_transformation/bipartite_tearing/modia_tearing.jl

@@ -48,9 +48,8 @@ function tearEquations!(ict::IncrementalCycleTracker, Gsolvable, es::Vector{Int}
     return ict
 end
 
-function tear_graph_block_modia!(var_eq_matching, vargraph, solvable_graph, eqs, vars,
+function tear_graph_block_modia!(var_eq_matching, ict, solvable_graph, eqs, vars,
                                  isder::F) where {F}
-    ict = IncrementalCycleTracker(vargraph; dir = :in)
     tearEquations!(ict, solvable_graph.fadjlist, eqs, vars, isder)
     for var in vars
         var_eq_matching[var] = ict.graph.matching[var]
@@ -77,6 +76,7 @@ function tear_graph_modia(structure::SystemStructure, isder::F = nothing,
     var_eq_matching = complete(maximal_matching(graph, eqfilter, varfilter, U))
     var_sccs::Vector{Union{Vector{Int}, Int}} = find_var_sccs(graph, var_eq_matching)
     vargraph = DiCMOBiGraph{true}(graph)
+    ict = IncrementalCycleTracker(vargraph; dir = :in)
 
     ieqs = Int[]
     filtered_vars = BitSet()
@@ -90,7 +90,7 @@ function tear_graph_modia(structure::SystemStructure, isder::F = nothing,
             end
             var_eq_matching[var] = unassigned
         end
-        tear_graph_block_modia!(var_eq_matching, vargraph, solvable_graph, ieqs,
+        tear_graph_block_modia!(var_eq_matching, ict, solvable_graph, ieqs,
                                 filtered_vars,
                                 isder)