Commit fc2ab1f1 authored by Kilian Weishaupt's avatar Kilian Weishaupt
Browse files

[cleanup] Use conv functions in 2pdfm

Use volvars.saturation(phaseIdx) and volvars.density(phaseIdx)
instead of volvars.fluidstate().x(phaseIdx)
parent 583e1bf7
......@@ -227,10 +227,10 @@ private:
// calculate the phase density at the integration point. we
// only do this if the wetting phase is present in both cells
Scalar SI = elemVolVars[this->face().i].fluidState().saturation(phaseIdx);
Scalar SJ = elemVolVars[this->face().j].fluidState().saturation(phaseIdx);
Scalar rhoI = elemVolVars[this->face().i].fluidState().density(phaseIdx);
Scalar rhoJ = elemVolVars[this->face().j].fluidState().density(phaseIdx);
Scalar SI = elemVolVars[this->face().i].saturation(phaseIdx);
Scalar SJ = elemVolVars[this->face().j].saturation(phaseIdx);
Scalar rhoI = elemVolVars[this->face().i].density(phaseIdx);
Scalar rhoJ = elemVolVars[this->face().j].density(phaseIdx);
Scalar fI = std::max(0.0, std::min(SI/1e-5, 0.5));
Scalar fJ = std::max(0.0, std::min(SJ/1e-5, 0.5));
if (Dune::FloatCmp::eq<Scalar, Dune::FloatCmp::absolute>(fI + fJ, 0.0, 1.0e-30))
......
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