diff --git a/base/serial/impl/psb_c_coo_impl.f90 b/base/serial/impl/psb_c_coo_impl.f90
index dc840176..9646a91a 100644
--- a/base/serial/impl/psb_c_coo_impl.f90
+++ b/base/serial/impl/psb_c_coo_impl.f90
@@ -451,6 +451,8 @@ subroutine psb_c_coo_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_c_csc_impl.f90 b/base/serial/impl/psb_c_csc_impl.f90
index 1f8b732b..270a4f76 100644
--- a/base/serial/impl/psb_c_csc_impl.f90
+++ b/base/serial/impl/psb_c_csc_impl.f90
@@ -2782,6 +2782,8 @@ subroutine psb_c_csc_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_c_csr_impl.f90 b/base/serial/impl/psb_c_csr_impl.f90
index 51d45825..c91f7cc8 100644
--- a/base/serial/impl/psb_c_csr_impl.f90
+++ b/base/serial/impl/psb_c_csr_impl.f90
@@ -2609,6 +2609,8 @@ subroutine psb_c_csr_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_d_coo_impl.f90 b/base/serial/impl/psb_d_coo_impl.f90
index 1c257972..f923822b 100644
--- a/base/serial/impl/psb_d_coo_impl.f90
+++ b/base/serial/impl/psb_d_coo_impl.f90
@@ -451,6 +451,8 @@ subroutine psb_d_coo_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_d_csc_impl.f90 b/base/serial/impl/psb_d_csc_impl.f90
index 04572bf7..33e511b5 100644
--- a/base/serial/impl/psb_d_csc_impl.f90
+++ b/base/serial/impl/psb_d_csc_impl.f90
@@ -2782,6 +2782,8 @@ subroutine psb_d_csc_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_d_csr_impl.f90 b/base/serial/impl/psb_d_csr_impl.f90
index 094f141e..4a5841b3 100644
--- a/base/serial/impl/psb_d_csr_impl.f90
+++ b/base/serial/impl/psb_d_csr_impl.f90
@@ -2609,6 +2609,8 @@ subroutine psb_d_csr_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_s_coo_impl.f90 b/base/serial/impl/psb_s_coo_impl.f90
index 4ad5bf23..452d8c15 100644
--- a/base/serial/impl/psb_s_coo_impl.f90
+++ b/base/serial/impl/psb_s_coo_impl.f90
@@ -451,6 +451,8 @@ subroutine psb_s_coo_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_s_csc_impl.f90 b/base/serial/impl/psb_s_csc_impl.f90
index aed2ada9..78ece378 100644
--- a/base/serial/impl/psb_s_csc_impl.f90
+++ b/base/serial/impl/psb_s_csc_impl.f90
@@ -2782,6 +2782,8 @@ subroutine psb_s_csc_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_s_csr_impl.f90 b/base/serial/impl/psb_s_csr_impl.f90
index 7fe08ee0..5117279a 100644
--- a/base/serial/impl/psb_s_csr_impl.f90
+++ b/base/serial/impl/psb_s_csr_impl.f90
@@ -2609,6 +2609,8 @@ subroutine psb_s_csr_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_z_coo_impl.f90 b/base/serial/impl/psb_z_coo_impl.f90
index f625f23f..fb698514 100644
--- a/base/serial/impl/psb_z_coo_impl.f90
+++ b/base/serial/impl/psb_z_coo_impl.f90
@@ -451,6 +451,8 @@ subroutine psb_z_coo_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_z_csc_impl.f90 b/base/serial/impl/psb_z_csc_impl.f90
index 0354bef8..944a7fca 100644
--- a/base/serial/impl/psb_z_csc_impl.f90
+++ b/base/serial/impl/psb_z_csc_impl.f90
@@ -2782,6 +2782,8 @@ subroutine psb_z_csc_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then 
diff --git a/base/serial/impl/psb_z_csr_impl.f90 b/base/serial/impl/psb_z_csr_impl.f90
index 4713da34..43b5e849 100644
--- a/base/serial/impl/psb_z_csr_impl.f90
+++ b/base/serial/impl/psb_z_csr_impl.f90
@@ -2609,6 +2609,8 @@ subroutine psb_z_csr_print(iout,a,iv,head,ivr,ivc)
   nc = a%get_ncols()
   nz = a%get_nzeros()
   nmx = max(nr,nc,1)
+  if (present(ivr)) nmx = max(nmx,maxval(abs(ivr)))
+  if (present(ivc)) nmx = max(nmx,maxval(abs(ivc)))
   ni  = floor(log10(1.0*nmx)) + 1
 
   if (datatype=='real') then