Merge branch 'final_testing' of https://github.com/Neuroblox/Neuroblo…

…x.jl into final_testing

src/Neuroblox.jl

@@ -181,5 +181,6 @@ export create_adjacency_edges!
 export get_namespaced_sys, nameof
 export run_experiment!, run_trial!
 export addnontunableparams, get_hemodynamic_observers
+export get_weights
 
 end

src/blox/blox_utilities.jl

@@ -235,6 +235,22 @@ function get_event_time(kwargs, name_blox1, name_blox2)
     end
 end
 
+function get_weights(agent::Agent, blox_out, blox_in)
+    ps = parameters(agent.odesystem)
+    pv = agent.problem.p
+    map_idxs = Int.(ModelingToolkit.varmap_to_vars([ps[i] => i for i in eachindex(ps)], ps))
+
+    name_out = String(namespaced_nameof(blox_out))
+    name_in = String(namespaced_nameof(blox_in))
+
+    idxs_weight = findall(ps) do p
+        n = String(Symbol(p))
+        occursin("w_", n) && occursin(name_out, n) && occursin(name_in, n)
+    end
+
+    return pv[map_idxs[idxs_weight]]
+end
+
 function find_spikes(x::AbstractVector{T}; minprom=zero(T), maxprom=nothing, minheight=zero(T), maxheight=nothing) where {T}
     spikes, _ = argmaxima(x)
     peakproms!(spikes, x; minprom, maxheight)