From 71511757ce60220214e4a054aedf88851a5fe11b Mon Sep 17 00:00:00 2001 From: Harry Moffat Date: Mon, 14 Aug 2006 20:48:36 +0000 Subject: [PATCH] Fixed a roundoff issue --- test_problems/cathermo/HMW_graph_VvT/HMW_graph_VvT.cpp | 10 +++++++--- 1 file changed, 7 insertions(+), 3 deletions(-) diff --git a/test_problems/cathermo/HMW_graph_VvT/HMW_graph_VvT.cpp b/test_problems/cathermo/HMW_graph_VvT/HMW_graph_VvT.cpp index fd1053aa6..00757bd07 100755 --- a/test_problems/cathermo/HMW_graph_VvT/HMW_graph_VvT.cpp +++ b/test_problems/cathermo/HMW_graph_VvT/HMW_graph_VvT.cpp @@ -151,8 +151,8 @@ int main(int argc, char **argv) */ double T; - double V0_NaCl, V0_Naplus, V0_Clminus, Delta_V0s, V0_H2O; - double V_NaCl, V_Naplus, V_Clminus, V_H2O; + double V0_NaCl = 0.0, V0_Naplus = 0.0, V0_Clminus = 0.0, Delta_V0s = 0.0, V0_H2O = 0.0; + double V_NaCl = 0.0, V_Naplus = 0.0, V_Clminus = 0.0, V_H2O = 0.0; double molarV0; #ifdef DEBUG_HKM FILE *ttt = fopen("table.csv","w"); @@ -325,8 +325,12 @@ int main(int argc, char **argv) V0_Clminus*1.E3, V_Clminus*1.E3, (V_Clminus - V0_Clminus)*1.E3); printf(" NaCl(s)"); + double dd = V_NaCl*1.E3 - V0_NaCl*1.E3; + if (fabs(dd) < 1.0E-12) { + dd = 0.0; + } printf("%13g %13g %13g %13g\n", 1.0, - V0_NaCl*1.E3 , V_NaCl*1.E3, V_NaCl*1.E3 - V0_NaCl*1.E3); + V0_NaCl*1.E3 , V_NaCl*1.E3, dd); delete HMW;