delete comments

Author: vyudu

src/structural_transformation/utils.jl

@@ -271,16 +271,13 @@ end
 
 function find_solvables!(state::TearingState; kwargs...)
     @assert state.structure.solvable_graph === nothing
-    println("in find_solvables")
-    @show eqs
     eqs = equations(state)
     graph = state.structure.graph
     state.structure.solvable_graph = BipartiteGraph(nsrcs(graph), ndsts(graph))
     to_rm = Int[]
     for ieq in 1:length(eqs)
         find_eq_solvables!(state, ieq, to_rm; kwargs...)
     end
-    @show eqs
     return nothing
 end
 

src/systems/systems.jl

@@ -41,10 +41,10 @@ function structural_simplify(
     end
     if newsys isa DiscreteSystem &&
        any(eq -> symbolic_type(eq.lhs) == NotSymbolic(), equations(newsys))
-        # error("""
-        #     Encountered algebraic equations when simplifying discrete system. This is \
-        #     not yet supported.
-        # """)
+        error("""
+            Encountered algebraic equations when simplifying discrete system. This is \
+            not yet supported.
+        """)
     end
     for pass in additional_passes
         newsys = pass(newsys)

test/structural_transformation/utils.jl

@@ -162,19 +162,3 @@ end
     structural_simplify(sys; additional_passes = [pass])
     @test value[] == 1
 end
-
-@testset "Shift simplification" begin
-    @variables x(t) y(t) z(t)
-    @parameters a b c
-    
-    # Expand shifts
-    @test isequal(ST.expand_shifts(Shift(t, -3)(x)), Shift(t, -1)(Shift(t, -1)(Shift(t, -1)(x))))
-    expr = a * Shift(t, -2)(x) + Shift(t, 2)(y) + b
-    @test isequal(ST.expand_shifts(expr), 
-                    a * Shift(t, -1)(Shift(t, -1)(x)) + Shift(t, 1)(Shift(t, 1)(y)) + b)
-    @test isequal(ST.expand_shifts(Shift(t, 2)(Shift(t, 1)(a))), a)
-
-
-    # Distribute shifts
-
-end