-
Notifications
You must be signed in to change notification settings - Fork 8
/
incidence.jl
85 lines (70 loc) · 2.16 KB
/
incidence.jl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
"""
Maps from equation number to a Dict: variable -> Set(time_periods)
"""
struct IncidenceTable
by_eq::Dict{Int,Dict{Symbol,Set{Int}}}
by_var::Dict{Symbol,Set{Int}}
by_date::Dict{Int,Set{Symbol}}
end
IncidenceTable() = IncidenceTable(Dict(), Dict(), Dict())
function IncidenceTable(eqs::AbstractVector, skip::AbstractVector=Symbol[])
# create incidence
it = IncidenceTable()
for (i, eq) in enumerate(eqs)
visit!(it, eq, i, 0, skip)
# make sure the Dict by_eq has an entry for this equation, even if it
# is empty.
get!(it.by_eq, i, Dict{Symbol,Set{Int}}())
end
it
end
function IncidenceTable(eq::Expr, skip::AbstractVector=Symbol[])
it = IncidenceTable([eq], skip)
end
Base.getindex(it::IncidenceTable, i::Int) = it.by_date[i]
Base.getindex(it::IncidenceTable, s::Symbol) = it.by_var[s]
function visit!(it::IncidenceTable, s::Symbol, n::Int, shift::Int,
skip::AbstractVector=Symbol[])
if s in skip
return nothing
end
for_eq = get!(it.by_eq, n, Dict{Symbol,Set{Int}}())
for_sym = get!(for_eq, s, Set{Int}())
push!(for_sym, shift)
# update by_var
push!(get!(it.by_var, s, Set{Int}()), shift)
# upate by_date
push!(get!(it.by_date, shift, Set{Symbol}()), s)
nothing
end
# don't worry about anything else
function visit!(it::IncidenceTable, s::Any, n::Int, shift::Int,
skip::AbstractVector=Symbol[])
nothing
end
function visit!(it::IncidenceTable, ex::Expr, n::Int, shift::Int,
skip::AbstractVector=Symbol[])
if ex.head == :line
# this is a line number node in a quoted expression. There is nothing
# for us to do here.
return nothing
end
if is_time_shift(ex)
var = ex.args[1]
i = ex.args[2]
visit!(it, var, n, i+shift, skip)
return nothing
end
# don't visit the function name
if ex.head == :call
for an_arg in ex.args[2:end]
visit!(it, an_arg, n, shift, skip)
end
return nothing
end
# Otherwise just visit everything
for an_arg in ex.args
visit!(it, an_arg, n, shift, skip)
end
nothing
end