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

[vangenuchten] Improve doc

parent 7e3e2ad7
......@@ -356,15 +356,14 @@ public:
// phase from the parameters
const Scalar swThHigh = params.krwHighSw();
// derivative of the regualarization
if (swe < 0) {
// the slope is zero
// derivative of the regularization
// the slope is zero below sw=0.0 and above sw=1.0
if (swe < 0)
return 0.0;
}
else if (swe > 1 - std::numeric_limits<Scalar>::epsilon()) {
// the slope is zero
else if (swe > 1 - std::numeric_limits<Scalar>::epsilon())
return 0.0;
}
// for high sw we need the slope of the interpolation spline
else if (swe > swThHigh) {
using Spline = Dumux::Spline<Scalar>;
Spline sp(swThHigh, 1.0, // x1, x2
......
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