diff --git a/src/Fortran_libraries/SOLVER_src/direct_solver/m_ludcmp_band.f90 b/src/Fortran_libraries/SOLVER_src/direct_solver/m_ludcmp_band.f90
index f257cf25..b4c2b915 100644
--- a/src/Fortran_libraries/SOLVER_src/direct_solver/m_ludcmp_band.f90
+++ b/src/Fortran_libraries/SOLVER_src/direct_solver/m_ludcmp_band.f90
@@ -115,7 +115,8 @@ subroutine ludcmp_band(n, nband, band_a, band_lu ,indx, d)
           aamax = max(aamax, abs(band_a(k,j)))
 !          aamax = max(aamax, abs(a(i,j)))
         end do
-        if (aamax.eq.0.0d0) pause 'singular matrix in ludcmp_band'
+        if (aamax.eq.0.0d0) write(*,*)                                  &
+     &                    '!!!!!!  singular matrix in ludcmp  !!!!!!'
         vv(i)=1.0d0 / aamax
         idx_org(i) = i
       end do