add the general framework for specifying auxiliary equations.

this does not disrupt the block nature of the linearized matrix
(i.e. Dune::BCRSMatrix is still used), but if the number of auxiliary
equations is smaller than that of the "main" discretization, the
superfluous equations are padded. if the number of additional
equations are larger than that of the equation, additional DOFs are
added.
This commit is contained in:
Andreas Lauser 2014-10-06 16:23:21 +02:00
parent 1d18473404
commit 2050913aef

View File

@ -351,7 +351,7 @@ public:
// calculate the intersection index
const auto &scvf = context.stencil(timeIdx).interiorFace(localIntersectionIdx);
int numElements = this->model().numDof();
int numElements = this->model().numGridDof();
size_t interiorElemIdx = context.globalSpaceIndex(scvf.interiorIndex(), timeIdx);
size_t exteriorElemIdx = context.globalSpaceIndex(scvf.exteriorIndex(), timeIdx);
@ -501,7 +501,7 @@ private:
auto eclipseState = this->simulator().gridManager().eclipseState();
const auto &grid = this->simulator().gridManager().grid();
size_t numDof = this->model().numDof();
size_t numDof = this->model().numGridDof();
intrinsicPermeability_.resize(numDof);
porosity_.resize(numDof);
@ -718,7 +718,7 @@ private:
OPM_THROW(std::runtime_error,
"The Eclipse input file requires the RV keyword to be non-present");
size_t numDof = this->model().numDof();
size_t numDof = this->model().numGridDof();
initialFluidStates_.resize(numDof);