diff --git a/examples/pdegen/mld_dexample_1lev.f90 b/examples/pdegen/mld_dexample_1lev.f90
index c7681181..fc32206d 100644
--- a/examples/pdegen/mld_dexample_1lev.f90
+++ b/examples/pdegen/mld_dexample_1lev.f90
@@ -90,6 +90,7 @@ program mld_dexample_1lev
   ! other variables
   integer            :: i,info,j
   integer(psb_epk_) :: amatsize, precsize, descsize
+  integer(psb_epk_) :: system_size
   integer            :: idim, nlev, ierr, ircode
   real(psb_dpk_) :: resmx, resmxp
   real(psb_dpk_) :: t1, t2, tprec
@@ -189,6 +190,7 @@ program mld_dexample_1lev
   amatsize = a%sizeof()
   descsize = desc_a%sizeof()
   precsize = p%sizeof()
+  system_size = desc_a%get_global_rows()
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
@@ -199,6 +201,7 @@ program mld_dexample_1lev
     write(*,'(" ")')
     write(*,'("Matrix from PDE example")')
     write(*,'("Computed solution on ",i8," processors")')np
+    write(*,'("Linear system size        : ",i12)') system_size
     write(*,'("Krylov method             : ",a)') kmethod
     write(*,'("Iterations to convergence : ",i6)')iter
     write(*,'("Error estimate on exit    : ",es12.5)')err
diff --git a/examples/pdegen/mld_dexample_ml.f90 b/examples/pdegen/mld_dexample_ml.f90
index 5a5b53c2..549e598e 100644
--- a/examples/pdegen/mld_dexample_ml.f90
+++ b/examples/pdegen/mld_dexample_ml.f90
@@ -109,6 +109,7 @@ program mld_dexample_ml
   integer            :: choice       
   integer            :: i,info,j
   integer(psb_epk_) :: amatsize, precsize, descsize
+  integer(psb_epk_) :: system_size
   integer            :: idim, ierr, ircode
   real(psb_dpk_)     :: resmx, resmxp
   real(psb_dpk_)     :: t1, t2, tprec
@@ -240,6 +241,7 @@ program mld_dexample_ml
   amatsize = a%sizeof()
   descsize = desc_a%sizeof()
   precsize = p%sizeof()
+  system_size = desc_a%get_global_rows()
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
@@ -250,6 +252,7 @@ program mld_dexample_ml
     write(*,'(" ")')
     write(*,'("Matrix from PDE example")')
     write(*,'("Computed solution on ",i8," processors")')np
+    write(*,'("Linear system size        : ",i12)') system_size
     write(*,'("Krylov method             : ",a)') kmethod
     write(*,'("Iterations to convergence : ",i6)')iter
     write(*,'("Error estimate on exit    : ",es12.5)')err
diff --git a/examples/pdegen/mld_sexample_1lev.f90 b/examples/pdegen/mld_sexample_1lev.f90
index 024f698b..866fe11a 100644
--- a/examples/pdegen/mld_sexample_1lev.f90
+++ b/examples/pdegen/mld_sexample_1lev.f90
@@ -90,6 +90,7 @@ program mld_sexample_1lev
   ! other variables
   integer            :: i,info,j
   integer(psb_epk_) :: amatsize, precsize, descsize
+  integer(psb_epk_) :: system_size
   integer            :: idim, nlev, ierr, ircode
   real(psb_spk_) :: resmx, resmxp
   real(psb_dpk_) :: t1, t2, tprec
@@ -189,6 +190,7 @@ program mld_sexample_1lev
   amatsize = a%sizeof()
   descsize = desc_a%sizeof()
   precsize = p%sizeof()
+  system_size = desc_a%get_global_rows()
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
@@ -199,6 +201,7 @@ program mld_sexample_1lev
     write(*,'(" ")')
     write(*,'("Matrix from PDE example")')
     write(*,'("Computed solution on ",i8," processors")')np
+    write(*,'("Linear system size        : ",i12)') system_size
     write(*,'("Krylov method             : ",a)') kmethod
     write(*,'("Iterations to convergence : ",i6)')iter
     write(*,'("Error estimate on exit    : ",es12.5)')err
diff --git a/examples/pdegen/mld_sexample_ml.f90 b/examples/pdegen/mld_sexample_ml.f90
index 22a04d1d..f32e8eb8 100644
--- a/examples/pdegen/mld_sexample_ml.f90
+++ b/examples/pdegen/mld_sexample_ml.f90
@@ -109,6 +109,7 @@ program mld_sexample_ml
   integer            :: choice       
   integer            :: i,info,j
   integer(psb_epk_) :: amatsize, precsize, descsize
+  integer(psb_epk_) :: system_size
   integer            :: idim, ierr, ircode
   real(psb_spk_)     :: resmx, resmxp
   real(psb_dpk_)     :: t1, t2, tprec
@@ -240,6 +241,7 @@ program mld_sexample_ml
   amatsize = a%sizeof()
   descsize = desc_a%sizeof()
   precsize = p%sizeof()
+  system_size = desc_a%get_global_rows()
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
@@ -250,6 +252,7 @@ program mld_sexample_ml
     write(*,'(" ")')
     write(*,'("Matrix from PDE example")')
     write(*,'("Computed solution on ",i8," processors")')np
+    write(*,'("Linear system size        : ",i12)') system_size
     write(*,'("Krylov method             : ",a)') kmethod
     write(*,'("Iterations to convergence : ",i6)')iter
     write(*,'("Error estimate on exit    : ",es12.5)')err
