Commit da420aed authored by Lukas Riedel's avatar Lukas Riedel

Create FV solver instance in KnoFu backend test

parent ef36a574
......@@ -27,7 +27,9 @@
using Grid = Dune::YaspGrid<2>;
using BaseTraits = Dune::Dorie::BaseTraits<
Grid, Dune::GeometryType::BasicType::cube>;
using Traits = Dune::Dorie::RichardsSimulationTraits<BaseTraits, 1>;
template<int order>
using Traits = Dune::Dorie::RichardsSimulationTraits<BaseTraits, order>;
/// Assert that inserting and retrieving does not alter the vector
/** Draw a random state vector within given limits, store it into the model
......@@ -182,20 +184,21 @@ int main(int argc, char** argv)
spdlog::level::from_str(log_level));
log->info("Performing solution set/get test for matric head");
Dune::Dorie::KnoFuInterface<Traits> sim_1(inifile, gc, helper,
Dune::Dorie::KnoFuInterface<Traits<1>> sim_1(inifile, gc, helper,
Dune::Dorie::FilterSolutionType::matric_head);
check_set_get_solution(sim_1, -100.0, 100.0, 0.0);
log->info("Performing solution set/get test for water content");
Dune::Dorie::KnoFuInterface<Traits> sim_2(inifile, gc, helper,
Dune::Dorie::KnoFuInterface<Traits<0>> sim_2(inifile, gc, helper,
Dune::Dorie::FilterSolutionType::water_content);
check_set_get_solution(sim_2, v_min, v_max);
log->info("Testing the state position vector");
verify_position_vector(sim_1, grid);
verify_position_vector(sim_2, grid);
log->info("Performing the time step test");
Dune::Dorie::KnoFuInterface<Traits> sim_3(inifile, gc, helper,
Dune::Dorie::KnoFuInterface<Traits<1>> sim_3(inifile, gc, helper,
Dune::Dorie::FilterSolutionType::matric_head);
test_sim_time(sim_3);
......
......@@ -7,7 +7,6 @@ dimensions = 2
initialLevel = 0
gridType = rectangular
cells = 50 50
FEorder = 1
extensions = 1 1
[grid.mapping]
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment