Skip to content
Snippets Groups Projects
Commit eb47a676 authored by Stefan Girke's avatar Stefan Girke
Browse files

[bugfix] pressure is problem dependent

parent 367e0537
No related branches found
No related tags found
No related merge requests found
......@@ -200,8 +200,8 @@ class NSModel : public DefaultModel < NSModelTraits< GridPartType, ProblemImp >
const RangeType& uRight,
RangeType& jump) const
{
const double pl = pressure( uLeft );
const double pr = pressure( uRight );
const double pl = problem_.pressure( uLeft );
const double pr = problem_.pressure( uRight );
// take pressure as shock detection values
jump = (pl-pr)/(0.5*(pl+pr));
......@@ -232,7 +232,7 @@ class NSModel : public DefaultModel < NSModelTraits< GridPartType, ProblemImp >
return false;
else
{
const double p = pressure( u );
const double p = problem_.pressure( u );
return p > 1e-8;
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment