FlapJackApps/subsurface
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Don't do Boyle compensation when still on first stop

Browse Source 
There is nothing to compensate but we would hit a 0/0 numerical
instability there.

Signed-off-by: Robert C. Helling <helling@atdotde.de>
Signed-off-by: Dirk Hohndel <dirk@hohndel.org>
This commit is contained in:
Robert C. Helling 2015-08-21 17:19:32 +02:00 committed by Dirk Hohndel
parent 0f096d0e53
commit d9b4f2b7e3
1 changed files with 4 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
deco.c
View File

@ -362,6 +362,10 @@ void boyles_law(double next_stop_pressure)
{
int ci;
// This should be a tautology but prevents a numerical instability.
if (IS_FP_SAME(next_stop_pressure, first_stop_pressure))
return;
if (!first_stop_pressure)
return;
for (ci = 0; ci < 16; ++ci) {