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

mld2p4-2

Browse Source 
mlprec/mld_c_prec_type.f90
 mlprec/mld_d_prec_type.f90
 mlprec/mld_s_prec_type.f90
 mlprec/mld_z_prec_type.f90

Add optional argument to PRECDESCR to allow printing from any
process.
stopcriterion
Salvatore Filippone 12 years ago
parent 71dd991b9b
commit 1d40ab3303
4 changed files with 52 additions and 9 deletions
Show Stats Download Patch File Download Diff File

15
mlprec/mld_c_prec_type.f90
Unescape Escape View File

@ -401,19 +401,24 @@ contains
! The id of the file where the preconditioner description
! will be printed. If iout is not present, then the standard
! output is condidered.
! root - integer, input, optional.
! The id of the process printing the message; -1 acts as a wildcard.
! Default is psb_root_
!
subroutine mld_cfile_prec_descr(p,info,iout)
subroutine mld_cfile_prec_descr(p,info,iout,root)
implicit none
! Arguments
type(mld_cprec_type), intent(in) :: p
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(in), optional :: iout
integer(psb_ipk_), intent(in), optional :: root
! Local variables
integer(psb_ipk_) :: ilev, nlev
integer(psb_ipk_) :: ictxt, me, np
character(len=20), parameter :: name='mld_file_prec_descr'
integer(psb_ipk_) :: iout_
integer(psb_ipk_) :: root_
info = psb_success_
if (present(iout)) then
@ -428,6 +433,12 @@ contains
if (allocated(p%precv)) then
call psb_info(ictxt,me,np)
if (present(root)) then
root_ = root
else
root_ = psb_root_
end if
if (root_ == -1) root_ = me
!
! The preconditioner description is printed by processor psb_root_.
@ -435,7 +446,7 @@ contains
! preconditioner have the same values on all the procs (this is
! ensured by mld_precbld).
!
if (me == psb_root_) then
if (me == root_) then
nlev = size(p%precv)
do ilev = 1, nlev
if (.not.allocated(p%precv(ilev)%sm)) then

16
mlprec/mld_d_prec_type.f90
Unescape Escape View File

@ -401,19 +401,24 @@ contains
! The id of the file where the preconditioner description
! will be printed. If iout is not present, then the standard
! output is condidered.
! root - integer, input, optional.
! The id of the process printing the message; -1 acts as a wildcard.
! Default is psb_root_
!
subroutine mld_dfile_prec_descr(p,info,iout)
subroutine mld_dfile_prec_descr(p,info,iout,root)
implicit none
! Arguments
type(mld_dprec_type), intent(in) :: p
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(in), optional :: iout
integer(psb_ipk_), intent(in), optional :: root
! Local variables
integer(psb_ipk_) :: ilev, nlev
integer(psb_ipk_) :: ictxt, me, np
character(len=20), parameter :: name='mld_file_prec_descr'
integer(psb_ipk_) :: iout_
integer(psb_ipk_) :: root_
info = psb_success_
if (present(iout)) then
@ -428,14 +433,19 @@ contains
if (allocated(p%precv)) then
call psb_info(ictxt,me,np)
if (present(root)) then
root_ = root
else
root_ = psb_root_
end if
if (root_ == -1) root_ = me
!
! The preconditioner description is printed by processor psb_root_.
! This agrees with the fact that all the parameters defining the
! preconditioner have the same values on all the procs (this is
! ensured by mld_precbld).
!
if (me == psb_root_) then
if (me == root_) then
nlev = size(p%precv)
do ilev = 1, nlev
if (.not.allocated(p%precv(ilev)%sm)) then

15
mlprec/mld_s_prec_type.f90
Unescape Escape View File

@ -401,19 +401,24 @@ contains
! The id of the file where the preconditioner description
! will be printed. If iout is not present, then the standard
! output is condidered.
! root - integer, input, optional.
! The id of the process printing the message; -1 acts as a wildcard.
! Default is psb_root_
!
subroutine mld_sfile_prec_descr(p,info,iout)
subroutine mld_sfile_prec_descr(p,info,iout,root)
implicit none
! Arguments
type(mld_sprec_type), intent(in) :: p
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(in), optional :: iout
integer(psb_ipk_), intent(in), optional :: root
! Local variables
integer(psb_ipk_) :: ilev, nlev
integer(psb_ipk_) :: ictxt, me, np
character(len=20), parameter :: name='mld_file_prec_descr'
integer(psb_ipk_) :: iout_
integer(psb_ipk_) :: root_
info = psb_success_
if (present(iout)) then
@ -428,6 +433,12 @@ contains
if (allocated(p%precv)) then
call psb_info(ictxt,me,np)
if (present(root)) then
root_ = root
else
root_ = psb_root_
end if
if (root_ == -1) root_ = me
!
! The preconditioner description is printed by processor psb_root_.
@ -435,7 +446,7 @@ contains
! preconditioner have the same values on all the procs (this is
! ensured by mld_precbld).
!
if (me == psb_root_) then
if (me == root_) then
nlev = size(p%precv)
do ilev = 1, nlev
if (.not.allocated(p%precv(ilev)%sm)) then

15
mlprec/mld_z_prec_type.f90
Unescape Escape View File

@ -401,19 +401,24 @@ contains
! The id of the file where the preconditioner description
! will be printed. If iout is not present, then the standard
! output is condidered.
! root - integer, input, optional.
! The id of the process printing the message; -1 acts as a wildcard.
! Default is psb_root_
!
subroutine mld_zfile_prec_descr(p,info,iout)
subroutine mld_zfile_prec_descr(p,info,iout,root)
implicit none
! Arguments
type(mld_zprec_type), intent(in) :: p
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(in), optional :: iout
integer(psb_ipk_), intent(in), optional :: root
! Local variables
integer(psb_ipk_) :: ilev, nlev
integer(psb_ipk_) :: ictxt, me, np
character(len=20), parameter :: name='mld_file_prec_descr'
integer(psb_ipk_) :: iout_
integer(psb_ipk_) :: root_
info = psb_success_
if (present(iout)) then
@ -428,6 +433,12 @@ contains
if (allocated(p%precv)) then
call psb_info(ictxt,me,np)
if (present(root)) then
root_ = root
else
root_ = psb_root_
end if
if (root_ == -1) root_ = me
!
! The preconditioner description is printed by processor psb_root_.
@ -435,7 +446,7 @@ contains
! preconditioner have the same values on all the procs (this is
! ensured by mld_precbld).
!
if (me == psb_root_) then
if (me == root_) then
nlev = size(p%precv)
do ilev = 1, nlev
if (.not.allocated(p%precv(ilev)%sm)) then

Write Preview
Loading…
Cancel
Save