diff --git a/test/pargen/psb_d_pde2d.f90 b/test/pargen/psb_d_pde2d.f90
index 8440281b..94bd77c4 100644
--- a/test/pargen/psb_d_pde2d.f90
+++ b/test/pargen/psb_d_pde2d.f90
@@ -537,7 +537,7 @@ program psb_d_pde2d
   ! input parameters
   character(len=20) :: kmethd, ptype
   character(len=5)  :: afmt
-  integer(psb_ipk_) :: idim
+  integer(psb_ipk_) :: idim, system_size
 
   ! miscellaneous 
   real(psb_dpk_), parameter :: one = done
@@ -654,10 +654,11 @@ program psb_d_pde2d
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
-
+  system_size = desc_a%get_global_rows()
   if (iam == psb_root_) then
     write(psb_out_unit,'(" ")')
     write(psb_out_unit,'("Number of processes           : ",i0)')np
+    write(psb_out_unit,'("Linear system size            : ",i0)') system_size
     write(psb_out_unit,'("Time to solve system          : ",es12.5)')t2
     write(psb_out_unit,'("Time per iteration            : ",es12.5)')t2/iter
     write(psb_out_unit,'("Number of iterations          : ",i0)')iter
diff --git a/test/pargen/psb_d_pde3d.f90 b/test/pargen/psb_d_pde3d.f90
index a34f462e..5f723f32 100644
--- a/test/pargen/psb_d_pde3d.f90
+++ b/test/pargen/psb_d_pde3d.f90
@@ -577,7 +577,7 @@ program psb_d_pde3d
   ! input parameters
   character(len=20) :: kmethd, ptype
   character(len=5)  :: afmt
-  integer(psb_ipk_) :: idim
+  integer(psb_ipk_) :: idim, system_size
 
   ! miscellaneous 
   real(psb_dpk_), parameter :: one = done
@@ -695,10 +695,11 @@ program psb_d_pde3d
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
-
+  system_size = desc_a%get_global_rows()
   if (iam == psb_root_) then
     write(psb_out_unit,'(" ")')
     write(psb_out_unit,'("Number of processes           : ",i0)')np
+    write(psb_out_unit,'("Linear system size            : ",i0)') system_size
     write(psb_out_unit,'("Time to solve system          : ",es12.5)')t2
     write(psb_out_unit,'("Time per iteration            : ",es12.5)')t2/iter
     write(psb_out_unit,'("Number of iterations          : ",i0)')iter
diff --git a/test/pargen/psb_s_pde2d.f90 b/test/pargen/psb_s_pde2d.f90
index ab33a8ba..1f6ea8eb 100644
--- a/test/pargen/psb_s_pde2d.f90
+++ b/test/pargen/psb_s_pde2d.f90
@@ -537,7 +537,7 @@ program psb_s_pde2d
   ! input parameters
   character(len=20) :: kmethd, ptype
   character(len=5)  :: afmt
-  integer(psb_ipk_) :: idim
+  integer(psb_ipk_) :: idim, system_size
 
   ! miscellaneous 
   real(psb_spk_), parameter :: one = sone
@@ -654,10 +654,11 @@ program psb_s_pde2d
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
-
+  system_size = desc_a%get_global_rows()
   if (iam == psb_root_) then
     write(psb_out_unit,'(" ")')
     write(psb_out_unit,'("Number of processes           : ",i0)')np
+    write(psb_out_unit,'("Linear system size            : ",i0)') system_size
     write(psb_out_unit,'("Time to solve system          : ",es12.5)')t2
     write(psb_out_unit,'("Time per iteration            : ",es12.5)')t2/iter
     write(psb_out_unit,'("Number of iterations          : ",i0)')iter
diff --git a/test/pargen/psb_s_pde3d.f90 b/test/pargen/psb_s_pde3d.f90
index 32c47593..a74d414b 100644
--- a/test/pargen/psb_s_pde3d.f90
+++ b/test/pargen/psb_s_pde3d.f90
@@ -577,7 +577,7 @@ program psb_s_pde3d
   ! input parameters
   character(len=20) :: kmethd, ptype
   character(len=5)  :: afmt
-  integer(psb_ipk_) :: idim
+  integer(psb_ipk_) :: idim, system_size
 
   ! miscellaneous 
   real(psb_spk_), parameter :: one = sone
@@ -695,10 +695,11 @@ program psb_s_pde3d
   call psb_sum(ictxt,amatsize)
   call psb_sum(ictxt,descsize)
   call psb_sum(ictxt,precsize)
-
+  system_size = desc_a%get_global_rows()
   if (iam == psb_root_) then
     write(psb_out_unit,'(" ")')
     write(psb_out_unit,'("Number of processes           : ",i0)')np
+    write(psb_out_unit,'("Linear system size            : ",i0)') system_size
     write(psb_out_unit,'("Time to solve system          : ",es12.5)')t2
     write(psb_out_unit,'("Time per iteration            : ",es12.5)')t2/iter
     write(psb_out_unit,'("Number of iterations          : ",i0)')iter