|
|
|
@ -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)
|
|
|
|
|