Some cleaning.

src/nullrange.jl

@@ -213,7 +213,7 @@ function grange(sys::DescriptorStateSpace{T}; zeros = "none", inner::Bool = fals
       # move unstable zeros to stable positions and compress the reduced system to one
       # with full row rank transfer function matrix
       if disc
-         X, _,F, = gared(A,E,B,D'*D,C'*C,C'*D); 
+         # X, _,F, = gared(A,E,B,D'*D,C'*C,C'*D); 
          X, _,F, = try 
             gared(A,E,B,D'*D,C'*C,C'*D) 
          catch err

src/types/DescriptorStateSpace.jl

@@ -92,23 +92,6 @@ index2range(ind::Colon) = ind
 function Base.lastindex(sys::DST, dim::Int) where DST <: DescriptorStateSpace
     lastindex(sys.D,dim)
 end
-# Base.getindex(sys::Vector{DST}, ind) where DST <: DescriptorStateSpace = getindex(sys,ind)
-# function Base.getindex(sys::Vector{DST}, inds1, inds2) where DST <: DescriptorStateSpace
-#     println("inds1 = $inds1 inds2 = $inds2")
-#     size(inds2, 1) != 1 &&
-#         error("Must specify 2 indices to index descriptor state-space model")
-#     # rows, cols = index2range(inds...) 
-#     # rows, cols = index2range(inds...) 
-#     rows = inds1; cols = inds2
-#     println("rows = $rows cols = $cols")
-#     N = length(sys)
-#     sysR = similar(sys,N)
-#     for i = 1:N
-#         #sysR[i] = DescriptorStateSpace{eltype(sys[i])}(copy(sys[i].A), copy(sys[i].E), sys[i].B[:, cols], sys[i].C[rows, :], sys[i].D[rows, cols], sys[i].Ts)
-#         sysR[i] = sys[i][rows, cols]
-#     end  
-#     return sysR
-# end
 
 # Basic Operations
 function ==(sys1::DST1, sys2::DST2) where {DST1<:DescriptorStateSpace, DST2<:DescriptorStateSpace}