You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
228 lines
15 KiB
228 lines
15 KiB
#include "gtest/gtest.h"
|
|
#include "storm-config.h"
|
|
#include "src/settings/SettingMemento.h"
|
|
#include "src/parser/PrismParser.h"
|
|
#include "src/parser/FormulaParser.h"
|
|
#include "src/logic/Formulas.h"
|
|
#include "src/builder/ExplicitPrismModelBuilder.h"
|
|
|
|
#include "src/utility/solver.h"
|
|
#include "src/modelchecker/csl/SparseCtmcCslModelChecker.h"
|
|
#include "src/modelchecker/results/ExplicitQuantitativeCheckResult.h"
|
|
|
|
#include "src/settings/SettingsManager.h"
|
|
|
|
TEST(SparseCtmcCslModelCheckerTest, Cluster) {
|
|
// Set the PRISM compatibility mode temporarily. It is set to its old value once the returned object is destructed.
|
|
std::unique_ptr<storm::settings::SettingMemento> enablePrismCompatibility = storm::settings::mutableGeneralSettings().overridePrismCompatibilityMode(true);
|
|
|
|
// Parse the model description.
|
|
storm::prism::Program program = storm::parser::PrismParser::parse(STORM_CPP_TESTS_BASE_PATH "/functional/builder/cluster2.sm");
|
|
storm::parser::FormulaParser formulaParser(program.getManager().getSharedPointer());
|
|
std::shared_ptr<storm::logic::Formula> formula(nullptr);
|
|
|
|
// Build the model.
|
|
typename storm::builder::ExplicitPrismModelBuilder<double>::Options options;
|
|
options.buildRewards = true;
|
|
options.rewardModelName = "num_repairs";
|
|
std::shared_ptr<storm::models::sparse::Model<double>> model = storm::builder::ExplicitPrismModelBuilder<double>::translateProgram(program, options);
|
|
ASSERT_EQ(storm::models::ModelType::Ctmc, model->getType());
|
|
std::shared_ptr<storm::models::sparse::Ctmc<double>> ctmc = model->as<storm::models::sparse::Ctmc<double>>();
|
|
uint_fast64_t initialState = *ctmc->getInitialStates().begin();
|
|
|
|
// Create model checker.
|
|
storm::modelchecker::SparseCtmcCslModelChecker<double> modelchecker(*ctmc, std::unique_ptr<storm::utility::solver::LinearEquationSolverFactory<double>>(new storm::utility::solver::NativeLinearEquationSolverFactory<double>()));
|
|
|
|
// Start checking properties.
|
|
formula = formulaParser.parseFromString("P=? [ F<=100 !\"minimum\"]");
|
|
std::unique_ptr<storm::modelchecker::CheckResult> checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult1 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(5.5461254704419085E-5, quantitativeCheckResult1[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("P=? [ F[100,100] !\"minimum\"]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult2 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(2.3397873548343415E-6, quantitativeCheckResult2[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("P=? [ \"minimum\" U<=10 \"premium\"]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult3 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(1, quantitativeCheckResult3[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("P=? [ !\"minimum\" U[1,inf] \"minimum\"]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult4 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(0, quantitativeCheckResult4[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("R=? [C<=100]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult5 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(0.8602815057967503, quantitativeCheckResult5[initialState], storm::settings::generalSettings().getPrecision());
|
|
}
|
|
|
|
TEST(SparseCtmcCslModelCheckerTest, Embedded) {
|
|
// Set the PRISM compatibility mode temporarily. It is set to its old value once the returned object is destructed.
|
|
std::unique_ptr<storm::settings::SettingMemento> enablePrismCompatibility = storm::settings::mutableGeneralSettings().overridePrismCompatibilityMode(true);
|
|
|
|
// Parse the model description.
|
|
storm::prism::Program program = storm::parser::PrismParser::parse(STORM_CPP_TESTS_BASE_PATH "/functional/builder/embedded2.sm");
|
|
storm::parser::FormulaParser formulaParser(program.getManager().getSharedPointer());
|
|
std::shared_ptr<storm::logic::Formula> formula(nullptr);
|
|
|
|
// Build the model.
|
|
typename storm::builder::ExplicitPrismModelBuilder<double>::Options options;
|
|
options.buildRewards = true;
|
|
options.rewardModelName = "up";
|
|
std::shared_ptr<storm::models::sparse::Model<double>> model = storm::builder::ExplicitPrismModelBuilder<double>::translateProgram(program, options);
|
|
ASSERT_EQ(storm::models::ModelType::Ctmc, model->getType());
|
|
std::shared_ptr<storm::models::sparse::Ctmc<double>> ctmc = model->as<storm::models::sparse::Ctmc<double>>();
|
|
uint_fast64_t initialState = *ctmc->getInitialStates().begin();
|
|
|
|
// Create model checker.
|
|
storm::modelchecker::SparseCtmcCslModelChecker<double> modelchecker(*ctmc, std::unique_ptr<storm::utility::solver::LinearEquationSolverFactory<double>>(new storm::utility::solver::NativeLinearEquationSolverFactory<double>()));
|
|
|
|
// Start checking properties.
|
|
formula = formulaParser.parseFromString("P=? [ F<=10000 \"down\"]");
|
|
std::unique_ptr<storm::modelchecker::CheckResult> checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult1 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(0.0019216435246119591, quantitativeCheckResult1[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("P=? [ !\"down\" U<=10000 \"fail_actuators\"]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult2 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(3.7079151806696567E-6, quantitativeCheckResult2[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("P=? [ !\"down\" U<=10000 \"fail_io\"]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult3 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(0.001556839327673734, quantitativeCheckResult3[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("P=? [ !\"down\" U<=10000 \"fail_sensors\"]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult4 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(4.429620626755424E-5, quantitativeCheckResult4[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("R=? [C<=10000]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult5 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(2.7720429852369946, quantitativeCheckResult5[initialState], storm::settings::generalSettings().getPrecision());
|
|
}
|
|
|
|
TEST(SparseCtmcCslModelCheckerTest, Polling) {
|
|
// Set the PRISM compatibility mode temporarily. It is set to its old value once the returned object is destructed.
|
|
std::unique_ptr<storm::settings::SettingMemento> enablePrismCompatibility = storm::settings::mutableGeneralSettings().overridePrismCompatibilityMode(true);
|
|
|
|
// Parse the model description.
|
|
storm::prism::Program program = storm::parser::PrismParser::parse(STORM_CPP_TESTS_BASE_PATH "/functional/builder/polling2.sm");
|
|
storm::parser::FormulaParser formulaParser(program.getManager().getSharedPointer());
|
|
std::shared_ptr<storm::logic::Formula> formula(nullptr);
|
|
|
|
// Build the model.
|
|
std::shared_ptr<storm::models::sparse::Model<double>> model = storm::builder::ExplicitPrismModelBuilder<double>::translateProgram(program);
|
|
ASSERT_EQ(storm::models::ModelType::Ctmc, model->getType());
|
|
std::shared_ptr<storm::models::sparse::Ctmc<double>> ctmc = model->as<storm::models::sparse::Ctmc<double>>();
|
|
uint_fast64_t initialState = *ctmc->getInitialStates().begin();
|
|
|
|
// Create model checker.
|
|
storm::modelchecker::SparseCtmcCslModelChecker<double> modelchecker(*ctmc, std::unique_ptr<storm::utility::solver::LinearEquationSolverFactory<double>>(new storm::utility::solver::NativeLinearEquationSolverFactory<double>()));
|
|
|
|
// Start checking properties.
|
|
formula = formulaParser.parseFromString("P=?[ F<=10 \"target\"]");
|
|
std::unique_ptr<storm::modelchecker::CheckResult> checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult1 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(1, quantitativeCheckResult1[initialState], storm::settings::generalSettings().getPrecision());
|
|
}
|
|
|
|
TEST(SparseCtmcCslModelCheckerTest, Fms) {
|
|
// Set the PRISM compatibility mode temporarily. It is set to its old value once the returned object is destructed.
|
|
std::unique_ptr<storm::settings::SettingMemento> enablePrismCompatibility = storm::settings::mutableGeneralSettings().overridePrismCompatibilityMode(true);
|
|
|
|
// No properties to check at this point.
|
|
}
|
|
|
|
TEST(SparseCtmcCslModelCheckerTest, Tandem) {
|
|
// Set the PRISM compatibility mode temporarily. It is set to its old value once the returned object is destructed.
|
|
std::unique_ptr<storm::settings::SettingMemento> enablePrismCompatibility = storm::settings::mutableGeneralSettings().overridePrismCompatibilityMode(true);
|
|
|
|
// Parse the model description.
|
|
storm::prism::Program program = storm::parser::PrismParser::parse(STORM_CPP_TESTS_BASE_PATH "/functional/builder/tandem5.sm");
|
|
storm::parser::FormulaParser formulaParser(program.getManager().getSharedPointer());
|
|
std::shared_ptr<storm::logic::Formula> formula(nullptr);
|
|
|
|
// Build the model with the customers reward structure.
|
|
typename storm::builder::ExplicitPrismModelBuilder<double>::Options options;
|
|
options.buildRewards = true;
|
|
options.rewardModelName = "customers";
|
|
std::shared_ptr<storm::models::sparse::Model<double>> model = storm::builder::ExplicitPrismModelBuilder<double>::translateProgram(program, options);
|
|
ASSERT_EQ(storm::models::ModelType::Ctmc, model->getType());
|
|
std::shared_ptr<storm::models::sparse::Ctmc<double>> ctmc = model->as<storm::models::sparse::Ctmc<double>>();
|
|
uint_fast64_t initialState = *ctmc->getInitialStates().begin();
|
|
|
|
// Create model checker.
|
|
storm::modelchecker::SparseCtmcCslModelChecker<double> modelchecker(*ctmc, std::unique_ptr<storm::utility::solver::LinearEquationSolverFactory<double>>(new storm::utility::solver::NativeLinearEquationSolverFactory<double>()));
|
|
|
|
// Start checking properties.
|
|
formula = formulaParser.parseFromString("P=? [ F<=10 \"network_full\" ]");
|
|
std::unique_ptr<storm::modelchecker::CheckResult> checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult1 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(0.015446370562428037, quantitativeCheckResult1[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("P=? [ F<=10 \"first_queue_full\" ]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult2 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(0.999999837225515, quantitativeCheckResult2[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("P=? [\"second_queue_full\" U<=1 !\"second_queue_full\"]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult3 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(1, quantitativeCheckResult3[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("R=? [I=10]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult4 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(5.679243850315877, quantitativeCheckResult4[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("R=? [C<=10]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult5 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(55.44792186036232, quantitativeCheckResult5[initialState], storm::settings::generalSettings().getPrecision());
|
|
|
|
formula = formulaParser.parseFromString("R=? [F \"first_queue_full\"&\"second_queue_full\"]");
|
|
checkResult = modelchecker.check(*formula);
|
|
|
|
ASSERT_TRUE(checkResult->isExplicitQuantitativeCheckResult());
|
|
storm::modelchecker::ExplicitQuantitativeCheckResult<double> quantitativeCheckResult6 = checkResult->asExplicitQuantitativeCheckResult<double>();
|
|
EXPECT_NEAR(262.78584491454814, quantitativeCheckResult6[initialState], storm::settings::generalSettings().getPrecision());
|
|
}
|