OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-22 07:23:27 -06:00
Code Issues Packages Projects Releases Wiki Activity

Added testcase for TUNING/TSINIT and NEXTSTEP - depends on https://github.com/OPM/opm-tests/pull/1028

Browse Source
This commit is contained in:
Vegard Kippe 2023-08-25 13:59:11 +02:00
parent f5c31588cb
commit 4ff08ffd2d
3 changed files with 67 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
CMakeLists.txt
View File

@ -335,6 +335,13 @@ opm_add_test(test_tuning_xxxMBE
${Boost_UNIT_TEST_FRAMEWORK_LIBRARY}
ONLY_COMPILE)
opm_add_test(test_tuning_tsinit_nextstep
SOURCES
tests/test_tuning_TSINIT_NEXTSTEP.cpp
LIBRARIES
${Boost_UNIT_TEST_FRAMEWORK_LIBRARY} opmcommon
ONLY_COMPILE)
if (HAVE_OPM_TESTS)
include (${CMAKE_CURRENT_SOURCE_DIR}/compareECLFiles.cmake)
endif()

8
compareECLFiles.cmake
View File

@ -312,6 +312,14 @@ add_test_runSimulator(CASENAME notuning_xxxmbe
set_tests_properties(runSimulator/notuning_xxxmbe PROPERTIES WILL_FAIL TRUE)
add_test_runSimulator(CASENAME tuning_tsinit_nextstep
FILENAME 02_TUNING_TSINIT_NEXTSTEP
SIMULATOR flow
DIR tuning
TEST_ARGS --enable-tuning=true
POST_COMMAND $<TARGET_FILE:test_tuning_tsinit_nextstep>)
include (${CMAKE_CURRENT_SOURCE_DIR}/regressionTests.cmake)
include (${CMAKE_CURRENT_SOURCE_DIR}/restartTests.cmake)

52
tests/test_tuning_TSINIT_NEXTSTEP.cpp Normal file
View File

@ -0,0 +1,52 @@
#include <opm/io/eclipse/ESmry.hpp>
#include <string>
#include <cmath>
#include <iostream>
#define BOOST_TEST_MODULE TestTuningTSINIT
#include <boost/test/unit_test.hpp>
constexpr double TOLERANCE = 1.0e-10;
inline bool is_close(double a, double b, double tol = TOLERANCE) {
return std::abs(a-b) <= tol;
}
BOOST_AUTO_TEST_CASE(CheckTSINITAndNEXTSTEP)
{
//std::string case_name(boost::unit_test::framework::master_test_suite().argv[1]);
std::string case_name("02_TUNING_TSINIT_NEXTSTEP");
BOOST_TEST_MESSAGE("---------------------------------------------------------------------------");
BOOST_TEST_MESSAGE("Checking TSINIT and NEXTSTEP, see file " + case_name + ".DATA");
BOOST_TEST_MESSAGE("---------------------------------------------------------------------------");
Opm::EclIO::ESmry smry(case_name, false);
smry.loadData({"TIME"});
const auto& time = smry.get("TIME");
// First time step 1 day
BOOST_CHECK_CLOSE(time[0], 1.0, TOLERANCE);
for (size_t i=0; i<time.size(); ++i) {
std::cout << "######################################" << std::endl;
std::cout << time[i] << std::endl;
// Max time step 3 days initially
if (time[i] < 14.0) BOOST_CHECK(time[i+1] <= (time[i]+3.0));
// No short next step
if (is_close(time[i], 14.0)) BOOST_CHECK(time[i+1] > 15.1);
// Persistent NEXTSTEP=0.5
if (is_close(time[i], 31.0)) BOOST_CHECK_CLOSE(time[i+1], 31.5, TOLERANCE);
if (is_close(time[i], 45.0)) BOOST_CHECK_CLOSE(time[i+1], 45.5, TOLERANCE);
// Non-persistent NEXTSTEP=1.0
if (is_close(time[i], 60.0)) BOOST_CHECK_CLOSE(time[i+1], 61.0, TOLERANCE);
if (is_close(time[i], 74.0)) BOOST_CHECK(time[i+1] > 75.1);
// TSINIT=0.5
if (is_close(time[i], 91.0)) BOOST_CHECK_CLOSE(time[i+1], 91.5, TOLERANCE);
if (is_close(time[i], 105.0)) BOOST_CHECK(time[i+1] > 105.6);
}
}