fix: handle residuals array is empty

Signed-off-by: Endre Moen <[email protected]>
py-why · Jan 3, 2025 · 2ecc5d3 · 2ecc5d3
1 parent 8c72b52
commit 2ecc5d3
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/causallearn/search/ScoreBased/ExactSearch.py b/causallearn/search/ScoreBased/ExactSearch.py
@@ -365,6 +365,8 @@ def bic_score_node(X, i, structure):
                                             b=X[:, i],
                                             rcond=None)
     bic = n * np.log(residual / n) + len(structure) * np.log(n)
+    if bic.size == 0:
+        return NEGINF   # Return negative infinity if bic is empty
     return bic.item()