diff --git a/tests/pdegen/mld_d_pde2d.f90 b/tests/pdegen/mld_d_pde2d.f90
index 908fada1..2f8bd87b 100644
--- a/tests/pdegen/mld_d_pde2d.f90
+++ b/tests/pdegen/mld_d_pde2d.f90
@@ -548,6 +548,7 @@ program mld_d_pde2d
   character(len=20) :: kmethd, ptype
   character(len=5)  :: afmt
   integer(psb_ipk_) :: idim
+  integer(psb_epk_) :: system_size
 
   ! miscellaneous 
   real(psb_dpk_) :: t1, t2, tprec, thier, tslv
@@ -845,12 +846,14 @@ program mld_d_pde2d
   amatsize = a%sizeof()
   descsize = desc_a%sizeof()
   precsize = prec%sizeof()
+  system_size = desc_a%get_global_rows()
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
   call prec%descr(iout=psb_out_unit)
   if (iam == psb_root_) then 
     write(psb_out_unit,'("Computed solution on ",i8," processors")')   np
+    write(psb_out_unit,'("Linear system size                 : ",i12)') system_size
     write(psb_out_unit,'("Krylov method                      : ",a)')  trim(s_choice%kmethd)
     write(psb_out_unit,'("Preconditioner                     : ",a)')  trim(p_choice%descr)
     write(psb_out_unit,'("Iterations to convergence          : ",i12)')   iter
diff --git a/tests/pdegen/mld_d_pde3d.f90 b/tests/pdegen/mld_d_pde3d.f90
index 47a5da7b..9f38ee50 100644
--- a/tests/pdegen/mld_d_pde3d.f90
+++ b/tests/pdegen/mld_d_pde3d.f90
@@ -588,6 +588,7 @@ program mld_d_pde3d
   character(len=20) :: kmethd, ptype
   character(len=5)  :: afmt
   integer(psb_ipk_) :: idim
+  integer(psb_epk_) :: system_size
 
   ! miscellaneous 
   real(psb_dpk_) :: t1, t2, tprec, thier, tslv
@@ -886,12 +887,14 @@ program mld_d_pde3d
   amatsize = a%sizeof()
   descsize = desc_a%sizeof()
   precsize = prec%sizeof()
+  system_size = desc_a%get_global_rows()
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
   call prec%descr(iout=psb_out_unit)
   if (iam == psb_root_) then 
     write(psb_out_unit,'("Computed solution on ",i8," processors")')  np
+    write(psb_out_unit,'("Linear system size                 : ",i12)') system_size
     write(psb_out_unit,'("Krylov method                      : ",a)') trim(s_choice%kmethd)
     write(psb_out_unit,'("Preconditioner                     : ",a)') trim(p_choice%descr)
     write(psb_out_unit,'("Iterations to convergence          : ",i12)')    iter
diff --git a/tests/pdegen/mld_s_pde2d.f90 b/tests/pdegen/mld_s_pde2d.f90
index 0d825b67..13aa5904 100644
--- a/tests/pdegen/mld_s_pde2d.f90
+++ b/tests/pdegen/mld_s_pde2d.f90
@@ -548,6 +548,7 @@ program mld_s_pde2d
   character(len=20) :: kmethd, ptype
   character(len=5)  :: afmt
   integer(psb_ipk_) :: idim
+  integer(psb_epk_) :: system_size
 
   ! miscellaneous 
   real(psb_dpk_) :: t1, t2, tprec, thier, tslv
@@ -845,12 +846,14 @@ program mld_s_pde2d
   amatsize = a%sizeof()
   descsize = desc_a%sizeof()
   precsize = prec%sizeof()
+  system_size = desc_a%get_global_rows()
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
   call prec%descr(iout=psb_out_unit)
   if (iam == psb_root_) then 
     write(psb_out_unit,'("Computed solution on ",i8," processors")')   np
+    write(psb_out_unit,'("Linear system size                 : ",i12)') system_size
     write(psb_out_unit,'("Krylov method                      : ",a)')  trim(s_choice%kmethd)
     write(psb_out_unit,'("Preconditioner                     : ",a)')  trim(p_choice%descr)
     write(psb_out_unit,'("Iterations to convergence          : ",i12)')   iter
diff --git a/tests/pdegen/mld_s_pde3d.f90 b/tests/pdegen/mld_s_pde3d.f90
index 1ac05b44..5e41c784 100644
--- a/tests/pdegen/mld_s_pde3d.f90
+++ b/tests/pdegen/mld_s_pde3d.f90
@@ -588,6 +588,7 @@ program mld_s_pde3d
   character(len=20) :: kmethd, ptype
   character(len=5)  :: afmt
   integer(psb_ipk_) :: idim
+  integer(psb_epk_) :: system_size
 
   ! miscellaneous 
   real(psb_dpk_) :: t1, t2, tprec, thier, tslv
@@ -886,12 +887,14 @@ program mld_s_pde3d
   amatsize = a%sizeof()
   descsize = desc_a%sizeof()
   precsize = prec%sizeof()
+  system_size = desc_a%get_global_rows()
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
   call prec%descr(iout=psb_out_unit)
   if (iam == psb_root_) then 
     write(psb_out_unit,'("Computed solution on ",i8," processors")')  np
+    write(psb_out_unit,'("Linear system size                 : ",i12)') system_size
     write(psb_out_unit,'("Krylov method                      : ",a)') trim(s_choice%kmethd)
     write(psb_out_unit,'("Preconditioner                     : ",a)') trim(p_choice%descr)
     write(psb_out_unit,'("Iterations to convergence          : ",i12)')    iter