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

psblas3:

Browse Source 
base/tools/psb_ccdbldext.F90
 base/tools/psb_dcdbldext.F90
 base/tools/psb_scdbldext.F90
 base/tools/psb_zcdbldext.F90

Allow for empyt input matrices in BLDEXT
psblas-3.4-maint
Salvatore Filippone 10 years ago
parent 992a760fe8
commit 003a1096ce
4 changed files with 8 additions and 20 deletions
Show Stats Download Patch File Download Diff File

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

@ -190,15 +190,12 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)
! LOVR= (NNZ/NROW)*N_HALO*NOVR This assumes that the local average ! LOVR= (NNZ/NROW)*N_HALO*NOVR This assumes that the local average
! nonzeros per row is the same as the global. ! nonzeros per row is the same as the global.
! !
nztot = a%get_nzeros() ! Allow for empty matrices.
nztot = max(ione,a%get_nzeros())
if (nztot>0) then if (nztot>0) then
lovr = ((nztot+m-1)/m)*nhalo*novr lovr = ((nztot+m-1)/m)*nhalo*novr
lworks = ((nztot+m-1)/m)*nhalo lworks = ((nztot+m-1)/m)*nhalo
lworkr = ((nztot+m-1)/m)*nhalo lworkr = ((nztot+m-1)/m)*nhalo
else
info = -1
call psb_errpush(info,name)
goto 9999
endif endif
If (debug_level >= psb_debug_outer_)& If (debug_level >= psb_debug_outer_)&
& Write(debug_unit,*) me,' ',trim(name),':ovr_est done',novr,lovr & Write(debug_unit,*) me,' ',trim(name),':ovr_est done',novr,lovr

7
base/tools/psb_dcdbldext.F90
Unescape Escape View File

@ -190,15 +190,12 @@ Subroutine psb_dcdbldext(a,desc_a,novr,desc_ov,info, extype)
! LOVR= (NNZ/NROW)*N_HALO*NOVR This assumes that the local average ! LOVR= (NNZ/NROW)*N_HALO*NOVR This assumes that the local average
! nonzeros per row is the same as the global. ! nonzeros per row is the same as the global.
! !
nztot = a%get_nzeros() ! Allow for empty matrices.
nztot = max(ione,a%get_nzeros())
if (nztot>0) then if (nztot>0) then
lovr = ((nztot+m-1)/m)*nhalo*novr lovr = ((nztot+m-1)/m)*nhalo*novr
lworks = ((nztot+m-1)/m)*nhalo lworks = ((nztot+m-1)/m)*nhalo
lworkr = ((nztot+m-1)/m)*nhalo lworkr = ((nztot+m-1)/m)*nhalo
else
info = -1
call psb_errpush(info,name)
goto 9999
endif endif
If (debug_level >= psb_debug_outer_)& If (debug_level >= psb_debug_outer_)&
& Write(debug_unit,*) me,' ',trim(name),':ovr_est done',novr,lovr & Write(debug_unit,*) me,' ',trim(name),':ovr_est done',novr,lovr

7
base/tools/psb_scdbldext.F90
Unescape Escape View File

@ -190,15 +190,12 @@ Subroutine psb_scdbldext(a,desc_a,novr,desc_ov,info, extype)
! LOVR= (NNZ/NROW)*N_HALO*NOVR This assumes that the local average ! LOVR= (NNZ/NROW)*N_HALO*NOVR This assumes that the local average
! nonzeros per row is the same as the global. ! nonzeros per row is the same as the global.
! !
nztot = a%get_nzeros() ! Allow for empty matrices.
nztot = max(ione,a%get_nzeros())
if (nztot>0) then if (nztot>0) then
lovr = ((nztot+m-1)/m)*nhalo*novr lovr = ((nztot+m-1)/m)*nhalo*novr
lworks = ((nztot+m-1)/m)*nhalo lworks = ((nztot+m-1)/m)*nhalo
lworkr = ((nztot+m-1)/m)*nhalo lworkr = ((nztot+m-1)/m)*nhalo
else
info = -1
call psb_errpush(info,name)
goto 9999
endif endif
If (debug_level >= psb_debug_outer_)& If (debug_level >= psb_debug_outer_)&
& Write(debug_unit,*) me,' ',trim(name),':ovr_est done',novr,lovr & Write(debug_unit,*) me,' ',trim(name),':ovr_est done',novr,lovr

7
base/tools/psb_zcdbldext.F90
Unescape Escape View File

@ -190,15 +190,12 @@ Subroutine psb_zcdbldext(a,desc_a,novr,desc_ov,info, extype)
! LOVR= (NNZ/NROW)*N_HALO*NOVR This assumes that the local average ! LOVR= (NNZ/NROW)*N_HALO*NOVR This assumes that the local average
! nonzeros per row is the same as the global. ! nonzeros per row is the same as the global.
! !
nztot = a%get_nzeros() ! Allow for empty matrices.
nztot = max(ione,a%get_nzeros())
if (nztot>0) then if (nztot>0) then
lovr = ((nztot+m-1)/m)*nhalo*novr lovr = ((nztot+m-1)/m)*nhalo*novr
lworks = ((nztot+m-1)/m)*nhalo lworks = ((nztot+m-1)/m)*nhalo
lworkr = ((nztot+m-1)/m)*nhalo lworkr = ((nztot+m-1)/m)*nhalo
else
info = -1
call psb_errpush(info,name)
goto 9999
endif endif
If (debug_level >= psb_debug_outer_)& If (debug_level >= psb_debug_outer_)&
& Write(debug_unit,*) me,' ',trim(name),':ovr_est done',novr,lovr & Write(debug_unit,*) me,' ',trim(name),':ovr_est done',novr,lovr

Write Preview
Loading…
Cancel
Save