fix2

src/component_functions.jl

@@ -964,7 +964,7 @@ function _fill_defaults(T, @nospecialize(kwargs))
     #### External Inputs
     ####
     if haskey(dict, :extin)
-        @assert dict[:extin] isa Vector{<:SymbolicIndex}
+        @assert dict[:extin] isa Union{Nothing, Vector{<:SymbolicIndex}}
     else
         dict[:extin] = nothing
     end

src/external_inputs.jl

@@ -13,11 +13,13 @@ function ExtMap(im::IndexManager)
     isempty(map) && error("There are no external inputs.")
 
     for (vi, vm) in pairs(im.vertexm)
+        has_external_input(vm) || continue
         for (i, si) in pairs(vm.extin)
             map[im.v_ext[vi][i]] = _symidx_to_extidx(im, si)
         end
     end
     for (ei, em) in pairs(im.edgem)
+        has_external_input(em) || continue
         for (i, si) in pairs(em.extin)
             map[im.e_ext[ei][i]] = _symidx_to_extidx(im, si)
         end

src/symbolicindexing.jl

@@ -507,14 +507,14 @@ function _get_observed_f(nw::Network, cf::EdgeModel, eidx)
     ur    = nw.im.e_data[eidx]
     esrcr = nw.im.v_out[nw.im.edgevec[eidx].src]
     edstr = nw.im.v_out[nw.im.edgevec[eidx].dst]
-    extr = nw.im.v_out[eidx]
+    extr = nw.im.e_out[eidx]
     pr   =  nw.im.e_para[eidx]
     ret = Vector{Float64}(undef, N)
 
     @closure (u, outbuf, aggbuf, extbuf, p, t) -> begin
         ins = (view(outbuf, esrcr), view(outbuf, edstr))
         if has_external_input(cf)
-            (ins..., view(aggbuf, aggr), view(extbuf, extr))
+            (ins..., view(extbuf, extr))
         end
         cf.obsf(ret, view(u, ur), ins..., view(p, pr), t)
         ret