diff --git a/base/serial/psi_c_serial_impl.f90 b/base/serial/psi_c_serial_impl.f90
index 13795f4a..05924f59 100644
--- a/base/serial/psi_c_serial_impl.f90
+++ b/base/serial/psi_c_serial_impl.f90
@@ -113,7 +113,7 @@ subroutine psi_cgthmv(n,k,idx,alpha,x,beta,y)
   use psb_const_mod
   implicit none
 
-  integer(psb_ipk_) :: n, idx(:)
+  integer(psb_ipk_) :: n, k, idx(:)
   complex(psb_spk_) :: x(:,:), y(:),alpha,beta
 
   ! Locals
diff --git a/base/serial/psi_d_serial_impl.f90 b/base/serial/psi_d_serial_impl.f90
index 490ad385..5e911cbe 100644
--- a/base/serial/psi_d_serial_impl.f90
+++ b/base/serial/psi_d_serial_impl.f90
@@ -113,7 +113,7 @@ subroutine psi_dgthmv(n,k,idx,alpha,x,beta,y)
   use psb_const_mod
   implicit none
 
-  integer(psb_ipk_) :: n, idx(:)
+  integer(psb_ipk_) :: n, k, idx(:)
   real(psb_dpk_) :: x(:,:), y(:),alpha,beta
 
   ! Locals
diff --git a/base/serial/psi_i_serial_impl.f90 b/base/serial/psi_i_serial_impl.f90
index 3c3bc1c6..135c57cf 100644
--- a/base/serial/psi_i_serial_impl.f90
+++ b/base/serial/psi_i_serial_impl.f90
@@ -113,7 +113,7 @@ subroutine psi_igthmv(n,k,idx,alpha,x,beta,y)
   use psb_const_mod
   implicit none
 
-  integer(psb_ipk_) :: n, idx(:)
+  integer(psb_ipk_) :: n, k, idx(:)
   integer(psb_ipk_) :: x(:,:), y(:),alpha,beta
 
   ! Locals
diff --git a/base/serial/psi_s_serial_impl.f90 b/base/serial/psi_s_serial_impl.f90
index 287520a7..5e6fbac3 100644
--- a/base/serial/psi_s_serial_impl.f90
+++ b/base/serial/psi_s_serial_impl.f90
@@ -113,7 +113,7 @@ subroutine psi_sgthmv(n,k,idx,alpha,x,beta,y)
   use psb_const_mod
   implicit none
 
-  integer(psb_ipk_) :: n, idx(:)
+  integer(psb_ipk_) :: n, k, idx(:)
   real(psb_spk_) :: x(:,:), y(:),alpha,beta
 
   ! Locals
diff --git a/base/serial/psi_z_serial_impl.f90 b/base/serial/psi_z_serial_impl.f90
index 259c04aa..6b4c44ee 100644
--- a/base/serial/psi_z_serial_impl.f90
+++ b/base/serial/psi_z_serial_impl.f90
@@ -113,7 +113,7 @@ subroutine psi_zgthmv(n,k,idx,alpha,x,beta,y)
   use psb_const_mod
   implicit none
 
-  integer(psb_ipk_) :: n, idx(:)
+  integer(psb_ipk_) :: n, k, idx(:)
   complex(psb_dpk_) :: x(:,:), y(:),alpha,beta
 
   ! Locals