From 37e96e211c3ebf03360fcf68c10deef5b44112d5 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Erik=20Hide=20S=C3=A6ternes?= Date: Fri, 8 Sep 2023 14:57:57 +0200 Subject: [PATCH] Fixed indexing error in PID controller --- opm/simulators/timestepping/TimeStepControl.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/opm/simulators/timestepping/TimeStepControl.cpp b/opm/simulators/timestepping/TimeStepControl.cpp index d72251b1c..33c7b3665 100644 --- a/opm/simulators/timestepping/TimeStepControl.cpp +++ b/opm/simulators/timestepping/TimeStepControl.cpp @@ -206,7 +206,7 @@ namespace Opm const double kD = 0.01 ; const double newDt = (dt * std::pow( errors_[ 1 ] / errors_[ 2 ], kP ) * std::pow( tol_ / errors_[ 2 ], kI ) * - std::pow( errors_[0]*errors_[0]/errors_[ 1 ]/errors_[ 2 ], kD )); + std::pow( errors_[1]*errors_[1]/errors_[ 0 ]/errors_[ 2 ], kD )); if( verbose_ ) OpmLog::info(fmt::format("Computed step size (pow): {} days", unit::convert::to( newDt, unit::day ))); return newDt;