diff --git a/base/modules/psb_c_base_mat_mod.f90 b/base/modules/psb_c_base_mat_mod.f90
index b66a66bd..56049a1f 100644
--- a/base/modules/psb_c_base_mat_mod.f90
+++ b/base/modules/psb_c_base_mat_mod.f90
@@ -282,7 +282,7 @@ module psb_c_base_mat_mod
   !!  \param info  return code
   !!  \param jmin [1] minimum col index 
   !!  \param jmax [a\%get_ncols()] maximum col index 
-  !!  \param iren(:) [none] an array to return renumbered indices   (iren(ia(:)),iren(ja(:))
+  !!  \param iren(:) [none] an array to return renumbered indices (iren(ia(:)),iren(ja(:))
   !!  \param rscale [false] map [min(ia(:)):max(ia(:))] onto [1:max(ia(:))-min(ia(:))+1]
   !!  \param cscale [false] map [min(ja(:)):max(ja(:))] onto [1:max(ja(:))-min(ja(:))+1]
   !!          ( iren cannot be specified with rscale/cscale)
@@ -681,6 +681,10 @@ module psb_c_base_mat_mod
   !!
   !!        Compute
   !!           Y = alpha*op(A)*X + beta*Y
+  !!        Usually the unwrapping of the encapsulated vector is done
+  !!        here, so that all the derived classes need only the
+  !!        versions with the standard arrays. 
+  !!
   !!
   !! \param alpha  Scaling factor for Ax
   !! \param A      the input sparse matrix
@@ -704,13 +708,6 @@ module psb_c_base_mat_mod
     end subroutine psb_c_base_vect_mv
   end interface
   
-  !
-  ! Triangular system solve.
-  ! The CSSM/CSSV/VECT_SV outer methods are implemented at the base
-  ! level, and they take care of the SCALE and D control arguments.
-  ! So the derived classes need to override only the INNER_ methods.
-  !
-  
   !
   !> Function  cssm:
   !! \memberof  psb_c_base_sparse_mat
diff --git a/base/modules/psb_d_base_mat_mod.f90 b/base/modules/psb_d_base_mat_mod.f90
index 36738bbd..17ec36df 100644
--- a/base/modules/psb_d_base_mat_mod.f90
+++ b/base/modules/psb_d_base_mat_mod.f90
@@ -282,7 +282,7 @@ module psb_d_base_mat_mod
   !!  \param info  return code
   !!  \param jmin [1] minimum col index 
   !!  \param jmax [a\%get_ncols()] maximum col index 
-  !!  \param iren(:) [none] an array to return renumbered indices   (iren(ia(:)),iren(ja(:))
+  !!  \param iren(:) [none] an array to return renumbered indices (iren(ia(:)),iren(ja(:))
   !!  \param rscale [false] map [min(ia(:)):max(ia(:))] onto [1:max(ia(:))-min(ia(:))+1]
   !!  \param cscale [false] map [min(ja(:)):max(ja(:))] onto [1:max(ja(:))-min(ja(:))+1]
   !!          ( iren cannot be specified with rscale/cscale)
@@ -681,6 +681,10 @@ module psb_d_base_mat_mod
   !!
   !!        Compute
   !!           Y = alpha*op(A)*X + beta*Y
+  !!        Usually the unwrapping of the encapsulated vector is done
+  !!        here, so that all the derived classes need only the
+  !!        versions with the standard arrays. 
+  !!
   !!
   !! \param alpha  Scaling factor for Ax
   !! \param A      the input sparse matrix
@@ -704,13 +708,6 @@ module psb_d_base_mat_mod
     end subroutine psb_d_base_vect_mv
   end interface
   
-  !
-  ! Triangular system solve.
-  ! The CSSM/CSSV/VECT_SV outer methods are implemented at the base
-  ! level, and they take care of the SCALE and D control arguments.
-  ! So the derived classes need to override only the INNER_ methods.
-  !
-  
   !
   !> Function  cssm:
   !! \memberof  psb_d_base_sparse_mat
diff --git a/base/modules/psb_s_base_mat_mod.f90 b/base/modules/psb_s_base_mat_mod.f90
index 715d80f5..81b90e69 100644
--- a/base/modules/psb_s_base_mat_mod.f90
+++ b/base/modules/psb_s_base_mat_mod.f90
@@ -282,7 +282,7 @@ module psb_s_base_mat_mod
   !!  \param info  return code
   !!  \param jmin [1] minimum col index 
   !!  \param jmax [a\%get_ncols()] maximum col index 
-  !!  \param iren(:) [none] an array to return renumbered indices   (iren(ia(:)),iren(ja(:))
+  !!  \param iren(:) [none] an array to return renumbered indices (iren(ia(:)),iren(ja(:))
   !!  \param rscale [false] map [min(ia(:)):max(ia(:))] onto [1:max(ia(:))-min(ia(:))+1]
   !!  \param cscale [false] map [min(ja(:)):max(ja(:))] onto [1:max(ja(:))-min(ja(:))+1]
   !!          ( iren cannot be specified with rscale/cscale)
@@ -681,6 +681,10 @@ module psb_s_base_mat_mod
   !!
   !!        Compute
   !!           Y = alpha*op(A)*X + beta*Y
+  !!        Usually the unwrapping of the encapsulated vector is done
+  !!        here, so that all the derived classes need only the
+  !!        versions with the standard arrays. 
+  !!
   !!
   !! \param alpha  Scaling factor for Ax
   !! \param A      the input sparse matrix
@@ -704,13 +708,6 @@ module psb_s_base_mat_mod
     end subroutine psb_s_base_vect_mv
   end interface
   
-  !
-  ! Triangular system solve.
-  ! The CSSM/CSSV/VECT_SV outer methods are implemented at the base
-  ! level, and they take care of the SCALE and D control arguments.
-  ! So the derived classes need to override only the INNER_ methods.
-  !
-  
   !
   !> Function  cssm:
   !! \memberof  psb_s_base_sparse_mat
diff --git a/base/modules/psb_z_base_mat_mod.f90 b/base/modules/psb_z_base_mat_mod.f90
index 08b8c274..4ed514fb 100644
--- a/base/modules/psb_z_base_mat_mod.f90
+++ b/base/modules/psb_z_base_mat_mod.f90
@@ -282,7 +282,7 @@ module psb_z_base_mat_mod
   !!  \param info  return code
   !!  \param jmin [1] minimum col index 
   !!  \param jmax [a\%get_ncols()] maximum col index 
-  !!  \param iren(:) [none] an array to return renumbered indices   (iren(ia(:)),iren(ja(:))
+  !!  \param iren(:) [none] an array to return renumbered indices (iren(ia(:)),iren(ja(:))
   !!  \param rscale [false] map [min(ia(:)):max(ia(:))] onto [1:max(ia(:))-min(ia(:))+1]
   !!  \param cscale [false] map [min(ja(:)):max(ja(:))] onto [1:max(ja(:))-min(ja(:))+1]
   !!          ( iren cannot be specified with rscale/cscale)
@@ -681,6 +681,10 @@ module psb_z_base_mat_mod
   !!
   !!        Compute
   !!           Y = alpha*op(A)*X + beta*Y
+  !!        Usually the unwrapping of the encapsulated vector is done
+  !!        here, so that all the derived classes need only the
+  !!        versions with the standard arrays. 
+  !!
   !!
   !! \param alpha  Scaling factor for Ax
   !! \param A      the input sparse matrix
@@ -704,13 +708,6 @@ module psb_z_base_mat_mod
     end subroutine psb_z_base_vect_mv
   end interface
   
-  !
-  ! Triangular system solve.
-  ! The CSSM/CSSV/VECT_SV outer methods are implemented at the base
-  ! level, and they take care of the SCALE and D control arguments.
-  ! So the derived classes need to override only the INNER_ methods.
-  !
-  
   !
   !> Function  cssm:
   !! \memberof  psb_z_base_sparse_mat