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

Fix handling of large process count in filename generation.

Browse Source
stopcriterion
Salvatore Filippone 5 years ago
parent de25234293
commit ac1c81f112
4 changed files with 44 additions and 20 deletions
Show Stats Download Patch File Download Diff File

16
mlprec/impl/level/mld_c_base_onelev_dump.f90
Unescape Escape View File

@ -47,10 +47,10 @@ subroutine mld_c_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
character(len=*), intent(in), optional :: prefix, head character(len=*), intent(in), optional :: prefix, head
logical, optional, intent(in) :: ac, rp, smoother, solver, tprol, global_num logical, optional, intent(in) :: ac, rp, smoother, solver, tprol, global_num
! Local variables ! Local variables
integer(psb_ipk_) :: i, j, il1, iln, lname, lev integer(psb_ipk_) :: i, j, il1, iln, lname, lev, ni
integer(psb_ipk_) :: icontxt,iam, np integer(psb_ipk_) :: icontxt,iam, np
character(len=80) :: prefix_ character(len=80) :: prefix_, frmt
character(len=120) :: fname ! len should be at least 20 more than character(len=1024) :: fname
logical :: ac_, rp_, tprol_, global_num_ logical :: ac_, rp_, tprol_, global_num_
integer(psb_ipk_), allocatable :: ivr(:), ivc(:) integer(psb_ipk_), allocatable :: ivr(:), ivc(:)
@ -68,6 +68,7 @@ subroutine mld_c_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
else else
icontxt = -1 icontxt = -1
iam = -1 iam = -1
np = -1
end if end if
if (present(ac)) then if (present(ac)) then
ac_ = ac ac_ = ac
@ -91,8 +92,13 @@ subroutine mld_c_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
end if end if
lname = len_trim(prefix_) lname = len_trim(prefix_)
fname = trim(prefix_) fname = trim(prefix_)
write(fname(lname+1:lname+5),'(a,i3.3)') '_p',iam
lname = lname + 5 if (np > 0) then
ni = floor(log10(1.0*np)) + 1
write(frmt,'(a,i3.3,a,i3.3,a)') '(a,i',ni,'.',ni,')'
write(fname(lname+1:lname+ni+2),frmt) '_p',iam
lname = lname + ni + 2
end if
if (global_num_) then if (global_num_) then
if (level >= 2) then if (level >= 2) then

16
mlprec/impl/level/mld_d_base_onelev_dump.f90
Unescape Escape View File

@ -47,10 +47,10 @@ subroutine mld_d_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
character(len=*), intent(in), optional :: prefix, head character(len=*), intent(in), optional :: prefix, head
logical, optional, intent(in) :: ac, rp, smoother, solver, tprol, global_num logical, optional, intent(in) :: ac, rp, smoother, solver, tprol, global_num
! Local variables ! Local variables
integer(psb_ipk_) :: i, j, il1, iln, lname, lev integer(psb_ipk_) :: i, j, il1, iln, lname, lev, ni
integer(psb_ipk_) :: icontxt,iam, np integer(psb_ipk_) :: icontxt,iam, np
character(len=80) :: prefix_ character(len=80) :: prefix_, frmt
character(len=120) :: fname ! len should be at least 20 more than character(len=1024) :: fname
logical :: ac_, rp_, tprol_, global_num_ logical :: ac_, rp_, tprol_, global_num_
integer(psb_ipk_), allocatable :: ivr(:), ivc(:) integer(psb_ipk_), allocatable :: ivr(:), ivc(:)
@ -68,6 +68,7 @@ subroutine mld_d_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
else else
icontxt = -1 icontxt = -1
iam = -1 iam = -1
np = -1
end if end if
if (present(ac)) then if (present(ac)) then
ac_ = ac ac_ = ac
@ -91,8 +92,13 @@ subroutine mld_d_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
end if end if
lname = len_trim(prefix_) lname = len_trim(prefix_)
fname = trim(prefix_) fname = trim(prefix_)
write(fname(lname+1:lname+5),'(a,i3.3)') '_p',iam
lname = lname + 5 if (np > 0) then
ni = floor(log10(1.0*np)) + 1
write(frmt,'(a,i3.3,a,i3.3,a)') '(a,i',ni,'.',ni,')'
write(fname(lname+1:lname+ni+2),frmt) '_p',iam
lname = lname + ni + 2
end if
if (global_num_) then if (global_num_) then
if (level >= 2) then if (level >= 2) then

