diff --git a/base/modules/auxil/psi_i_serial_mod.f90 b/base/modules/auxil/psi_i_serial_mod.f90
index 894dcba9..b7802b00 100644
--- a/base/modules/auxil/psi_i_serial_mod.f90
+++ b/base/modules/auxil/psi_i_serial_mod.f90
@@ -30,12 +30,12 @@
 !   
 !    
 module psi_i_serial_mod
-  use psb_const_mod, only : psb_ipk_
+  use psb_const_mod, only : psb_ipk_, psb_lpk_
 
   interface psb_gelp
     ! 2-D version
     subroutine psb_igelp(trans,iperm,x,info)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none 
       integer(psb_ipk_), intent(inout)     ::  x(:,:)
       integer(psb_ipk_), intent(in)      ::  iperm(:)
@@ -43,7 +43,7 @@ module psi_i_serial_mod
       character, intent(in)              :: trans
     end subroutine psb_igelp
     subroutine psb_igelpv(trans,iperm,x,info)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none 
       integer(psb_ipk_), intent(inout)     ::  x(:)
       integer(psb_ipk_), intent(in)      ::  iperm(:)
@@ -54,7 +54,7 @@ module psi_i_serial_mod
   
   interface psb_geaxpby 
     subroutine psi_iaxpby(m,n,alpha, x, beta, y, info)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none 
       integer(psb_ipk_), intent(in)      :: m, n
       integer(psb_ipk_), intent (in)       ::  x(:,:)
@@ -63,7 +63,7 @@ module psi_i_serial_mod
       integer(psb_ipk_), intent(out)     :: info
     end subroutine psi_iaxpby
     subroutine psi_iaxpbyv(m,alpha, x, beta, y, info)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none 
       integer(psb_ipk_), intent(in)      :: m
       integer(psb_ipk_), intent (in)       ::  x(:)
@@ -75,33 +75,33 @@ module psi_i_serial_mod
 
   interface psi_gth
     subroutine psi_igthmv(n,k,idx,alpha,x,beta,y)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, k, idx(:)
       integer(psb_ipk_) :: x(:,:), y(:),alpha,beta
     end subroutine psi_igthmv
     subroutine psi_igthv(n,idx,alpha,x,beta,y)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, idx(:)
       integer(psb_ipk_) :: x(:), y(:),alpha,beta
     end subroutine psi_igthv
     subroutine psi_igthzmv(n,k,idx,x,y)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, k, idx(:)
       integer(psb_ipk_) :: x(:,:), y(:)
       
     end subroutine psi_igthzmv
     subroutine psi_igthzmm(n,k,idx,x,y)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, k, idx(:)
       integer(psb_ipk_) :: x(:,:), y(:,:)
       
     end subroutine psi_igthzmm
     subroutine psi_igthzv(n,idx,x,y)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none 
       integer(psb_ipk_) :: n, idx(:)
       integer(psb_ipk_) :: x(:), y(:)
@@ -110,19 +110,19 @@ module psi_i_serial_mod
 
   interface psi_sct
     subroutine psi_isctmm(n,k,idx,x,beta,y)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, k, idx(:)
       integer(psb_ipk_) :: beta, x(:,:), y(:,:)
     end subroutine psi_isctmm
     subroutine psi_isctmv(n,k,idx,x,beta,y)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, k, idx(:)
       integer(psb_ipk_) :: beta, x(:), y(:,:)
     end subroutine psi_isctmv
     subroutine psi_isctv(n,idx,x,beta,y)
-      import :: psb_ipk_
+      import :: psb_ipk_, psb_lpk_
       implicit none
       
       integer(psb_ipk_) :: n, idx(:)
diff --git a/base/modules/auxil/psi_l_serial_mod.f90 b/base/modules/auxil/psi_l_serial_mod.f90
index 7c35f309..e5e29c3f 100644
--- a/base/modules/auxil/psi_l_serial_mod.f90
+++ b/base/modules/auxil/psi_l_serial_mod.f90
@@ -30,13 +30,11 @@
 !   
 !    
 module psi_l_serial_mod
-  use psb_const_mod, only : psb_ipk_
   use psb_const_mod, only : psb_ipk_, psb_lpk_
 
   interface psb_gelp
     ! 2-D version
     subroutine psb_lgelp(trans,iperm,x,info)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none 
       integer(psb_lpk_), intent(inout)     ::  x(:,:)
@@ -45,7 +43,6 @@ module psi_l_serial_mod
       character, intent(in)              :: trans
     end subroutine psb_lgelp
     subroutine psb_lgelpv(trans,iperm,x,info)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none 
       integer(psb_lpk_), intent(inout)     ::  x(:)
@@ -57,7 +54,6 @@ module psi_l_serial_mod
   
   interface psb_geaxpby 
     subroutine psi_laxpby(m,n,alpha, x, beta, y, info)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none 
       integer(psb_ipk_), intent(in)      :: m, n
@@ -67,7 +63,6 @@ module psi_l_serial_mod
       integer(psb_ipk_), intent(out)     :: info
     end subroutine psi_laxpby
     subroutine psi_laxpbyv(m,alpha, x, beta, y, info)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none 
       integer(psb_ipk_), intent(in)      :: m
@@ -80,21 +75,18 @@ module psi_l_serial_mod
 
   interface psi_gth
     subroutine psi_lgthmv(n,k,idx,alpha,x,beta,y)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, k, idx(:)
       integer(psb_lpk_) :: x(:,:), y(:),alpha,beta
     end subroutine psi_lgthmv
     subroutine psi_lgthv(n,idx,alpha,x,beta,y)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, idx(:)
       integer(psb_lpk_) :: x(:), y(:),alpha,beta
     end subroutine psi_lgthv
     subroutine psi_lgthzmv(n,k,idx,x,y)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, k, idx(:)
@@ -102,7 +94,6 @@ module psi_l_serial_mod
       
     end subroutine psi_lgthzmv
     subroutine psi_lgthzmm(n,k,idx,x,y)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, k, idx(:)
@@ -110,7 +101,6 @@ module psi_l_serial_mod
       
     end subroutine psi_lgthzmm
     subroutine psi_lgthzv(n,idx,x,y)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none 
       integer(psb_ipk_) :: n, idx(:)
@@ -120,21 +110,18 @@ module psi_l_serial_mod
 
   interface psi_sct
     subroutine psi_lsctmm(n,k,idx,x,beta,y)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, k, idx(:)
       integer(psb_lpk_) :: beta, x(:,:), y(:,:)
     end subroutine psi_lsctmm
     subroutine psi_lsctmv(n,k,idx,x,beta,y)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none
       integer(psb_ipk_) :: n, k, idx(:)
       integer(psb_lpk_) :: beta, x(:), y(:,:)
     end subroutine psi_lsctmv
     subroutine psi_lsctv(n,idx,x,beta,y)
-      import :: psb_ipk_
       import :: psb_ipk_, psb_lpk_
       implicit none