Commit dea0b6a0 authored by Dennis Gläser's avatar Dennis Gläser Committed by Timo Koch
Browse files

[mpfa-o][iavol] use dimvector for the local omegas

If dimWorld > dim, using a GlobalPosition vector is wrong, we have to
use a vector of size=dim here.
parent 681c1e16
...@@ -116,6 +116,7 @@ class CCMpfaOInteractionVolume : public CCMpfaInteractionVolumeBase<TypeTag, Tra ...@@ -116,6 +116,7 @@ class CCMpfaOInteractionVolume : public CCMpfaInteractionVolumeBase<TypeTag, Tra
static const int dim = GridView::dimension; static const int dim = GridView::dimension;
using Element = typename GridView::template Codim<0>::Entity; using Element = typename GridView::template Codim<0>::Entity;
using DimVector = typename Traits::DimVector;
using GlobalPosition = typename Traits::GlobalPosition; using GlobalPosition = typename Traits::GlobalPosition;
using DynamicVector = typename Traits::Vector; using DynamicVector = typename Traits::Vector;
using DynamicMatrix = typename Traits::Matrix; using DynamicMatrix = typename Traits::Matrix;
...@@ -390,7 +391,7 @@ private: ...@@ -390,7 +391,7 @@ private:
B[idxInFluxFaces][posLocalScvIdx] += posWijk[localDir]; B[idxInFluxFaces][posLocalScvIdx] += posWijk[localDir];
} }
// If not on a boundary or interior dirichlet face, add entries for the "negative" scv // If face is not on a boundary or an interior dirichlet face, add entries for the "negative" scv
if (faceType == MpfaFaceTypes::interior) if (faceType == MpfaFaceTypes::interior)
{ {
const auto negLocalScvIdx = localScvf.outsideLocalScvIndex(); const auto negLocalScvIdx = localScvf.outsideLocalScvIndex();
...@@ -472,7 +473,7 @@ private: ...@@ -472,7 +473,7 @@ private:
} }
} }
GlobalPosition calculateOmegas_(const LocalScvType& localScv, DimVector calculateOmegas_(const LocalScvType& localScv,
const LocalScvfType& localScvf, const LocalScvfType& localScvf,
const Tensor& T) const const Tensor& T) const
{ {
...@@ -490,7 +491,7 @@ private: ...@@ -490,7 +491,7 @@ private:
return wijk; return wijk;
} }
GlobalPosition calculateOmegas_(const LocalScvType& localScv, DimVector calculateOmegas_(const LocalScvType& localScv,
const LocalScvfType& localScvf, const LocalScvfType& localScvf,
const Scalar t) const const Scalar t) const
{ {
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment