From d880650f6276f1440e7509079cbe053343ddec2d Mon Sep 17 00:00:00 2001
From: Mahmoud Bentriou <mahmoud.bentriou@centralesupelec.fr>
Date: Sun, 10 Jan 2021 13:24:16 +0100
Subject: [PATCH] implementation of probability_var_value_lha
---
core/MarkovProcesses.jl | 2 +-
core/model.jl | 15 +++++++++++++++
2 files changed, 16 insertions(+), 1 deletion(-)
diff --git a/core/MarkovProcesses.jl b/core/MarkovProcesses.jl
index d971b80..f6e5d04 100644
--- a/core/MarkovProcesses.jl
+++ b/core/MarkovProcesses.jl
@@ -28,7 +28,7 @@ export get_index, get_value, length_var, isaccepted
# Model related methods
export simulate, volatile_simulate
-export distribute_mean_value_lha, mean_value_lha, distribute_prob_accept_lha
+export distribute_mean_value_lha, mean_value_lha, distribute_prob_accept_lha, probability_var_value_lha
export set_param!, set_x0!, set_time_bound!, set_observed_var!, observe_all!
export get_param, get_x0, getproperty, get_proba_model, get_observed_var
export isbounded, isaccepted, check_consistency
diff --git a/core/model.jl b/core/model.jl
index 7857f93..431e528 100644
--- a/core/model.jl
+++ b/core/model.jl
@@ -382,6 +382,21 @@ function mean_value_lha(sm::SynchronizedModel, sym_var::VariableAutomaton, nbr_s
end
return sum_val / nbr_sim
end
+"""
+ `distribute_var_value_lha(sm::SynchronizedModel, nbr_sim::Int, value = 0, sym_var = :d)
+
+Compute the probability that the variable `sym_var` is equal to `value`
+of a LHA over `nbr_sim` simulations of the model.
+"""
+function probability_var_value_lha(sm::SynchronizedModel, nbr_sim::Int;
+ value::Float64 = 0.0, sym_var::VariableAutomaton = :d)
+ sum_val = @distributed (+) for i = 1:nbr_sim
+ S = volatile_simulate(sm)
+ S[sym_var] == value
+ end
+ return sum_val / nbr_sim
+end
+
function distribute_prob_accept_lha(sm::SynchronizedModel, nbr_sim::Int)
sum_val = @distributed (+) for i = 1:nbr_sim
--
GitLab