diff --git a/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal3dpressure2p.hh b/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal3dpressure2p.hh index 91c56e8dfb980fab70732ac6716ce9e2c7bae1ce..7dd3e72d94729f5bb07f59883d4addcc4ef83176 100644 --- a/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal3dpressure2p.hh +++ b/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal3dpressure2p.hh @@ -506,11 +506,11 @@ void FvMpfaL3dPressure2p<TypeTag>::initializeMatrixRowSize() std::set<int> neighborIndices; - int numCorners = eIt->geometry().corners(); + int numVertices = eIt->geometry().corners(); - for (int cornerNum = 0; cornerNum < numCorners; cornerNum++) + for (int vIdx = 0; vIdx < numVertices; vIdx++) { - int globalVertIdx = problem_.variables().vertexMapper().map(*eIt, cornerNum, dim); + int globalVertIdx = problem_.variables().vertexMapper().map(*eIt, vIdx, dim); InteractionVolume& interactionVolume = interactionVolumes_.interactionVolume(globalVertIdx); @@ -546,11 +546,11 @@ void FvMpfaL3dPressure2p<TypeTag>::initializeMatrixIndices() // add diagonal index this->A_.addindex(globalIdxI, globalIdxI); - int numCorners = eIt->geometry().corners(); + int numVertices = eIt->geometry().corners(); - for (int cornerNum = 0; cornerNum < numCorners; cornerNum++) + for (int vIdx = 0; vIdx < numVertices; vIdx++) { - int globalVertIdx = problem_.variables().vertexMapper().map(*eIt, cornerNum, dim); + int globalVertIdx = problem_.variables().vertexMapper().map(*eIt, vIdx, dim); InteractionVolume& interactionVolume = interactionVolumes_.interactionVolume(globalVertIdx); for (int subVolumeIdx = 0; subVolumeIdx < InteractionVolume::subVolumeTotalNum; subVolumeIdx++) diff --git a/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal3dpressure2padaptive.hh b/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal3dpressure2padaptive.hh index fffc4234b7ab9f0f8b7498fab42a04f410f59df9..e300fcd601586e70c876dd0be2ed126a4c538355 100644 --- a/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal3dpressure2padaptive.hh +++ b/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal3dpressure2padaptive.hh @@ -263,11 +263,11 @@ void FvMpfaL3dPressure2pAdaptive<TypeTag>::initializeMatrixRowSize() std::set<int> neighborIndices; - int numCorners = eIt->geometry().corners(); + int numVertices = eIt->geometry().corners(); - for (int cornerNum = 0; cornerNum < numCorners; cornerNum++) + for (int vIdx = 0; vIdx < numVertices; vIdx++) { - int globalVertIdx = problem_.variables().vertexMapper().map(*eIt, cornerNum, dim); + int globalVertIdx = problem_.variables().vertexMapper().map(*eIt, vIdx, dim); InteractionVolume& interactionVolume = this->interactionVolumes_.interactionVolume(globalVertIdx); @@ -284,7 +284,7 @@ void FvMpfaL3dPressure2pAdaptive<TypeTag>::initializeMatrixRowSize() { if (neighbor->level() == levelI + 2) { - for (int vIdx = 0; vIdx < numCorners; vIdx++) + for (int vIdx = 0; vIdx < numVertices; vIdx++) { int globalVertIdxJ = problem_.variables().vertexMapper().map(*neighbor, vIdx, dim); @@ -351,11 +351,11 @@ void FvMpfaL3dPressure2pAdaptive<TypeTag>::initializeMatrixIndices() // add diagonal index this->A_.addindex(globalIdxI, globalIdxI); - int numCorners = eIt->geometry().corners(); + int numVertices = eIt->geometry().corners(); - for (int cornerNum = 0; cornerNum < numCorners; cornerNum++) + for (int vIdx = 0; vIdx < numVertices; vIdx++) { - int globalVertIdx = problem_.variables().vertexMapper().map(*eIt, cornerNum, dim); + int globalVertIdx = problem_.variables().vertexMapper().map(*eIt, vIdx, dim); InteractionVolume& interactionVolume = this->interactionVolumes_.interactionVolume(globalVertIdx); for (int subVolumeIdx = 0; subVolumeIdx < InteractionVolume::subVolumeTotalNum; subVolumeIdx++)