From 48e3a01b4cbb4f3cf4eac4afb6e54d84556816a8 Mon Sep 17 00:00:00 2001 From: agchesebro <76024790+agchesebro@users.noreply.github.com> Date: Thu, 21 Nov 2024 09:18:27 -0500 Subject: [PATCH 1/3] Fix warning --- src/blox/neural_mass.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/blox/neural_mass.jl b/src/blox/neural_mass.jl index e98a0605..06ed6bcf 100644 --- a/src/blox/neural_mass.jl +++ b/src/blox/neural_mass.jl @@ -528,7 +528,7 @@ struct PYR_Izh <: NeuralMassBlox ω=0.0) p = paramscoping(Δ=Δ, α=α, gₛ=gₛ, η̄=η̄, I_ext=I_ext, eᵣ=eᵣ, a=a, b=b, wⱼ=wⱼ, sⱼ=sⱼ, κ=κ) Δ, α, gₛ, η̄, I_ext, eᵣ, a, b, wⱼ, sⱼ, κ = p - sts = @variables r(t)=0.0 v(t)=0.0 w(t)=0.0 s(t)=0.0 [output=true] jcn(t)=0.0 [input=true] + sts = @variables r(t)=0.0 v(t)=0.0 w(t)=0.0 s(t)=0.0 [output=true] jcn(t) [input=true] eqs = [ D(r) ~ Δ/π + 2*r*v - (α+gₛ*s)*r, D(v) ~ v^2 - α*v - w + η̄ + I_ext*sin(ω*t) + gₛ*s*κ*(eᵣ - v) + jcn - (π*r)^2, D(w) ~ a*(b*v - w) + wⱼ*r, From 6ebc513795d945be6d9dd4adf104468e732f8607 Mon Sep 17 00:00:00 2001 From: Mason Protter Date: Thu, 21 Nov 2024 16:16:33 +0100 Subject: [PATCH 2/3] set PING jcn to zero --- src/blox/ping_neuron_examples.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/blox/ping_neuron_examples.jl b/src/blox/ping_neuron_examples.jl index 6aa0277a..330fe296 100644 --- a/src/blox/ping_neuron_examples.jl +++ b/src/blox/ping_neuron_examples.jl @@ -56,7 +56,7 @@ struct PINGNeuronExci <: AbstractPINGNeuron τ_D=2.0) p = paramscoping(C=C, g_Na=g_Na, V_Na=V_Na, g_K=g_K, V_K=V_K, g_L=g_L, V_L=V_L, I_ext=I_ext, τ_R=τ_R, τ_D=τ_D) C, g_Na, V_Na, g_K, V_K, g_L, V_L, I_ext, τ_R, τ_D = p - sts = @variables V(t)=0.0 [output=true] n(t)=0.0 h(t)=0.0 s(t)=0.0 jcn(t)=0.0 [input=true] + sts = @variables V(t)=0.0 [output=true] n(t)=0.0 h(t)=0.0 s(t)=0.0 jcn(t) [input=true] a_m(v) = 0.32*(v+54.0)/(1.0 - exp(-(v+54.0)/4.0)) b_m(v) = 0.28*(v+27.0)/(exp((v+27.0)/5.0) - 1.0) @@ -131,7 +131,7 @@ struct PINGNeuronInhib <: AbstractPINGNeuron τ_D=10.0) p = paramscoping(C=C, g_Na=g_Na, V_Na=V_Na, g_K=g_K, V_K=V_K, g_L=g_L, V_L=V_L, I_ext=I_ext, τ_R=τ_R, τ_D=τ_D) C, g_Na, V_Na, g_K, V_K, g_L, V_L, I_ext, τ_R, τ_D = p - sts = @variables V(t)=0.0 [output=true] n(t)=0.0 h(t)=0.0 s(t)=0.0 jcn(t)=0.0 [input=true] + sts = @variables V(t)=0.0 [output=true] n(t)=0.0 h(t)=0.0 s(t)=0.0 jcn(t) [input=true] a_m(v) = 0.1*(v+35.0)/(1.0 - exp(-(v+35.0)/10.0)) b_m(v) = 4*exp(-(v+60.0)/18.0) @@ -191,4 +191,4 @@ function (bc::BloxConnector)( eq = sys_in.jcn ~ w*s*(V_I-v_in) accumulate_equation!(bc, eq) -end \ No newline at end of file +end From 9f927299357d665140d95d88e49397310f793503 Mon Sep 17 00:00:00 2001 From: agchesebro <76024790+agchesebro@users.noreply.github.com> Date: Thu, 21 Nov 2024 10:25:43 -0500 Subject: [PATCH 3/3] Removing jcn(t)=0 from another neural mass --- src/blox/neural_mass.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/blox/neural_mass.jl b/src/blox/neural_mass.jl index 06ed6bcf..f879a06b 100644 --- a/src/blox/neural_mass.jl +++ b/src/blox/neural_mass.jl @@ -558,7 +558,7 @@ struct QIF_PING_NGNMM <: NeuralMassBlox A=0.0) p = paramscoping(Δ=Δ, τₘ=τₘ, H=H, I_ext=I_ext, J_internal=J_internal) Δ, τₘ, H, I_ext, J_internal = p - sts = @variables r(t)=0.0 [output=true] v(t)=0.0 jcn(t)=0.0 [input=true] + sts = @variables r(t)=0.0 [output=true] v(t)=0.0 jcn(t) [input=true] @brownian ξ eqs = [D(r) ~ Δ/(π*τₘ^2) + 2*r*v/τₘ, D(v) ~ (v^2 + H + I_ext*sin(ω*t))/τₘ - τₘ*(π*r)^2 + J_internal*r + A*ξ + jcn]