16
mlprec/impl/level/mld_s_base_onelev_dump.f90
Unescape Escape View File

@ -47,10 +47,10 @@ subroutine mld_s_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
character(len=*), intent(in), optional :: prefix, head character(len=*), intent(in), optional :: prefix, head
logical, optional, intent(in) :: ac, rp, smoother, solver, tprol, global_num logical, optional, intent(in) :: ac, rp, smoother, solver, tprol, global_num
! Local variables ! Local variables
integer(psb_ipk_) :: i, j, il1, iln, lname, lev integer(psb_ipk_) :: i, j, il1, iln, lname, lev, ni
integer(psb_ipk_) :: icontxt,iam, np integer(psb_ipk_) :: icontxt,iam, np
character(len=80) :: prefix_ character(len=80) :: prefix_, frmt
character(len=120) :: fname ! len should be at least 20 more than character(len=1024) :: fname
logical :: ac_, rp_, tprol_, global_num_ logical :: ac_, rp_, tprol_, global_num_
integer(psb_ipk_), allocatable :: ivr(:), ivc(:) integer(psb_ipk_), allocatable :: ivr(:), ivc(:)
@ -68,6 +68,7 @@ subroutine mld_s_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
else else
icontxt = -1 icontxt = -1
iam = -1 iam = -1
np = -1
end if end if
if (present(ac)) then if (present(ac)) then
ac_ = ac ac_ = ac
@ -91,8 +92,13 @@ subroutine mld_s_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
end if end if
lname = len_trim(prefix_) lname = len_trim(prefix_)
fname = trim(prefix_) fname = trim(prefix_)
write(fname(lname+1:lname+5),'(a,i3.3)') '_p',iam
lname = lname + 5 if (np > 0) then
ni = floor(log10(1.0*np)) + 1
write(frmt,'(a,i3.3,a,i3.3,a)') '(a,i',ni,'.',ni,')'
write(fname(lname+1:lname+ni+2),frmt) '_p',iam
lname = lname + ni + 2
end if
if (global_num_) then if (global_num_) then
if (level >= 2) then if (level >= 2) then

16
mlprec/impl/level/mld_z_base_onelev_dump.f90
Unescape Escape View File

@ -47,10 +47,10 @@ subroutine mld_z_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
character(len=*), intent(in), optional :: prefix, head character(len=*), intent(in), optional :: prefix, head
logical, optional, intent(in) :: ac, rp, smoother, solver, tprol, global_num logical, optional, intent(in) :: ac, rp, smoother, solver, tprol, global_num
! Local variables ! Local variables
integer(psb_ipk_) :: i, j, il1, iln, lname, lev integer(psb_ipk_) :: i, j, il1, iln, lname, lev, ni
integer(psb_ipk_) :: icontxt,iam, np integer(psb_ipk_) :: icontxt,iam, np
character(len=80) :: prefix_ character(len=80) :: prefix_, frmt
character(len=120) :: fname ! len should be at least 20 more than character(len=1024) :: fname
logical :: ac_, rp_, tprol_, global_num_ logical :: ac_, rp_, tprol_, global_num_
integer(psb_ipk_), allocatable :: ivr(:), ivc(:) integer(psb_ipk_), allocatable :: ivr(:), ivc(:)
@ -68,6 +68,7 @@ subroutine mld_z_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
else else
icontxt = -1 icontxt = -1
iam = -1 iam = -1
np = -1
end if end if
if (present(ac)) then if (present(ac)) then
ac_ = ac ac_ = ac
@ -91,8 +92,13 @@ subroutine mld_z_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
end if end if
lname = len_trim(prefix_) lname = len_trim(prefix_)
fname = trim(prefix_) fname = trim(prefix_)
write(fname(lname+1:lname+5),'(a,i3.3)') '_p',iam
lname = lname + 5 if (np > 0) then
ni = floor(log10(1.0*np)) + 1
write(frmt,'(a,i3.3,a,i3.3,a)') '(a,i',ni,'.',ni,')'
write(fname(lname+1:lname+ni+2),frmt) '_p',iam
lname = lname + ni + 2
end if
if (global_num_) then if (global_num_) then
if (level >= 2) then if (level >= 2) then

Write Preview
Loading…
Cancel
Save