Commit 04d44c23 authored by lindnemi's avatar lindnemi
Browse files

removing redundant computations

parent ea0ede9e
......@@ -71,9 +71,9 @@ for sig in [10,5,2,1.0,0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2,0.1,0.05,0]
@inbounds for e in e_s # edges for which vertex is source
[1] += e[1]
end
@inbounds for e in e_d
[1] += e[2] # edges for which vertex is destination
end
#@inbounds for e in e_d
# dξ[1] += e[2] # edges for which vertex is destination
#end
nothing
end
......@@ -81,11 +81,11 @@ for sig in [10,5,2,1.0,0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2,0.1,0.05,0]
# The coupling is not symmetric wrt. change of source and destination
# Hence we compute the flow in both directions
e[1] = v_d[3] * v_d[4] * (v_d[1] - v_s[1]) * p[1]
e[2] = v_s[3] * v_s[4] * (v_s[1] - v_d[1]) * p[2]
#e[2] = v_s[3] * v_s[4] * (v_s[1] - v_d[1]) * p[2]
nothing
end
electricaledge = StaticEdge(f! = electrical_edge!, dim = 2)
electricaledge = StaticEdge(f! = electrical_edge!, dim = 1)
odevertex = ODEVertex(f! = fhn_stp!, dim = 4, sym = [:u, :v, :x, :y])
fhn_network! = network_dynamics(odevertex, electricaledge, g);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment