JIANG/EPANET
1
0
Fork 0
Code Issues Pull Requests Actions 4 Packages Projects Releases Wiki Activity

Avoid possible NAN in leakage.c

Browse Source
This commit is contained in:
Lew Rossman
2025-09-10 11:06:05 -04:00
parent b1630f3091
commit d345541af7
1 changed files with 11 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

17
src/leakage.c
View File

@@ -7,7 +7,7 @@
Authors: see AUTHORS Authors: see AUTHORS
Copyright: see AUTHORS Copyright: see AUTHORS
License: see LICENSE License: see LICENSE
Last Updated: 06/14/2024 Last Updated: 09/10/2025
****************************************************************************** ******************************************************************************
*/ */
/* /*
@@ -432,12 +432,15 @@ int leakagehasconverged(Project *pr)
// Directly compute a reference leakage at this pressure head // Directly compute a reference leakage at this pressure head
qref = 0.0; qref = 0.0;
// Contribution from pipes with fixed area leaks if ( h > 0.0)
if (hyd->Leakage[i].cfa > 0.0) {
qref = sqrt(h / hyd->Leakage[i].cfa); // Contribution from pipes with fixed area leaks
// Contribution from pipes with variable area leaks if (hyd->Leakage[i].cfa > 0.0)
if (hyd->Leakage[i].cva > 0.0) qref = sqrt(h / hyd->Leakage[i].cfa);
qref += pow((h / hyd->Leakage[i].cva), 1.5); // Contribution from pipes with variable area leaks
if (hyd->Leakage[i].cva > 0.0)
qref += pow((h / hyd->Leakage[i].cva), 1.5);
}
// Compare reference leakage to solution leakage // Compare reference leakage to solution leakage
qtest = hyd->Leakage[i].qfa + hyd->Leakage[i].qva; qtest = hyd->Leakage[i].qfa + hyd->Leakage[i].qva;