psblas3:

base/Makefile base/comm/psb_cgather.f90 base/comm/psb_chalo.f90 base/comm/psb_covrl.f90 base/comm/psb_cscatter.F90 base/comm/psb_dgather.f90 base/comm/psb_dhalo.f90 base/comm/psb_dovrl.f90 base/comm/psb_dscatter.F90 base/comm/psb_igather.f90 base/comm/psb_ihalo.f90 base/comm/psb_iovrl.f90 base/comm/psb_iscatter.F90 base/comm/psb_sgather.f90 base/comm/psb_shalo.f90 base/comm/psb_sovrl.f90 base/comm/psb_sscatter.F90 base/comm/psb_zgather.f90 base/comm/psb_zhalo.f90 base/comm/psb_zovrl.f90 base/comm/psb_zscatter.F90 base/modules/Makefile base/modules/psb_serial_mod.f90 base/modules/psb_sparse_mod.f90 base/psblas/psb_camax.f90 base/psblas/psb_casum.f90 base/psblas/psb_caxpby.f90 base/psblas/psb_cdot.f90 base/psblas/psb_cnrmi.f90 base/psblas/psb_cspmm.f90 base/psblas/psb_cspsm.f90 base/psblas/psb_damax.f90 base/psblas/psb_dasum.f90 base/psblas/psb_daxpby.f90 base/psblas/psb_ddot.f90 base/psblas/psb_dspmm.f90 base/psblas/psb_dspnrm1.f90 base/psblas/psb_dspsm.f90 base/psblas/psb_samax.f90 base/psblas/psb_sasum.f90 base/psblas/psb_saxpby.f90 base/psblas/psb_sdot.f90 base/psblas/psb_sspmm.f90 base/psblas/psb_sspsm.f90 base/psblas/psb_sxdot.f90 base/psblas/psb_zamax.f90 base/psblas/psb_zasum.f90 base/psblas/psb_zaxpby.f90 base/psblas/psb_zdot.f90 base/psblas/psb_zspmm.f90 base/psblas/psb_zspsm.f90 base/serial/psb_cnumbmm.f90 base/serial/psb_crwextd.f90 base/serial/psb_csymbmm.f90 base/serial/psb_dnumbmm.f90 base/serial/psb_drwextd.f90 base/serial/psb_dsymbmm.f90 base/serial/psb_snumbmm.f90 base/serial/psb_srwextd.f90 base/serial/psb_ssymbmm.f90 base/serial/psb_znumbmm.f90 base/serial/psb_zrwextd.f90 base/serial/psb_zsymbmm.f90 base/tools/psb_callc.f90 base/tools/psb_casb.f90 base/tools/psb_ccdbldext.F90 base/tools/psb_cd_inloc.f90 base/tools/psb_cd_lstext.f90 base/tools/psb_cd_reinit.f90 base/tools/psb_cd_set_bld.f90 base/tools/psb_cd_switch_ovl_indxmap.f90 base/tools/psb_cdals.f90 base/tools/psb_cdalv.f90 base/tools/psb_cdcpy.F90 base/tools/psb_cdins.f90 base/tools/psb_cdprt.f90 base/tools/psb_cdren.f90 base/tools/psb_cdrep.f90 base/tools/psb_cfree.f90 base/tools/psb_cins.f90 base/tools/psb_cspalloc.f90 base/tools/psb_cspasb.f90 base/tools/psb_cspfree.f90 base/tools/psb_csphalo.F90 base/tools/psb_cspins.f90 base/tools/psb_csprn.f90 base/tools/psb_dallc.f90 base/tools/psb_dasb.f90 base/tools/psb_dcdbldext.F90 base/tools/psb_dfree.f90 base/tools/psb_dins.f90 base/tools/psb_dspalloc.f90 base/tools/psb_dspasb.f90 base/tools/psb_dspfree.f90 base/tools/psb_dsphalo.F90 base/tools/psb_dspins.f90 base/tools/psb_dsprn.f90 base/tools/psb_get_overlap.f90 base/tools/psb_glob_to_loc.f90 base/tools/psb_ialloc.f90 base/tools/psb_iasb.f90 base/tools/psb_icdasb.F90 base/tools/psb_ifree.f90 base/tools/psb_iins.f90 base/tools/psb_linmap.f90 base/tools/psb_loc_to_glob.f90 base/tools/psb_map.f90 base/tools/psb_sallc.f90 base/tools/psb_sasb.f90 base/tools/psb_scdbldext.F90 base/tools/psb_sfree.f90 base/tools/psb_sins.f90 base/tools/psb_sspalloc.f90 base/tools/psb_sspasb.f90 base/tools/psb_sspfree.f90 base/tools/psb_ssphalo.F90 base/tools/psb_sspins.f90 base/tools/psb_ssprn.f90 base/tools/psb_zallc.f90 base/tools/psb_zasb.f90 base/tools/psb_zcdbldext.F90 base/tools/psb_zfree.f90 base/tools/psb_zins.f90 base/tools/psb_zspalloc.f90 base/tools/psb_zspasb.f90 base/tools/psb_zspfree.f90 base/tools/psb_zsphalo.F90 base/tools/psb_zspins.f90 base/tools/psb_zsprn.f90 krylov/Makefile krylov/psb_base_inner_krylov_mod.f90 krylov/psb_c_inner_krylov_mod.f90 krylov/psb_cbicg.f90 krylov/psb_ccg.f90 krylov/psb_ccgs.f90 krylov/psb_ccgstab.f90 krylov/psb_ccgstabl.f90 krylov/psb_ckrylov.f90 krylov/psb_crgmres.f90 krylov/psb_d_inner_krylov_mod.f90 krylov/psb_dbicg.f90 krylov/psb_dcg.F90 krylov/psb_dcgs.f90 krylov/psb_dcgstab.F90 krylov/psb_dcgstabl.f90 krylov/psb_dkrylov.f90 krylov/psb_drgmres.f90 krylov/psb_krylov_mod.f90 krylov/psb_s_inner_krylov_mod.f90 krylov/psb_sbicg.f90 krylov/psb_scg.F90 krylov/psb_scgs.f90 krylov/psb_scgstab.F90 krylov/psb_scgstabl.f90 krylov/psb_skrylov.f90 krylov/psb_srgmres.f90 krylov/psb_z_inner_krylov_mod.f90 krylov/psb_zbicg.f90 krylov/psb_zcg.F90 krylov/psb_zcgs.f90 krylov/psb_zcgstab.f90 krylov/psb_zcgstabl.f90 krylov/psb_zkrylov.f90 krylov/psb_zrgmres.f90 opt/psb_d_ell_impl.f90 opt/psb_d_rsb_mat_mod.F90 prec/Makefile prec/psb_c_base_prec_mod.f90 prec/psb_c_bjacprec.f90 prec/psb_c_diagprec.f90 prec/psb_c_nullprec.f90 prec/psb_c_prec_type.f90 prec/psb_cilu_fct.f90 prec/psb_cprecbld.f90 prec/psb_cprecinit.f90 prec/psb_cprecset.f90 prec/psb_d_base_prec_mod.f90 prec/psb_d_bjacprec.f90 prec/psb_d_diagprec.f90 prec/psb_d_nullprec.f90 prec/psb_d_prec_type.f90 prec/psb_dilu_fct.f90 prec/psb_dprecbld.f90 prec/psb_dprecinit.f90 prec/psb_dprecset.f90 prec/psb_prec_const_mod.f90 prec/psb_prec_mod.f90 prec/psb_s_base_prec_mod.f90 prec/psb_s_bjacprec.f90 prec/psb_s_diagprec.f90 prec/psb_s_nullprec.f90 prec/psb_s_prec_type.f90 prec/psb_silu_fct.f90 prec/psb_sprecbld.f90 prec/psb_sprecinit.f90 prec/psb_sprecset.f90 prec/psb_z_base_prec_mod.f90 prec/psb_z_bjacprec.f90 prec/psb_z_diagprec.f90 prec/psb_z_nullprec.f90 prec/psb_z_prec_type.f90 prec/psb_zilu_fct.f90 prec/psb_zprecbld.f90 prec/psb_zprecinit.f90 prec/psb_zprecset.f90 test/fileread/cf_sample.f90 test/fileread/df_sample.f90 test/fileread/getp.f90 test/fileread/sf_sample.f90 test/fileread/zf_sample.f90 test/newfmt/ppde.F90 test/newfmt/spde.f90 test/pargen/ppde.f90 test/pargen/spde.f90 test/serial/d_coo_matgen.f90 test/serial/d_matgen.f90 test/serial/psb_d_cxx_impl.f90 test/serial/psb_d_cyy_impl.f90 test/torture/psb_c_mvsv_tester.f90 test/torture/psb_d_mvsv_tester.f90 test/torture/psb_s_mvsv_tester.f90 test/torture/psb_z_mvsv_tester.f90 test/torture/psbtf.f90 test/util/dhb2mm.f90 test/util/dmm2hb.f90 test/util/zhb2mm.f90 test/util/zmm2hb.f90 util/psb_hbio_impl.f90 util/psb_hbio_mod.f90 util/psb_mat_dist_impl.f90 util/psb_mat_dist_mod.f90 util/psb_metispart_mod.F90 util/psb_mmio_impl.f90 util/psb_mmio_mod.f90 Changed psb_sparse_mod into psb_base_mod, added psb_csprt. This allows the cbind stuff to be shared between version 2 and version 3.
14 years ago · 324830ef7b
parent 9993b41288
commit 324830ef7b
234 changed files with 1031 additions and 923 deletions
--- a/base/Makefile
+++ b/base/Makefile
@ -3,7 +3,7 @@ include ../Make.inc
 HERE=.
 LIBDIR=../lib
 LIBNAME=$(BASELIBNAME)
-LIBMOD=psb_sparse_mod$(.mod)
+LIBMOD=psb_base_mod$(.mod)
 lib: mods sr cm in pb tl
 	/bin/cp -p $(HERE)/$(LIBNAME) $(LIBDIR)
 	/bin/cp -p $(LIBMOD) *$(.mod) $(LIBDIR)
