Another crack at debug output for Cholesky factorization errors in WSMP

src/solve_with_pwssmp.f90

@@ -211,12 +211,25 @@ call pwssmp(n_iproc,ia,ja,avals,diag,perm,invp,b,ldb,nrhs,aux,naux,mrp,iparm,dpa
 time2=mpi_wtime()
 
 if (iproc.eq.0) then
-   if (iparm(64) .ne. 0) then
+   if (iparm(64) .gt. 0) then
       call write_streepjes(6,1)
       write (*,*) 'Cholesky factorization problem at location ',iparm(64),' in WSMP'
       call write_streepjes(6,1)
+   elseif (iparm(64).lt.0) then
+      call write_streepjes(6,1)
+      write (*,*) 'Cholesky factorization problem ',iparm(64),' in WSMP'
+      call write_streepjes(6,1)
+   endif
+endif
+
+if (iparm(64).gt.0) then
+   if (iparm(64).ge.n_iproc_st .and. iparm(64).le.n_iproc_end) then
+      badproc=iproc
+      write (*,*) 'WSMP failed on process ',badproc,'.'
+   endif
+   if (iproc.eq.badproc) then
       if (params%debug.gt.0) then
-        write (*,*) 'Debug output:'
+        write (*,*) 'Debug output from process ',badproc,':'
         write (*,*) 'Bad main diagonal value (avals): ',avals(ia(iparm(64)))
         idof=mod(iparm(64),ndof)
         if (idof.eq.0) idof=3
@@ -254,6 +267,12 @@ if (iproc.eq.0) then
         enddo
       endif
       call write_streepjes(6,1)
+      deallocate(bad_elems)
+   endif
+endif
+
+if (iproc.eq.0) then
+   if (iparm(64).ne.0) then
       !print *,'WSMP-> Cholesky factorisation pb',iparm(64)
       stop
    else