OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Make it compile with dune-fem

Browse Source
This commit is contained in:
Tor Harald Sandve 2020-11-19 11:13:55 +01:00
parent 0ac2f922d3
commit c4c0b7a756
1 changed files with 2 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
opm/simulators/linalg/ISTLSolverEbos.hpp
View File

@ -90,6 +90,7 @@ namespace Opm
using FlexibleSolverType = Dune::FlexibleSolver<Matrix, Vector>; using FlexibleSolverType = Dune::FlexibleSolver<Matrix, Vector>;
using AbstractOperatorType = Dune::AssembledLinearOperator<Matrix, Vector, Vector>; using AbstractOperatorType = Dune::AssembledLinearOperator<Matrix, Vector, Vector>;
using WellModelOperator = WellModelAsLinearOperator<WellModel, Vector, Vector>; using WellModelOperator = WellModelAsLinearOperator<WellModel, Vector, Vector>;
using ElementMapper = GetPropType<TypeTag, Properties::ElementMapper>;
#if HAVE_CUDA || HAVE_OPENCL #if HAVE_CUDA || HAVE_OPENCL
static const unsigned int block_size = Matrix::block_type::rows; static const unsigned int block_size = Matrix::block_type::rows;
@ -152,9 +153,7 @@ namespace Opm
// For some reason simulator_.model().elementMapper() is not initialized at this stage // For some reason simulator_.model().elementMapper() is not initialized at this stage
// Hence const auto& elemMapper = simulator_.model().elementMapper(); does not work. // Hence const auto& elemMapper = simulator_.model().elementMapper(); does not work.
// Set it up manually // Set it up manually
using ElementMapper = ElementMapper elemMapper(simulator_.vanguard().gridView(), Dune::mcmgElementLayout());
Dune::MultipleCodimMultipleGeomTypeMapper<GridView>;
ElementMapper elemMapper(simulator_.vanguard().grid().leafGridView(), Dune::mcmgElementLayout());
detail::findOverlapAndInterior(simulator_.vanguard().grid(), elemMapper, overlapRows_, interiorRows_); detail::findOverlapAndInterior(simulator_.vanguard().grid(), elemMapper, overlapRows_, interiorRows_);
useWellConn_ = EWOMS_GET_PARAM(TypeTag, bool, MatrixAddWellContributions); useWellConn_ = EWOMS_GET_PARAM(TypeTag, bool, MatrixAddWellContributions);