Fix divide-by-zero in aromaticity index calculation

robertyoung3 · robertyoung3 · commit 615108ab570c · 2026-03-03T11:24:50.000+10:30
diff --git a/corems/molecular_formula/calc/MolecularFormulaCalc.py b/corems/molecular_formula/calc/MolecularFormulaCalc.py
@@ -560,7 +560,7 @@ def _calc_aromaticity_index_mod(self):
         )
         ai_d = ai_es["C"] - (0.5 * ai_es["O"]) - ai_es["N"] - ai_es["S"]
 
-        ai = ai_n / ai_d
+        ai = ai_n / ai_d if ai_d != 0 else 0
 
         if ai < 0:
             ai = 0
@@ -601,7 +601,7 @@ def _calc_aromaticity_index(self):
                 )
         ai_d = ai_es["C"] - (ai_es["O"]) - ai_es["N"] - ai_es["S"]
 
-        ai = ai_n / ai_d
+        ai = ai_n / ai_d if ai_d != 0 else 0
 
         if ai < 0:
             ai = 0

Original file line number	Diff line number	Diff line change
`@@ -560,7 +560,7 @@ def _calc_aromaticity_index_mod(self):`
`560`	`560`	`)`
`561`	`561`	`ai_d = ai_es["C"] - (0.5 * ai_es["O"]) - ai_es["N"] - ai_es["S"]`
`562`	`562`
`563`		`- ai = ai_n / ai_d`
	`563`	`+ ai = ai_n / ai_d if ai_d != 0 else 0`
`564`	`564`
`565`	`565`	`if ai < 0:`
`566`	`566`	`ai = 0`
`@@ -601,7 +601,7 @@ def _calc_aromaticity_index(self):`
`601`	`601`	`)`
`602`	`602`	`ai_d = ai_es["C"] - (ai_es["O"]) - ai_es["N"] - ai_es["S"]`
`603`	`603`
`604`		`- ai = ai_n / ai_d`
	`604`	`+ ai = ai_n / ai_d if ai_d != 0 else 0`
`605`	`605`
`606`	`606`	`if ai < 0:`
`607`	`607`	`ai = 0`