Commit ba62ac27 authored by Lukas Riedel's avatar Lukas Riedel

Fix bug in call to 'adapt_grid

For the correct query of the boundary conditions, the time argument must
be 'time+dt', not the new time.
parent 517c5a39
......@@ -160,20 +160,20 @@ bool RichardsSimulation<Traits>::compute_time_step ()
template<typename Traits>
void RichardsSimulation<Traits>::run ()
{
auto time = controller->getTime();
const auto time = controller->getTime();
const RF dt = controller->getDT();
write_data(time);
while(controller->doStep()) {
if(!compute_time_step()){
continue;
}
time = controller->getTime();
write_data(time);
write_data(controller->getTime());
if(controller->doStep()
&& adaptivity->adapt_grid(*grid, gv, *gfs, *param, *fboundary,
time, *u))
time+dt, *u))
{
operator_setup(); // reset operators if grid changes
}
......
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