diff --git a/base/serial/aux/dasr.f90 b/base/serial/aux/dasr.f90
index 77502684..ea4c8ab8 100644
--- a/base/serial/aux/dasr.f90
+++ b/base/serial/aux/dasr.f90
@@ -134,8 +134,8 @@ subroutine dasr(n,x,dir)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='dasr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -258,8 +258,8 @@ subroutine dasr(n,x,dir)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='dasr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/dasrx.f90 b/base/serial/aux/dasrx.f90
index 35e0c912..b408bacc 100644
--- a/base/serial/aux/dasrx.f90
+++ b/base/serial/aux/dasrx.f90
@@ -160,8 +160,8 @@ subroutine dasrx(n,x,indx,dir,flag)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='dasrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -298,8 +298,8 @@ subroutine dasrx(n,x,indx,dir,flag)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='dasrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/dsr.f90 b/base/serial/aux/dsr.f90
index c1885cf5..614f57e7 100644
--- a/base/serial/aux/dsr.f90
+++ b/base/serial/aux/dsr.f90
@@ -134,8 +134,8 @@ subroutine dsr(n,x,dir)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='dsr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -258,8 +258,8 @@ subroutine dsr(n,x,dir)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='dsr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/dsrx.f90 b/base/serial/aux/dsrx.f90
index 82d11d08..5c1d2091 100644
--- a/base/serial/aux/dsrx.f90
+++ b/base/serial/aux/dsrx.f90
@@ -161,8 +161,8 @@ subroutine dsrx(n,x,indx,dir,flag)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='dsrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -300,8 +300,8 @@ subroutine dsrx(n,x,indx,dir,flag)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='dsrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/iasr.f90 b/base/serial/aux/iasr.f90
index ba719305..ddf7e6b9 100644
--- a/base/serial/aux/iasr.f90
+++ b/base/serial/aux/iasr.f90
@@ -134,8 +134,8 @@ subroutine iasr(n,x,dir)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='iasr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -258,8 +258,8 @@ subroutine iasr(n,x,dir)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='iasr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/iasrx.f90 b/base/serial/aux/iasrx.f90
index b86cff83..ca98864f 100644
--- a/base/serial/aux/iasrx.f90
+++ b/base/serial/aux/iasrx.f90
@@ -160,8 +160,8 @@ subroutine iasrx(n,x,indx,dir,flag)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='iasrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -298,8 +298,8 @@ subroutine iasrx(n,x,indx,dir,flag)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='iasrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/isr.f90 b/base/serial/aux/isr.f90
index f0319999..fed4d9ea 100644
--- a/base/serial/aux/isr.f90
+++ b/base/serial/aux/isr.f90
@@ -134,8 +134,8 @@ subroutine isr(n,x,dir)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='isr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -258,8 +258,8 @@ subroutine isr(n,x,dir)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='isr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/isrx.f90 b/base/serial/aux/isrx.f90
index 325a7927..12b4e8cd 100644
--- a/base/serial/aux/isrx.f90
+++ b/base/serial/aux/isrx.f90
@@ -159,9 +159,9 @@ subroutine isrx(n,x,indx,dir,flag)
           end if
         end do outer_up
         if (i == ilx) then 
-          if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+          if (x(i) /= piv) then 
+            call psb_errpush(4001,r_name='isrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -300,8 +300,8 @@ subroutine isrx(n,x,indx,dir,flag)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='isrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/zalsr.f90 b/base/serial/aux/zalsr.f90
index 08faca08..89364edc 100644
--- a/base/serial/aux/zalsr.f90
+++ b/base/serial/aux/zalsr.f90
@@ -135,8 +135,8 @@ subroutine zalsr(n,x,dir)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zalsr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -259,8 +259,8 @@ subroutine zalsr(n,x,dir)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zalsr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/zalsrx.f90 b/base/serial/aux/zalsrx.f90
index e0811f81..4895ea06 100644
--- a/base/serial/aux/zalsrx.f90
+++ b/base/serial/aux/zalsrx.f90
@@ -162,8 +162,8 @@ subroutine zalsrx(n,x,indx,dir,flag)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zalsrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -302,8 +302,8 @@ subroutine zalsrx(n,x,indx,dir,flag)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zalsrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/zasr.f90 b/base/serial/aux/zasr.f90
index ca391153..4aacf1e9 100644
--- a/base/serial/aux/zasr.f90
+++ b/base/serial/aux/zasr.f90
@@ -135,8 +135,8 @@ subroutine zasr(n,x,dir)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zasr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -259,8 +259,8 @@ subroutine zasr(n,x,dir)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zasr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/zasrx.f90 b/base/serial/aux/zasrx.f90
index de79b048..3bc516bf 100644
--- a/base/serial/aux/zasrx.f90
+++ b/base/serial/aux/zasrx.f90
@@ -162,8 +162,8 @@ subroutine zasrx(n,x,indx,dir,flag)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zasrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -302,8 +302,8 @@ subroutine zasrx(n,x,indx,dir,flag)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zasrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/zlsr.f90 b/base/serial/aux/zlsr.f90
index d7c12abd..1b3b40be 100644
--- a/base/serial/aux/zlsr.f90
+++ b/base/serial/aux/zlsr.f90
@@ -135,8 +135,8 @@ subroutine zlsr(n,x,dir)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zlsr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -259,8 +259,8 @@ subroutine zlsr(n,x,dir)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zlsr',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
diff --git a/base/serial/aux/zlsrx.f90 b/base/serial/aux/zlsrx.f90
index 33c0b9c5..76c2c37a 100644
--- a/base/serial/aux/zlsrx.f90
+++ b/base/serial/aux/zlsrx.f90
@@ -162,8 +162,8 @@ subroutine zlsrx(n,x,indx,dir,flag)
         end do outer_up
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zlsrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif
@@ -302,8 +302,8 @@ subroutine zlsrx(n,x,indx,dir,flag)
         end do outer_dw
         if (i == ilx) then 
           if (x(i) /= piv) then
-            write(0,*) 'Should never ever get here????!!!!'
-            stop
+            call psb_errpush(4001,r_name='zlsrx',a_err='impossible pivot condition')
+            call psb_error()
           endif
           i = i + 1 
         endif