From a86bd1c5856302ef052a3e828b31e01de74338e4 Mon Sep 17 00:00:00 2001 From: hnil Date: Tue, 1 Nov 2022 21:39:55 +0100 Subject: [PATCH] removed assert in cprw for strictly positive diagonal well coupling --- opm/simulators/wells/MultisegmentWell_impl.hpp | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/opm/simulators/wells/MultisegmentWell_impl.hpp b/opm/simulators/wells/MultisegmentWell_impl.hpp index ac69f0ff9..76ce35d7c 100644 --- a/opm/simulators/wells/MultisegmentWell_impl.hpp +++ b/opm/simulators/wells/MultisegmentWell_impl.hpp @@ -815,7 +815,11 @@ namespace Opm } } - assert(diag_ell > 0.0); + if(not(diag_ell > 0.0)){ + std::cout << "Diagonal element for cprw on " + << this->name() + << " is " << diag_ell << std::endl; + } jacobian[welldof_ind][welldof_ind] = diag_ell; }else{ jacobian[welldof_ind][welldof_ind] = 1.0; // maybe we could have used diag_ell if calculated