Don't do Boyle compensation when still on first stop

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>

1 files changed, 4 insertions, 0 deletions

diff --git a/deco.c b/deco.c
index 7b8d7c2be..16089e32b 100644
--- a/deco.c
+++ b/deco.c
@@ -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) {