--- a/base/comm/psb_cgather.f90
+++ b/base/comm/psb_cgather.f90
@ -46,7 +46,7 @@
 !                                              the processes will have a copy.
 !
 subroutine  psb_cgatherm(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_cgatherm
+  use psb_base_mod, psb_protect_name => psb_cgatherm
  implicit none

  complex(psb_spk_), intent(in)    :: locx(:,:)
@ -219,7 +219,7 @@ end subroutine psb_cgatherm
 !                                              default: -1
 !
 subroutine  psb_cgatherv(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_cgatherv
+  use psb_base_mod, psb_protect_name => psb_cgatherv
  implicit none

  complex(psb_spk_), intent(in)    :: locx(:)
--- a/base/comm/psb_chalo.f90
+++ b/base/comm/psb_chalo.f90
@ -54,7 +54,7 @@
 !
 !
 subroutine  psb_chalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
-  use psb_sparse_mod, psb_protect_name => psb_chalom
+  use psb_base_mod, psb_protect_name => psb_chalom
  use psi_mod
  implicit none

@ -276,7 +276,7 @@ end subroutine psb_chalom
 !
 !
 subroutine  psb_chalov(x,desc_a,info,alpha,work,tran,mode,data)
-  use psb_sparse_mod, psb_protect_name => psb_chalov
+  use psb_base_mod, psb_protect_name => psb_chalov
  use psi_mod
  implicit none

--- a/base/comm/psb_covrl.f90
+++ b/base/comm/psb_covrl.f90
@ -65,7 +65,7 @@
 !
 !
 subroutine  psb_covrlm(x,desc_a,info,jx,ik,work,update,mode)
-  use psb_sparse_mod, psb_protect_name => psb_covrlm
+  use psb_base_mod, psb_protect_name => psb_covrlm
  use psi_mod
  implicit none

@ -264,7 +264,7 @@ end subroutine psb_covrlm
 !
 !
 subroutine  psb_covrlv(x,desc_a,info,work,update,mode)
-  use psb_sparse_mod, psb_protect_name => psb_covrlv
+  use psb_base_mod, psb_protect_name => psb_covrlv
  use psi_mod
  implicit none

--- a/base/comm/psb_cscatter.F90
+++ b/base/comm/psb_cscatter.F90
@ -45,7 +45,7 @@
 !                                              Default -1
 subroutine  psb_cscatterm(globx, locx, desc_a, info, iroot)

-  use psb_sparse_mod, psb_protect_name => psb_cscatterm
+  use psb_base_mod, psb_protect_name => psb_cscatterm
 #ifdef MPI_MOD
  use mpi
 #endif
@ -278,7 +278,7 @@ end subroutine psb_cscatterm
 !                                              the processes have a copy.
 !
 subroutine  psb_cscatterv(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_cscatterv
+  use psb_base_mod, psb_protect_name => psb_cscatterv
 #ifdef MPI_MOD
  use mpi
 #endif
--- a/base/comm/psb_dgather.f90
+++ b/base/comm/psb_dgather.f90
@ -47,7 +47,7 @@
 !                                              Default: -1. 
 !
 subroutine  psb_dgatherm(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_dgatherm
+  use psb_base_mod, psb_protect_name => psb_dgatherm
  implicit none

  real(psb_dpk_), intent(in)    :: locx(:,:)
@ -218,7 +218,7 @@ end subroutine psb_dgatherm
 !                                              the processes will have a copy.
 !
 subroutine  psb_dgatherv(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_dgatherv
+  use psb_base_mod, psb_protect_name => psb_dgatherv
  implicit none

  real(psb_dpk_), intent(in)    :: locx(:)
--- a/base/comm/psb_dhalo.f90
+++ b/base/comm/psb_dhalo.f90
@ -54,7 +54,7 @@
 !
 !
 subroutine  psb_dhalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
-  use psb_sparse_mod, psb_protect_name => psb_dhalom
+  use psb_base_mod, psb_protect_name => psb_dhalom
  use psi_mod
  implicit none

@ -275,7 +275,7 @@ end subroutine psb_dhalom
 !
 !
 subroutine  psb_dhalov(x,desc_a,info,alpha,work,tran,mode,data)
-  use psb_sparse_mod, psb_protect_name => psb_dhalov
+  use psb_base_mod, psb_protect_name => psb_dhalov
  use psi_mod
  implicit none

--- a/base/comm/psb_dovrl.f90
+++ b/base/comm/psb_dovrl.f90
@ -64,7 +64,7 @@
 !
 !
 subroutine  psb_dovrlm(x,desc_a,info,jx,ik,work,update,mode)
-  use psb_sparse_mod, psb_protect_name => psb_dovrlm
+  use psb_base_mod, psb_protect_name => psb_dovrlm
  use psi_mod
  implicit none

@ -266,7 +266,7 @@ end subroutine psb_dovrlm
 !
 !
 subroutine  psb_dovrlv(x,desc_a,info,work,update,mode)
-  use psb_sparse_mod, psb_protect_name => psb_dovrlv
+  use psb_base_mod, psb_protect_name => psb_dovrlv
  use psi_mod
  implicit none

--- a/base/comm/psb_dscatter.F90
+++ b/base/comm/psb_dscatter.F90
@ -45,7 +45,7 @@
 !
 subroutine  psb_dscatterm(globx, locx, desc_a, info, iroot)

-  use psb_sparse_mod, psb_protect_name => psb_dscatterm
+  use psb_base_mod, psb_protect_name => psb_dscatterm
 #ifdef MPI_MOD
  use mpi
 #endif
@ -278,7 +278,7 @@ end subroutine psb_dscatterm
 !                                              the processes have a copy.
 !
 subroutine  psb_dscatterv(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_dscatterv
+  use psb_base_mod, psb_protect_name => psb_dscatterv

 #ifdef MPI_MOD
  use mpi
--- a/base/comm/psb_igather.f90
+++ b/base/comm/psb_igather.f90
@ -47,7 +47,7 @@
 !                                              Default: -1. 
 !
 subroutine  psb_igatherm(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_igatherm
+  use psb_base_mod, psb_protect_name => psb_igatherm
  implicit none

  integer, intent(in)    :: locx(:,:)
@ -218,7 +218,7 @@ end subroutine psb_igatherm
 !                                              the processes will have a copy.
 !
 subroutine  psb_igatherv(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_igatherv
+  use psb_base_mod, psb_protect_name => psb_igatherv
  implicit none

  integer, intent(in)    :: locx(:)
--- a/base/comm/psb_ihalo.f90
+++ b/base/comm/psb_ihalo.f90
@ -55,7 +55,7 @@
 !
 !
 subroutine  psb_ihalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
-  use psb_sparse_mod, psb_protect_name => psb_ihalom
+  use psb_base_mod, psb_protect_name => psb_ihalom
  use psi_mod
  implicit none

@ -279,7 +279,7 @@ end subroutine psb_ihalom
 !
 !
 subroutine  psb_ihalov(x,desc_a,info,alpha,work,tran,mode,data)
-  use psb_sparse_mod, psb_protect_name => psb_ihalov
+  use psb_base_mod, psb_protect_name => psb_ihalov
  use psi_mod
  implicit none

--- a/base/comm/psb_iovrl.f90
+++ b/base/comm/psb_iovrl.f90
@ -64,7 +64,7 @@
 !
 !
 subroutine  psb_iovrlm(x,desc_a,info,jx,ik,work,update,mode)
-  use psb_sparse_mod, psb_protect_name => psb_iovrlm
+  use psb_base_mod, psb_protect_name => psb_iovrlm
  use psi_mod
  implicit none

@ -265,7 +265,7 @@ end subroutine psb_iovrlm
 !
 !
 subroutine  psb_iovrlv(x,desc_a,info,work,update,mode)
-  use psb_sparse_mod, psb_protect_name => psb_iovrlv
+  use psb_base_mod, psb_protect_name => psb_iovrlv
  use psi_mod
  implicit none

--- a/base/comm/psb_iscatter.F90
+++ b/base/comm/psb_iscatter.F90
@ -44,7 +44,7 @@
 !                                              the processes have a copy.
 !
 subroutine  psb_iscatterm(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_iscatterm
+  use psb_base_mod, psb_protect_name => psb_iscatterm
 #ifdef MPI_MOD
  use mpi
 #endif
@ -276,7 +276,7 @@ end subroutine psb_iscatterm
 !                                              the processes have a copy.
 !
 subroutine  psb_iscatterv(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_iscatterv
+  use psb_base_mod, psb_protect_name => psb_iscatterv
 #ifdef MPI_MOD
  use mpi
 #endif
--- a/base/comm/psb_sgather.f90
+++ b/base/comm/psb_sgather.f90
@ -47,7 +47,7 @@
 !                                              Default: -1. 
 !
 subroutine  psb_sgatherm(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_sgatherm
+  use psb_base_mod, psb_protect_name => psb_sgatherm
  implicit none

  real(psb_spk_), intent(in)    :: locx(:,:)
@ -218,7 +218,7 @@ end subroutine psb_sgatherm
 !                                              the processes will have a copy.
 !
 subroutine  psb_sgatherv(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_sgatherv
+  use psb_base_mod, psb_protect_name => psb_sgatherv
  implicit none

  real(psb_spk_), intent(in)    :: locx(:)
--- a/base/comm/psb_shalo.f90
+++ b/base/comm/psb_shalo.f90
@ -54,7 +54,7 @@
 !
 !
 subroutine  psb_shalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
-  use psb_sparse_mod, psb_protect_name => psb_shalom
+  use psb_base_mod, psb_protect_name => psb_shalom
  use psi_mod
  implicit none

@ -275,7 +275,7 @@ end subroutine psb_shalom
 !
 !
 subroutine  psb_shalov(x,desc_a,info,alpha,work,tran,mode,data)
-  use psb_sparse_mod, psb_protect_name => psb_shalov
+  use psb_base_mod, psb_protect_name => psb_shalov
  use psi_mod
  implicit none

--- a/base/comm/psb_sovrl.f90
+++ b/base/comm/psb_sovrl.f90
@ -64,7 +64,7 @@
 !
 !
 subroutine  psb_sovrlm(x,desc_a,info,jx,ik,work,update,mode)
-  use psb_sparse_mod, psb_protect_name => psb_sovrlm
+  use psb_base_mod, psb_protect_name => psb_sovrlm
  use psi_mod
  implicit none

@ -266,7 +266,7 @@ end subroutine psb_sovrlm
 !
 !
 subroutine  psb_sovrlv(x,desc_a,info,work,update,mode)
-  use psb_sparse_mod, psb_protect_name => psb_sovrlv
+  use psb_base_mod, psb_protect_name => psb_sovrlv
  use psi_mod
  implicit none

--- a/base/comm/psb_sscatter.F90
+++ b/base/comm/psb_sscatter.F90
@ -44,7 +44,7 @@
 !                                              the processes have a copy. Default -1.
 !
 subroutine  psb_sscatterm(globx, locx, desc_a, info, iroot) 
-  use psb_sparse_mod, psb_protect_name => psb_sscatterm
+  use psb_base_mod, psb_protect_name => psb_sscatterm

 #ifdef MPI_MOD
  use mpi
@ -278,7 +278,7 @@ end subroutine psb_sscatterm
 !                                              the processes have a copy.
 !
 subroutine  psb_sscatterv(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_sscatterv
+  use psb_base_mod, psb_protect_name => psb_sscatterv
 #ifdef MPI_MOD
  use mpi
 #endif
--- a/base/comm/psb_zgather.f90
+++ b/base/comm/psb_zgather.f90
@ -46,7 +46,7 @@
 !                                              the processes will have a copy.
 !
 subroutine  psb_zgatherm(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_zgatherm
+  use psb_base_mod, psb_protect_name => psb_zgatherm
  implicit none

  complex(psb_dpk_), intent(in)    :: locx(:,:)
@ -219,7 +219,7 @@ end subroutine psb_zgatherm
 !                                              default: -1
 !
 subroutine  psb_zgatherv(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_zgatherv
+  use psb_base_mod, psb_protect_name => psb_zgatherv
  implicit none

  complex(psb_dpk_), intent(in)    :: locx(:)
--- a/base/comm/psb_zhalo.f90
+++ b/base/comm/psb_zhalo.f90
@ -54,7 +54,7 @@
 !
 !
 subroutine  psb_zhalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
-  use psb_sparse_mod, psb_protect_name => psb_zhalom
+  use psb_base_mod, psb_protect_name => psb_zhalom
  use psi_mod
  implicit none

@ -276,7 +276,7 @@ end subroutine psb_zhalom
 !
 !
 subroutine  psb_zhalov(x,desc_a,info,alpha,work,tran,mode,data)
-  use psb_sparse_mod, psb_protect_name => psb_zhalov
+  use psb_base_mod, psb_protect_name => psb_zhalov
  use psi_mod
  implicit none

--- a/base/comm/psb_zovrl.f90
+++ b/base/comm/psb_zovrl.f90
@ -65,7 +65,7 @@
 !
 !
 subroutine  psb_zovrlm(x,desc_a,info,jx,ik,work,update,mode)
-  use psb_sparse_mod, psb_protect_name => psb_zovrlm
+  use psb_base_mod, psb_protect_name => psb_zovrlm
  use psi_mod
  implicit none

@ -264,7 +264,7 @@ end subroutine psb_zovrlm
 !
 !
 subroutine  psb_zovrlv(x,desc_a,info,work,update,mode)
-  use psb_sparse_mod, psb_protect_name => psb_zovrlv
+  use psb_base_mod, psb_protect_name => psb_zovrlv
  use psi_mod
  implicit none

--- a/base/comm/psb_zscatter.F90
+++ b/base/comm/psb_zscatter.F90
@ -44,7 +44,7 @@
 !                                              If -1 all the processes have a copy. 
 !                                              Default -1
 subroutine  psb_zscatterm(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_zscatterm
+  use psb_base_mod, psb_protect_name => psb_zscatterm
 #ifdef MPI_MOD
  use mpi
 #endif
@ -277,7 +277,7 @@ end subroutine psb_zscatterm
 !                                              the processes have a copy.
 !
 subroutine  psb_zscatterv(globx, locx, desc_a, info, iroot)
-  use psb_sparse_mod, psb_protect_name => psb_zscatterv
+  use psb_base_mod, psb_protect_name => psb_zscatterv
 #ifdef MPI_MOD
  use mpi
 #endif
--- a/base/modules/Makefile
+++ b/base/modules/Makefile
@ -25,7 +25,7 @@ UTIL_MODS = psb_string_mod.o psb_desc_const_mod.o psb_indx_map_mod.o\

 MODULES=$(BASIC_MODS) $(UTIL_MODS)

-LIBMOD=psb_sparse_mod$(.mod)
+LIBMOD=psb_base_mod$(.mod)
 OBJS = error.o psb_sparse_mod.o $(EXTRA_COBJS) cutil.o
 LIBDIR=..
 CINCLUDES=-I. 
--- a/base/modules/psb_serial_mod.f90
+++ b/base/modules/psb_serial_mod.f90
@ -98,7 +98,7 @@ module psb_serial_mod
      type(psb_z_csr_sparse_mat), intent(out)  :: c
      integer, intent(out)                     :: info
    end subroutine psb_zbase_symbmm
-  end interface
+  end interface psb_symbmm

  interface psb_numbmm
    subroutine psb_snumbmm(a,b,c)
@ -149,7 +149,7 @@ module psb_serial_mod
      class(psb_z_base_sparse_mat), intent(in) :: a,b
      type(psb_z_csr_sparse_mat), intent(inout)  :: c
    end subroutine psb_zbase_numbmm
-  end interface
+  end interface psb_numbmm

  interface psb_rwextd
    subroutine psb_srwextd(nr,a,info,b,rowscale)
@ -224,7 +224,7 @@ module psb_serial_mod
      class(psb_z_base_sparse_mat), intent(in), optional :: b
      logical,intent(in), optional                       :: rowscale
    end subroutine psb_zbase_rwextd
-  end interface
+  end interface psb_rwextd


  interface psb_geprt
@ -324,9 +324,12 @@ module psb_serial_mod
      complex(psb_dpk_), intent(in)  :: a(:)
      character(len=*), optional     :: head
    end subroutine psb_zgeprt1
-  end interface
-  
+  end interface psb_geprt

+  interface psb_csprt
+    module procedure psb_scsprt, psb_scsprtn, psb_dcsprt, psb_dcsprtn, &
+         & psb_ccsprt, psb_ccsprtn, psb_zcsprt, psb_zcsprtn
+  end interface psb_csprt

  interface psb_spdot_srtd
    function psb_s_spdot_srtd(nv1,iv1,v1,nv2,iv2,v2) result(dot) 
@ -360,7 +363,7 @@ module psb_serial_mod
      complex(psb_dpk_), intent(in) :: v1(*),v2(*)
      complex(psb_dpk_)      :: dot
    end function psb_z_spdot_srtd
-  end interface
+  end interface psb_spdot_srtd


  interface psb_spge_dot
@ -395,7 +398,7 @@ module psb_serial_mod
      complex(psb_dpk_), intent(in) :: v1(*),v2(*)
      complex(psb_dpk_)      :: dot
    end function psb_z_spge_dot
-  end interface
+  end interface psb_spge_dot


  interface psb_nspaxpby
@ -410,7 +413,7 @@ module psb_serial_mod
      real(psb_dpk_), intent (in)       :: alpha, beta
      integer, intent(out)              :: info
    end subroutine psb_d_nspaxpby
-  end interface
+  end interface psb_nspaxpby

  interface psb_aspxpby
    subroutine psb_s_aspxpby(alpha, nx, ix, x, beta, y, info)
@ -449,8 +452,113 @@ module psb_serial_mod
      complex(psb_dpk_), intent (in)    :: alpha, beta
      integer, intent(out)              :: info
    end subroutine psb_z_aspxpby
-  end interface
+  end interface psb_aspxpby
+
+contains
+
+  subroutine psb_scsprt(iout,a,iv,irs,ics,head,ivr,ivc)
+    use psb_mat_mod
+    integer, intent(in)       :: iout
+    type(psb_sspmat_type), intent(in) :: a
+    integer, intent(in), optional :: iv(:)
+    integer, intent(in), optional :: irs,ics
+    character(len=*), optional    :: head
+    integer, intent(in), optional :: ivr(:),ivc(:)
+
+    call a%print(iout,iv,irs,ics,head,ivr,ivc)
+
+  end subroutine psb_scsprt
+
+  subroutine psb_scsprtn(fname,a,iv,irs,ics,head,ivr,ivc)
+    use psb_mat_mod
+    character(len=*), intent(in)  :: fname   
+    type(psb_sspmat_type), intent(in) :: a
+    integer, intent(in), optional :: iv(:)
+    integer, intent(in), optional :: irs,ics
+    character(len=*), optional    :: head
+    integer, intent(in), optional :: ivr(:),ivc(:)
+
+    call a%print(fname,iv,irs,ics,head,ivr,ivc)
+
+  end subroutine psb_scsprtn
+
+  subroutine psb_dcsprt(iout,a,iv,irs,ics,head,ivr,ivc)
+    use psb_mat_mod
+    integer, intent(in)       :: iout
+    type(psb_dspmat_type), intent(in) :: a
+    integer, intent(in), optional :: iv(:)
+    integer, intent(in), optional :: irs,ics
+    character(len=*), optional    :: head
+    integer, intent(in), optional :: ivr(:),ivc(:)
+
+    call a%print(iout,iv,irs,ics,head,ivr,ivc)
+
+  end subroutine psb_dcsprt
+
+  subroutine psb_dcsprtn(fname,a,iv,irs,ics,head,ivr,ivc)
+    use psb_mat_mod
+    character(len=*), intent(in)  :: fname   
+    type(psb_dspmat_type), intent(in) :: a
+    integer, intent(in), optional :: iv(:)
+    integer, intent(in), optional :: irs,ics
+    character(len=*), optional    :: head
+    integer, intent(in), optional :: ivr(:),ivc(:)
+
+    call a%print(fname,iv,irs,ics,head,ivr,ivc)
+
+  end subroutine psb_dcsprtn
+
+  subroutine psb_ccsprt(iout,a,iv,irs,ics,head,ivr,ivc)
+    use psb_mat_mod
+    integer, intent(in)       :: iout
+    type(psb_cspmat_type), intent(in) :: a
+    integer, intent(in), optional :: iv(:)
+    integer, intent(in), optional :: irs,ics
+    character(len=*), optional    :: head
+    integer, intent(in), optional :: ivr(:),ivc(:)
+
+    call a%print(iout,iv,irs,ics,head,ivr,ivc)
+
+  end subroutine psb_ccsprt
+
+  subroutine psb_ccsprtn(fname,a,iv,irs,ics,head,ivr,ivc)
+    use psb_mat_mod
+    character(len=*), intent(in)  :: fname   
+    type(psb_cspmat_type), intent(in) :: a
+    integer, intent(in), optional :: iv(:)
+    integer, intent(in), optional :: irs,ics
+    character(len=*), optional    :: head
+    integer, intent(in), optional :: ivr(:),ivc(:)
+
+    call a%print(fname,iv,irs,ics,head,ivr,ivc)
+
+  end subroutine psb_ccsprtn
+
+  subroutine psb_zcsprt(iout,a,iv,irs,ics,head,ivr,ivc)
+    use psb_mat_mod
+    integer, intent(in)       :: iout
+    type(psb_zspmat_type), intent(in) :: a
+    integer, intent(in), optional :: iv(:)
+    integer, intent(in), optional :: irs,ics
+    character(len=*), optional    :: head
+    integer, intent(in), optional :: ivr(:),ivc(:)
+
+    call a%print(iout,iv,irs,ics,head,ivr,ivc)
+
+  end subroutine psb_zcsprt
+
+  subroutine psb_zcsprtn(fname,a,iv,irs,ics,head,ivr,ivc)
+    use psb_mat_mod
+    character(len=*), intent(in)  :: fname   
+    type(psb_zspmat_type), intent(in) :: a
+    integer, intent(in), optional :: iv(:)
+    integer, intent(in), optional :: irs,ics
+    character(len=*), optional    :: head
+    integer, intent(in), optional :: ivr(:),ivc(:)
+
+    call a%print(fname,iv,irs,ics,head,ivr,ivc)

+  end subroutine psb_zcsprtn

 end module psb_serial_mod

--- a/base/modules/psb_sparse_mod.f90
+++ b/base/modules/psb_sparse_mod.f90
@ -29,7 +29,7 @@
 !!$  POSSIBILITY OF SUCH DAMAGE.
 !!$ 
 !!$  
-module psb_sparse_mod
+module psb_base_mod
  use psb_string_mod
  use psb_error_mod
  use psb_penv_mod
@ -42,4 +42,4 @@ module psb_sparse_mod
  use psb_psblas_mod
  use psb_gps_mod
  use psb_tools_mod
-end module psb_sparse_mod
+end module psb_base_mod
--- a/base/psblas/psb_camax.f90
+++ b/base/psblas/psb_camax.f90
@ -45,7 +45,7 @@
 !    jx     -  integer(optional).    The column offset.
 !
 function psb_camax(x,desc_a, info, jx)
-  use psb_sparse_mod, psb_protect_name => psb_camax
+  use psb_base_mod, psb_protect_name => psb_camax

  implicit none

@ -175,7 +175,7 @@ end function psb_camax
 !    info   -  integer.              Return code
 !
 function psb_camaxv (x,desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_camaxv
+  use psb_base_mod, psb_protect_name => psb_camaxv

  implicit none

@ -303,7 +303,7 @@ end function psb_camaxv
 !    jx     -  integer(optional).   The column offset.
 !
 subroutine psb_camaxvs(res,x,desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_camaxvs
+  use psb_base_mod, psb_protect_name => psb_camaxvs

  implicit none

@ -428,7 +428,7 @@ end subroutine psb_camaxvs
 !    info   -  integer.             Return code
 !
 subroutine psb_cmamaxs(res,x,desc_a, info,jx)
-  use psb_sparse_mod, psb_protect_name => psb_cmamaxs
+  use psb_base_mod, psb_protect_name => psb_cmamaxs

  implicit none

--- a/base/psblas/psb_casum.f90
+++ b/base/psblas/psb_casum.f90
@ -45,7 +45,7 @@
 !    jx     -  integer(optional).     The column offset.
 !
 function psb_casum (x,desc_a, info, jx)
-  use psb_sparse_mod, psb_protect_name => psb_casum
+  use psb_base_mod, psb_protect_name => psb_casum

  implicit none

@ -187,7 +187,7 @@ end function psb_casum
 !    info   -  integer.              Return code
 !
 function psb_casumv(x,desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_casumv
+  use psb_base_mod, psb_protect_name => psb_casumv

  implicit none

@ -325,7 +325,7 @@ end function psb_casumv
 !    jx     -  integer(optional).    The column offset.
 !
 subroutine psb_casumvs(res,x,desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_casumvs
+  use psb_base_mod, psb_protect_name => psb_casumvs

  implicit none

--- a/base/psblas/psb_caxpby.f90
+++ b/base/psblas/psb_caxpby.f90
@ -51,7 +51,7 @@
 !    jy     -  integer(optional)    The column offset for Y 
 !
 subroutine  psb_caxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
-  use psb_sparse_mod, psb_protect_name => psb_caxpby
+  use psb_base_mod, psb_protect_name => psb_caxpby

  implicit none                    

@ -202,7 +202,7 @@ end subroutine psb_caxpby
 !
 !
 subroutine  psb_caxpbyv(alpha, x, beta,y,desc_a,info)
-  use psb_sparse_mod, psb_protect_name => psb_caxpbyv
+  use psb_base_mod, psb_protect_name => psb_caxpbyv
  implicit none                    

  integer, intent(out)            :: info
--- a/base/psblas/psb_cdot.f90
+++ b/base/psblas/psb_cdot.f90
@ -49,7 +49,7 @@
 !    jy     -  integer(optional).    The column offset for sub( Y ).
 !
 function psb_cdot(x, y,desc_a, info, jx, jy)  
-  use psb_sparse_mod, psb_protect_name => psb_cdot
+  use psb_base_mod, psb_protect_name => psb_cdot
  implicit none

  complex(psb_spk_), intent(in)     :: x(:,:), y(:,:)
@ -199,7 +199,7 @@ end function psb_cdot
 !    info   -  integer.              Return code
 !
 function psb_cdotv(x, y,desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_cdotv
+  use psb_base_mod, psb_protect_name => psb_cdotv
  implicit none

  complex(psb_spk_), intent(in)  :: x(:), y(:)
@ -333,7 +333,7 @@ end function psb_cdotv
 !    info   -  integer.             Return code
 !
 subroutine psb_cdotvs(res, x, y,desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_cdotvs
+  use psb_base_mod, psb_protect_name => psb_cdotvs
  implicit none

  complex(psb_spk_), intent(in)     :: x(:), y(:)
@ -465,7 +465,7 @@ end subroutine psb_cdotvs
 !    info   -  integer.             Return code
 !
 subroutine psb_cmdots(res, x, y, desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_cmdots
+  use psb_base_mod, psb_protect_name => psb_cmdots
  implicit none

  complex(psb_spk_), intent(in)     :: x(:,:), y(:,:)
--- a/base/psblas/psb_cnrmi.f90
+++ b/base/psblas/psb_cnrmi.f90
@ -42,7 +42,7 @@
 !    info   -  integer.                   Return code
 !
 function psb_cnrmi(a,desc_a,info)  
-  use psb_sparse_mod, psb_protect_name => psb_cnrmi
+  use psb_base_mod, psb_protect_name => psb_cnrmi
  implicit none

  type(psb_cspmat_type), intent(in)   :: a
--- a/base/psblas/psb_cspmm.f90
+++ b/base/psblas/psb_cspmm.f90
@ -64,7 +64,7 @@
 ! 
 subroutine  psb_cspmm(alpha,a,x,beta,y,desc_a,info,&
     & trans, k, jx, jy, work, doswap)   
-  use psb_sparse_mod, psb_protect_name => psb_cspmm
+  use psb_base_mod, psb_protect_name => psb_cspmm
  use psi_mod
  implicit none

@ -416,7 +416,7 @@ end subroutine psb_cspmm
 ! 
 subroutine  psb_cspmv(alpha,a,x,beta,y,desc_a,info,&
     & trans, work, doswap)   
-  use psb_sparse_mod, psb_protect_name => psb_cspmv
+  use psb_base_mod, psb_protect_name => psb_cspmv
  use psi_mod
  implicit none

--- a/base/psblas/psb_cspsm.f90
+++ b/base/psblas/psb_cspsm.f90
@ -75,7 +75,7 @@
 ! 
 subroutine  psb_cspsm(alpha,a,x,beta,y,desc_a,info,&
     & trans, scale, choice, diag, k, jx, jy, work)   
-  use psb_sparse_mod, psb_protect_name => psb_cspsm
+  use psb_base_mod, psb_protect_name => psb_cspsm
  use psi_mod
  implicit none

@ -354,7 +354,7 @@ end subroutine psb_cspsm
 ! 
 subroutine  psb_cspsv(alpha,a,x,beta,y,desc_a,info,&
     & trans, scale, choice, diag, work)   
-  use psb_sparse_mod, psb_protect_name => psb_cspsv
+  use psb_base_mod, psb_protect_name => psb_cspsv
  use psi_mod
  implicit none 

--- a/base/psblas/psb_damax.f90
+++ b/base/psblas/psb_damax.f90
@ -301,7 +301,7 @@ end function psb_damaxv
 !    jx     -  integer(optional).         The column offset.
 !
 subroutine psb_damaxvs (res,x,desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_damaxvs 
+  use psb_base_mod, psb_protect_name => psb_damaxvs 
  implicit none

  real(psb_dpk_), intent(in)      :: x(:)
@ -421,7 +421,7 @@ end subroutine psb_damaxvs
 !    info   -  integer.                 Return code
 !
 subroutine psb_dmamaxs (res,x,desc_a, info,jx)
-  use psb_sparse_mod, psb_protect_name => psb_dmamaxs 
+  use psb_base_mod, psb_protect_name => psb_dmamaxs 
  implicit none

  real(psb_dpk_), intent(in)      :: x(:,:)
--- a/base/psblas/psb_dasum.f90
+++ b/base/psblas/psb_dasum.f90
@ -327,7 +327,7 @@ end function psb_dasumv
 !    jx     -  integer(optional).         The column offset.
 !
 subroutine psb_dasumvs(res,x,desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_dasumvs
+  use psb_base_mod, psb_protect_name => psb_dasumvs
  implicit none

  real(psb_dpk_), intent(in)      :: x(:)
--- a/base/psblas/psb_daxpby.f90
+++ b/base/psblas/psb_daxpby.f90
@ -51,7 +51,7 @@
 !    jy     -  integer(optional).   The column offset for sub( Y ).
 !
 subroutine  psb_daxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
-  use psb_sparse_mod, psb_protect_name => psb_daxpby
+  use psb_base_mod, psb_protect_name => psb_daxpby
  implicit none                    

  integer, intent(in), optional   :: n, jx, jy
@ -202,7 +202,7 @@ end subroutine psb_daxpby
 !
 !
 subroutine  psb_daxpbyv(alpha, x, beta,y,desc_a,info)
-  use psb_sparse_mod, psb_protect_name => psb_daxpbyv
+  use psb_base_mod, psb_protect_name => psb_daxpbyv
  implicit none                    

  integer, intent(out)            :: info
--- a/base/psblas/psb_ddot.f90
+++ b/base/psblas/psb_ddot.f90
@ -340,7 +340,7 @@ end function psb_ddotv
 !    info   -  integer.                 Return code
 !
 subroutine psb_ddotvs(res, x, y,desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_ddotvs
+  use psb_base_mod, psb_protect_name => psb_ddotvs
  implicit none

  real(psb_dpk_), intent(in)     :: x(:), y(:)
@ -472,7 +472,7 @@ end subroutine psb_ddotvs
 !    info   -  integer.             Return code
 !
 subroutine psb_dmdots(res, x, y, desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_dmdots
+  use psb_base_mod, psb_protect_name => psb_dmdots
  implicit none

  real(psb_dpk_), intent(in)     :: x(:,:), y(:,:)
--- a/base/psblas/psb_dspmm.f90
+++ b/base/psblas/psb_dspmm.f90
@ -64,7 +64,7 @@
 ! 
 subroutine  psb_dspmm(alpha,a,x,beta,y,desc_a,info,&
     & trans, k, jx, jy, work, doswap)   
-  use psb_sparse_mod, psb_protect_name => psb_dspmm
+  use psb_base_mod, psb_protect_name => psb_dspmm
  use psi_mod
  implicit none

@ -416,7 +416,7 @@ end subroutine psb_dspmm
 ! 
 subroutine  psb_dspmv(alpha,a,x,beta,y,desc_a,info,&
     & trans, work, doswap)   
-  use psb_sparse_mod, psb_protect_name => psb_dspmv
+  use psb_base_mod, psb_protect_name => psb_dspmv
  use psi_mod
  implicit none

--- a/base/psblas/psb_dspnrm1.f90
+++ b/base/psblas/psb_dspnrm1.f90
@ -49,7 +49,7 @@ function psb_dspnrm1(a,desc_a,info)
 !!$  use psb_penv_mod
 !!$  use psb_mat_mod
 !!$  use psb_tools_mod
-  use psb_sparse_mod, psb_protect_name => psb_dspnrm1
+  use psb_base_mod, psb_protect_name => psb_dspnrm1
  implicit none

  type(psb_dspmat_type), intent(in) :: a
--- a/base/psblas/psb_dspsm.f90
+++ b/base/psblas/psb_dspsm.f90
@ -76,7 +76,7 @@
 ! 
 subroutine  psb_dspsm(alpha,a,x,beta,y,desc_a,info,&
     & trans, scale, choice, diag, k, jx, jy, work)   
-  use psb_sparse_mod, psb_protect_name => psb_dspsm
+  use psb_base_mod, psb_protect_name => psb_dspsm
  use psi_mod
  implicit none

@ -355,7 +355,7 @@ end subroutine psb_dspsm
 ! 
 subroutine  psb_dspsv(alpha,a,x,beta,y,desc_a,info,&
     & trans, scale, choice, diag, work)   
-  use psb_sparse_mod, psb_protect_name => psb_dspsv
+  use psb_base_mod, psb_protect_name => psb_dspsv
  use psi_mod
  implicit none 

--- a/base/psblas/psb_samax.f90
+++ b/base/psblas/psb_samax.f90
@ -301,7 +301,7 @@ end function psb_samaxv
 !    jx     -  integer(optional).         The column offset.
 !
 subroutine psb_samaxvs (res,x,desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_samaxvs 
+  use psb_base_mod, psb_protect_name => psb_samaxvs 
  implicit none

  real(psb_spk_), intent(in)      :: x(:)
@ -421,7 +421,7 @@ end subroutine psb_samaxvs
 !    info   -  integer.                 Return code
 !
 subroutine psb_smamaxs (res,x,desc_a, info,jx)
-  use psb_sparse_mod, psb_protect_name => psb_smamaxs 
+  use psb_base_mod, psb_protect_name => psb_smamaxs 
  implicit none

  real(psb_spk_), intent(in)      :: x(:,:)
--- a/base/psblas/psb_sasum.f90
+++ b/base/psblas/psb_sasum.f90
@ -327,7 +327,7 @@ end function psb_sasumv
 !    jx     -  integer(optional).         The column offset.
 !
 subroutine psb_sasumvs(res,x,desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_sasumvs
+  use psb_base_mod, psb_protect_name => psb_sasumvs
  implicit none

  real(psb_spk_), intent(in)      :: x(:)
--- a/base/psblas/psb_saxpby.f90
+++ b/base/psblas/psb_saxpby.f90
@ -51,7 +51,7 @@
 !    jy     -  integer(optional).   The column offset for sub( Y ).
 !
 subroutine  psb_saxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
-  use psb_sparse_mod, psb_protect_name => psb_saxpby
+  use psb_base_mod, psb_protect_name => psb_saxpby
  implicit none                    

  integer, intent(in), optional   :: n, jx, jy
@ -202,7 +202,7 @@ end subroutine psb_saxpby
 !
 !
 subroutine  psb_saxpbyv(alpha, x, beta,y,desc_a,info)
-  use psb_sparse_mod, psb_protect_name => psb_saxpbyv
+  use psb_base_mod, psb_protect_name => psb_saxpbyv
  implicit none                    

  integer, intent(out)            :: info
--- a/base/psblas/psb_sdot.f90
+++ b/base/psblas/psb_sdot.f90
@ -339,7 +339,7 @@ end function psb_sdotv
 !    info   -  integer.                 Return code
 !
 subroutine psb_sdotvs(res, x, y,desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_sdotvs
+  use psb_base_mod, psb_protect_name => psb_sdotvs
  implicit none

  real(psb_spk_), intent(in)     :: x(:), y(:)
@ -471,7 +471,7 @@ end subroutine psb_sdotvs
 !    info   -  integer.             Return code
 !
 subroutine psb_smdots(res, x, y, desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_smdots
+  use psb_base_mod, psb_protect_name => psb_smdots
  implicit none

  real(psb_spk_), intent(in)     :: x(:,:), y(:,:)
--- a/base/psblas/psb_sspmm.f90
+++ b/base/psblas/psb_sspmm.f90
@ -64,7 +64,7 @@
 ! 
 subroutine  psb_sspmm(alpha,a,x,beta,y,desc_a,info,&
     & trans, k, jx, jy, work, doswap)   
-  use psb_sparse_mod, psb_protect_name => psb_sspmm
+  use psb_base_mod, psb_protect_name => psb_sspmm
  use psi_mod
  implicit none

@ -416,7 +416,7 @@ end subroutine psb_sspmm
 ! 
 subroutine  psb_sspmv(alpha,a,x,beta,y,desc_a,info,&
     & trans, work, doswap)   
-  use psb_sparse_mod, psb_protect_name => psb_sspmv
+  use psb_base_mod, psb_protect_name => psb_sspmv
  use psi_mod
  implicit none

--- a/base/psblas/psb_sspsm.f90
+++ b/base/psblas/psb_sspsm.f90
@ -76,7 +76,7 @@
 ! 
 subroutine  psb_sspsm(alpha,a,x,beta,y,desc_a,info,&
     & trans, scale, choice, diag, k, jx, jy, work)   
-  use psb_sparse_mod, psb_protect_name => psb_sspsm
+  use psb_base_mod, psb_protect_name => psb_sspsm
  use psi_mod
  implicit none

@ -355,7 +355,7 @@ end subroutine psb_sspsm
 ! 
 subroutine  psb_sspsv(alpha,a,x,beta,y,desc_a,info,&
     & trans, scale, choice, diag, work)   
-  use psb_sparse_mod, psb_protect_name => psb_sspsv
+  use psb_base_mod, psb_protect_name => psb_sspsv
  use psi_mod
  implicit none 

--- a/base/psblas/psb_sxdot.f90
+++ b/base/psblas/psb_sxdot.f90
@ -360,7 +360,7 @@ end function sxdot
 !    info   -  integer.                 Return code
 !
 subroutine psb_sxdotvs(res, x, y,desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_sxdotvs
+  use psb_base_mod, psb_protect_name => psb_sxdotvs
  implicit none

  real(psb_spk_), intent(in)     :: x(:), y(:)
@ -492,7 +492,7 @@ end subroutine psb_sxdotvs
 !    info   -  integer.             Return code
 !
 subroutine psb_sxmdots(res, x, y, desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_sxmdots
+  use psb_base_mod, psb_protect_name => psb_sxmdots
  implicit none

  real(psb_spk_), intent(in)     :: x(:,:), y(:,:)
--- a/base/psblas/psb_zamax.f90
+++ b/base/psblas/psb_zamax.f90
@ -309,7 +309,7 @@ end function psb_zamaxv
 !    jx     -  integer(optional).   The column offset.
 !
 subroutine psb_zamaxvs(res,x,desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_zamaxvs
+  use psb_base_mod, psb_protect_name => psb_zamaxvs
  implicit none

  complex(psb_dpk_), intent(in)   :: x(:)
@ -433,7 +433,7 @@ end subroutine psb_zamaxvs
 !    info   -  integer.             Return code
 !
 subroutine psb_zmamaxs(res,x,desc_a, info,jx)
-  use psb_sparse_mod, psb_protect_name => psb_zmamaxs
+  use psb_base_mod, psb_protect_name => psb_zmamaxs
  implicit none

  complex(psb_dpk_), intent(in)   :: x(:,:)
--- a/base/psblas/psb_zasum.f90
+++ b/base/psblas/psb_zasum.f90
@ -337,7 +337,7 @@ end function psb_zasumv
 !    jx     -  integer(optional).    The column offset.
 !
 subroutine psb_zasumvs(res,x,desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_zasumvs
+  use psb_base_mod, psb_protect_name => psb_zasumvs
  implicit none

  complex(psb_dpk_), intent(in)   :: x(:)
--- a/base/psblas/psb_zaxpby.f90
+++ b/base/psblas/psb_zaxpby.f90
@ -51,7 +51,7 @@
 !    jy     -  integer(optional)    The column offset for Y 
 !
 subroutine  psb_zaxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
-  use psb_sparse_mod, psb_protect_name => psb_zaxpby
+  use psb_base_mod, psb_protect_name => psb_zaxpby
  implicit none                    

  integer, intent(in), optional   :: n, jx, jy
@ -201,7 +201,7 @@ end subroutine psb_zaxpby
 !
 !
 subroutine  psb_zaxpbyv(alpha, x, beta,y,desc_a,info)
-  use psb_sparse_mod, psb_protect_name => psb_zaxpbyv
+  use psb_base_mod, psb_protect_name => psb_zaxpbyv
  implicit none                    

  integer, intent(out)            :: info
--- a/base/psblas/psb_zdot.f90
+++ b/base/psblas/psb_zdot.f90
@ -339,7 +339,7 @@ end function psb_zdotv
 !    info   -  integer.             Return code
 !
 subroutine psb_zdotvs(res, x, y,desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_zdotvs
+  use psb_base_mod, psb_protect_name => psb_zdotvs
  implicit none

  complex(psb_dpk_), intent(in)     :: x(:), y(:)
@ -471,7 +471,7 @@ end subroutine psb_zdotvs
 !    info   -  integer.             Return code
 !
 subroutine psb_zmdots(res, x, y, desc_a, info)  
-  use psb_sparse_mod, psb_protect_name => psb_zmdots
+  use psb_base_mod, psb_protect_name => psb_zmdots
  implicit none

  complex(psb_dpk_), intent(in)     :: x(:,:), y(:,:)
--- a/base/psblas/psb_zspmm.f90
+++ b/base/psblas/psb_zspmm.f90
@ -64,7 +64,7 @@
 ! 
 subroutine  psb_zspmm(alpha,a,x,beta,y,desc_a,info,&
     & trans, k, jx, jy, work, doswap)   
-  use psb_sparse_mod, psb_protect_name => psb_zspmm
+  use psb_base_mod, psb_protect_name => psb_zspmm
  use psi_mod
  implicit none

@ -416,7 +416,7 @@ end subroutine psb_zspmm
 ! 
 subroutine  psb_zspmv(alpha,a,x,beta,y,desc_a,info,&
     & trans, work, doswap)   
-  use psb_sparse_mod, psb_protect_name => psb_zspmv
+  use psb_base_mod, psb_protect_name => psb_zspmv
  use psi_mod
  implicit none

--- a/base/psblas/psb_zspsm.f90
+++ b/base/psblas/psb_zspsm.f90
@ -75,7 +75,7 @@
 ! 
 subroutine  psb_zspsm(alpha,a,x,beta,y,desc_a,info,&
     & trans, scale, choice, diag, k, jx, jy, work)   
-  use psb_sparse_mod, psb_protect_name => psb_zspsm
+  use psb_base_mod, psb_protect_name => psb_zspsm
  use psi_mod
  implicit none

@ -354,7 +354,7 @@ end subroutine psb_zspsm
 ! 
 subroutine  psb_zspsv(alpha,a,x,beta,y,desc_a,info,&
     & trans, scale, choice, diag, work)   
-  use psb_sparse_mod, psb_protect_name => psb_zspsv
+  use psb_base_mod, psb_protect_name => psb_zspsv
  use psi_mod
  implicit none 

--- a/base/serial/psb_cnumbmm.f90
+++ b/base/serial/psb_cnumbmm.f90
@ -41,7 +41,7 @@
 !

 subroutine psb_cnumbmm(a,b,c)
-  use psb_sparse_mod, psb_protect_name => psb_cnumbmm
+  use psb_base_mod, psb_protect_name => psb_cnumbmm
  implicit none 

  type(psb_cspmat_type), intent(in) :: a,b
--- a/base/serial/psb_crwextd.f90
+++ b/base/serial/psb_crwextd.f90
@ -39,7 +39,7 @@
 !
 !
 subroutine psb_crwextd(nr,a,info,b,rowscale)
-  use psb_sparse_mod, psb_protect_name => psb_crwextd
+  use psb_base_mod, psb_protect_name => psb_crwextd
  implicit none

  ! Extend matrix A up to NR rows with empty ones (i.e.: all zeroes)
@ -99,7 +99,7 @@ subroutine psb_crwextd(nr,a,info,b,rowscale)

 end subroutine psb_crwextd
 subroutine psb_cbase_rwextd(nr,a,info,b,rowscale)
-  use psb_sparse_mod, psb_protect_name => psb_cbase_rwextd
+  use psb_base_mod, psb_protect_name => psb_cbase_rwextd
  implicit none

  ! Extend matrix A up to NR rows with empty ones (i.e.: all zeroes)
--- a/base/serial/psb_csymbmm.f90
+++ b/base/serial/psb_csymbmm.f90
@ -40,7 +40,7 @@
 !

 subroutine psb_csymbmm(a,b,c,info)
-  use psb_sparse_mod, psb_protect_name => psb_csymbmm
+  use psb_base_mod, psb_protect_name => psb_csymbmm
  implicit none 

  type(psb_cspmat_type), intent(in)    :: a,b
--- a/base/serial/psb_dnumbmm.f90
+++ b/base/serial/psb_dnumbmm.f90
@ -41,7 +41,7 @@
 !

 subroutine psb_dnumbmm(a,b,c)
-  use psb_sparse_mod, psb_protect_name => psb_dnumbmm
+  use psb_base_mod, psb_protect_name => psb_dnumbmm
  implicit none 

  type(psb_dspmat_type), intent(in) :: a,b
--- a/base/serial/psb_drwextd.f90
+++ b/base/serial/psb_drwextd.f90
@ -39,7 +39,7 @@
 !
 !
 subroutine psb_drwextd(nr,a,info,b,rowscale)
-  use psb_sparse_mod, psb_protect_name => psb_drwextd
+  use psb_base_mod, psb_protect_name => psb_drwextd
  implicit none

  ! Extend matrix A up to NR rows with empty ones (i.e.: all zeroes)
@ -99,7 +99,7 @@ subroutine psb_drwextd(nr,a,info,b,rowscale)

 end subroutine psb_drwextd
 subroutine psb_dbase_rwextd(nr,a,info,b,rowscale)
-  use psb_sparse_mod, psb_protect_name => psb_dbase_rwextd
+  use psb_base_mod, psb_protect_name => psb_dbase_rwextd
  implicit none

  ! Extend matrix A up to NR rows with empty ones (i.e.: all zeroes)
--- a/base/serial/psb_dsymbmm.f90
+++ b/base/serial/psb_dsymbmm.f90
@ -40,7 +40,7 @@
 !

 subroutine psb_dsymbmm(a,b,c,info)
-  use psb_sparse_mod, psb_protect_name => psb_dsymbmm
+  use psb_base_mod, psb_protect_name => psb_dsymbmm
  implicit none 

  type(psb_dspmat_type), intent(in)    :: a,b
--- a/base/serial/psb_snumbmm.f90
+++ b/base/serial/psb_snumbmm.f90
@ -41,7 +41,7 @@
 !

 subroutine psb_snumbmm(a,b,c)
-  use psb_sparse_mod, psb_protect_name => psb_snumbmm
+  use psb_base_mod, psb_protect_name => psb_snumbmm
  implicit none 

  type(psb_sspmat_type), intent(in) :: a,b
--- a/base/serial/psb_srwextd.f90
+++ b/base/serial/psb_srwextd.f90
@ -39,7 +39,7 @@
 !
 !
 subroutine psb_srwextd(nr,a,info,b,rowscale)
-  use psb_sparse_mod, psb_protect_name => psb_srwextd
+  use psb_base_mod, psb_protect_name => psb_srwextd
  implicit none

  ! Extend matrix A up to NR rows with empty ones (i.e.: all zeroes)
@ -99,7 +99,7 @@ subroutine psb_srwextd(nr,a,info,b,rowscale)

 end subroutine psb_srwextd
 subroutine psb_sbase_rwextd(nr,a,info,b,rowscale)
-  use psb_sparse_mod, psb_protect_name => psb_sbase_rwextd
+  use psb_base_mod, psb_protect_name => psb_sbase_rwextd
  implicit none

  ! Extend matrix A up to NR rows with empty ones (i.e.: all zeroes)
--- a/base/serial/psb_ssymbmm.f90
+++ b/base/serial/psb_ssymbmm.f90
@ -40,7 +40,7 @@
 !

 subroutine psb_ssymbmm(a,b,c,info)
-  use psb_sparse_mod, psb_protect_name => psb_ssymbmm
+  use psb_base_mod, psb_protect_name => psb_ssymbmm
  implicit none 

  type(psb_sspmat_type), intent(in)    :: a,b
--- a/base/serial/psb_znumbmm.f90
+++ b/base/serial/psb_znumbmm.f90
@ -41,7 +41,7 @@
 !

 subroutine psb_znumbmm(a,b,c)
-  use psb_sparse_mod, psb_protect_name => psb_znumbmm
+  use psb_base_mod, psb_protect_name => psb_znumbmm
  implicit none 

  type(psb_zspmat_type), intent(in) :: a,b
--- a/base/serial/psb_zrwextd.f90
+++ b/base/serial/psb_zrwextd.f90
@ -39,7 +39,7 @@
 !
 !
 subroutine psb_zrwextd(nr,a,info,b,rowscale)
-  use psb_sparse_mod, psb_protect_name => psb_zrwextd
+  use psb_base_mod, psb_protect_name => psb_zrwextd
  implicit none

  ! Extend matrix A up to NR rows with empty ones (i.e.: all zeroes)
@ -99,7 +99,7 @@ subroutine psb_zrwextd(nr,a,info,b,rowscale)

 end subroutine psb_zrwextd
 subroutine psb_zbase_rwextd(nr,a,info,b,rowscale)
-  use psb_sparse_mod, psb_protect_name => psb_zbase_rwextd
+  use psb_base_mod, psb_protect_name => psb_zbase_rwextd
  implicit none

  ! Extend matrix A up to NR rows with empty ones (i.e.: all zeroes)
--- a/base/serial/psb_zsymbmm.f90
+++ b/base/serial/psb_zsymbmm.f90
@ -40,7 +40,7 @@
 !

 subroutine psb_zsymbmm(a,b,c,info)
-  use psb_sparse_mod, psb_protect_name => psb_zsymbmm
+  use psb_base_mod, psb_protect_name => psb_zsymbmm
  implicit none 

  type(psb_zspmat_type), intent(in)    :: a,b
--- a/base/tools/psb_callc.f90
+++ b/base/tools/psb_callc.f90
@ -43,7 +43,7 @@
 !    n      - optional number of columns.
 !    lb     - optional lower bound on column indices
 subroutine psb_calloc(x, desc_a, info, n, lb)
-  use psb_sparse_mod, psb_protect_name => psb_calloc
+  use psb_base_mod, psb_protect_name => psb_calloc
  use psi_mod
  implicit none

@ -176,7 +176,7 @@ end subroutine psb_calloc
 !    desc_a - the communication descriptor.
 !    info   - return code
 subroutine psb_callocv(x, desc_a,info,n)
-  use psb_sparse_mod, psb_protect_name => psb_callocv
+  use psb_base_mod, psb_protect_name => psb_callocv
  use psi_mod
  implicit none

--- a/base/tools/psb_casb.f90
+++ b/base/tools/psb_casb.f90
@ -43,7 +43,7 @@
 !    desc_a  - type(psb_desc_type).  The communication descriptor.
 !    info    - integer.                return code
 subroutine psb_casb(x, desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_casb
+  use psb_base_mod, psb_protect_name => psb_casb
  implicit none

  type(psb_desc_type), intent(in) ::  desc_a
@ -176,7 +176,7 @@ end subroutine psb_casb
 !    desc_a  - type(psb_desc_type).  The communication descriptor.
 !    info    - integer.                Return  code
 subroutine psb_casbv(x, desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_casbv
+  use psb_base_mod, psb_protect_name => psb_casbv
  implicit none

  type(psb_desc_type), intent(in)                 ::  desc_a
--- a/base/tools/psb_ccdbldext.F90
+++ b/base/tools/psb_ccdbldext.F90
@ -61,7 +61,7 @@
 !
 Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)

-  use psb_sparse_mod, psb_protect_name => psb_ccdbldext
+  use psb_base_mod, psb_protect_name => psb_ccdbldext
  use psi_mod

 #ifdef MPI_MOD
--- a/base/tools/psb_cd_inloc.f90
+++ b/base/tools/psb_cd_inloc.f90
@ -42,7 +42,7 @@
 !    desc  - type(psb_desc_type).         The communication descriptor.
 !    info    - integer.                       Eventually returns an error code
 subroutine psb_cd_inloc(v, ictxt, desc, info, globalcheck)
-  use psb_sparse_mod
+  use psb_base_mod
  use psi_mod
  use psb_repl_map_mod
  use psb_list_map_mod
--- a/base/tools/psb_cd_lstext.f90
+++ b/base/tools/psb_cd_lstext.f90
@ -31,7 +31,7 @@
 !!$  
 Subroutine psb_cd_lstext(desc_a,in_list,desc_ov,info, mask,extype)

-  use psb_sparse_mod, psb_protect_name => psb_cd_lstext
+  use psb_base_mod, psb_protect_name => psb_cd_lstext
 !!$  use psi_mod

  Implicit None
--- a/base/tools/psb_cd_reinit.f90
+++ b/base/tools/psb_cd_reinit.f90
@ -32,7 +32,7 @@

 Subroutine psb_cd_reinit(desc,info)

-  use psb_sparse_mod, psb_protect_name => psb_cd_reinit
+  use psb_base_mod, psb_protect_name => psb_cd_reinit
  use psi_mod

  Implicit None
--- a/base/tools/psb_cd_set_bld.f90
+++ b/base/tools/psb_cd_set_bld.f90
@ -30,7 +30,7 @@
 !!$ 
 !!$  
 subroutine psb_cd_set_ovl_bld(desc,info)
-  use psb_sparse_mod, psb_protect_name => psb_cd_set_ovl_bld
+  use psb_base_mod, psb_protect_name => psb_cd_set_ovl_bld
  implicit none
  type(psb_desc_type), intent(inout) :: desc
  integer                            :: info
@ -48,7 +48,7 @@ subroutine psb_cd_set_ovl_bld(desc,info)
 end subroutine psb_cd_set_ovl_bld

 subroutine psb_cd_set_bld(desc,info)
-  use psb_sparse_mod, psb_protect_name => psb_cd_set_bld
+  use psb_base_mod, psb_protect_name => psb_cd_set_bld
  use psi_mod
  implicit none
  type(psb_desc_type), intent(inout) :: desc
--- a/base/tools/psb_cd_switch_ovl_indxmap.f90
+++ b/base/tools/psb_cd_switch_ovl_indxmap.f90
@ -34,7 +34,7 @@
 !
 Subroutine psb_cd_switch_ovl_indxmap(desc,info)

-  use psb_sparse_mod, psb_protect_name => psb_cd_switch_ovl_indxmap
+  use psb_base_mod, psb_protect_name => psb_cd_switch_ovl_indxmap
  use psi_mod


--- a/base/tools/psb_cdals.f90
+++ b/base/tools/psb_cdals.f90
@ -44,7 +44,7 @@
 !    desc  - type(psb_desc_type).         The communication descriptor.
 !    info    - integer.                       Error code (if any).
 subroutine psb_cdals(m, n, parts, ictxt, desc, info)
-  use psb_sparse_mod
+  use psb_base_mod
  use psi_mod
  use psb_repl_map_mod
  use psb_list_map_mod
--- a/base/tools/psb_cdalv.f90
+++ b/base/tools/psb_cdalv.f90
@ -44,7 +44,7 @@
 !    info    - integer.                       Return code
 !    flag    - integer.                       Are V's contents 0- or 1-based?
 subroutine psb_cdalv(v, ictxt, desc, info, flag)
-  use psb_sparse_mod
+  use psb_base_mod
  use psi_mod
  use psb_repl_map_mod
  use psb_glist_map_mod
--- a/base/tools/psb_cdcpy.F90
+++ b/base/tools/psb_cdcpy.F90
@ -39,7 +39,7 @@
  !    info     - integer.                       Return code.
 subroutine psb_cdcpy(desc_in, desc_out, info)

-  use psb_sparse_mod, psb_protect_name => psb_cdcpy
+  use psb_base_mod, psb_protect_name => psb_cdcpy

  implicit none
  !....parameters...
--- a/base/tools/psb_cdins.f90
+++ b/base/tools/psb_cdins.f90
@ -46,7 +46,7 @@
 !    jla(:)   - integer, optional              The col indices in local numbering
 !
 subroutine psb_cdinsrc(nz,ia,ja,desc_a,info,ila,jla)
-  use psb_sparse_mod, psb_protect_name => psb_cdinsrc
+  use psb_base_mod, psb_protect_name => psb_cdinsrc
  use psi_mod
  implicit none

@ -168,7 +168,7 @@ end subroutine psb_cdinsrc
 !    mask(:)  - logical, optional, target
 !
 subroutine psb_cdinsc(nz,ja,desc,info,jla,mask)
-  use psb_sparse_mod, psb_protect_name => psb_cdinsc
+  use psb_base_mod, psb_protect_name => psb_cdinsc
  use psi_mod
  implicit none

--- a/base/tools/psb_cdprt.f90
+++ b/base/tools/psb_cdprt.f90
@ -41,7 +41,7 @@
 !    short         - logical(optional).      Used to choose a verbose output.
 !
 subroutine psb_cdprt(iout,desc_p,glob,short)
-  use psb_sparse_mod, psb_protect_name => psb_cdprt
+  use psb_base_mod, psb_protect_name => psb_cdprt
  implicit none 
  type(psb_desc_type), intent(in)    :: desc_p
  integer, intent(in)                :: iout
--- a/base/tools/psb_cdren.f90
+++ b/base/tools/psb_cdren.f90
@ -45,7 +45,7 @@
 !    info     - integer.                       Return code
 !
 subroutine psb_cdren(trans,iperm,desc_a,info)
-  use psb_sparse_mod, psb_protect_name => psb_cdren
+  use psb_base_mod, psb_protect_name => psb_cdren
  use psi_mod
  implicit none

--- a/base/tools/psb_cdrep.f90
+++ b/base/tools/psb_cdrep.f90
@ -102,7 +102,7 @@
  !
  !
 subroutine psb_cdrep(m, ictxt, desc, info)
-  use psb_sparse_mod
+  use psb_base_mod
  use psi_mod
  use psb_repl_map_mod
  implicit None
--- a/base/tools/psb_cfree.f90
+++ b/base/tools/psb_cfree.f90
@ -39,7 +39,7 @@
 !    desc_a   - type(psb_desc_type).        The communication descriptor.
 !    info     - integer.                      Return code
 subroutine psb_cfree(x, desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_cfree
+  use psb_base_mod, psb_protect_name => psb_cfree
  implicit none

  !....parameters...
@ -110,7 +110,7 @@ end subroutine psb_cfree
 !    desc_a   - type(psb_desc_type).      The communication descriptor.
 !    info     - integer.                    Return code
 subroutine psb_cfreev(x, desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_cfreev
+  use psb_base_mod, psb_protect_name => psb_cfreev
  implicit none
  !....parameters...
  complex(psb_spk_),allocatable, intent(inout)    :: x(:)
--- a/base/tools/psb_cins.f90
+++ b/base/tools/psb_cins.f90
@ -46,7 +46,7 @@
 !                                     psb_dupl_ovwrt_    overwrite
 !                                     psb_dupl_add_      add         
 subroutine psb_cinsvi(m, irw, val, x, desc_a, info, dupl)
-  use psb_sparse_mod, psb_protect_name => psb_cinsvi
+  use psb_base_mod, psb_protect_name => psb_cinsvi
  use psi_mod
  implicit none

@ -227,7 +227,7 @@ end subroutine psb_cinsvi
 !                                     psb_dupl_ovwrt_    overwrite
 !                                     psb_dupl_add_      add         
 subroutine psb_cinsi(m, irw, val, x, desc_a, info, dupl)
-  use psb_sparse_mod, psb_protect_name => psb_cinsi
+  use psb_base_mod, psb_protect_name => psb_cinsi
  use psi_mod
  implicit none

--- a/base/tools/psb_cspalloc.f90
+++ b/base/tools/psb_cspalloc.f90
@ -42,7 +42,7 @@
 !                                              (local, user estimate)
 !
 subroutine psb_cspalloc(a, desc_a, info, nnz)
-  use psb_sparse_mod, psb_protect_name => psb_cspalloc
+  use psb_base_mod, psb_protect_name => psb_cspalloc
  implicit none

  !....parameters...
--- a/base/tools/psb_cspasb.f90
+++ b/base/tools/psb_cspasb.f90
@ -49,7 +49,7 @@
 ! 
 !
 subroutine psb_cspasb(a,desc_a, info, afmt, upd, dupl,mold)
-  use psb_sparse_mod, psb_protect_name => psb_cspasb
+  use psb_base_mod, psb_protect_name => psb_cspasb
  use psi_mod
  implicit none

--- a/base/tools/psb_cspfree.f90
+++ b/base/tools/psb_cspfree.f90
@ -40,7 +40,7 @@
 !    info     - integer.                          return code.
 !
 subroutine psb_cspfree(a, desc_a,info)
-  use psb_sparse_mod, psb_protect_name => psb_cspfree
+  use psb_base_mod, psb_protect_name => psb_cspfree
  implicit none

  !....parameters...
--- a/base/tools/psb_csphalo.F90
+++ b/base/tools/psb_csphalo.F90
@ -58,7 +58,7 @@
 Subroutine psb_csphalo(a,desc_a,blk,info,rowcnv,colcnv,&
     &  rowscale,colscale,outfmt,data)

-  use psb_sparse_mod, psb_protect_name => psb_csphalo
+  use psb_base_mod, psb_protect_name => psb_csphalo

 #ifdef MPI_MOD
  use mpi
--- a/base/tools/psb_cspins.f90
+++ b/base/tools/psb_cspins.f90
@ -49,7 +49,7 @@
 !                                           certain circumstances.
 !
 subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild)
-  use psb_sparse_mod, psb_protect_name => psb_cspins
+  use psb_base_mod, psb_protect_name => psb_cspins
  use psi_mod
  implicit none

@ -195,7 +195,7 @@ end subroutine psb_cspins


 subroutine psb_cspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
-  use psb_sparse_mod, psb_protect_name => psb_cspins_2desc
+  use psb_base_mod, psb_protect_name => psb_cspins_2desc
  use psi_mod
  implicit none

--- a/base/tools/psb_csprn.f90
+++ b/base/tools/psb_csprn.f90
@ -43,7 +43,7 @@
 !                                               default .true.          
 !
 Subroutine psb_csprn(a, desc_a,info,clear)
-  use psb_sparse_mod, psb_protect_name => psb_csprn
+  use psb_base_mod, psb_protect_name => psb_csprn
  Implicit None

  !....Parameters...
--- a/base/tools/psb_dallc.f90
+++ b/base/tools/psb_dallc.f90
@ -43,7 +43,7 @@
 !    n      - optional number of columns.
 !    lb     - optional lower bound on column indices
 subroutine psb_dalloc(x, desc_a, info, n, lb)
-  use psb_sparse_mod, psb_protect_name => psb_dalloc
+  use psb_base_mod, psb_protect_name => psb_dalloc
  use psi_mod
  implicit none

@ -177,7 +177,7 @@ end subroutine psb_dalloc
 !    desc_a - the communication descriptor.
 !    info   - return code
 subroutine psb_dallocv(x, desc_a,info,n)
-  use psb_sparse_mod, psb_protect_name => psb_dallocv
+  use psb_base_mod, psb_protect_name => psb_dallocv
  use psi_mod
  implicit none

--- a/base/tools/psb_dasb.f90
+++ b/base/tools/psb_dasb.f90
@ -43,7 +43,7 @@
 !    desc_a  - type(psb_desc_type).  The communication descriptor.
 !    info    - integer.                return code
 subroutine psb_dasb(x, desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_dasb
+  use psb_base_mod, psb_protect_name => psb_dasb
  implicit none

  type(psb_desc_type), intent(in) ::  desc_a
@ -176,7 +176,7 @@ end subroutine psb_dasb
 !    desc_a  - type(psb_desc_type).  The communication descriptor.
 !    info    - integer.                Return code
 subroutine psb_dasbv(x, desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_dasbv
+  use psb_base_mod, psb_protect_name => psb_dasbv
  implicit none

  type(psb_desc_type), intent(in) ::  desc_a
--- a/base/tools/psb_dcdbldext.F90
+++ b/base/tools/psb_dcdbldext.F90
@ -61,7 +61,7 @@
 !
 Subroutine psb_dcdbldext(a,desc_a,novr,desc_ov,info, extype)

-  use psb_sparse_mod, psb_protect_name => psb_dcdbldext
+  use psb_base_mod, psb_protect_name => psb_dcdbldext
  use psi_mod

 #ifdef MPI_MOD
--- a/base/tools/psb_dfree.f90
+++ b/base/tools/psb_dfree.f90
@ -39,7 +39,7 @@
 !    desc_a   - type(psb_desc_type).      The communication descriptor.
 !    info     - integer.                    Return code
 subroutine psb_dfree(x, desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_dfree
+  use psb_base_mod, psb_protect_name => psb_dfree
  implicit none

  !....parameters...
@ -110,7 +110,7 @@ end subroutine psb_dfree
 !    desc_a   - type(psb_desc_type).         The communication descriptor.
 !    info     - integer.                       Return code
 subroutine psb_dfreev(x, desc_a, info)
-  use psb_sparse_mod, psb_protect_name => psb_dfreev
+  use psb_base_mod, psb_protect_name => psb_dfreev
  implicit none
  !....parameters...
  real(psb_dpk_),allocatable, intent(inout)    :: x(:)
--- a/base/tools/psb_dins.f90
+++ b/base/tools/psb_dins.f90
@ -46,7 +46,7 @@
 !                                     psb_dupl_ovwrt_    overwrite
 !                                     psb_dupl_add_      add         
 subroutine psb_dinsvi(m, irw, val, x, desc_a, info, dupl)
-  use psb_sparse_mod, psb_protect_name => psb_dinsvi
+  use psb_base_mod, psb_protect_name => psb_dinsvi
  use psi_mod
  implicit none

@ -227,7 +227,7 @@ end subroutine psb_dinsvi
 !                                     psb_dupl_ovwrt_    overwrite
 !                                     psb_dupl_add_      add         
 subroutine psb_dinsi(m, irw, val, x, desc_a, info, dupl)
-  use psb_sparse_mod, psb_protect_name => psb_dinsi
+  use psb_base_mod, psb_protect_name => psb_dinsi
  use psi_mod
  implicit none

--- a/base/tools/psb_dspalloc.f90
+++ b/base/tools/psb_dspalloc.f90
@ -42,7 +42,7 @@
 !                                              (local, user estimate)
 !
 subroutine psb_dspalloc(a, desc_a, info, nnz)
-  use psb_sparse_mod, psb_protect_name => psb_dspalloc
+  use psb_base_mod, psb_protect_name => psb_dspalloc
  implicit none

  !....parameters...
--- a/base/tools/psb_dspasb.f90
+++ b/base/tools/psb_dspasb.f90
@ -49,7 +49,7 @@
 ! 
 !
 subroutine psb_dspasb(a,desc_a, info, afmt, upd, dupl, mold)
-  use psb_sparse_mod, psb_protect_name => psb_dspasb
+  use psb_base_mod, psb_protect_name => psb_dspasb
  use psi_mod
  implicit none

--- a/base/tools/psb_dspfree.f90
+++ b/base/tools/psb_dspfree.f90
@ -40,7 +40,7 @@
 !    info     - integer.                          return code.
 !
 subroutine psb_dspfree(a, desc_a,info)
-  use psb_sparse_mod, psb_protect_name => psb_dspfree
+  use psb_base_mod, psb_protect_name => psb_dspfree
  implicit none

  !....parameters...
--- a/base/tools/psb_dsphalo.F90
+++ b/base/tools/psb_dsphalo.F90
@ -58,7 +58,7 @@
 !
 Subroutine psb_dsphalo(a,desc_a,blk,info,rowcnv,colcnv,&
     &  rowscale,colscale,outfmt,data)
-  use psb_sparse_mod, psb_protect_name => psb_dsphalo
+  use psb_base_mod, psb_protect_name => psb_dsphalo

 #ifdef MPI_MOD
  use mpi
--- a/base/tools/psb_dspins.f90
+++ b/base/tools/psb_dspins.f90
@ -49,7 +49,7 @@
 !                                           certain circumstances.
 !
 subroutine psb_dspins(nz,ia,ja,val,a,desc_a,info,rebuild)
-  use psb_sparse_mod, psb_protect_name => psb_dspins
+  use psb_base_mod, psb_protect_name => psb_dspins
  implicit none

  !....parameters...
@ -195,7 +195,7 @@ end subroutine psb_dspins


 subroutine psb_dspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
-  use psb_sparse_mod, psb_protect_name => psb_dspins_2desc
+  use psb_base_mod, psb_protect_name => psb_dspins_2desc
  implicit none

  !....parameters...
--- a/base/tools/psb_dsprn.f90
+++ b/base/tools/psb_dsprn.f90
@ -42,7 +42,7 @@
 !    clear    - logical, optional               Whether the coefficients should be zeroed
 !                                               default .true.          
 Subroutine psb_dsprn(a, desc_a,info,clear)
-  use psb_sparse_mod, psb_protect_name => psb_dsprn
+  use psb_base_mod, psb_protect_name => psb_dsprn
  Implicit None

  !....Parameters...
--- a/base/tools/psb_get_overlap.f90
+++ b/base/tools/psb_get_overlap.f90
@ -43,7 +43,7 @@
 !    info     - integer.                  return code.
 !
 subroutine psb_get_ovrlap(ovrel,desc,info)
-  use psb_sparse_mod, psb_protect_name => psb_get_ovrlap
+  use psb_base_mod, psb_protect_name => psb_get_ovrlap
  implicit none 
  integer, allocatable, intent(out) :: ovrel(:)
  type(psb_desc_type), intent(in) :: desc
--- a/base/tools/psb_glob_to_loc.f90
+++ b/base/tools/psb_glob_to_loc.f90
@ -47,7 +47,7 @@
 !                                         owned by the process, i.e. excludes halo.
 !
 subroutine psb_glob_to_loc2(x,y,desc_a,info,iact,owned)
-  use psb_sparse_mod, psb_protect_name => psb_glob_to_loc2
+  use psb_base_mod, psb_protect_name => psb_glob_to_loc2
  use psi_mod
  implicit none

@ -166,7 +166,7 @@ end subroutine psb_glob_to_loc2
 !                                         owned by the process, i.e. excludes halo.
 !
 subroutine psb_glob_to_loc(x,desc_a,info,iact,owned)
-  use psb_sparse_mod, psb_protect_name => psb_glob_to_loc
+  use psb_base_mod, psb_protect_name => psb_glob_to_loc
  use psi_mod
  implicit none

@ -241,7 +241,7 @@ subroutine psb_glob_to_loc(x,desc_a,info,iact,owned)
 end subroutine psb_glob_to_loc

 subroutine psb_glob_to_loc2s(x,y,desc_a,info,iact,owned)
-  use psb_sparse_mod, psb_protect_name => psb_glob_to_loc2s
+  use psb_base_mod, psb_protect_name => psb_glob_to_loc2s
  implicit none 
  type(psb_desc_type), intent(in)    ::  desc_a
  integer,intent(in)                 ::  x
@ -258,7 +258,7 @@ subroutine psb_glob_to_loc2s(x,y,desc_a,info,iact,owned)
 end subroutine psb_glob_to_loc2s

 subroutine psb_glob_to_locs(x,desc_a,info,iact,owned)
-  use psb_sparse_mod, psb_protect_name => psb_glob_to_locs
+  use psb_base_mod, psb_protect_name => psb_glob_to_locs
  implicit none 
  type(psb_desc_type), intent(in)    ::  desc_a
  integer,intent(inout)              ::  x  
--- a/base/tools/psb_ialloc.f90
+++ b/base/tools/psb_ialloc.f90
@ -42,7 +42,7 @@
 !    n      - optional number of columns.
 !    lb     - optional lower bound on column indices
 subroutine psb_ialloc(x, desc_a, info, n, lb)
-  use psb_sparse_mod, psb_protect_name => psb_ialloc
+  use psb_base_mod, psb_protect_name => psb_ialloc
  implicit none
  
  !....parameters...
@ -175,7 +175,7 @@ end subroutine psb_ialloc
 !    desc_a - the communication descriptor.
 !    info   - return code
 subroutine psb_iallocv(x, desc_a, info,n)
-  use psb_sparse_mod, psb_protect_name => psb_iallocv
+  use psb_base_mod, psb_protect_name => psb_iallocv
  implicit none

  !....parameters...
--- a/Show More
+++ b/Show More