From 8e4a100ad926b27b1ab58d0c67af6ffbdec5103b Mon Sep 17 00:00:00 2001 From: Jesse Perla Date: Thu, 25 Aug 2022 12:11:19 -0700 Subject: [PATCH] Bug fix for verify_steady_state (#162) --- Project.toml | 2 +- src/generate_perturbation.jl | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/Project.toml b/Project.toml index ae4e8a6..b151eff 100644 --- a/Project.toml +++ b/Project.toml @@ -1,7 +1,7 @@ name = "DifferentiableStateSpaceModels" uuid = "beacd9db-9e5e-4956-9b09-459a4b2028df" authors = ["Jesse Perla and contributors"] -version = "0.5.1" +version = "0.5.2" [deps] ChainRulesCore = "d360d2e6-b24c-11e9-a2a3-2a2ae2dbcce4" diff --git a/src/generate_perturbation.jl b/src/generate_perturbation.jl index 99797ff..b6d5860 100644 --- a/src/generate_perturbation.jl +++ b/src/generate_perturbation.jl @@ -28,7 +28,7 @@ function verify_steady_state(m, p_d, p_f; atol = 1e-8, args...) p_d_symbols = collect(Symbol.(keys(p_d))) p = order_vector_by_symbols(merge(p_d, p_f), m.mod.m.p_symbols) - w = hcat(sol.y, sol.x) # get a vector for the proposed steadystate + w = vcat(sol.y, sol.x) # get a vector for the proposed steadystate H = Vector{Float64}(undef, length(w)) # allocate it, but leave undef to make sure we can see if it goes to 0 or not m.mod.m.H̄!(H, w, p) # evaluate in place return (norm(H) < atol)