diff --git a/src/utils/linear_algebra.irp.f b/src/utils/linear_algebra.irp.f
index 26d0b78c..d858341d 100644
--- a/src/utils/linear_algebra.irp.f
+++ b/src/utils/linear_algebra.irp.f
@@ -1550,6 +1550,7 @@ subroutine nullify_small_elements(m,n,A,LDA,thresh)
       amax = max(dabs(A(i,j)), amax)
     enddo
   enddo
+  if (amax == 0.d0) return
   amax = 1.d0/amax
 
   ! Remove tiny elements