From 497c81cbb10c105c499417df51b8f1208e16a6bf Mon Sep 17 00:00:00 2001 From: timmyfaraday Date: Mon, 22 Jan 2024 12:39:23 +0100 Subject: [PATCH] final commit before merge into main --- Project.toml | 10 +++++----- src/form/semi_markov_process.jl | 5 +++-- 2 files changed, 8 insertions(+), 7 deletions(-) diff --git a/Project.toml b/Project.toml index d259bc6..4284516 100644 --- a/Project.toml +++ b/Project.toml @@ -18,14 +18,14 @@ Unitful = "1986cc42-f94f-5a68-af5c-568840ba703d" [compat] AdditionalUnits = "0.1.4" -Graphs = "1.8.0" +Graphs = "1.9.0" Interpolations = "0.14.7" -Measurements = "2.9.0" +Measurements = "2.11.0" Multigraphs = "0.3.0" -OrdinaryDiffEq = "6.55.0" -QuadGK = "2.8.2" +OrdinaryDiffEq = "6.58.2" +QuadGK = "2.9.4" SpecialFunctions = "2.3.1" -Unitful = "1.17.0" +Unitful = "1.19.0" julia = "1" [extras] diff --git a/src/form/semi_markov_process.jl b/src/form/semi_markov_process.jl index 2c2f5d3..b79f0bd 100644 --- a/src/form/semi_markov_process.jl +++ b/src/form/semi_markov_process.jl @@ -106,6 +106,7 @@ function _set_p!(std::AbstractSTD, st::Int, init::Real, tol::Real, DST::Vector, _fill_p!(init, dst, t, h, p) end + set_prop!(std, st, :h, h) set_prop!(std, st, :prob, p) end function _fill_p!(init::Real, dst::AbstractDistribution, t::StepRangeLen, @@ -136,8 +137,8 @@ function solve!(std::AbstractSTD, cls::AbstractSemiMarkovProcess; t = zero(dt):dt:tsim # solve the problem - U = get_U(std,t,tol) - A = ustrip(get_A(std,t,tol)) + U = get_U(std, t, tol) + A = ustrip(get_A(std, t, tol)) H = U \ A * unit(1/dt) set_p!(std, t, H, tol)