From cd3e8d874178332c6636fc58e83c27fba81c8460 Mon Sep 17 00:00:00 2001 From: Harry Moffat Date: Mon, 31 Mar 2008 23:51:41 +0000 Subject: [PATCH] Replaced a naked floating point equality check with a fuzzy check. This was causing problems in the optimized, DEBUG_MODE code. --- Cantera/src/equil/vcs_solve_TP.cpp | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/Cantera/src/equil/vcs_solve_TP.cpp b/Cantera/src/equil/vcs_solve_TP.cpp index 464ab401c..0b7dd6b88 100644 --- a/Cantera/src/equil/vcs_solve_TP.cpp +++ b/Cantera/src/equil/vcs_solve_TP.cpp @@ -948,7 +948,8 @@ namespace VCSnonideal { * This should keep the amount of material constant. */ #ifdef DEBUG_MODE - if (ds[kspec] != dx) { + if (fabs(ds[kspec] -dx) > 1.0E-14*(fabs(ds[kspec]) + fabs(dx) + 1.0E-32)) { + plogf(" ds[kspec] = %20.16g dx = %20.16g , kspec = %d\n", ds[kspec], dx, kspec); plogf("we have a problem!\n"); exit(-1); }