diff --git a/test/fileread/cf_sample.f90 b/test/fileread/cf_sample.f90
index 625873f1..2a777924 100644
--- a/test/fileread/cf_sample.f90
+++ b/test/fileread/cf_sample.f90
@@ -273,7 +273,7 @@ program cf_sample
     write(psb_out_unit,'("Iterations to convergence: ",i6)')iter
     write(psb_out_unit,'("Error estimate on exit   : ",es12.5)') err
     write(psb_out_unit,'("Time to buil prec.       : ",es12.5)')tprec
-    write(psb_out_unit,'("Time to solve matrix     : ",es12.5)')t2
+    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,'("Total time               : ",es12.5)')t2+tprec
     write(psb_out_unit,'("Residual norm 2          : ",es12.5)')resmx
diff --git a/test/fileread/df_sample.f90 b/test/fileread/df_sample.f90
index 19461566..bfab6928 100644
--- a/test/fileread/df_sample.f90
+++ b/test/fileread/df_sample.f90
@@ -279,7 +279,7 @@ program df_sample
     write(psb_out_unit,'("Iterations to convergence: ",i6)')iter
     write(psb_out_unit,'("Error estimate on exit   : ",es12.5)') err
     write(psb_out_unit,'("Time to buil prec.       : ",es12.5)')tprec
-    write(psb_out_unit,'("Time to solve matrix     : ",es12.5)')t2
+    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,'("Total time               : ",es12.5)')t2+tprec
     write(psb_out_unit,'("Residual norm 2          : ",es12.5)')resmx
diff --git a/test/fileread/sf_sample.f90 b/test/fileread/sf_sample.f90
index aff706d5..39fba0d7 100644
--- a/test/fileread/sf_sample.f90
+++ b/test/fileread/sf_sample.f90
@@ -276,7 +276,7 @@ program sf_sample
     write(psb_out_unit,'("Iterations to convergence: ",i6)')iter
     write(psb_out_unit,'("Error estimate on exit   : ",es12.5)') err
     write(psb_out_unit,'("Time to buil prec.       : ",es12.5)')tprec
-    write(psb_out_unit,'("Time to solve matrix     : ",es12.5)')t2
+    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,'("Total time               : ",es12.5)')t2+tprec
     write(psb_out_unit,'("Residual norm 2          : ",es12.5)')resmx
diff --git a/test/fileread/zf_sample.f90 b/test/fileread/zf_sample.f90
index c8b312ba..994a7b2e 100644
--- a/test/fileread/zf_sample.f90
+++ b/test/fileread/zf_sample.f90
@@ -273,7 +273,7 @@ program zf_sample
     write(psb_out_unit,'("Iterations to convergence: ",i6)')iter
     write(psb_out_unit,'("Error estimate on exit   : ",es12.5)') err
     write(psb_out_unit,'("Time to buil prec.       : ",es12.5)')tprec
-    write(psb_out_unit,'("Time to solve matrix     : ",es12.5)')t2
+    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,'("Total time               : ",es12.5)')t2+tprec
     write(psb_out_unit,'("Residual norm 2          : ",es12.5)')resmx
diff --git a/test/pargen/ppde2d.f90 b/test/pargen/ppde2d.f90
index 922c4925..f0feccbb 100644
--- a/test/pargen/ppde2d.f90
+++ b/test/pargen/ppde2d.f90
@@ -235,7 +235,7 @@ program ppde2d
 
   if (iam == psb_root_) then
     write(psb_out_unit,'(" ")')
-    write(psb_out_unit,'("Time to solve matrix          : ",es12.5)')t2
+    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
     write(psb_out_unit,'("Convergence indicator on exit : ",es12.5)')err
diff --git a/test/pargen/ppde3d.f90 b/test/pargen/ppde3d.f90
index a8ad8f1d..2250a704 100644
--- a/test/pargen/ppde3d.f90
+++ b/test/pargen/ppde3d.f90
@@ -247,7 +247,7 @@ program ppde3d
 
   if (iam == psb_root_) then
     write(psb_out_unit,'(" ")')
-    write(psb_out_unit,'("Time to solve matrix          : ",es12.5)')t2
+    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
     write(psb_out_unit,'("Convergence indicator on exit : ",es12.5)')err
diff --git a/test/pargen/spde2d.f90 b/test/pargen/spde2d.f90
index ffa5d60d..c777f153 100644
--- a/test/pargen/spde2d.f90
+++ b/test/pargen/spde2d.f90
@@ -234,7 +234,7 @@ program spde2d
 
   if (iam == psb_root_) then
     write(psb_out_unit,'(" ")')
-    write(psb_out_unit,'("Time to solve matrix          : ",es12.5)')t2
+    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
     write(psb_out_unit,'("Convergence indicator on exit : ",es12.5)')err
diff --git a/test/pargen/spde3d.f90 b/test/pargen/spde3d.f90
index 554ca997..0f04aeca 100644
--- a/test/pargen/spde3d.f90
+++ b/test/pargen/spde3d.f90
@@ -248,7 +248,7 @@ program spde3d
 
   if (iam == psb_root_) then
     write(psb_out_unit,'(" ")')
-    write(psb_out_unit,'("Time to solve matrix          : ",es12.5)')t2
+    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
     write(psb_out_unit,'("Convergence indicator on exit : ",es12.5)')err