From b89f53f40ba57799186663e4c4500ac293266ea0 Mon Sep 17 00:00:00 2001 From: Anthony Scemama Date: Tue, 8 Dec 2020 22:05:23 +0100 Subject: [PATCH] Remove possible division by zero --- src/utils/linear_algebra.irp.f | 1 + 1 file changed, 1 insertion(+) 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