diff --git a/test/kernel/d_file_spmv.f90 b/test/kernel/d_file_spmv.f90 index ec48a68d..40a79982 100644 --- a/test/kernel/d_file_spmv.f90 +++ b/test/kernel/d_file_spmv.f90 @@ -39,7 +39,8 @@ program d_file_spmv character(len=512) :: mtrx_file, rhs_file ! sparse matrices - type(psb_dspmat_type) :: a, aux_a + type(psb_dspmat_type) :: a + type(psb_ldspmat_type) :: aux_a ! dense matrices real(psb_dpk_), allocatable, target :: aux_b(:,:), d(:) @@ -57,6 +58,7 @@ program d_file_spmv ! solver paramters integer(psb_ipk_) :: iter, itmax, ierr, itrace, ircode, ipart,& & methd, istopc, irst, nr + integer(psb_lpk_) :: lnp integer(psb_epk_) :: amatsize, descsize, annz, nbytes real(psb_dpk_) :: err, eps,cond @@ -79,7 +81,7 @@ program d_file_spmv call psb_init(ctxt) call psb_info(ctxt,iam,np) - + lnp = np if (iam < 0) then ! This should not happen, but just in case call psb_exit(ctxt) @@ -184,7 +186,7 @@ program d_file_spmv write(psb_out_unit,'("Partition type: graph")') write(psb_out_unit,'(" ")') ! write(psb_err_unit,'("Build type: graph")') - call build_mtpart(aux_a,np) + call build_mtpart(aux_a,lnp) endif call psb_barrier(ctxt) diff --git a/test/kernel/s_file_spmv.f90 b/test/kernel/s_file_spmv.f90 index 0f75518e..ebbb85cd 100644 --- a/test/kernel/s_file_spmv.f90 +++ b/test/kernel/s_file_spmv.f90 @@ -38,7 +38,8 @@ program s_file_spmv character(len=40) :: kmethd, ptype, mtrx_file, rhs_file ! sparse matrices - type(psb_sspmat_type) :: a, aux_a + type(psb_sspmat_type) :: a + type(psb_lsspmat_type) :: aux_a ! dense matrices real(psb_spk_), allocatable, target :: aux_b(:,:), d(:) @@ -56,6 +57,7 @@ program s_file_spmv ! solver paramters integer(psb_ipk_) :: iter, itmax, ierr, itrace, ircode, ipart,& & methd, istopc, irst, nr + integer(psb_lpk_) :: lnp integer(psb_epk_) :: amatsize, descsize, annz, nbytes real(psb_spk_) :: err, eps,cond @@ -183,7 +185,7 @@ program s_file_spmv write(psb_out_unit,'("Partition type: graph")') write(psb_out_unit,'(" ")') ! write(psb_err_unit,'("Build type: graph")') - call build_mtpart(aux_a,np) + call build_mtpart(aux_a,lnp) endif call psb_barrier(ctxt)