fdurastante/psblas3
1
0
Fork 0
mirror of https://github.com/sfilippone/psblas3.git synced 2026-10-18 11:52:11 +00:00
Code Issues Packages Projects Releases Wiki Activity

New fnd_owner.

Browse Source
This commit is contained in:
Salvatore Filippone
2022-03-31 17:16:08 +02:00
parent d6c232c676
commit b4c27ec4eb
8 changed files with 28 additions and 48 deletions
Download Patch File Download Diff File Expand all files Collapse all files
base/tools/psb_c_remote_mat.F90
+7 -11
View File
@@ -95,7 +95,7 @@ Subroutine psb_lc_remote_mat(a,desc_a,b,info)
integer(psb_lpk_) :: r, k, irmin, irmax, icmin, icmax, iszs, iszr, &
& lidx, l1, lnr, lnc, lnnz, idx, ngtz, tot_elem
integer(psb_lpk_) :: nz,nouth
integer(psb_ipk_) :: nnp, nrcvs, nsnds
integer(psb_ipk_) :: nrcvs, nsnds
integer(psb_mpk_) :: icomm, minfo
integer(psb_mpk_), allocatable :: brvindx(:), &
& rvsz(:), bsdindx(:),sdsz(:), sdsi(:), rvsi(:)
@@ -103,7 +103,7 @@ Subroutine psb_lc_remote_mat(a,desc_a,b,info)
complex(psb_spk_), allocatable :: valsnd(:)
type(psb_lc_coo_sparse_mat), allocatable :: acoo
class(psb_i_base_vect_type), pointer :: pdxv
integer(psb_ipk_), allocatable :: ladj(:), ila(:), iprc(:)
integer(psb_ipk_), allocatable :: ila(:), iprc(:)
logical :: rowcnv_,colcnv_,rowscale_,colscale_
character(len=5) :: outfmt_
integer(psb_ipk_) :: debug_level, debug_unit, err_act
@@ -149,9 +149,7 @@ Subroutine psb_lc_remote_mat(a,desc_a,b,info)
if ((nouth/=0).and.(me==0)) &
& write(0,*) 'Warning: would require reinit of DESC_A'
call psi_graph_fnd_owner(a%ia(1:nz),iprc,ladj,desc_a%indxmap,info)
call psb_msort_unique(ladj,nnp)
!write(0,*) me,name,' Processes:',ladj(1:nnp)
call desc_a%indxmap%fnd_owner(a%ia(1:nz),iprc,info)
icomm = desc_a%get_mpic()
sdsz(:)=0
@@ -304,13 +302,13 @@ subroutine psb_c_remote_vect(v,desc_a, info, dupl)
integer(psb_lpk_) :: r, k, irmin, irmax, icmin, icmax, iszs, iszr, &
& lidx, l1, lnr, lnc, lnnz, idx, ngtz, tot_elem
integer(psb_lpk_) :: nz,nouth
integer(psb_ipk_) :: nnp, nrcvs, nsnds
integer(psb_ipk_) :: nrcvs, nsnds
integer(psb_mpk_) :: icomm, minfo
integer(psb_mpk_), allocatable :: brvindx(:), &
& rvsz(:), bsdindx(:),sdsz(:), sdsi(:), rvsi(:)
integer(psb_lpk_), allocatable :: iasnd(:), jasnd(:)
complex(psb_spk_), allocatable :: valsnd(:)
integer(psb_ipk_), allocatable :: ladj(:), ila(:), iprc(:)
integer(psb_ipk_), allocatable :: ila(:), iprc(:)
logical :: rowcnv_,colcnv_,rowscale_,colscale_
character(len=5) :: outfmt_
integer(psb_ipk_) :: debug_level, debug_unit, err_act
@@ -331,7 +329,7 @@ subroutine psb_c_remote_vect(v,desc_a, info, dupl)
if (v%is_remote_build()) then
dupl_ = psb_dupl_add_
else
dupl_ = psb_dupl_ovwrt_
dupl_ = psb_dupl_def_
end if
endif
@@ -366,9 +364,7 @@ subroutine psb_c_remote_vect(v,desc_a, info, dupl)
!!$ if ((nouth/=0).and.(me==0)) &
!!$ & write(0,*) 'Warning: would require reinit of DESC_A'
!!$
!!$ call psi_graph_fnd_owner(a%ia(1:nz),iprc,ladj,desc_a%indxmap,info)
!!$ call psb_msort_unique(ladj,nnp)
!!$ !write(0,*) me,name,' Processes:',ladj(1:nnp)
!!$ call desc_a%indxmap%fnd_owner(a%ia(1:nz),iprc,info)
!!$
!!$ icomm = desc_a%get_mpic()
!!$ sdsz(:)=0
base/tools/psb_cspasb.f90
-1
View File
@@ -94,7 +94,6 @@ subroutine psb_cspasb(a,desc_a, info, afmt, upd, dupl, mold)
goto 9999
end if
if (debug_level >= psb_debug_ext_)&
& write(debug_unit, *) me,' ',trim(name),&
& ' Begin matrix assembly...'
base/tools/psb_d_remote_mat.F90
+7 -11
View File
@@ -95,7 +95,7 @@ Subroutine psb_ld_remote_mat(a,desc_a,b,info)
integer(psb_lpk_) :: r, k, irmin, irmax, icmin, icmax, iszs, iszr, &
& lidx, l1, lnr, lnc, lnnz, idx, ngtz, tot_elem
integer(psb_lpk_) :: nz,nouth
integer(psb_ipk_) :: nnp, nrcvs, nsnds
integer(psb_ipk_) :: nrcvs, nsnds
integer(psb_mpk_) :: icomm, minfo
integer(psb_mpk_), allocatable :: brvindx(:), &
& rvsz(:), bsdindx(:),sdsz(:), sdsi(:), rvsi(:)
@@ -103,7 +103,7 @@ Subroutine psb_ld_remote_mat(a,desc_a,b,info)
real(psb_dpk_), allocatable :: valsnd(:)
type(psb_ld_coo_sparse_mat), allocatable :: acoo
class(psb_i_base_vect_type), pointer :: pdxv
integer(psb_ipk_), allocatable :: ladj(:), ila(:), iprc(:)
integer(psb_ipk_), allocatable :: ila(:), iprc(:)
logical :: rowcnv_,colcnv_,rowscale_,colscale_
character(len=5) :: outfmt_
integer(psb_ipk_) :: debug_level, debug_unit, err_act
@@ -149,9 +149,7 @@ Subroutine psb_ld_remote_mat(a,desc_a,b,info)
if ((nouth/=0).and.(me==0)) &
& write(0,*) 'Warning: would require reinit of DESC_A'
call psi_graph_fnd_owner(a%ia(1:nz),iprc,ladj,desc_a%indxmap,info)
call psb_msort_unique(ladj,nnp)
!write(0,*) me,name,' Processes:',ladj(1:nnp)
call desc_a%indxmap%fnd_owner(a%ia(1:nz),iprc,info)
icomm = desc_a%get_mpic()
sdsz(:)=0
@@ -304,13 +302,13 @@ subroutine psb_d_remote_vect(v,desc_a, info, dupl)
integer(psb_lpk_) :: r, k, irmin, irmax, icmin, icmax, iszs, iszr, &
& lidx, l1, lnr, lnc, lnnz, idx, ngtz, tot_elem
integer(psb_lpk_) :: nz,nouth
integer(psb_ipk_) :: nnp, nrcvs, nsnds
integer(psb_ipk_) :: nrcvs, nsnds
integer(psb_mpk_) :: icomm, minfo
integer(psb_mpk_), allocatable :: brvindx(:), &
& rvsz(:), bsdindx(:),sdsz(:), sdsi(:), rvsi(:)
integer(psb_lpk_), allocatable :: iasnd(:), jasnd(:)
real(psb_dpk_), allocatable :: valsnd(:)
integer(psb_ipk_), allocatable :: ladj(:), ila(:), iprc(:)
integer(psb_ipk_), allocatable :: ila(:), iprc(:)
logical :: rowcnv_,colcnv_,rowscale_,colscale_
character(len=5) :: outfmt_
integer(psb_ipk_) :: debug_level, debug_unit, err_act
@@ -331,7 +329,7 @@ subroutine psb_d_remote_vect(v,desc_a, info, dupl)
if (v%is_remote_build()) then
dupl_ = psb_dupl_add_
else
dupl_ = psb_dupl_ovwrt_
dupl_ = psb_dupl_def_
end if
endif
@@ -366,9 +364,7 @@ subroutine psb_d_remote_vect(v,desc_a, info, dupl)
!!$ if ((nouth/=0).and.(me==0)) &
!!$ & write(0,*) 'Warning: would require reinit of DESC_A'
!!$
!!$ call psi_graph_fnd_owner(a%ia(1:nz),iprc,ladj,desc_a%indxmap,info)
!!$ call psb_msort_unique(ladj,nnp)
!!$ !write(0,*) me,name,' Processes:',ladj(1:nnp)
!!$ call desc_a%indxmap%fnd_owner(a%ia(1:nz),iprc,info)
!!$
!!$ icomm = desc_a%get_mpic()
!!$ sdsz(:)=0
base/tools/psb_dspasb.f90
-1
View File
@@ -94,7 +94,6 @@ subroutine psb_dspasb(a,desc_a, info, afmt, upd, dupl, mold)
goto 9999
end if
if (debug_level >= psb_debug_ext_)&
& write(debug_unit, *) me,' ',trim(name),&
& ' Begin matrix assembly...'
base/tools/psb_s_remote_mat.F90
+7 -11
View File
@@ -95,7 +95,7 @@ Subroutine psb_ls_remote_mat(a,desc_a,b,info)
integer(psb_lpk_) :: r, k, irmin, irmax, icmin, icmax, iszs, iszr, &
& lidx, l1, lnr, lnc, lnnz, idx, ngtz, tot_elem
integer(psb_lpk_) :: nz,nouth
integer(psb_ipk_) :: nnp, nrcvs, nsnds
integer(psb_ipk_) :: nrcvs, nsnds
integer(psb_mpk_) :: icomm, minfo
integer(psb_mpk_), allocatable :: brvindx(:), &
& rvsz(:), bsdindx(:),sdsz(:), sdsi(:), rvsi(:)
@@ -103,7 +103,7 @@ Subroutine psb_ls_remote_mat(a,desc_a,b,info)
real(psb_spk_), allocatable :: valsnd(:)
type(psb_ls_coo_sparse_mat), allocatable :: acoo
class(psb_i_base_vect_type), pointer :: pdxv
integer(psb_ipk_), allocatable :: ladj(:), ila(:), iprc(:)
integer(psb_ipk_), allocatable :: ila(:), iprc(:)
logical :: rowcnv_,colcnv_,rowscale_,colscale_
character(len=5) :: outfmt_
integer(psb_ipk_) :: debug_level, debug_unit, err_act
@@ -149,9 +149,7 @@ Subroutine psb_ls_remote_mat(a,desc_a,b,info)
if ((nouth/=0).and.(me==0)) &
& write(0,*) 'Warning: would require reinit of DESC_A'
call psi_graph_fnd_owner(a%ia(1:nz),iprc,ladj,desc_a%indxmap,info)
call psb_msort_unique(ladj,nnp)
!write(0,*) me,name,' Processes:',ladj(1:nnp)
call desc_a%indxmap%fnd_owner(a%ia(1:nz),iprc,info)
icomm = desc_a%get_mpic()
sdsz(:)=0
@@ -304,13 +302,13 @@ subroutine psb_s_remote_vect(v,desc_a, info, dupl)
integer(psb_lpk_) :: r, k, irmin, irmax, icmin, icmax, iszs, iszr, &
& lidx, l1, lnr, lnc, lnnz, idx, ngtz, tot_elem
integer(psb_lpk_) :: nz,nouth
integer(psb_ipk_) :: nnp, nrcvs, nsnds
integer(psb_ipk_) :: nrcvs, nsnds
integer(psb_mpk_) :: icomm, minfo
integer(psb_mpk_), allocatable :: brvindx(:), &
& rvsz(:), bsdindx(:),sdsz(:), sdsi(:), rvsi(:)
integer(psb_lpk_), allocatable :: iasnd(:), jasnd(:)
real(psb_spk_), allocatable :: valsnd(:)
integer(psb_ipk_), allocatable :: ladj(:), ila(:), iprc(:)
integer(psb_ipk_), allocatable :: ila(:), iprc(:)
logical :: rowcnv_,colcnv_,rowscale_,colscale_
character(len=5) :: outfmt_
integer(psb_ipk_) :: debug_level, debug_unit, err_act
@@ -331,7 +329,7 @@ subroutine psb_s_remote_vect(v,desc_a, info, dupl)
if (v%is_remote_build()) then
dupl_ = psb_dupl_add_
else
dupl_ = psb_dupl_ovwrt_
dupl_ = psb_dupl_def_
end if
endif
@@ -366,9 +364,7 @@ subroutine psb_s_remote_vect(v,desc_a, info, dupl)
!!$ if ((nouth/=0).and.(me==0)) &
!!$ & write(0,*) 'Warning: would require reinit of DESC_A'
!!$
!!$ call psi_graph_fnd_owner(a%ia(1:nz),iprc,ladj,desc_a%indxmap,info)
!!$ call psb_msort_unique(ladj,nnp)
!!$ !write(0,*) me,name,' Processes:',ladj(1:nnp)
!!$ call desc_a%indxmap%fnd_owner(a%ia(1:nz),iprc,info)
!!$
!!$ icomm = desc_a%get_mpic()
!!$ sdsz(:)=0
base/tools/psb_sspasb.f90
-1
View File
@@ -94,7 +94,6 @@ subroutine psb_sspasb(a,desc_a, info, afmt, upd, dupl, mold)
goto 9999
end if
if (debug_level >= psb_debug_ext_)&
& write(debug_unit, *) me,' ',trim(name),&
& ' Begin matrix assembly...'
base/tools/psb_z_remote_mat.F90
+7 -11
View File
@@ -95,7 +95,7 @@ Subroutine psb_lz_remote_mat(a,desc_a,b,info)
integer(psb_lpk_) :: r, k, irmin, irmax, icmin, icmax, iszs, iszr, &
& lidx, l1, lnr, lnc, lnnz, idx, ngtz, tot_elem
integer(psb_lpk_) :: nz,nouth
integer(psb_ipk_) :: nnp, nrcvs, nsnds
integer(psb_ipk_) :: nrcvs, nsnds
integer(psb_mpk_) :: icomm, minfo
integer(psb_mpk_), allocatable :: brvindx(:), &
& rvsz(:), bsdindx(:),sdsz(:), sdsi(:), rvsi(:)
@@ -103,7 +103,7 @@ Subroutine psb_lz_remote_mat(a,desc_a,b,info)
complex(psb_dpk_), allocatable :: valsnd(:)
type(psb_lz_coo_sparse_mat), allocatable :: acoo
class(psb_i_base_vect_type), pointer :: pdxv
integer(psb_ipk_), allocatable :: ladj(:), ila(:), iprc(:)
integer(psb_ipk_), allocatable :: ila(:), iprc(:)
logical :: rowcnv_,colcnv_,rowscale_,colscale_
character(len=5) :: outfmt_
integer(psb_ipk_) :: debug_level, debug_unit, err_act
@@ -149,9 +149,7 @@ Subroutine psb_lz_remote_mat(a,desc_a,b,info)
if ((nouth/=0).and.(me==0)) &
& write(0,*) 'Warning: would require reinit of DESC_A'
call psi_graph_fnd_owner(a%ia(1:nz),iprc,ladj,desc_a%indxmap,info)
call psb_msort_unique(ladj,nnp)
!write(0,*) me,name,' Processes:',ladj(1:nnp)
call desc_a%indxmap%fnd_owner(a%ia(1:nz),iprc,info)
icomm = desc_a%get_mpic()
sdsz(:)=0
@@ -304,13 +302,13 @@ subroutine psb_z_remote_vect(v,desc_a, info, dupl)
integer(psb_lpk_) :: r, k, irmin, irmax, icmin, icmax, iszs, iszr, &
& lidx, l1, lnr, lnc, lnnz, idx, ngtz, tot_elem
integer(psb_lpk_) :: nz,nouth
integer(psb_ipk_) :: nnp, nrcvs, nsnds
integer(psb_ipk_) :: nrcvs, nsnds
integer(psb_mpk_) :: icomm, minfo
integer(psb_mpk_), allocatable :: brvindx(:), &
& rvsz(:), bsdindx(:),sdsz(:), sdsi(:), rvsi(:)
integer(psb_lpk_), allocatable :: iasnd(:), jasnd(:)
complex(psb_dpk_), allocatable :: valsnd(:)
integer(psb_ipk_), allocatable :: ladj(:), ila(:), iprc(:)
integer(psb_ipk_), allocatable :: ila(:), iprc(:)
logical :: rowcnv_,colcnv_,rowscale_,colscale_
character(len=5) :: outfmt_
integer(psb_ipk_) :: debug_level, debug_unit, err_act
@@ -331,7 +329,7 @@ subroutine psb_z_remote_vect(v,desc_a, info, dupl)
if (v%is_remote_build()) then
dupl_ = psb_dupl_add_
else
dupl_ = psb_dupl_ovwrt_
dupl_ = psb_dupl_def_
end if
endif
@@ -366,9 +364,7 @@ subroutine psb_z_remote_vect(v,desc_a, info, dupl)
!!$ if ((nouth/=0).and.(me==0)) &
!!$ & write(0,*) 'Warning: would require reinit of DESC_A'
!!$
!!$ call psi_graph_fnd_owner(a%ia(1:nz),iprc,ladj,desc_a%indxmap,info)
!!$ call psb_msort_unique(ladj,nnp)
!!$ !write(0,*) me,name,' Processes:',ladj(1:nnp)
!!$ call desc_a%indxmap%fnd_owner(a%ia(1:nz),iprc,info)
!!$
!!$ icomm = desc_a%get_mpic()
!!$ sdsz(:)=0
base/tools/psb_zspasb.f90
-1
View File
@@ -94,7 +94,6 @@ subroutine psb_zspasb(a,desc_a, info, afmt, upd, dupl, mold)
goto 9999
end if
if (debug_level >= psb_debug_ext_)&
& write(debug_unit, *) me,' ',trim(name),&
& ' Begin matrix assembly...'