fdurastante
/
psblas3
mirror of https://github.com/sfilippone/psblas3.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

psblas3:

Browse Source 
fixes for:
1. Interface problems uncovered by Cray compiler;
2. Use desc_data% methods instead of direct calls.
psblas3-type-indexed
Salvatore Filippone 14 years ago
parent c05ed63026
commit f9e16f8340
195 changed files with 1224 additions and 1230 deletions
Show Stats Download Patch File Download Diff File

16
base/comm/psb_cgather.f90
Unescape Escape View File

@ -68,7 +68,7 @@ subroutine psb_cgatherm(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -103,8 +103,8 @@ subroutine psb_cgatherm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx,1)
lda_locx = size(locx, 1)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
lock=size(locx,2)-jlocx+1
globk=size(globx,2)-jglobx+1
@ -134,7 +134,7 @@ subroutine psb_cgatherm(globx, locx, desc_a, info, iroot)
globx(:,:)=0.d0
do j=1,k
do i=1,psb_cd_get_local_rows(desc_a)
do i=1,desc_a%get_local_rows()
call psb_loc_to_glob(i,idx,desc_a,info)
globx(idx,jglobx+j-1) = locx(i,jlx+j-1)
end do
@ -241,7 +241,7 @@ subroutine psb_cgatherv(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -271,8 +271,8 @@ subroutine psb_cgatherv(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
k = 1
@ -297,7 +297,7 @@ subroutine psb_cgatherv(globx, locx, desc_a, info, iroot)
globx(:)=0.d0
do i=1,psb_cd_get_local_rows(desc_a)
do i=1,desc_a%get_local_rows()
call psb_loc_to_glob(i,idx,desc_a,info)
globx(idx) = locx(i)
end do

16
base/comm/psb_chalo.f90
Unescape Escape View File

@ -80,7 +80,7 @@ subroutine psb_chalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -97,9 +97,9 @@ subroutine psb_chalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
maxk=size(x,2)-ijx+1
@ -301,7 +301,7 @@ subroutine psb_chalov(x,desc_a,info,alpha,work,tran,mode,data)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -314,9 +314,9 @@ subroutine psb_chalov(x,desc_a,info,alpha,work,tran,mode,data)
ix = 1
ijx = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
if (present(tran)) then
tran_ = psb_toupper(tran)

20
base/comm/psb_covrl.f90
Unescape Escape View File

@ -89,7 +89,7 @@ subroutine psb_covrlm(x,desc_a,info,jx,ik,work,update,mode)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -106,10 +106,10 @@ subroutine psb_covrlm(x,desc_a,info,jx,ik,work,update,mode)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
maxk=size(x,2)-ijx+1
@ -288,7 +288,7 @@ subroutine psb_covrlv(x,desc_a,info,work,update,mode)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -301,10 +301,10 @@ subroutine psb_covrlv(x,desc_a,info,work,update,mode)
ix = 1
ijx = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
k = 1

20
base/comm/psb_cscatter.F90
Unescape Escape View File

@ -75,7 +75,7 @@ subroutine psb_cscatterm(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -107,8 +107,8 @@ subroutine psb_cscatterm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx,1)
lda_locx = size(locx, 1)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
lock=size(locx,2)-jlocx+1
globk=size(globx,2)-jglobx+1
@ -121,8 +121,8 @@ subroutine psb_cscatterm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
call psb_bcast(ictxt,k,root=iiroot)
@ -143,7 +143,7 @@ subroutine psb_cscatterm(globx, locx, desc_a, info, iroot)
goto 9999
end if
nrow=psb_cd_get_local_rows(desc_a)
nrow=desc_a%get_local_rows()
if ((root == -1).or.(np == 1)) then
! extract my chunk
@ -308,7 +308,7 @@ subroutine psb_cscatterv(globx, locx, desc_a, info, iroot)
if (psb_get_errstatus() /= 0) return
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
@ -343,8 +343,8 @@ subroutine psb_cscatterv(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
k = 1
! there should be a global check on k here!!!
@ -365,7 +365,7 @@ subroutine psb_cscatterv(globx, locx, desc_a, info, iroot)
goto 9999
end if
nrow = psb_cd_get_local_rows(desc_a)
nrow = desc_a%get_local_rows()
if ((root == -1).or.(np == 1)) then
! extract my chunk

8
base/comm/psb_cspgather.F90
Unescape Escape View File

@ -30,8 +30,8 @@ subroutine psb_csp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keep
info=psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt, me, np)
if (present(keepnum)) then
@ -47,8 +47,8 @@ subroutine psb_csp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keep
call globa%free()
if (keepnum_) then
nrg = psb_cd_get_global_rows(desc_a)
ncg = psb_cd_get_global_rows(desc_a)
nrg = desc_a%get_global_rows()
ncg = desc_a%get_global_rows()
allocate(nzbr(np), idisp(np),stat=info)
if (info /= psb_success_) then

16
base/comm/psb_dgather.f90
Unescape Escape View File

@ -68,7 +68,7 @@ subroutine psb_dgatherm(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -103,8 +103,8 @@ subroutine psb_dgatherm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx,1)
lda_locx = size(locx, 1)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
lock=size(locx,2)-jlocx+1
globk=size(globx,2)-jglobx+1
@ -134,7 +134,7 @@ subroutine psb_dgatherm(globx, locx, desc_a, info, iroot)
globx(:,:)=0.d0
do j=1,k
do i=1,psb_cd_get_local_rows(desc_a)
do i=1,desc_a%get_local_rows()
call psb_loc_to_glob(i,idx,desc_a,info)
globx(idx,jglobx+j-1) = locx(i,jlx+j-1)
end do
@ -240,7 +240,7 @@ subroutine psb_dgatherv(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -270,8 +270,8 @@ subroutine psb_dgatherv(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
k = 1
@ -296,7 +296,7 @@ subroutine psb_dgatherv(globx, locx, desc_a, info, iroot)
globx(:)=0.d0
do i=1,psb_cd_get_local_rows(desc_a)
do i=1,desc_a%get_local_rows()
call psb_loc_to_glob(i,idx,desc_a,info)
globx(idx) = locx(i)
end do

16
base/comm/psb_dhalo.f90
Unescape Escape View File

@ -80,7 +80,7 @@ subroutine psb_dhalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -97,9 +97,9 @@ subroutine psb_dhalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
maxk=size(x,2)-ijx+1
@ -301,7 +301,7 @@ subroutine psb_dhalov(x,desc_a,info,alpha,work,tran,mode,data)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -314,9 +314,9 @@ subroutine psb_dhalov(x,desc_a,info,alpha,work,tran,mode,data)
ix = 1
ijx = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
if (present(tran)) then
tran_ = psb_toupper(tran)

20
base/comm/psb_dovrl.f90
Unescape Escape View File

@ -88,7 +88,7 @@ subroutine psb_dovrlm(x,desc_a,info,jx,ik,work,update,mode)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -105,10 +105,10 @@ subroutine psb_dovrlm(x,desc_a,info,jx,ik,work,update,mode)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
maxk=size(x,2)-ijx+1
@ -290,7 +290,7 @@ subroutine psb_dovrlv(x,desc_a,info,work,update,mode)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -303,10 +303,10 @@ subroutine psb_dovrlv(x,desc_a,info,work,update,mode)
ix = 1
ijx = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
k = 1

20
base/comm/psb_dscatter.F90
Unescape Escape View File

@ -75,7 +75,7 @@ subroutine psb_dscatterm(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -107,8 +107,8 @@ subroutine psb_dscatterm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx,1)
lda_locx = size(locx, 1)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
lock=size(locx,2)-jlocx+1
globk=size(globx,2)-jglobx+1
@ -121,8 +121,8 @@ subroutine psb_dscatterm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
call psb_bcast(ictxt,k,root=iiroot)
@ -143,7 +143,7 @@ subroutine psb_dscatterm(globx, locx, desc_a, info, iroot)
goto 9999
end if
nrow=psb_cd_get_local_rows(desc_a)
nrow=desc_a%get_local_rows()
if ((root == -1).or.(np == 1)) then
! extract my chunk
@ -309,7 +309,7 @@ subroutine psb_dscatterv(globx, locx, desc_a, info, iroot)
if (psb_get_errstatus() /= 0) return
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
@ -344,8 +344,8 @@ subroutine psb_dscatterv(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
k = 1
! there should be a global check on k here!!!
@ -366,7 +366,7 @@ subroutine psb_dscatterv(globx, locx, desc_a, info, iroot)
goto 9999
end if
nrow = psb_cd_get_local_rows(desc_a)
nrow = desc_a%get_local_rows()
if ((root == -1).or.(np == 1)) then
! extract my chunk

8
base/comm/psb_dspgather.F90
Unescape Escape View File

@ -30,8 +30,8 @@ subroutine psb_dsp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keep
info=psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt, me, np)
if (present(keepnum)) then
@ -47,8 +47,8 @@ subroutine psb_dsp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keep
call globa%free()
if (keepnum_) then
nrg = psb_cd_get_global_rows(desc_a)
ncg = psb_cd_get_global_rows(desc_a)
nrg = desc_a%get_global_rows()
ncg = desc_a%get_global_rows()
allocate(nzbr(np), idisp(np),stat=info)
if (info /= psb_success_) then

16
base/comm/psb_igather.f90
Unescape Escape View File

@ -68,7 +68,7 @@ subroutine psb_igatherm(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -103,8 +103,8 @@ subroutine psb_igatherm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx,1)
lda_locx = size(locx, 1)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
lock=size(locx,2)-jlocx+1
globk=size(globx,2)-jglobx+1
@ -134,7 +134,7 @@ subroutine psb_igatherm(globx, locx, desc_a, info, iroot)
globx(:,:)=0
do j=1,k
do i=1,psb_cd_get_local_rows(desc_a)
do i=1,desc_a%get_local_rows()
call psb_loc_to_glob(i,idx,desc_a,info)
globx(idx,jglobx+j-1) = locx(i,jlx+j-1)
end do
@ -240,7 +240,7 @@ subroutine psb_igatherv(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -270,8 +270,8 @@ subroutine psb_igatherv(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
k = 1
@ -296,7 +296,7 @@ subroutine psb_igatherv(globx, locx, desc_a, info, iroot)
globx(:)=0
do i=1,psb_cd_get_local_rows(desc_a)
do i=1,desc_a%get_local_rows()
call psb_loc_to_glob(i,idx,desc_a,info)
globx(idx) = locx(i)
end do

18
base/comm/psb_ihalo.f90
Unescape Escape View File

@ -81,7 +81,7 @@ subroutine psb_ihalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -98,9 +98,9 @@ subroutine psb_ihalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
maxk=size(x,2)-ijx+1
@ -305,7 +305,7 @@ subroutine psb_ihalov(x,desc_a,info,alpha,work,tran,mode,data)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -318,10 +318,10 @@ subroutine psb_ihalov(x,desc_a,info,alpha,work,tran,mode,data)
ix = 1
ijx = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
! ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
! ncol = desc_a%get_local_cols()
if (present(tran)) then

20
base/comm/psb_iovrl.f90
Unescape Escape View File

@ -88,7 +88,7 @@ subroutine psb_iovrlm(x,desc_a,info,jx,ik,work,update,mode)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -105,10 +105,10 @@ subroutine psb_iovrlm(x,desc_a,info,jx,ik,work,update,mode)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
maxk=size(x,2)-ijx+1
@ -289,7 +289,7 @@ subroutine psb_iovrlv(x,desc_a,info,work,update,mode)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -302,10 +302,10 @@ subroutine psb_iovrlv(x,desc_a,info,work,update,mode)
ix = 1
ijx = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
k = 1

20
base/comm/psb_iscatter.F90
Unescape Escape View File

@ -73,7 +73,7 @@ subroutine psb_iscatterm(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -105,8 +105,8 @@ subroutine psb_iscatterm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx,1)
lda_locx = size(locx, 1)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
lock=size(locx,2)-jlocx+1
globk=size(globx,2)-jglobx+1
@ -119,8 +119,8 @@ subroutine psb_iscatterm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
call psb_bcast(ictxt,k,root=iiroot)
@ -141,7 +141,7 @@ subroutine psb_iscatterm(globx, locx, desc_a, info, iroot)
goto 9999
end if
nrow=psb_cd_get_local_rows(desc_a)
nrow=desc_a%get_local_rows()
if ((root == -1).or.(np == 1)) then
! extract my chunk
@ -306,7 +306,7 @@ subroutine psb_iscatterv(globx, locx, desc_a, info, iroot)
if (psb_get_errstatus() /= 0) return
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
@ -341,8 +341,8 @@ subroutine psb_iscatterv(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
k = 1
! there should be a global check on k here!!!
@ -363,7 +363,7 @@ subroutine psb_iscatterv(globx, locx, desc_a, info, iroot)
goto 9999
end if
nrow = psb_cd_get_local_rows(desc_a)
nrow = desc_a%get_local_rows()
if ((root == -1).or.(np == 1)) then
! extract my chunk

16
base/comm/psb_sgather.f90
Unescape Escape View File

@ -68,7 +68,7 @@ subroutine psb_sgatherm(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -103,8 +103,8 @@ subroutine psb_sgatherm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx,1)
lda_locx = size(locx, 1)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
lock=size(locx,2)-jlocx+1
globk=size(globx,2)-jglobx+1
@ -134,7 +134,7 @@ subroutine psb_sgatherm(globx, locx, desc_a, info, iroot)
globx(:,:)=0.d0
do j=1,k
do i=1,psb_cd_get_local_rows(desc_a)
do i=1,desc_a%get_local_rows()
call psb_loc_to_glob(i,idx,desc_a,info)
globx(idx,jglobx+j-1) = locx(i,jlx+j-1)
end do
@ -240,7 +240,7 @@ subroutine psb_sgatherv(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -270,8 +270,8 @@ subroutine psb_sgatherv(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
k = 1
@ -296,7 +296,7 @@ subroutine psb_sgatherv(globx, locx, desc_a, info, iroot)
globx(:)=0.d0
do i=1,psb_cd_get_local_rows(desc_a)
do i=1,desc_a%get_local_rows()
call psb_loc_to_glob(i,idx,desc_a,info)
globx(idx) = locx(i)
end do

16
base/comm/psb_shalo.f90
Unescape Escape View File

@ -80,7 +80,7 @@ subroutine psb_shalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -97,9 +97,9 @@ subroutine psb_shalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
maxk=size(x,2)-ijx+1
@ -301,7 +301,7 @@ subroutine psb_shalov(x,desc_a,info,alpha,work,tran,mode,data)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -314,9 +314,9 @@ subroutine psb_shalov(x,desc_a,info,alpha,work,tran,mode,data)
ix = 1
ijx = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
if (present(tran)) then
tran_ = psb_toupper(tran)

20
base/comm/psb_sovrl.f90
Unescape Escape View File

@ -88,7 +88,7 @@ subroutine psb_sovrlm(x,desc_a,info,jx,ik,work,update,mode)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -105,10 +105,10 @@ subroutine psb_sovrlm(x,desc_a,info,jx,ik,work,update,mode)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
maxk=size(x,2)-ijx+1
@ -290,7 +290,7 @@ subroutine psb_sovrlv(x,desc_a,info,work,update,mode)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -303,10 +303,10 @@ subroutine psb_sovrlv(x,desc_a,info,work,update,mode)
ix = 1
ijx = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
k = 1

20
base/comm/psb_sscatter.F90
Unescape Escape View File

@ -75,7 +75,7 @@ subroutine psb_sscatterm(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -107,8 +107,8 @@ subroutine psb_sscatterm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx,1)
lda_locx = size(locx, 1)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
lock=size(locx,2)-jlocx+1
globk=size(globx,2)-jglobx+1
@ -121,8 +121,8 @@ subroutine psb_sscatterm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
call psb_bcast(ictxt,k,root=iiroot)
@ -143,7 +143,7 @@ subroutine psb_sscatterm(globx, locx, desc_a, info, iroot)
goto 9999
end if
nrow=psb_cd_get_local_rows(desc_a)
nrow=desc_a%get_local_rows()
if ((root == -1).or.(np == 1)) then
! extract my chunk
@ -308,7 +308,7 @@ subroutine psb_sscatterv(globx, locx, desc_a, info, iroot)
if (psb_get_errstatus() /= 0) return
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
@ -343,8 +343,8 @@ subroutine psb_sscatterv(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
k = 1
! there should be a global check on k here!!!
@ -365,7 +365,7 @@ subroutine psb_sscatterv(globx, locx, desc_a, info, iroot)
goto 9999
end if
nrow = psb_cd_get_local_rows(desc_a)
nrow = desc_a%get_local_rows()
if ((root == -1).or.(np == 1)) then
! extract my chunk

8
base/comm/psb_sspgather.F90
Unescape Escape View File

@ -30,8 +30,8 @@ subroutine psb_ssp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keep
info=psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt, me, np)
if (present(keepnum)) then
@ -47,8 +47,8 @@ subroutine psb_ssp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keep
call globa%free()
if (keepnum_) then
nrg = psb_cd_get_global_rows(desc_a)
ncg = psb_cd_get_global_rows(desc_a)
nrg = desc_a%get_global_rows()
ncg = desc_a%get_global_rows()
allocate(nzbr(np), idisp(np),stat=info)
if (info /= psb_success_) then

16
base/comm/psb_zgather.f90
Unescape Escape View File

@ -68,7 +68,7 @@ subroutine psb_zgatherm(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -103,8 +103,8 @@ subroutine psb_zgatherm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx,1)
lda_locx = size(locx, 1)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
lock=size(locx,2)-jlocx+1
globk=size(globx,2)-jglobx+1
@ -134,7 +134,7 @@ subroutine psb_zgatherm(globx, locx, desc_a, info, iroot)
globx(:,:)=0.d0
do j=1,k
do i=1,psb_cd_get_local_rows(desc_a)
do i=1,desc_a%get_local_rows()
call psb_loc_to_glob(i,idx,desc_a,info)
globx(idx,jglobx+j-1) = locx(i,jlx+j-1)
end do
@ -241,7 +241,7 @@ subroutine psb_zgatherv(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -271,8 +271,8 @@ subroutine psb_zgatherv(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
k = 1
@ -297,7 +297,7 @@ subroutine psb_zgatherv(globx, locx, desc_a, info, iroot)
globx(:)=0.d0
do i=1,psb_cd_get_local_rows(desc_a)
do i=1,desc_a%get_local_rows()
call psb_loc_to_glob(i,idx,desc_a,info)
globx(idx) = locx(i)
end do

16
base/comm/psb_zhalo.f90
Unescape Escape View File

@ -80,7 +80,7 @@ subroutine psb_zhalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -97,9 +97,9 @@ subroutine psb_zhalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
maxk=size(x,2)-ijx+1
@ -301,7 +301,7 @@ subroutine psb_zhalov(x,desc_a,info,alpha,work,tran,mode,data)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -314,9 +314,9 @@ subroutine psb_zhalov(x,desc_a,info,alpha,work,tran,mode,data)
ix = 1
ijx = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
if (present(tran)) then
tran_ = psb_toupper(tran)

20
base/comm/psb_zovrl.f90
Unescape Escape View File

@ -89,7 +89,7 @@ subroutine psb_zovrlm(x,desc_a,info,jx,ik,work,update,mode)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -106,10 +106,10 @@ subroutine psb_zovrlm(x,desc_a,info,jx,ik,work,update,mode)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
maxk=size(x,2)-ijx+1
@ -288,7 +288,7 @@ subroutine psb_zovrlv(x,desc_a,info,work,update,mode)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -301,10 +301,10 @@ subroutine psb_zovrlv(x,desc_a,info,work,update,mode)
ix = 1
ijx = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
k = 1

20
base/comm/psb_zscatter.F90
Unescape Escape View File

@ -74,7 +74,7 @@ subroutine psb_zscatterm(globx, locx, desc_a, info, iroot)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)
@ -106,8 +106,8 @@ subroutine psb_zscatterm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx,1)
lda_locx = size(locx, 1)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
lock=size(locx,2)-jlocx+1
globk=size(globx,2)-jglobx+1
@ -120,8 +120,8 @@ subroutine psb_zscatterm(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
call psb_bcast(ictxt,k,root=iiroot)
@ -142,7 +142,7 @@ subroutine psb_zscatterm(globx, locx, desc_a, info, iroot)
goto 9999
end if
nrow=psb_cd_get_local_rows(desc_a)
nrow=desc_a%get_local_rows()
if ((root == -1).or.(np == 1)) then
! extract my chunk
@ -307,7 +307,7 @@ subroutine psb_zscatterv(globx, locx, desc_a, info, iroot)
if (psb_get_errstatus() /= 0) return
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
@ -342,8 +342,8 @@ subroutine psb_zscatterv(globx, locx, desc_a, info, iroot)
lda_globx = size(globx)
lda_locx = size(locx)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
k = 1
! there should be a global check on k here!!!
@ -364,7 +364,7 @@ subroutine psb_zscatterv(globx, locx, desc_a, info, iroot)
goto 9999
end if
nrow = psb_cd_get_local_rows(desc_a)
nrow = desc_a%get_local_rows()
if ((root == -1).or.(np == 1)) then
! extract my chunk

8
base/comm/psb_zspgather.F90
Unescape Escape View File

@ -30,8 +30,8 @@ subroutine psb_zsp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keep
info=psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt, me, np)
if (present(keepnum)) then
@ -47,8 +47,8 @@ subroutine psb_zsp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keep
call globa%free()
if (keepnum_) then
nrg = psb_cd_get_global_rows(desc_a)
ncg = psb_cd_get_global_rows(desc_a)
nrg = desc_a%get_global_rows()
ncg = desc_a%get_global_rows()
allocate(nzbr(np), idisp(np),stat=info)
if (info /= psb_success_) then

8
base/internals/psi_bld_tmphalo.f90
Unescape Escape View File

@ -68,10 +68,10 @@ subroutine psi_bld_tmphalo(desc,info)
name = 'psi_bld_tmphalo'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc)
icomm = psb_cd_get_mpic(desc)
n_row = psb_cd_get_local_rows(desc)
n_col = psb_cd_get_local_cols(desc)
ictxt = desc%get_context()
icomm = desc%get_mpic()
n_row = desc%get_local_rows()
n_col = desc%get_local_cols()
! check on blacs grid
call psb_info(ictxt, me, np)

4
base/internals/psi_bld_tmpovrl.f90
Unescape Escape View File

@ -76,8 +76,8 @@ subroutine psi_bld_tmpovrl(iv,desc,info)
call psb_erractionsave(err_act)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt = psb_cd_get_context(desc)
icomm = psb_cd_get_mpic(desc)
ictxt = desc%get_context()
icomm = desc%get_mpic()
! check on blacs grid
call psb_info(ictxt, me, np)

8
base/internals/psi_cswapdata.F90
Unescape Escape View File

@ -111,8 +111,8 @@ subroutine psi_cswapdatam(flag,n,beta,y,desc_a,work,info,data)
name='psi_swap_data'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_
@ -601,7 +601,7 @@ subroutine psi_cswapdatav(flag,beta,y,desc_a,work,info,data)
name='psi_swap_datav'
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_
@ -615,7 +615,7 @@ subroutine psi_cswapdatav(flag,beta,y,desc_a,work,info,data)
goto 9999
endif
icomm = psb_cd_get_mpic(desc_a)
icomm = desc_a%get_mpic()
if(present(data)) then
data_ = data

8
base/internals/psi_cswaptran.F90
Unescape Escape View File

@ -116,8 +116,8 @@ subroutine psi_cswaptranm(flag,n,beta,y,desc_a,work,info,data)
name='psi_swap_tran'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
@ -605,8 +605,8 @@ subroutine psi_cswaptranv(flag,beta,y,desc_a,work,info,data)
name='psi_swap_tranv'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_

2
base/internals/psi_desc_impl.f90
Unescape Escape View File

@ -88,7 +88,7 @@ subroutine psi_cnv_dsc(halo_in,ovrlap_in,ext_in,cdesc, info)
debug_unit = psb_get_debug_unit()
info = psb_success_
ictxt = psb_cd_get_context(cdesc)
ictxt = cdesc%get_context()
call psb_info(ictxt,me,np)
if (np == -1) then

8
base/internals/psi_dswapdata.F90
Unescape Escape View File

@ -112,8 +112,8 @@ subroutine psi_dswapdatam(flag,n,beta,y,desc_a,work,info,data)
name='psi_swap_data'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_
@ -601,7 +601,7 @@ subroutine psi_dswapdatav(flag,beta,y,desc_a,work,info,data)
name='psi_swap_datav'
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_
@ -615,7 +615,7 @@ subroutine psi_dswapdatav(flag,beta,y,desc_a,work,info,data)
goto 9999
endif
icomm = psb_cd_get_mpic(desc_a)
icomm = desc_a%get_mpic()
if(present(data)) then
data_ = data

8
base/internals/psi_dswaptran.F90
Unescape Escape View File

@ -116,8 +116,8 @@ subroutine psi_dswaptranm(flag,n,beta,y,desc_a,work,info,data)
name='psi_swap_tran'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
@ -605,8 +605,8 @@ subroutine psi_dswaptranv(flag,beta,y,desc_a,work,info,data)
name='psi_swap_tranv'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_

8
base/internals/psi_fnd_owner.F90
Unescape Escape View File

@ -82,10 +82,10 @@ subroutine psi_fnd_owner(nv,idx,iprc,desc,info)
name = 'psi_fnd_owner'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc)
icomm = psb_cd_get_mpic(desc)
n_row = psb_cd_get_local_rows(desc)
n_col = psb_cd_get_local_cols(desc)
ictxt = desc%get_context()
icomm = desc%get_mpic()
n_row = desc%get_local_rows()
n_col = desc%get_local_cols()
! check on blacs grid

20
base/internals/psi_idx_cnv.f90
Unescape Escape View File

@ -70,11 +70,11 @@ subroutine psi_idx_cnv1(nv,idxin,desc,info,mask,owned)
name = 'psb_idx_cnv'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc)
mglob = psb_cd_get_global_rows(desc)
nglob = psb_cd_get_global_cols(desc)
nrow = psb_cd_get_local_rows(desc)
ncol = psb_cd_get_local_cols(desc)
ictxt = desc%get_context()
mglob = desc%get_global_rows()
nglob = desc%get_global_cols()
nrow = desc%get_local_rows()
ncol = desc%get_local_cols()
call psb_info(ictxt, me, np)
@ -206,11 +206,11 @@ subroutine psi_idx_cnv2(nv,idxin,idxout,desc,info,mask,owned)
name = 'psb_idx_cnv'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc)
mglob = psb_cd_get_global_rows(desc)
nglob = psb_cd_get_global_cols(desc)
nrow = psb_cd_get_local_rows(desc)
ncol = psb_cd_get_local_cols(desc)
ictxt = desc%get_context()
mglob = desc%get_global_rows()
nglob = desc%get_global_cols()
nrow = desc%get_local_rows()
ncol = desc%get_local_cols()
call psb_info(ictxt, me, np)

20
base/internals/psi_idx_ins_cnv.f90
Unescape Escape View File

@ -73,11 +73,11 @@ subroutine psi_idx_ins_cnv1(nv,idxin,desc,info,mask)
name = 'psb_idx_ins_cnv'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc)
mglob = psb_cd_get_global_rows(desc)
nglob = psb_cd_get_global_cols(desc)
nrow = psb_cd_get_local_rows(desc)
ncol = psb_cd_get_local_cols(desc)
ictxt = desc%get_context()
mglob = desc%get_global_rows()
nglob = desc%get_global_cols()
nrow = desc%get_local_rows()
ncol = desc%get_local_cols()
call psb_info(ictxt, me, np)
@ -209,11 +209,11 @@ subroutine psi_idx_ins_cnv2(nv,idxin,idxout,desc,info,mask)
name = 'psb_idx_ins_cnv'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc)
mglob = psb_cd_get_global_rows(desc)
nglob = psb_cd_get_global_cols(desc)
nrow = psb_cd_get_local_rows(desc)
ncol = psb_cd_get_local_cols(desc)
ictxt = desc%get_context()
mglob = desc%get_global_rows()
nglob = desc%get_global_cols()
nrow = desc%get_local_rows()
ncol = desc%get_local_cols()
call psb_info(ictxt, me, np)

8
base/internals/psi_iswapdata.F90
Unescape Escape View File

@ -111,8 +111,8 @@ subroutine psi_iswapdatam(flag,n,beta,y,desc_a,work,info,data)
name='psi_swap_data'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_
@ -601,7 +601,7 @@ subroutine psi_iswapdatav(flag,beta,y,desc_a,work,info,data)
name='psi_swap_datav'
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_
@ -615,7 +615,7 @@ subroutine psi_iswapdatav(flag,beta,y,desc_a,work,info,data)
goto 9999
endif
icomm = psb_cd_get_mpic(desc_a)
icomm = desc_a%get_mpic()
if(present(data)) then
data_ = data

8
base/internals/psi_iswaptran.F90
Unescape Escape View File

@ -116,8 +116,8 @@ subroutine psi_iswaptranm(flag,n,beta,y,desc_a,work,info,data)
name='psi_swap_tran'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
@ -605,8 +605,8 @@ subroutine psi_iswaptranv(flag,beta,y,desc_a,work,info,data)
name='psi_swap_tranv'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_

8
base/internals/psi_ldsc_pre_halo.f90
Unescape Escape View File

@ -67,10 +67,10 @@ subroutine psi_ldsc_pre_halo(desc,ext_hv,info)
name = 'psi_ldsc_pre_halo'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc)
icomm = psb_cd_get_mpic(desc)
n_row = psb_cd_get_local_rows(desc)
n_col = psb_cd_get_local_cols(desc)
ictxt = desc%get_context()
icomm = desc%get_mpic()
n_row = desc%get_local_rows()
n_col = desc%get_local_cols()
! check on blacs grid
call psb_info(ictxt, me, np)

20
base/internals/psi_ovrl_restr.f90
Unescape Escape View File

@ -47,7 +47,7 @@ subroutine psi_sovrl_restrr1(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -93,7 +93,7 @@ subroutine psi_sovrl_restrr2(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -147,7 +147,7 @@ subroutine psi_dovrl_restrr1(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -194,7 +194,7 @@ subroutine psi_dovrl_restrr2(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -248,7 +248,7 @@ subroutine psi_covrl_restrr1(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -294,7 +294,7 @@ subroutine psi_covrl_restrr2(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -348,7 +348,7 @@ subroutine psi_zovrl_restrr1(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -394,7 +394,7 @@ subroutine psi_zovrl_restrr2(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -448,7 +448,7 @@ subroutine psi_iovrl_restrr1(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -494,7 +494,7 @@ subroutine psi_iovrl_restrr2(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_

20
base/internals/psi_ovrl_save.f90
Unescape Escape View File

@ -49,7 +49,7 @@ subroutine psi_sovrl_saver1(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -102,7 +102,7 @@ subroutine psi_sovrl_saver2(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -157,7 +157,7 @@ subroutine psi_dovrl_saver1(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -211,7 +211,7 @@ subroutine psi_dovrl_saver2(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -265,7 +265,7 @@ subroutine psi_covrl_saver1(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -319,7 +319,7 @@ subroutine psi_covrl_saver2(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -375,7 +375,7 @@ subroutine psi_zovrl_saver1(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -430,7 +430,7 @@ subroutine psi_zovrl_saver2(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -486,7 +486,7 @@ subroutine psi_iovrl_saver1(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -542,7 +542,7 @@ subroutine psi_iovrl_saver2(x,xs,desc_a,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_

20
base/internals/psi_ovrl_upd.f90
Unescape Escape View File

@ -48,7 +48,7 @@ subroutine psi_sovrl_updr1(x,desc_a,update,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -118,7 +118,7 @@ subroutine psi_sovrl_updr2(x,desc_a,update,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -187,7 +187,7 @@ subroutine psi_dovrl_updr1(x,desc_a,update,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -257,7 +257,7 @@ subroutine psi_dovrl_updr2(x,desc_a,update,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -326,7 +326,7 @@ subroutine psi_covrl_updr1(x,desc_a,update,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -396,7 +396,7 @@ subroutine psi_covrl_updr2(x,desc_a,update,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -465,7 +465,7 @@ subroutine psi_zovrl_updr1(x,desc_a,update,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -535,7 +535,7 @@ subroutine psi_zovrl_updr2(x,desc_a,update,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -604,7 +604,7 @@ subroutine psi_iovrl_updr1(x,desc_a,update,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -675,7 +675,7 @@ subroutine psi_iovrl_updr2(x,desc_a,update,info)
if (psb_get_errstatus() /= 0) return
info = psb_success_
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_

8
base/internals/psi_sswapdata.F90
Unescape Escape View File

@ -112,8 +112,8 @@ subroutine psi_sswapdatam(flag,n,beta,y,desc_a,work,info,data)
name='psi_swap_data'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_
@ -601,7 +601,7 @@ subroutine psi_sswapdatav(flag,beta,y,desc_a,work,info,data)
name='psi_swap_datav'
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_
@ -615,7 +615,7 @@ subroutine psi_sswapdatav(flag,beta,y,desc_a,work,info,data)
goto 9999
endif
icomm = psb_cd_get_mpic(desc_a)
icomm = desc_a%get_mpic()
if(present(data)) then
data_ = data

8
base/internals/psi_sswaptran.F90
Unescape Escape View File

@ -116,8 +116,8 @@ subroutine psi_sswaptranm(flag,n,beta,y,desc_a,work,info,data)
name='psi_swap_tran'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
@ -605,8 +605,8 @@ subroutine psi_sswaptranv(flag,beta,y,desc_a,work,info,data)
name='psi_swap_tranv'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_

8
base/internals/psi_zswapdata.F90
Unescape Escape View File

@ -111,8 +111,8 @@ subroutine psi_zswapdatam(flag,n,beta,y,desc_a,work,info,data)
name='psi_swap_data'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_
@ -601,7 +601,7 @@ subroutine psi_zswapdatav(flag,beta,y,desc_a,work,info,data)
name='psi_swap_datav'
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_
@ -615,7 +615,7 @@ subroutine psi_zswapdatav(flag,beta,y,desc_a,work,info,data)
goto 9999
endif
icomm = psb_cd_get_mpic(desc_a)
icomm = desc_a%get_mpic()
if(present(data)) then
data_ = data

8
base/internals/psi_zswaptran.F90
Unescape Escape View File

@ -116,8 +116,8 @@ subroutine psi_zswaptranm(flag,n,beta,y,desc_a,work,info,data)
name='psi_swap_tran'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
@ -605,8 +605,8 @@ subroutine psi_zswaptranv(flag,beta,y,desc_a,work,info,data)
name='psi_swap_tranv'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
call psb_info(ictxt,me,np)
if (np == -1) then
info=psb_err_context_error_

10
base/modules/psb_comm_mod.f90
Unescape Escape View File

@ -42,7 +42,7 @@ module psb_comm_mod
end subroutine psb_sovrlm
subroutine psb_sovrlv(x,desc_a,info,work,update,mode)
use psb_descriptor_type
real(psb_spk_), intent(inout) :: x(:)
real(psb_spk_), intent(inout), target :: x(:)
type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info
real(psb_spk_), intent(inout), optional, target :: work(:)
@ -58,7 +58,7 @@ module psb_comm_mod
end subroutine psb_dovrlm
subroutine psb_dovrlv(x,desc_a,info,work,update,mode)
use psb_descriptor_type
real(psb_dpk_), intent(inout) :: x(:)
real(psb_dpk_), intent(inout), target :: x(:)
type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info
real(psb_dpk_), intent(inout), optional, target :: work(:)
@ -74,7 +74,7 @@ module psb_comm_mod
end subroutine psb_iovrlm
subroutine psb_iovrlv(x,desc_a,info,work,update,mode)
use psb_descriptor_type
integer, intent(inout) :: x(:)
integer, intent(inout), target :: x(:)
type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info
integer, intent(inout), optional, target :: work(:)
@ -90,7 +90,7 @@ module psb_comm_mod
end subroutine psb_covrlm
subroutine psb_covrlv(x,desc_a,info,work,update,mode)
use psb_descriptor_type
complex(psb_spk_), intent(inout) :: x(:)
complex(psb_spk_), intent(inout), target :: x(:)
type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info
complex(psb_spk_), intent(inout), optional, target :: work(:)
@ -106,7 +106,7 @@ module psb_comm_mod
end subroutine psb_zovrlm
subroutine psb_zovrlv(x,desc_a,info,work,update,mode)
use psb_descriptor_type
complex(psb_dpk_), intent(inout) :: x(:)
complex(psb_dpk_), intent(inout), target :: x(:)
type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info
complex(psb_dpk_), intent(inout), optional, target :: work(:)

6
base/modules/psb_desc_type.f90
Unescape Escape View File

@ -217,6 +217,8 @@ module psb_descriptor_type
procedure, pass(desc) :: is_asb => psb_is_asb_desc
procedure, pass(desc) :: is_ovl => psb_is_ovl_desc
procedure, pass(desc) :: is_repl => psb_is_repl_desc
procedure, pass(desc) :: get_mpic => psb_cd_get_mpic
procedure, pass(desc) :: get_dectype => psb_cd_get_dectype
procedure, pass(desc) :: get_context => psb_cd_get_context
procedure, pass(desc) :: get_local_rows => psb_cd_get_local_rows
procedure, pass(desc) :: get_local_cols => psb_cd_get_local_cols
@ -456,7 +458,7 @@ contains
integer function psb_cd_get_dectype(desc)
use psb_error_mod
type(psb_desc_type), intent(in) :: desc
class(psb_desc_type), intent(in) :: desc
if (allocated(desc%indxmap)) then
psb_cd_get_dectype = desc%indxmap%get_state()
@ -470,7 +472,7 @@ contains
integer function psb_cd_get_mpic(desc)
use psb_error_mod
type(psb_desc_type), intent(in) :: desc
class(psb_desc_type), intent(in) :: desc
if (allocated(desc%indxmap)) then
psb_cd_get_mpic = desc%indxmap%get_mpic()

32
base/psblas/psb_camax.f90
Unescape Escape View File

@ -71,7 +71,7 @@ function psb_camax(x,desc_a, info, jx)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -87,7 +87,7 @@ function psb_camax(x,desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -104,8 +104,8 @@ function psb_camax(x,desc_a, info, jx)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=icamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix,jjx),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=icamax(desc_a%get_local_rows()-iix+1,x(iix,jjx),1)
amax=cabs1(x(iix+imax-1,jjx))
else
amax = szero
@ -201,7 +201,7 @@ function psb_camaxv (x,desc_a, info)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -213,7 +213,7 @@ function psb_camaxv (x,desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -230,8 +230,8 @@ function psb_camaxv (x,desc_a, info)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=icamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=icamax(desc_a%get_local_rows()-iix+1,x(iix),1)
cmax=(x(iix+imax-1))
amax=cabs1(cmax)
else
@ -329,7 +329,7 @@ subroutine psb_camaxvs(res,x,desc_a, info)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -341,7 +341,7 @@ subroutine psb_camaxvs(res,x,desc_a, info)
ix = 1
ijx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
@ -357,8 +357,8 @@ subroutine psb_camaxvs(res,x,desc_a, info)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=icamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=icamax(desc_a%get_local_rows()-iix+1,x(iix),1)
cmax=(x(iix+imax-1))
amax=cabs1(cmax)
else
@ -455,7 +455,7 @@ subroutine psb_cmamaxs(res,x,desc_a, info,jx)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -471,7 +471,7 @@ subroutine psb_cmamaxs(res,x,desc_a, info,jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
k = min(size(x,2),size(res,1))
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -489,9 +489,9 @@ subroutine psb_cmamaxs(res,x,desc_a, info,jx)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
do i=1,k
imax=icamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix,jjx+i-1),1)
imax=icamax(desc_a%get_local_rows()-iix+1,x(iix,jjx+i-1),1)
cmax=(x(iix+imax-1,jjx+i-1))
res(i)=cabs1(cmax)
end do

24
base/psblas/psb_casum.f90
Unescape Escape View File

@ -72,7 +72,7 @@ function psb_casum (x,desc_a, info, jx)
asum=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -88,7 +88,7 @@ function psb_casum (x,desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -107,8 +107,8 @@ function psb_casum (x,desc_a, info, jx)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=scasum(psb_cd_get_local_rows(desc_a)-iix+1,x(iix:,jjx),ione)
if(desc_a%get_local_rows() > 0) then
asum=scasum(desc_a%get_local_rows()-iix+1,x(iix:,jjx),ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -213,7 +213,7 @@ function psb_casumv(x,desc_a, info)
asum=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -225,7 +225,7 @@ function psb_casumv(x,desc_a, info)
ix = 1
jx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
@ -244,8 +244,8 @@ function psb_casumv(x,desc_a, info)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=scasum(psb_cd_get_local_rows(desc_a),x,ione)
if(desc_a%get_local_rows() > 0) then
asum=scasum(desc_a%get_local_rows(),x,ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -351,7 +351,7 @@ subroutine psb_casumvs(res,x,desc_a, info)
asum=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -363,7 +363,7 @@ subroutine psb_casumvs(res,x,desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
@ -382,8 +382,8 @@ subroutine psb_casumvs(res,x,desc_a, info)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=scasum(psb_cd_get_local_rows(desc_a),x,ione)
if(desc_a%get_local_rows() > 0) then
asum=scasum(desc_a%get_local_rows(),x,ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)

16
base/psblas/psb_caxpby.f90
Unescape Escape View File

@ -72,7 +72,7 @@ subroutine psb_caxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
info = psb_err_context_error_
@ -111,7 +111,7 @@ subroutine psb_caxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -131,8 +131,8 @@ subroutine psb_caxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
end if
if ((in /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
call caxpby(psb_cd_get_local_cols(desc_a),in,&
if(desc_a%get_local_rows() > 0) then
call caxpby(desc_a%get_local_cols(),in,&
& alpha,x(iix:,jjx),size(x,1),beta,&
& y(iiy:,jjy),size(y,1),info)
end if
@ -222,7 +222,7 @@ subroutine psb_caxpbyv(alpha, x, beta,y,desc_a,info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -234,7 +234,7 @@ subroutine psb_caxpbyv(alpha, x, beta,y,desc_a,info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x),ix,ione,desc_a,info,iix,jjx)
@ -257,8 +257,8 @@ subroutine psb_caxpbyv(alpha, x, beta,y,desc_a,info)
call psb_errpush(info,name)
end if
if(psb_cd_get_local_rows(desc_a) > 0) then
call caxpby(psb_cd_get_local_cols(desc_a),ione,&
if(desc_a%get_local_rows() > 0) then
call caxpby(desc_a%get_local_cols(),ione,&
& alpha,x,size(x),beta,&
& y,size(y),info)
end if

24
base/psblas/psb_cdot.f90
Unescape Escape View File

@ -70,7 +70,7 @@ function psb_cdot(x, y,desc_a, info, jx, jy)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
info = psb_err_context_error_
@ -98,7 +98,7 @@ function psb_cdot(x, y,desc_a, info, jx, jy)
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -118,7 +118,7 @@ function psb_cdot(x, y,desc_a, info, jx, jy)
end if
if(m /= 0) then
nr = psb_cd_get_local_rows(desc_a)
nr = desc_a%get_local_rows()
if(nr > 0) then
dot_local = cdotc(nr, x(iix:,jjx),ione,y(iiy:,jjy),ione)
! adjust dot_local because overlapped elements are computed more than once
@ -219,7 +219,7 @@ function psb_cdotv(x, y,desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -232,7 +232,7 @@ function psb_cdotv(x, y,desc_a, info)
iy = ione
jx = ione
jy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,jx,desc_a,info,iix,jjx)
@ -252,7 +252,7 @@ function psb_cdotv(x, y,desc_a, info)
end if
if(m /= 0) then
nr = psb_cd_get_local_rows(desc_a)
nr = desc_a%get_local_rows()
if(nr > 0) then
dot_local = cdotc(nr, x,ione,y,ione)
! adjust dot_local because overlapped elements are computed more than once
@ -353,7 +353,7 @@ subroutine psb_cdotvs(res, x, y,desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -364,7 +364,7 @@ subroutine psb_cdotvs(res, x, y,desc_a, info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ix,desc_a,info,iix,jjx)
if (info == psb_success_) &
@ -383,7 +383,7 @@ subroutine psb_cdotvs(res, x, y,desc_a, info)
end if
if(m /= 0) then
nr = psb_cd_get_local_rows(desc_a)
nr = desc_a%get_local_rows()
if(nr > 0) then
dot_local = cdotc(nr, x,ione,y,ione)
! adjust dot_local because overlapped elements are computed more than once
@ -485,7 +485,7 @@ subroutine psb_cmdots(res, x, y, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -497,7 +497,7 @@ subroutine psb_cmdots(res, x, y, desc_a, info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ix,desc_a,info,iix,jjx)
@ -525,7 +525,7 @@ subroutine psb_cmdots(res, x, y, desc_a, info)
allocate(dot_local(k))
if(m /= 0) then
nr = psb_cd_get_local_rows(desc_a)
nr = desc_a%get_local_rows()
if(nr > 0) then
do j=1,k
dot_local(j) = cdotc(nr, x(1:,j),ione,y(1:,j),ione)

24
base/psblas/psb_cnrm2.f90
Unescape Escape View File

@ -70,7 +70,7 @@ function psb_cnrm2(x, desc_a, info, jx)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -86,7 +86,7 @@ function psb_cnrm2(x, desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
@ -101,8 +101,8 @@ function psb_cnrm2(x, desc_a, info, jx)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = scnrm2( ndim, x(iix:,jjx), ione )
! adjust because overlapped elements are computed more than once
@ -206,7 +206,7 @@ function psb_cnrm2v(x, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -217,7 +217,7 @@ function psb_cnrm2v(x, desc_a, info)
ix = 1
jx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -233,8 +233,8 @@ function psb_cnrm2v(x, desc_a, info)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = scnrm2( ndim, x, ione )
! adjust because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -339,7 +339,7 @@ subroutine psb_cnrm2vs(res, x, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -350,7 +350,7 @@ subroutine psb_cnrm2vs(res, x, desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -366,8 +366,8 @@ subroutine psb_cnrm2vs(res, x, desc_a, info)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = scnrm2( ndim, x, ione )
! adjust because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)

6
base/psblas/psb_cnrmi.f90
Unescape Escape View File

@ -61,7 +61,7 @@ function psb_cnrmi(a,desc_a,info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -72,8 +72,8 @@ function psb_cnrmi(a,desc_a,info)
ia = 1
ja = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
call psb_chkmat(m,n,ia,ja,desc_a,info,iia,jja)
if(info /= psb_success_) then

20
base/psblas/psb_cspmm.f90
Unescape Escape View File

@ -99,7 +99,7 @@ subroutine psb_cspmm(alpha,a,x,beta,y,desc_a,info,&
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -151,10 +151,10 @@ subroutine psb_cspmm(alpha,a,x,beta,y,desc_a,info,&
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x,1)
lldy = size(y,1)
@ -450,7 +450,7 @@ subroutine psb_cspmv(alpha,a,x,beta,y,desc_a,info,&
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -486,10 +486,10 @@ subroutine psb_cspmv(alpha,a,x,beta,y,desc_a,info,&
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x)
lldy = size(y)

16
base/psblas/psb_cspsm.f90
Unescape Escape View File

@ -109,7 +109,7 @@ subroutine psb_cspsm(alpha,a,x,beta,y,desc_a,info,&
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -168,9 +168,9 @@ subroutine psb_cspsm(alpha,a,x,beta,y,desc_a,info,&
itrans = 'N'
endif
m = psb_cd_get_global_rows(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x,1)
lldy = size(y,1)
@ -387,7 +387,7 @@ subroutine psb_cspsv(alpha,a,x,beta,y,desc_a,info,&
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -430,9 +430,9 @@ subroutine psb_cspsv(alpha,a,x,beta,y,desc_a,info,&
itrans = 'N'
endif
m = psb_cd_get_global_rows(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x)
lldy = size(y)

32
base/psblas/psb_damax.f90
Unescape Escape View File

@ -71,7 +71,7 @@ function psb_damax (x,desc_a, info, jx)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -87,7 +87,7 @@ function psb_damax (x,desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -104,8 +104,8 @@ function psb_damax (x,desc_a, info, jx)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=idamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix,jjx),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=idamax(desc_a%get_local_rows()-iix+1,x(iix,jjx),1)
amax=abs(x(iix+imax-1,jjx))
else
amax = dzero
@ -200,7 +200,7 @@ function psb_damaxv (x,desc_a, info)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -212,7 +212,7 @@ function psb_damaxv (x,desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -229,8 +229,8 @@ function psb_damaxv (x,desc_a, info)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=idamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=idamax(desc_a%get_local_rows()-iix+1,x(iix),1)
amax=abs(x(iix+imax-1))
else
amax = dzero
@ -322,7 +322,7 @@ subroutine psb_damaxvs (res,x,desc_a, info)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -334,7 +334,7 @@ subroutine psb_damaxvs (res,x,desc_a, info)
ix = 1
ijx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -351,8 +351,8 @@ subroutine psb_damaxvs (res,x,desc_a, info)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=idamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=idamax(desc_a%get_local_rows()-iix+1,x(iix),1)
amax=abs(x(iix+imax-1))
else
amax = dzero
@ -443,7 +443,7 @@ subroutine psb_dmamaxs (res,x,desc_a, info,jx)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -459,7 +459,7 @@ subroutine psb_dmamaxs (res,x,desc_a, info,jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
k = min(size(x,2),size(res,1))
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -477,9 +477,9 @@ subroutine psb_dmamaxs (res,x,desc_a, info,jx)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
do i=1,k
imax=idamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix,jjx+i-1),1)
imax=idamax(desc_a%get_local_rows()-iix+1,x(iix,jjx+i-1),1)
res(i)=abs(x(iix+imax-1,jjx+i-1))
end do
else

24
base/psblas/psb_dasum.f90
Unescape Escape View File

@ -72,7 +72,7 @@ function psb_dasum (x,desc_a, info, jx)
asum=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -88,7 +88,7 @@ function psb_dasum (x,desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -107,8 +107,8 @@ function psb_dasum (x,desc_a, info, jx)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=dasum(psb_cd_get_local_rows(desc_a)-iix+1,x(iix:,jjx),ione)
if(desc_a%get_local_rows() > 0) then
asum=dasum(desc_a%get_local_rows()-iix+1,x(iix:,jjx),ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -214,7 +214,7 @@ function psb_dasumv (x,desc_a, info)
asum=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -226,7 +226,7 @@ function psb_dasumv (x,desc_a, info)
ix = 1
jx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
@ -245,8 +245,8 @@ function psb_dasumv (x,desc_a, info)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=dasum(psb_cd_get_local_rows(desc_a),x,ione)
if(desc_a%get_local_rows() > 0) then
asum=dasum(desc_a%get_local_rows(),x,ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
idx = desc_a%ovrlap_elem(i,1)
@ -347,7 +347,7 @@ subroutine psb_dasumvs(res,x,desc_a, info)
asum=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -359,7 +359,7 @@ subroutine psb_dasumvs(res,x,desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
@ -378,8 +378,8 @@ subroutine psb_dasumvs(res,x,desc_a, info)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=dasum(psb_cd_get_local_rows(desc_a),x,ione)
if(desc_a%get_local_rows() > 0) then
asum=dasum(desc_a%get_local_rows(),x,ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)

16
base/psblas/psb_daxpby.f90
Unescape Escape View File

@ -71,7 +71,7 @@ subroutine psb_daxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -111,7 +111,7 @@ subroutine psb_daxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -131,8 +131,8 @@ subroutine psb_daxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
end if
if ((in /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
call daxpby(psb_cd_get_local_rows(desc_a),in,&
if(desc_a%get_local_rows() > 0) then
call daxpby(desc_a%get_local_rows(),in,&
& alpha,x(iix:,jjx),size(x,1),beta,&
& y(iiy:,jjy),size(y,1),info)
end if
@ -221,7 +221,7 @@ subroutine psb_daxpbyv(alpha, x, beta,y,desc_a,info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -233,7 +233,7 @@ subroutine psb_daxpbyv(alpha, x, beta,y,desc_a,info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x),ix,ione,desc_a,info,iix,jjx)
@ -256,8 +256,8 @@ subroutine psb_daxpbyv(alpha, x, beta,y,desc_a,info)
call psb_errpush(info,name)
end if
if(psb_cd_get_local_rows(desc_a) > 0) then
call daxpby(psb_cd_get_local_rows(desc_a),ione,&
if(desc_a%get_local_rows() > 0) then
call daxpby(desc_a%get_local_rows(),ione,&
& alpha,x,size(x),beta,&
& y,size(y),info)
end if

24
base/psblas/psb_ddot.f90
Unescape Escape View File

@ -73,7 +73,7 @@ function psb_ddot(x, y,desc_a, info, jx, jy)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
info = psb_err_context_error_
@ -101,7 +101,7 @@ function psb_ddot(x, y,desc_a, info, jx, jy)
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -121,7 +121,7 @@ function psb_ddot(x, y,desc_a, info, jx, jy)
end if
if(m /= 0) then
nr = psb_cd_get_local_rows(desc_a)
nr = desc_a%get_local_rows()
if(nr > 0) then
dot_local = ddot(nr, x(iix:,jjx),ione,y(iiy:,jjy),ione)
! adjust dot_local because overlapped elements are computed more than once
@ -225,7 +225,7 @@ function psb_ddotv(x, y,desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -238,7 +238,7 @@ function psb_ddotv(x, y,desc_a, info)
iy = ione
jx = ione
jy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,jx,desc_a,info,iix,jjx)
@ -258,7 +258,7 @@ function psb_ddotv(x, y,desc_a, info)
end if
if(m /= 0) then
nr = psb_cd_get_local_rows(desc_a)
nr = desc_a%get_local_rows()
if(nr > 0) then
dot_local = ddot(nr, x,ione,y,ione)
@ -360,7 +360,7 @@ subroutine psb_ddotvs(res, x, y,desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -371,7 +371,7 @@ subroutine psb_ddotvs(res, x, y,desc_a, info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ix,desc_a,info,iix,jjx)
if (info == psb_success_) &
@ -390,7 +390,7 @@ subroutine psb_ddotvs(res, x, y,desc_a, info)
end if
if(m /= 0) then
nr = psb_cd_get_local_rows(desc_a)
nr = desc_a%get_local_rows()
if(nr > 0) then
dot_local = ddot(nr, x,ione,y,ione)
! adjust dot_local because overlapped elements are computed more than once
@ -492,7 +492,7 @@ subroutine psb_dmdots(res, x, y, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -504,7 +504,7 @@ subroutine psb_dmdots(res, x, y, desc_a, info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ix,desc_a,info,iix,jjx)
@ -532,7 +532,7 @@ subroutine psb_dmdots(res, x, y, desc_a, info)
allocate(dot_local(k))
if(m /= 0) then
nr = psb_cd_get_local_rows(desc_a)
nr = desc_a%get_local_rows()
if(nr > 0) then
do j=1,k
dot_local(j) = ddot(nr,x(1:,j),ione,y(1:,j),ione)

24
base/psblas/psb_dnrm2.f90
Unescape Escape View File

@ -69,7 +69,7 @@ function psb_dnrm2(x, desc_a, info, jx)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -85,7 +85,7 @@ function psb_dnrm2(x, desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
@ -100,8 +100,8 @@ function psb_dnrm2(x, desc_a, info, jx)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = dnrm2( ndim, x(iix:,jjx), ione )
! adjust because overlapped elements are computed more than once
@ -204,7 +204,7 @@ function psb_dnrm2v(x, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -215,7 +215,7 @@ function psb_dnrm2v(x, desc_a, info)
ix = 1
jx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -231,8 +231,8 @@ function psb_dnrm2v(x, desc_a, info)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = dnrm2( ndim, x, ione )
! adjust because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -336,7 +336,7 @@ subroutine psb_dnrm2vs(res, x, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -347,7 +347,7 @@ subroutine psb_dnrm2vs(res, x, desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -363,8 +363,8 @@ subroutine psb_dnrm2vs(res, x, desc_a, info)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = dnrm2( ndim, x, ione )
! adjust because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)

6
base/psblas/psb_dnrmi.f90
Unescape Escape View File

@ -66,7 +66,7 @@ function psb_dnrmi(a,desc_a,info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -77,8 +77,8 @@ function psb_dnrmi(a,desc_a,info)
ia = 1
ja = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
call psb_chkmat(m,n,ia,ja,desc_a,info,iia,jja)
if(info /= psb_success_) then

20
base/psblas/psb_dspmm.f90
Unescape Escape View File

@ -99,7 +99,7 @@ subroutine psb_dspmm(alpha,a,x,beta,y,desc_a,info,&
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -151,10 +151,10 @@ subroutine psb_dspmm(alpha,a,x,beta,y,desc_a,info,&
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x,1)
lldy = size(y,1)
@ -450,7 +450,7 @@ subroutine psb_dspmv(alpha,a,x,beta,y,desc_a,info,&
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -486,10 +486,10 @@ subroutine psb_dspmv(alpha,a,x,beta,y,desc_a,info,&
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x)
lldy = size(y)

10
base/psblas/psb_dspnrm1.f90
Unescape Escape View File

@ -69,7 +69,7 @@ function psb_dspnrm1(a,desc_a,info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -80,10 +80,10 @@ function psb_dspnrm1(a,desc_a,info)
ia = 1
ja = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nr = psb_cd_get_local_rows(desc_a)
nc = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nr = desc_a%get_local_rows()
nc = desc_a%get_local_cols()
call psb_chkmat(m,n,ia,ja,desc_a,info,iia,jja)
if(info /= psb_success_) then

16
base/psblas/psb_dspsm.f90
Unescape Escape View File

@ -110,7 +110,7 @@ subroutine psb_dspsm(alpha,a,x,beta,y,desc_a,info,&
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -169,9 +169,9 @@ subroutine psb_dspsm(alpha,a,x,beta,y,desc_a,info,&
itrans = 'N'
endif
m = psb_cd_get_global_rows(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x,1)
lldy = size(y,1)
@ -388,7 +388,7 @@ subroutine psb_dspsv(alpha,a,x,beta,y,desc_a,info,&
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -431,9 +431,9 @@ subroutine psb_dspsv(alpha,a,x,beta,y,desc_a,info,&
itrans = 'N'
endif
m = psb_cd_get_global_rows(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x)
lldy = size(y)

32
base/psblas/psb_samax.f90
Unescape Escape View File

@ -71,7 +71,7 @@ function psb_samax (x,desc_a, info, jx)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -87,7 +87,7 @@ function psb_samax (x,desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -104,8 +104,8 @@ function psb_samax (x,desc_a, info, jx)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=isamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix,jjx),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=isamax(desc_a%get_local_rows()-iix+1,x(iix,jjx),1)
amax=abs(x(iix+imax-1,jjx))
else
amax = szero
@ -200,7 +200,7 @@ function psb_samaxv (x,desc_a, info)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -212,7 +212,7 @@ function psb_samaxv (x,desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -229,8 +229,8 @@ function psb_samaxv (x,desc_a, info)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=isamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=isamax(desc_a%get_local_rows()-iix+1,x(iix),1)
amax=abs(x(iix+imax-1))
else
amax = szero
@ -322,7 +322,7 @@ subroutine psb_samaxvs (res,x,desc_a, info)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -334,7 +334,7 @@ subroutine psb_samaxvs (res,x,desc_a, info)
ix = 1
ijx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -351,8 +351,8 @@ subroutine psb_samaxvs (res,x,desc_a, info)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=isamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=isamax(desc_a%get_local_rows()-iix+1,x(iix),1)
amax=abs(x(iix+imax-1))
else
amax = szero
@ -443,7 +443,7 @@ subroutine psb_smamaxs (res,x,desc_a, info,jx)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -459,7 +459,7 @@ subroutine psb_smamaxs (res,x,desc_a, info,jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
k = min(size(x,2),size(res,1))
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -477,9 +477,9 @@ subroutine psb_smamaxs (res,x,desc_a, info,jx)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
do i=1,k
imax=isamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix,jjx+i-1),1)
imax=isamax(desc_a%get_local_rows()-iix+1,x(iix,jjx+i-1),1)
res(i)=abs(x(iix+imax-1,jjx+i-1))
end do
else

24
base/psblas/psb_sasum.f90
Unescape Escape View File

@ -72,7 +72,7 @@ function psb_sasum (x,desc_a, info, jx)
asum=0.0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -88,7 +88,7 @@ function psb_sasum (x,desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -107,8 +107,8 @@ function psb_sasum (x,desc_a, info, jx)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=sasum(psb_cd_get_local_rows(desc_a)-iix+1,x(iix:,jjx),ione)
if(desc_a%get_local_rows() > 0) then
asum=sasum(desc_a%get_local_rows()-iix+1,x(iix:,jjx),ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -214,7 +214,7 @@ function psb_sasumv (x,desc_a, info)
asum=0.0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -226,7 +226,7 @@ function psb_sasumv (x,desc_a, info)
ix = 1
jx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
@ -245,8 +245,8 @@ function psb_sasumv (x,desc_a, info)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=sasum(psb_cd_get_local_rows(desc_a),x,ione)
if(desc_a%get_local_rows() > 0) then
asum=sasum(desc_a%get_local_rows(),x,ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
idx = desc_a%ovrlap_elem(i,1)
@ -347,7 +347,7 @@ subroutine psb_sasumvs(res,x,desc_a, info)
asum=0.0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -359,7 +359,7 @@ subroutine psb_sasumvs(res,x,desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
@ -378,8 +378,8 @@ subroutine psb_sasumvs(res,x,desc_a, info)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=sasum(psb_cd_get_local_rows(desc_a),x,ione)
if(desc_a%get_local_rows() > 0) then
asum=sasum(desc_a%get_local_rows(),x,ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)

16
base/psblas/psb_saxpby.f90
Unescape Escape View File

@ -71,7 +71,7 @@ subroutine psb_saxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -111,7 +111,7 @@ subroutine psb_saxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -131,8 +131,8 @@ subroutine psb_saxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
end if
if ((in /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
call saxpby(psb_cd_get_local_rows(desc_a),in,&
if(desc_a%get_local_rows() > 0) then
call saxpby(desc_a%get_local_rows(),in,&
& alpha,x(iix:,jjx),size(x,1),beta,&
& y(iiy:,jjy),size(y,1),info)
end if
@ -221,7 +221,7 @@ subroutine psb_saxpbyv(alpha, x, beta,y,desc_a,info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -233,7 +233,7 @@ subroutine psb_saxpbyv(alpha, x, beta,y,desc_a,info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x),ix,ione,desc_a,info,iix,jjx)
@ -256,8 +256,8 @@ subroutine psb_saxpbyv(alpha, x, beta,y,desc_a,info)
call psb_errpush(info,name)
end if
if(psb_cd_get_local_rows(desc_a) > 0) then
call saxpby(psb_cd_get_local_rows(desc_a),ione,&
if(desc_a%get_local_rows() > 0) then
call saxpby(desc_a%get_local_rows(),ione,&
& alpha,x,size(x),beta,&
& y,size(y),info)
end if

32
base/psblas/psb_sdot.f90
Unescape Escape View File

@ -73,7 +73,7 @@ function psb_sdot(x, y,desc_a, info, jx, jy)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
info = psb_err_context_error_
@ -101,7 +101,7 @@ function psb_sdot(x, y,desc_a, info, jx, jy)
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -121,8 +121,8 @@ function psb_sdot(x, y,desc_a, info, jx, jy)
end if
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
dot_local = sdot(psb_cd_get_local_rows(desc_a),&
if(desc_a%get_local_rows() > 0) then
dot_local = sdot(desc_a%get_local_rows(),&
& x(iix:,jjx),ione,y(iiy:,jjy),ione)
! adjust dot_local because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -225,7 +225,7 @@ function psb_sdotv(x, y,desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -238,7 +238,7 @@ function psb_sdotv(x, y,desc_a, info)
iy = ione
jx = ione
jy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,jx,desc_a,info,iix,jjx)
@ -258,8 +258,8 @@ function psb_sdotv(x, y,desc_a, info)
end if
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
dot_local = sdot(psb_cd_get_local_rows(desc_a),&
if(desc_a%get_local_rows() > 0) then
dot_local = sdot(desc_a%get_local_rows(),&
& x,ione,y,ione)
! adjust dot_local because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -359,7 +359,7 @@ subroutine psb_sdotvs(res, x, y,desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -370,7 +370,7 @@ subroutine psb_sdotvs(res, x, y,desc_a, info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ix,desc_a,info,iix,jjx)
if (info == psb_success_) &
@ -389,8 +389,8 @@ subroutine psb_sdotvs(res, x, y,desc_a, info)
end if
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
dot_local = sdot(psb_cd_get_local_rows(desc_a),&
if(desc_a%get_local_rows() > 0) then
dot_local = sdot(desc_a%get_local_rows(),&
& x,ione,y,ione)
! adjust dot_local because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -491,7 +491,7 @@ subroutine psb_smdots(res, x, y, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -503,7 +503,7 @@ subroutine psb_smdots(res, x, y, desc_a, info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ix,desc_a,info,iix,jjx)
@ -531,9 +531,9 @@ subroutine psb_smdots(res, x, y, desc_a, info)
allocate(dot_local(k))
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
if(desc_a%get_local_rows() > 0) then
do j=1,k
dot_local(j) = sdot(psb_cd_get_local_rows(desc_a),&
dot_local(j) = sdot(desc_a%get_local_rows(),&
& x(1:,j),ione,y(1:,j),ione)
! adjust dot_local because overlapped elements are computed more than once
end do

24
base/psblas/psb_snrm2.f90
Unescape Escape View File

@ -69,7 +69,7 @@ function psb_snrm2(x, desc_a, info, jx)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -85,7 +85,7 @@ function psb_snrm2(x, desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
@ -100,8 +100,8 @@ function psb_snrm2(x, desc_a, info, jx)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = snrm2( ndim, x(iix:,jjx), ione )
! adjust because overlapped elements are computed more than once
@ -204,7 +204,7 @@ function psb_snrm2v(x, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -215,7 +215,7 @@ function psb_snrm2v(x, desc_a, info)
ix = 1
jx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -231,8 +231,8 @@ function psb_snrm2v(x, desc_a, info)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = snrm2( ndim, x, ione )
! adjust because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -336,7 +336,7 @@ subroutine psb_snrm2vs(res, x, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -347,7 +347,7 @@ subroutine psb_snrm2vs(res, x, desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -363,8 +363,8 @@ subroutine psb_snrm2vs(res, x, desc_a, info)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = snrm2( ndim, x, ione )
! adjust because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)

6
base/psblas/psb_snrmi.f90
Unescape Escape View File

@ -66,7 +66,7 @@ function psb_snrmi(a,desc_a,info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -77,8 +77,8 @@ function psb_snrmi(a,desc_a,info)
ia = 1
ja = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
call psb_chkmat(m,n,ia,ja,desc_a,info,iia,jja)
if(info /= psb_success_) then

20
base/psblas/psb_sspmm.f90
Unescape Escape View File

@ -99,7 +99,7 @@ subroutine psb_sspmm(alpha,a,x,beta,y,desc_a,info,&
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -151,10 +151,10 @@ subroutine psb_sspmm(alpha,a,x,beta,y,desc_a,info,&
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x,1)
lldy = size(y,1)
@ -450,7 +450,7 @@ subroutine psb_sspmv(alpha,a,x,beta,y,desc_a,info,&
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -486,10 +486,10 @@ subroutine psb_sspmv(alpha,a,x,beta,y,desc_a,info,&
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x)
lldy = size(y)

16
base/psblas/psb_sspsm.f90
Unescape Escape View File

@ -110,7 +110,7 @@ subroutine psb_sspsm(alpha,a,x,beta,y,desc_a,info,&
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -169,9 +169,9 @@ subroutine psb_sspsm(alpha,a,x,beta,y,desc_a,info,&
itrans = 'N'
endif
m = psb_cd_get_global_rows(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x,1)
lldy = size(y,1)
@ -388,7 +388,7 @@ subroutine psb_sspsv(alpha,a,x,beta,y,desc_a,info,&
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -431,9 +431,9 @@ subroutine psb_sspsv(alpha,a,x,beta,y,desc_a,info,&
itrans = 'N'
endif
m = psb_cd_get_global_rows(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x)
lldy = size(y)

32
base/psblas/psb_sxdot.f90
Unescape Escape View File

@ -73,7 +73,7 @@ function psb_sxdot(x, y,desc_a, info, jx, jy)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
info = psb_err_context_error_
@ -101,7 +101,7 @@ function psb_sxdot(x, y,desc_a, info, jx, jy)
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -121,8 +121,8 @@ function psb_sxdot(x, y,desc_a, info, jx, jy)
end if
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
dot_local = sxdot(psb_cd_get_local_rows(desc_a),&
if(desc_a%get_local_rows() > 0) then
dot_local = sxdot(desc_a%get_local_rows(),&
& x(iix,jjx),ione,y(iiy,jjy),ione)
! adjust dot_local because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -225,7 +225,7 @@ function psb_sxdotv(x, y,desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -238,7 +238,7 @@ function psb_sxdotv(x, y,desc_a, info)
iy = ione
jx = ione
jy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,jx,desc_a,info,iix,jjx)
@ -258,8 +258,8 @@ function psb_sxdotv(x, y,desc_a, info)
end if
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
dot_local = sxdot(psb_cd_get_local_rows(desc_a),&
if(desc_a%get_local_rows() > 0) then
dot_local = sxdot(desc_a%get_local_rows(),&
& x,ione,y,ione)
! adjust dot_local because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -380,7 +380,7 @@ subroutine psb_sxdotvs(res, x, y,desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -391,7 +391,7 @@ subroutine psb_sxdotvs(res, x, y,desc_a, info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ix,desc_a,info,iix,jjx)
if (info == psb_success_) &
@ -410,8 +410,8 @@ subroutine psb_sxdotvs(res, x, y,desc_a, info)
end if
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
dot_local = sxdot(psb_cd_get_local_rows(desc_a),&
if(desc_a%get_local_rows() > 0) then
dot_local = sxdot(desc_a%get_local_rows(),&
& x,ione,y,ione)
! adjust dot_local because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -512,7 +512,7 @@ subroutine psb_sxmdots(res, x, y, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -524,7 +524,7 @@ subroutine psb_sxmdots(res, x, y, desc_a, info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ix,desc_a,info,iix,jjx)
@ -552,9 +552,9 @@ subroutine psb_sxmdots(res, x, y, desc_a, info)
allocate(dot_local(k))
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
if(desc_a%get_local_rows() > 0) then
do j=1,k
dot_local(j) = sxdot(psb_cd_get_local_rows(desc_a),&
dot_local(j) = sxdot(desc_a%get_local_rows(),&
& x(1,j),ione,y(1,j),ione)
! adjust dot_local because overlapped elements are computed more than once
end do

32
base/psblas/psb_zamax.f90
Unescape Escape View File

@ -74,7 +74,7 @@ function psb_zamax (x,desc_a, info, jx)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -90,7 +90,7 @@ function psb_zamax (x,desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -107,8 +107,8 @@ function psb_zamax (x,desc_a, info, jx)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=izamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix,jjx),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=izamax(desc_a%get_local_rows()-iix+1,x(iix,jjx),1)
amax=cabs1(x(iix+imax-1,jjx))
else
amax = dzero
@ -207,7 +207,7 @@ function psb_zamaxv (x,desc_a, info)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -219,7 +219,7 @@ function psb_zamaxv (x,desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -236,8 +236,8 @@ function psb_zamaxv (x,desc_a, info)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=izamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=izamax(desc_a%get_local_rows()-iix+1,x(iix),1)
cmax=(x(iix+imax-1))
amax=cabs1(cmax)
else
@ -334,7 +334,7 @@ subroutine psb_zamaxvs(res,x,desc_a, info)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -346,7 +346,7 @@ subroutine psb_zamaxvs(res,x,desc_a, info)
ix = 1
ijx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
@ -362,8 +362,8 @@ subroutine psb_zamaxvs(res,x,desc_a, info)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
imax=izamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix),1)
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
imax=izamax(desc_a%get_local_rows()-iix+1,x(iix),1)
cmax=(x(iix+imax-1))
amax=cabs1(cmax)
else
@ -459,7 +459,7 @@ subroutine psb_zmamaxs(res,x,desc_a, info,jx)
amax=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -475,7 +475,7 @@ subroutine psb_zmamaxs(res,x,desc_a, info,jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
k = min(size(x,2),size(res,1))
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -493,9 +493,9 @@ subroutine psb_zmamaxs(res,x,desc_a, info,jx)
end if
! compute local max
if ((psb_cd_get_local_rows(desc_a) > 0).and.(m /= 0)) then
if ((desc_a%get_local_rows() > 0).and.(m /= 0)) then
do i=1,k
imax=izamax(psb_cd_get_local_rows(desc_a)-iix+1,x(iix,jjx+i-1),1)
imax=izamax(desc_a%get_local_rows()-iix+1,x(iix,jjx+i-1),1)
cmax=(x(iix+imax-1,jjx+i-1))
res(i)=cabs1(cmax)
end do

24
base/psblas/psb_zasum.f90
Unescape Escape View File

@ -76,7 +76,7 @@ function psb_zasum (x,desc_a, info, jx)
asum=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -92,7 +92,7 @@ function psb_zasum (x,desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -111,8 +111,8 @@ function psb_zasum (x,desc_a, info, jx)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=dzasum(psb_cd_get_local_rows(desc_a)-iix+1,x(iix:,jjx),ione)
if(desc_a%get_local_rows() > 0) then
asum=dzasum(desc_a%get_local_rows()-iix+1,x(iix:,jjx),ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -223,7 +223,7 @@ function psb_zasumv(x,desc_a, info)
asum=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -235,7 +235,7 @@ function psb_zasumv(x,desc_a, info)
ix = 1
jx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
@ -254,8 +254,8 @@ function psb_zasumv(x,desc_a, info)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=dzasum(psb_cd_get_local_rows(desc_a),x,ione)
if(desc_a%get_local_rows() > 0) then
asum=dzasum(desc_a%get_local_rows(),x,ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -362,7 +362,7 @@ subroutine psb_zasumvs(res,x,desc_a, info)
asum=0.d0
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -374,7 +374,7 @@ subroutine psb_zasumvs(res,x,desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
@ -393,8 +393,8 @@ subroutine psb_zasumvs(res,x,desc_a, info)
! compute local max
if ((m /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
asum=dzasum(psb_cd_get_local_rows(desc_a),x,ione)
if(desc_a%get_local_rows() > 0) then
asum=dzasum(desc_a%get_local_rows(),x,ione)
! adjust asum because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)

16
base/psblas/psb_zaxpby.f90
Unescape Escape View File

@ -71,7 +71,7 @@ subroutine psb_zaxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
info = psb_err_context_error_
@ -110,7 +110,7 @@ subroutine psb_zaxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -130,8 +130,8 @@ subroutine psb_zaxpby(alpha, x, beta,y,desc_a,info, n, jx, jy)
end if
if ((in /= 0)) then
if(psb_cd_get_local_rows(desc_a) > 0) then
call zaxpby(psb_cd_get_local_cols(desc_a),in,&
if(desc_a%get_local_rows() > 0) then
call zaxpby(desc_a%get_local_cols(),in,&
& alpha,x(iix:,jjx),size(x,1),beta,&
& y(iiy:,jjy),size(y,1),info)
end if
@ -221,7 +221,7 @@ subroutine psb_zaxpbyv(alpha, x, beta,y,desc_a,info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -233,7 +233,7 @@ subroutine psb_zaxpbyv(alpha, x, beta,y,desc_a,info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x),ix,ione,desc_a,info,iix,jjx)
@ -256,8 +256,8 @@ subroutine psb_zaxpbyv(alpha, x, beta,y,desc_a,info)
call psb_errpush(info,name)
end if
if(psb_cd_get_local_rows(desc_a) > 0) then
call zaxpby(psb_cd_get_local_cols(desc_a),ione,&
if(desc_a%get_local_rows() > 0) then
call zaxpby(desc_a%get_local_cols(),ione,&
& alpha,x,size(x),beta,&
& y,size(y),info)
end if

32
base/psblas/psb_zdot.f90
Unescape Escape View File

@ -73,7 +73,7 @@ function psb_zdot(x, y,desc_a, info, jx, jy)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
info = psb_err_context_error_
@ -101,7 +101,7 @@ function psb_zdot(x, y,desc_a, info, jx, jy)
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ijx,desc_a,info,iix,jjx)
@ -121,8 +121,8 @@ function psb_zdot(x, y,desc_a, info, jx, jy)
end if
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
dot_local = zdotc(psb_cd_get_local_rows(desc_a),&
if(desc_a%get_local_rows() > 0) then
dot_local = zdotc(desc_a%get_local_rows(),&
& x(iix:,jjx),ione,y(iiy:,jjy),ione)
! adjust dot_local because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -225,7 +225,7 @@ function psb_zdotv(x, y,desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -238,7 +238,7 @@ function psb_zdotv(x, y,desc_a, info)
iy = ione
jx = ione
jy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,jx,desc_a,info,iix,jjx)
@ -258,8 +258,8 @@ function psb_zdotv(x, y,desc_a, info)
end if
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
dot_local = zdotc(psb_cd_get_local_rows(desc_a),&
if(desc_a%get_local_rows() > 0) then
dot_local = zdotc(desc_a%get_local_rows(),&
& x,ione,y,ione)
! adjust dot_local because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -359,7 +359,7 @@ subroutine psb_zdotvs(res, x, y,desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -370,7 +370,7 @@ subroutine psb_zdotvs(res, x, y,desc_a, info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ix,desc_a,info,iix,jjx)
if (info == psb_success_) &
@ -389,8 +389,8 @@ subroutine psb_zdotvs(res, x, y,desc_a, info)
end if
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
dot_local = zdotc(psb_cd_get_local_rows(desc_a),&
if(desc_a%get_local_rows() > 0) then
dot_local = zdotc(desc_a%get_local_rows(),&
& x,ione,y,ione)
! adjust dot_local because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -491,7 +491,7 @@ subroutine psb_zmdots(res, x, y, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -ione) then
@ -503,7 +503,7 @@ subroutine psb_zmdots(res, x, y, desc_a, info)
ix = ione
iy = ione
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
! check vector correctness
call psb_chkvect(m,ione,size(x,1),ix,ix,desc_a,info,iix,jjx)
@ -531,9 +531,9 @@ subroutine psb_zmdots(res, x, y, desc_a, info)
allocate(dot_local(k))
if(m /= 0) then
if(psb_cd_get_local_rows(desc_a) > 0) then
if(desc_a%get_local_rows() > 0) then
do j=1,k
dot_local(j) = zdotc(psb_cd_get_local_rows(desc_a),&
dot_local(j) = zdotc(desc_a%get_local_rows(),&
& x(1:,j),ione,y(1:,j),ione)
! adjust dot_local because overlapped elements are computed more than once
end do

24
base/psblas/psb_znrm2.f90
Unescape Escape View File

@ -70,7 +70,7 @@ function psb_znrm2(x, desc_a, info, jx)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -86,7 +86,7 @@ function psb_znrm2(x, desc_a, info, jx)
ijx = 1
endif
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x,1),ix,ijx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
@ -101,8 +101,8 @@ function psb_znrm2(x, desc_a, info, jx)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = dznrm2( ndim, x(iix:,jjx), ione )
! adjust because overlapped elements are computed more than once
@ -206,7 +206,7 @@ function psb_znrm2v(x, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -217,7 +217,7 @@ function psb_znrm2v(x, desc_a, info)
ix = 1
jx=1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -233,8 +233,8 @@ function psb_znrm2v(x, desc_a, info)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = dznrm2( ndim, x, ione )
! adjust because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)
@ -339,7 +339,7 @@ subroutine psb_znrm2vs(res, x, desc_a, info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -350,7 +350,7 @@ subroutine psb_znrm2vs(res, x, desc_a, info)
ix = 1
jx = 1
m = psb_cd_get_global_rows(desc_a)
m = desc_a%get_global_rows()
call psb_chkvect(m,1,size(x),ix,jx,desc_a,info,iix,jjx)
if(info /= psb_success_) then
@ -366,8 +366,8 @@ subroutine psb_znrm2vs(res, x, desc_a, info)
end if
if(m /= 0) then
if (psb_cd_get_local_rows(desc_a) > 0) then
ndim = psb_cd_get_local_rows(desc_a)
if (desc_a%get_local_rows() > 0) then
ndim = desc_a%get_local_rows()
nrm2 = dznrm2( ndim, x, ione )
! adjust because overlapped elements are computed more than once
do i=1,size(desc_a%ovrlap_elem,1)

6
base/psblas/psb_znrmi.f90
Unescape Escape View File

@ -66,7 +66,7 @@ function psb_znrmi(a,desc_a,info)
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -77,8 +77,8 @@ function psb_znrmi(a,desc_a,info)
ia = 1
ja = 1
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
call psb_chkmat(m,n,ia,ja,desc_a,info,iia,jja)
if(info /= psb_success_) then

20
base/psblas/psb_zspmm.f90
Unescape Escape View File

@ -99,7 +99,7 @@ subroutine psb_zspmm(alpha,a,x,beta,y,desc_a,info,&
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -151,10 +151,10 @@ subroutine psb_zspmm(alpha,a,x,beta,y,desc_a,info,&
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x,1)
lldy = size(y,1)
@ -450,7 +450,7 @@ subroutine psb_zspmv(alpha,a,x,beta,y,desc_a,info,&
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
info = psb_err_context_error_
@ -486,10 +486,10 @@ subroutine psb_zspmv(alpha,a,x,beta,y,desc_a,info,&
goto 9999
end if
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x)
lldy = size(y)

16
base/psblas/psb_zspsm.f90
Unescape Escape View File

@ -109,7 +109,7 @@ subroutine psb_zspsm(alpha,a,x,beta,y,desc_a,info,&
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -168,9 +168,9 @@ subroutine psb_zspsm(alpha,a,x,beta,y,desc_a,info,&
itrans = 'N'
endif
m = psb_cd_get_global_rows(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x,1)
lldy = size(y,1)
@ -387,7 +387,7 @@ subroutine psb_zspsv(alpha,a,x,beta,y,desc_a,info,&
info=psb_success_
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -430,9 +430,9 @@ subroutine psb_zspsv(alpha,a,x,beta,y,desc_a,info,&
itrans = 'N'
endif
m = psb_cd_get_global_rows(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
m = desc_a%get_global_rows()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
lldx = size(x)
lldy = size(y)

12
base/tools/psb_callc.f90
Unescape Escape View File

@ -66,7 +66,7 @@ subroutine psb_calloc(x, desc_a, info, n, lb)
int_err(1)=0
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -103,9 +103,9 @@ subroutine psb_calloc(x, desc_a, info, n, lb)
!....allocate x .....
if (psb_is_asb_desc(desc_a).or.psb_is_upd_desc(desc_a)) then
nr = max(1,psb_cd_get_local_cols(desc_a))
nr = max(1,desc_a%get_local_cols())
else if (psb_is_bld_desc(desc_a)) then
nr = max(1,psb_cd_get_local_rows(desc_a))
nr = max(1,desc_a%get_local_rows())
else
info = psb_err_internal_error_
call psb_errpush(info,name,int_err,a_err='Invalid desc_a')
@ -199,7 +199,7 @@ subroutine psb_callocv(x, desc_a,info,n)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
! ....verify blacs grid correctness..
@ -220,9 +220,9 @@ subroutine psb_callocv(x, desc_a,info,n)
!....allocate x .....
if (psb_is_asb_desc(desc_a).or.psb_is_upd_desc(desc_a)) then
nr = max(1,psb_cd_get_local_cols(desc_a))
nr = max(1,desc_a%get_local_cols())
else if (psb_is_bld_desc(desc_a)) then
nr = max(1,psb_cd_get_local_rows(desc_a))
nr = max(1,desc_a%get_local_rows())
else
info = psb_err_internal_error_
call psb_errpush(info,name,int_err,a_err='Invalid desc_a')

18
base/tools/psb_casb.f90
Unescape Escape View File

@ -68,14 +68,14 @@ subroutine psb_casb(x, desc_a, info)
call psb_errpush(info,name)
goto 9999
endif
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (debug_level >= psb_debug_ext_) &
& write(debug_unit,*) me,' ',trim(name),': start: ',np,&
& psb_cd_get_dectype(desc_a)
& desc_a%get_dectype()
! ....verify blacs grid correctness..
if (np == -1) then
info = psb_err_context_error_
@ -84,16 +84,16 @@ subroutine psb_casb(x, desc_a, info)
else if (.not.psb_is_asb_desc(desc_a)) then
if (debug_level >= psb_debug_ext_) &
& write(debug_unit,*) me,' ',trim(name),' error ',&
& psb_cd_get_dectype(desc_a)
& desc_a%get_dectype()
info = psb_err_input_matrix_unassembled_
call psb_errpush(info,name)
goto 9999
endif
! check size
ictxt = psb_cd_get_context(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
ictxt = desc_a%get_context()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
i1sz = size(x,dim=1)
i2sz = size(x,dim=2)
if (debug_level >= psb_debug_ext_) &
@ -193,7 +193,7 @@ subroutine psb_casbv(x, desc_a, info)
int_err(1) = 0
name = 'psb_cgeasb_v'
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
@ -210,8 +210,8 @@ subroutine psb_casbv(x, desc_a, info)
goto 9999
endif
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
if (debug_level >= psb_debug_ext_) &
& write(debug_unit,*) me,' ',trim(name),': sizes: ',nrow,ncol
i1sz = size(x)

20
base/tools/psb_ccdbldext.F90
Unescape Escape View File

@ -105,8 +105,8 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
Call psb_info(ictxt, me, np)
If (debug_level >= psb_debug_outer_) &
@ -118,9 +118,9 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)
else
extype_ = psb_ovt_xhal_
endif
m = psb_cd_get_local_rows(desc_a)
n_row = psb_cd_get_local_rows(desc_a)
n_col = psb_cd_get_local_cols(desc_a)
m = desc_a%get_local_rows()
n_row = desc_a%get_local_rows()
n_col = desc_a%get_local_cols()
nhalo = n_col-m
if (novr<0) then
@ -317,7 +317,7 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)
counter = 1
counter_t = 1
n_col_prev = psb_cd_get_local_cols(desc_ov)
n_col_prev = desc_ov%get_local_cols()
Do While (halo(counter) /= -1)
tot_elem=0
@ -543,7 +543,7 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)
& write(debug_unit,*) me,' ',trim(name),&
& ': going for fnd_owner', desc_ov%indxmap%get_state()
call psi_fnd_owner(iszs,works,temp,desc_a,info)
n_col = psb_cd_get_local_cols(desc_ov)
n_col = desc_ov%get_local_cols()
if (debug_level >= psb_debug_outer_) &
& write(debug_unit,*) me,' ',trim(name),&
@ -551,9 +551,9 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)
do i=1,iszs
idx = works(i)
n_col = psb_cd_get_local_cols(desc_ov)
n_col = desc_ov%get_local_cols()
call psi_idx_ins_cnv(idx,lidx,desc_ov,info)
if (psb_cd_get_local_cols(desc_ov) > n_col ) then
if (desc_ov%get_local_cols() > n_col ) then
!
! This is a new index. Assigning a local index as
! we receive them guarantees that all indices for HALO(I)
@ -575,7 +575,7 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)
counter_t = counter_t+3
endif
end Do
n_col = psb_cd_get_local_cols(desc_ov)
n_col = desc_ov%get_local_cols()
end if

10
base/tools/psb_cd_lstext.f90
Unescape Escape View File

@ -69,8 +69,8 @@ Subroutine psb_cd_lstext(desc_a,in_list,desc_ov,info, mask,extype)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
Call psb_info(ictxt, me, np)
If (debug_level >= psb_debug_outer_) &
@ -78,9 +78,9 @@ Subroutine psb_cd_lstext(desc_a,in_list,desc_ov,info, mask,extype)
m = psb_cd_get_local_rows(desc_a)
n_row = psb_cd_get_local_rows(desc_a)
n_col = psb_cd_get_local_cols(desc_a)
m = desc_a%get_local_rows()
n_row = desc_a%get_local_rows()
n_col = desc_a%get_local_cols()
nhalo = n_col-n_row
nl = size(in_list)

4
base/tools/psb_cd_reinit.f90
Unescape Escape View File

@ -55,8 +55,8 @@ Subroutine psb_cd_reinit(desc,info)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt = psb_cd_get_context(desc)
icomm = psb_cd_get_mpic(desc)
ictxt = desc%get_context()
icomm = desc%get_mpic()
Call psb_info(ictxt, me, np)
if (debug_level >= psb_debug_outer_) &
& write(debug_unit,*) me,' ',trim(name),': start'

2
base/tools/psb_cd_set_bld.f90
Unescape Escape View File

@ -63,7 +63,7 @@ subroutine psb_cd_set_bld(desc,info)
call psb_erractionsave(err_act)
name = 'psb_cd_set_bld'
ictxt = psb_cd_get_context(desc)
ictxt = desc%get_context()
if (debug) write(psb_err_unit,*)'Entered CDSETBLD',ictxt
! check on blacs grid

10
base/tools/psb_cd_switch_ovl_indxmap.f90
Unescape Escape View File

@ -58,17 +58,17 @@ Subroutine psb_cd_switch_ovl_indxmap(desc,info)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt = psb_cd_get_context(desc)
icomm = psb_cd_get_mpic(desc)
ictxt = desc%get_context()
icomm = desc%get_mpic()
Call psb_info(ictxt, me, np)
If (debug_level >= psb_debug_outer_) &
& Write(debug_unit,*) me,' ',trim(name),&
& ': start'
mglob = psb_cd_get_global_rows(desc)
n_row = psb_cd_get_local_rows(desc)
n_col = psb_cd_get_local_cols(desc)
mglob = desc%get_global_rows()
n_row = desc%get_local_rows()
n_col = desc%get_local_cols()
Allocate(vl(n_col),stat=info)
if (info /= psb_success_) then

2
base/tools/psb_cdcpy.F90
Unescape Escape View File

@ -61,7 +61,7 @@ subroutine psb_cdcpy(desc_in, desc_out, info)
call psb_erractionsave(err_act)
name = 'psb_cdcpy'
ictxt = psb_cd_get_context(desc_in)
ictxt = desc_in%get_context()
! check on blacs grid
call psb_info(ictxt, me, np)

24
base/tools/psb_cdins.f90
Unescape Escape View File

@ -70,12 +70,12 @@ subroutine psb_cdinsrc(nz,ia,ja,desc_a,info,ila,jla)
name = 'psb_cdins'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
dectype = psb_cd_get_dectype(desc_a)
mglob = psb_cd_get_global_rows(desc_a)
nglob = psb_cd_get_global_cols(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
ictxt = desc_a%get_context()
dectype = desc_a%get_dectype()
mglob = desc_a%get_global_rows()
nglob = desc_a%get_global_cols()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
call psb_info(ictxt, me, np)
@ -195,12 +195,12 @@ subroutine psb_cdinsc(nz,ja,desc,info,jla,mask)
name = 'psb_cdins'
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc)
dectype = psb_cd_get_dectype(desc)
mglob = psb_cd_get_global_rows(desc)
nglob = psb_cd_get_global_cols(desc)
nrow = psb_cd_get_local_rows(desc)
ncol = psb_cd_get_local_cols(desc)
ictxt = desc%get_context()
dectype = desc%get_dectype()
mglob = desc%get_global_rows()
nglob = desc%get_global_cols()
nrow = desc%get_local_rows()
ncol = desc%get_local_cols()
call psb_info(ictxt, me, np)

8
base/tools/psb_cdren.f90
Unescape Escape View File

@ -69,10 +69,10 @@ subroutine psb_cdren(trans,iperm,desc_a,info)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt = psb_cd_get_context(desc_a)
dectype = psb_cd_get_dectype(desc_a)
n_row = psb_cd_get_local_rows(desc_a)
n_col = psb_cd_get_local_cols(desc_a)
ictxt = desc_a%get_context()
dectype = desc_a%get_dectype()
n_row = desc_a%get_local_rows()
n_col = desc_a%get_local_cols()
! check on blacs grid
call psb_info(ictxt, me, np)

4
base/tools/psb_cfree.f90
Unescape Escape View File

@ -62,7 +62,7 @@ subroutine psb_cfree(x, desc_a, info)
return
end if
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
! ....verify blacs grid correctness..
@ -133,7 +133,7 @@ subroutine psb_cfreev(x, desc_a, info)
call psb_errpush(info,name)
goto 9999
end if
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then

24
base/tools/psb_cins.f90
Unescape Escape View File

@ -82,7 +82,7 @@ subroutine psb_cinsvi(m, irw, val, x, desc_a, info, dupl)
return
end if
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -100,10 +100,10 @@ subroutine psb_cinsvi(m, irw, val, x, desc_a, info, dupl)
goto 9999
else if (.not.psb_is_ok_desc(desc_a)) then
info = psb_err_input_matrix_unassembled_
int_err(1) = psb_cd_get_dectype(desc_a)
int_err(1) = desc_a%get_dectype()
call psb_errpush(info,name,int_err)
goto 9999
else if (size(x, dim=1) < psb_cd_get_local_rows(desc_a)) then
else if (size(x, dim=1) < desc_a%get_local_rows()) then
info = 310
int_err(1) = 5
int_err(2) = 4
@ -112,9 +112,9 @@ subroutine psb_cinsvi(m, irw, val, x, desc_a, info, dupl)
endif
if (m == 0) return
loc_rows = psb_cd_get_local_rows(desc_a)
loc_cols = psb_cd_get_local_cols(desc_a)
mglob = psb_cd_get_global_rows(desc_a)
loc_rows = desc_a%get_local_rows()
loc_cols = desc_a%get_local_cols()
mglob = desc_a%get_global_rows()
allocate(irl(m),stat=info)
if (info /= psb_success_) then
@ -263,7 +263,7 @@ subroutine psb_cinsi(m, irw, val, x, desc_a, info, dupl)
return
end if
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -281,10 +281,10 @@ subroutine psb_cinsi(m, irw, val, x, desc_a, info, dupl)
goto 9999
else if (.not.psb_is_ok_desc(desc_a)) then
info = psb_err_input_matrix_unassembled_
int_err(1) = psb_cd_get_dectype(desc_a)
int_err(1) = desc_a%get_dectype()
call psb_errpush(info,name,int_err)
goto 9999
else if (size(x, dim=1) < psb_cd_get_local_rows(desc_a)) then
else if (size(x, dim=1) < desc_a%get_local_rows()) then
info = 310
int_err(1) = 5
int_err(2) = 4
@ -293,9 +293,9 @@ subroutine psb_cinsi(m, irw, val, x, desc_a, info, dupl)
endif
if (m == 0) return
loc_rows = psb_cd_get_local_rows(desc_a)
loc_cols = psb_cd_get_local_cols(desc_a)
mglob = psb_cd_get_global_rows(desc_a)
loc_rows = desc_a%get_local_rows()
loc_cols = desc_a%get_local_cols()
mglob = desc_a%get_global_rows()
n = min(size(val,2),size(x,2))

14
base/tools/psb_cspalloc.f90
Unescape Escape View File

@ -66,8 +66,8 @@ subroutine psb_cspalloc(a, desc_a, info, nnz)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt = psb_cd_get_context(desc_a)
dectype = psb_cd_get_dectype(desc_a)
ictxt = desc_a%get_context()
dectype = desc_a%get_dectype()
call psb_info(ictxt, me, np)
! ....verify blacs grid correctness..
@ -78,10 +78,10 @@ subroutine psb_cspalloc(a, desc_a, info, nnz)
endif
loc_row = psb_cd_get_local_rows(desc_a)
loc_col = psb_cd_get_local_cols(desc_a)
m = psb_cd_get_global_rows(desc_a)
n = psb_cd_get_global_cols(desc_a)
loc_row = desc_a%get_local_rows()
loc_col = desc_a%get_local_cols()
m = desc_a%get_global_rows()
n = desc_a%get_global_cols()
!...allocate matrix data...
if (present(nnz))then
@ -113,7 +113,7 @@ subroutine psb_cspalloc(a, desc_a, info, nnz)
if (debug_level >= psb_debug_ext_) &
& write(debug_unit,*) me,' ',trim(name),': ', &
& psb_cd_get_dectype(desc_a),psb_desc_bld_
& desc_a%get_dectype(),psb_desc_bld_
return

12
base/tools/psb_cspasb.f90
Unescape Escape View File

@ -76,9 +76,9 @@ subroutine psb_cspasb(a,desc_a, info, afmt, upd, dupl,mold)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt = psb_cd_get_context(desc_a)
n_row = psb_cd_get_local_rows(desc_a)
n_col = psb_cd_get_local_cols(desc_a)
ictxt = desc_a%get_context()
n_row = desc_a%get_local_rows()
n_col = desc_a%get_local_cols()
! check on BLACS grid
call psb_info(ictxt, me, np)
@ -90,7 +90,7 @@ subroutine psb_cspasb(a,desc_a, info, afmt, upd, dupl,mold)
if (.not.psb_is_asb_desc(desc_a)) then
info = psb_err_spmat_invalid_state_
int_err(1) = psb_cd_get_dectype(desc_a)
int_err(1) = desc_a%get_dectype()
call psb_errpush(info,name)
goto 9999
endif
@ -107,8 +107,8 @@ subroutine psb_cspasb(a,desc_a, info, afmt, upd, dupl,mold)
! First case: we come from a fresh build.
!
n_row = psb_cd_get_local_rows(desc_a)
n_col = psb_cd_get_local_cols(desc_a)
n_row = desc_a%get_local_rows()
n_col = desc_a%get_local_cols()
call a%set_nrows(n_row)
call a%set_ncols(n_col)
end if

2
base/tools/psb_cspfree.f90
Unescape Escape View File

@ -61,7 +61,7 @@ subroutine psb_cspfree(a, desc_a,info)
call psb_errpush(info,name)
return
else
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
end if
!...deallocate a....

4
base/tools/psb_csphalo.F90
Unescape Escape View File

@ -133,8 +133,8 @@ Subroutine psb_csphalo(a,desc_a,blk,info,rowcnv,colcnv,&
outfmt_ = 'CSR'
endif
ictxt = psb_cd_get_context(desc_a)
icomm = psb_cd_get_mpic(desc_a)
ictxt = desc_a%get_context()
icomm = desc_a%get_mpic()
Call psb_info(ictxt, me, np)

16
base/tools/psb_cspins.f90
Unescape Escape View File

@ -75,7 +75,7 @@ subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild)
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
@ -131,8 +131,8 @@ subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild)
& a_err=ch_err,i_err=(/info,0,0,0,0/))
goto 9999
end if
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
if (a%is_bld()) then
call a%csput(nz,ila,jla,val,1,nrow,1,ncol,info)
@ -163,8 +163,8 @@ subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild)
jla(1:nz) = ja(1:nz)
call psb_glob_to_loc(ila(1:nz),desc_a,info,iact='I')
call psb_glob_to_loc(jla(1:nz),desc_a,info,iact='I')
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
call a%csput(nz,ila,jla,val,1,nrow,1,ncol,info)
if (info /= psb_success_) then
@ -220,7 +220,7 @@ subroutine psb_cspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
call psb_erractionsave(err_act)
ictxt = psb_cd_get_context(desc_ar)
ictxt = desc_ar%get_context()
call psb_info(ictxt, me, np)
@ -280,8 +280,8 @@ subroutine psb_cspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
goto 9999
end if
nrow = psb_cd_get_local_rows(desc_ar)
ncol = psb_cd_get_local_cols(desc_ac)
nrow = desc_ar%get_local_rows()
ncol = desc_ac%get_local_cols()
call a%csput(nz,ila,jla,val,1,nrow,1,ncol,info)
if (info /= psb_success_) then

2
base/tools/psb_csprn.f90
Unescape Escape View File

@ -67,7 +67,7 @@ Subroutine psb_csprn(a, desc_a,info,clear)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (debug_level >= psb_debug_outer_) &
& write(debug_unit,*) me,' ',trim(name),': start '

12
base/tools/psb_dallc.f90
Unescape Escape View File

@ -66,7 +66,7 @@ subroutine psb_dalloc(x, desc_a, info, n, lb)
int_err(1)=0
call psb_erractionsave(err_act)
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
if (np == -1) then
@ -103,9 +103,9 @@ subroutine psb_dalloc(x, desc_a, info, n, lb)
!....allocate x .....
if (psb_is_asb_desc(desc_a).or.psb_is_upd_desc(desc_a)) then
nr = max(1,psb_cd_get_local_cols(desc_a))
nr = max(1,desc_a%get_local_cols())
else if (psb_is_bld_desc(desc_a)) then
nr = max(1,psb_cd_get_local_rows(desc_a))
nr = max(1,desc_a%get_local_rows())
else
info = psb_err_internal_error_
call psb_errpush(info,name,int_err,a_err='Invalid desc_a')
@ -200,7 +200,7 @@ subroutine psb_dallocv(x, desc_a,info,n)
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
ictxt=psb_cd_get_context(desc_a)
ictxt=desc_a%get_context()
call psb_info(ictxt, me, np)
! ....verify blacs grid correctness..
@ -221,9 +221,9 @@ subroutine psb_dallocv(x, desc_a,info,n)
!....allocate x .....
if (psb_is_asb_desc(desc_a).or.psb_is_upd_desc(desc_a)) then
nr = max(1,psb_cd_get_local_cols(desc_a))
nr = max(1,desc_a%get_local_cols())
else if (psb_is_bld_desc(desc_a)) then
nr = max(1,psb_cd_get_local_rows(desc_a))
nr = max(1,desc_a%get_local_rows())
else
info = psb_err_internal_error_
call psb_errpush(info,name,int_err,a_err='Invalid desc_a')

18
base/tools/psb_dasb.f90
Unescape Escape View File

@ -68,14 +68,14 @@ subroutine psb_dasb(x, desc_a, info)
call psb_errpush(info,name)
goto 9999
endif
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
call psb_info(ictxt, me, np)
if (debug_level >= psb_debug_ext_) &
& write(debug_unit,*) me,' ',trim(name),': start: ',np,&
& psb_cd_get_dectype(desc_a)
& desc_a%get_dectype()
! ....verify blacs grid correctness..
if (np == -1) then
info = psb_err_context_error_
@ -84,16 +84,16 @@ subroutine psb_dasb(x, desc_a, info)
else if (.not.psb_is_asb_desc(desc_a)) then
if (debug_level >= psb_debug_ext_) &
& write(debug_unit,*) me,' ',trim(name),' error ',&
& psb_cd_get_dectype(desc_a)
& desc_a%get_dectype()
info = psb_err_input_matrix_unassembled_
call psb_errpush(info,name)
goto 9999
endif
! check size
ictxt = psb_cd_get_context(desc_a)
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
ictxt = desc_a%get_context()
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
i1sz = size(x,dim=1)
i2sz = size(x,dim=2)
if (debug_level >= psb_debug_ext_) &
@ -193,7 +193,7 @@ subroutine psb_dasbv(x, desc_a, info)
int_err(1) = 0
name = 'psb_dgeasb_v'
ictxt = psb_cd_get_context(desc_a)
ictxt = desc_a%get_context()
debug_unit = psb_get_debug_unit()
debug_level = psb_get_debug_level()
@ -210,8 +210,8 @@ subroutine psb_dasbv(x, desc_a, info)
goto 9999
endif
nrow = psb_cd_get_local_rows(desc_a)
ncol = psb_cd_get_local_cols(desc_a)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
if (debug_level >= psb_debug_ext_) &
& write(debug_unit,*) me,' ',trim(name),': sizes: ',nrow,ncol
i1sz = size(x)

Some files were not shown because too many files have changed in this diff Show More

Write Preview
Loading…
Cancel
Save