From 5f68b71d0fe248f34b3c734c3673986900851b1f Mon Sep 17 00:00:00 2001 From: andrthu Date: Wed, 4 Dec 2019 15:43:31 +0100 Subject: [PATCH] Remove makeOverlapInvalid --- opm/simulators/linalg/ISTLSolverEbos.hpp | 21 --------------------- 1 file changed, 21 deletions(-) diff --git a/opm/simulators/linalg/ISTLSolverEbos.hpp b/opm/simulators/linalg/ISTLSolverEbos.hpp index 74fbaaeeb..64b11c9a2 100644 --- a/opm/simulators/linalg/ISTLSolverEbos.hpp +++ b/opm/simulators/linalg/ISTLSolverEbos.hpp @@ -645,27 +645,6 @@ protected: #endif } - /// Zero out off-diagonal blocks on rows corresponding to overlap cells - /// Diagonal blocks on ovelap rows are set to diag(1.0). - void makeOverlapRowsInvalid(Matrix& ebosJacIgnoreOverlap) const - { - //value to set on diagonal - Dune::FieldMatrix diag_block(0.0); - for (int eq = 0; eq < numEq; ++eq) - diag_block[eq][eq] = 1.0; - - //loop over precalculated overlap rows and columns - for (auto row = overlapRows_.begin(); row != overlapRows_.end(); row++ ) - { - int lcell = *row; - // Zero out row. - ebosJacIgnoreOverlap[lcell] = 0.0; - - //diagonal block set to diag(1.0). - ebosJacIgnoreOverlap[lcell][lcell] = diag_block; - } - } - /// Create sparsity pattern of matrix without off-diagonal ghost entries. void noGhostAdjacency() {