psblas3:

configure.ac
 configure
 test/fileread/df_sample.f90
 util/psb_d_renum_impl.F90
 util/psb_renum_mod.f90

Further test machinery.

configure

@@ -9349,7 +9349,7 @@ UTILLIBNAME=libpsb_util.a
 
 
 PSBLASRULES='
-PSBLDLIBS=$(LAPACK) $(BLAS) $(METIS_LIB) $(AMD_LIBS) $(LIBS)
+PSBLDLIBS=$(LAPACK) $(BLAS) $(METIS_LIB) $(AMD_LIB) $(LIBS)
 CDEFINES=$(PSBCDEFINES)
 FDEFINES=$(PSBFDEFINES)
 

configure.ac

@@ -751,7 +751,7 @@ AC_SUBST(FINCLUDES)
 
 
 PSBLASRULES='
-PSBLDLIBS=$(LAPACK) $(BLAS) $(METIS_LIB) $(AMD_LIBS) $(LIBS)
+PSBLDLIBS=$(LAPACK) $(BLAS) $(METIS_LIB) $(AMD_LIB) $(LIBS)
 CDEFINES=$(PSBCDEFINES)
 FDEFINES=$(PSBFDEFINES)
 	

test/fileread/df_sample.f90

@@ -38,7 +38,7 @@ program df_sample
   implicit none
 
   ! input parameters
-  character(len=40) :: kmethd, ptype, mtrx_file, rhs_file
+  character(len=40) :: kmethd, ptype, mtrx_file, rhs_file,renum
 
   ! sparse matrices
   type(psb_dspmat_type) :: a, aux_a
@@ -140,7 +140,7 @@ program df_sample
     
     m_problem = aux_a%get_nrows()
     call psb_bcast(ictxt,m_problem)
-    call psb_mat_renum(psb_mat_renum_gps_,aux_a,info,perm) 
+    call psb_mat_renum(psb_mat_renum_amd_,aux_a,info,perm) 
 
     ! At this point aux_b may still be unallocated
     if (psb_size(aux_b,dim=1) == m_problem) then

util/psb_d_renum_impl.F90

@@ -1,3 +1,31 @@
+subroutine psb_d_mat_renums(alg,mat,info,perm)
+  use psb_base_mod
+  use psb_renum_mod, psb_protect_name => psb_d_mat_renums
+  implicit none 
+  character(len=*), intent(in) :: alg
+  type(psb_dspmat_type), intent(inout) :: mat
+  integer, intent(out) :: info
+  integer, allocatable, optional, intent(out) :: perm(:)
+  
+  integer            :: err_act, nr, nc, ialg
+  character(len=20)  :: name
+
+  info = psb_success_
+  name = 'mat_renum'
+  call psb_erractionsave(err_act)
+
+  info = psb_success_
+  select case (psb_toupper(alg))
+  case ('GPS')
+    ialg = psb_mat_renum_gps_
+  case ('AMD')
+    ialg = psb_mat_renum_amd_
+  case default
+    ialg = -1
+  end select
+  call psb_mat_renum(ialg,mat,info,perm)
+end subroutine psb_d_mat_renums
+  
 subroutine psb_d_mat_renum(alg,mat,info,perm)
   use psb_base_mod
   use psb_renum_mod, psb_protect_name => psb_d_mat_renum

util/psb_renum_mod.f90

@@ -6,6 +6,13 @@ module psb_renum_mod
 
   
   interface psb_mat_renum
+    subroutine psb_d_mat_renums(alg,mat,info,perm)
+      import  psb_dspmat_type
+      character(len=*), intent(in) :: alg
+      type(psb_dspmat_type), intent(inout) :: mat
+      integer, intent(out) :: info
+      integer, allocatable, optional, intent(out) :: perm(:)
+    end subroutine psb_d_mat_renums
     subroutine psb_d_mat_renum(alg,mat,info,perm)
       import  psb_dspmat_type
       integer, intent(in) :: alg