Skip to content

Commit

Permalink
still issues with #413
Browse files Browse the repository at this point in the history
  • Loading branch information
Jorrit Wronski committed Jan 20, 2015
1 parent 104175d commit 2a325ab
Showing 1 changed file with 14 additions and 3 deletions.
17 changes: 14 additions & 3 deletions src/Backends/Incompressible/IncompressibleFluid.cpp
Expand Up @@ -23,12 +23,11 @@ void IncompressibleFluid::set_reference_state(double T0, double p0, double x0, d
this->rhoref = rho(T0,p0,x0);
this->pref = p0;
this->href = h0;
this->sref = s0;
// Now we take care of the energy related values
this->uref = 0.0;
this->uref = u(T0,p0,x0) - h0; // (value without ref) - (desired ref)
//this->sref = 0.0;
//this->sref = s(T0,p0,x0) - s0; // (value without ref) - (desired ref)
this->sref = 0.0;
this->sref = s(T0,p0,x0) - s0; // (value without ref) - (desired ref)
}

void IncompressibleFluid::validate(){
Expand Down Expand Up @@ -665,8 +664,20 @@ TEST_CASE("Internal consistency checks and example use cases for the incompressi
// Compare reference state
{
res = XLT.h(Tref,pref,xref);
CAPTURE(Tref);
CAPTURE(pref);
CAPTURE(xref);
CAPTURE(href)
CAPTURE(res)
CHECK( check_abs(href,res,acc) );
}
{
res = XLT.s(Tref,pref,xref);
CAPTURE(Tref);
CAPTURE(pref);
CAPTURE(xref);
CAPTURE(sref)
CAPTURE(res)
CHECK( check_abs(sref,res,acc) );
}

Expand Down

0 comments on commit 2a325ab

Please sign in to comment.