Fix non-static zero entry

Author: janbruedigam

src/system/entry.jl

@@ -19,7 +19,10 @@ function Base.show(io::IO, mime::MIME{Symbol("text/plain")}, entry::Entry)
 end
 
 
-Base.zero(::Entry{ET}) where ET = Entry(zero(ET), false)
+function Base.zero(::Entry{ET}) where ET
+    dims = [ET.parameters[1].parameters...]
+    return Entry{ET.parameters[2]}(dims...)
+end
 Base.zero(::Type{Entry}) = Entry(nothing, false)
 
 Base.:+(E1::Entry, E2::Entry) = Entry(E1.value + E2.value, false)

test/algebra.jl

@@ -6,7 +6,7 @@ include("adjacency_matrix.jl")
 for i=1:5
     r1 = randn()
     r2 = randn()
-    system = System{Float64}(A, rand(0:3, size(A)[1]))
+    system = System{Float64}(A, rand(0:20, size(A)[1]))
     initialize!(system)
 
     F = full_matrix(system)