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/serial/impl/psb_c_csr_impl.f90
 base/serial/impl/psb_d_csr_impl.f90
 base/serial/impl/psb_s_csr_impl.f90
 base/serial/impl/psb_z_csr_impl.f90

Fix computation of output space for getrow & friends.
psblas-3.2.0
Salvatore Filippone 11 years ago
parent 82d292c8f6
commit fd82da8d8a
4 changed files with 24 additions and 8 deletions
Show Stats Download Patch File Download Diff File

8
base/serial/impl/psb_c_csr_impl.f90
Unescape Escape View File

@ -2029,7 +2029,7 @@ contains
logical, intent(in) :: append logical, intent(in) :: append
integer(psb_ipk_) :: info integer(psb_ipk_) :: info
integer(psb_ipk_), optional :: iren(:) integer(psb_ipk_), optional :: iren(:)
integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, nrd,ncd integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, icl,lcl,nrd,ncd
integer(psb_ipk_) :: debug_level, debug_unit integer(psb_ipk_) :: debug_level, debug_unit
character(len=20) :: name='csr_getptn' character(len=20) :: name='csr_getptn'
@ -2039,6 +2039,8 @@ contains
nza = a%get_nzeros() nza = a%get_nzeros()
irw = imin irw = imin
lrw = min(imax,a%get_nrows()) lrw = min(imax,a%get_nrows())
icl = jmin
lcl = min(jmax,a%get_ncols())
if (irw<0) then if (irw<0) then
info = psb_err_pivot_too_small_ info = psb_err_pivot_too_small_
return return
@ -2209,7 +2211,7 @@ contains
logical, intent(in) :: append logical, intent(in) :: append
integer(psb_ipk_) :: info integer(psb_ipk_) :: info
integer(psb_ipk_), optional :: iren(:) integer(psb_ipk_), optional :: iren(:)
integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, nrd, ncd integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, icl,lcl, nrd, ncd
integer(psb_ipk_) :: debug_level, debug_unit integer(psb_ipk_) :: debug_level, debug_unit
character(len=20) :: name='coo_getrow' character(len=20) :: name='coo_getrow'
@ -2219,6 +2221,8 @@ contains
nza = a%get_nzeros() nza = a%get_nzeros()
irw = imin irw = imin
lrw = min(imax,a%get_nrows()) lrw = min(imax,a%get_nrows())
icl = jmin
lcl = min(jmax,a%get_ncols())
if (irw<0) then if (irw<0) then
info = psb_err_pivot_too_small_ info = psb_err_pivot_too_small_
return return

8
base/serial/impl/psb_d_csr_impl.f90
Unescape Escape View File

@ -2029,7 +2029,7 @@ contains
logical, intent(in) :: append logical, intent(in) :: append
integer(psb_ipk_) :: info integer(psb_ipk_) :: info
integer(psb_ipk_), optional :: iren(:) integer(psb_ipk_), optional :: iren(:)
integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, nrd,ncd integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, icl,lcl,nrd,ncd
integer(psb_ipk_) :: debug_level, debug_unit integer(psb_ipk_) :: debug_level, debug_unit
character(len=20) :: name='csr_getptn' character(len=20) :: name='csr_getptn'
@ -2039,6 +2039,8 @@ contains
nza = a%get_nzeros() nza = a%get_nzeros()
irw = imin irw = imin
lrw = min(imax,a%get_nrows()) lrw = min(imax,a%get_nrows())
icl = jmin
lcl = min(jmax,a%get_ncols())
if (irw<0) then if (irw<0) then
info = psb_err_pivot_too_small_ info = psb_err_pivot_too_small_
return return
@ -2209,7 +2211,7 @@ contains
logical, intent(in) :: append logical, intent(in) :: append
integer(psb_ipk_) :: info integer(psb_ipk_) :: info
integer(psb_ipk_), optional :: iren(:) integer(psb_ipk_), optional :: iren(:)
integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, nrd, ncd integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, icl,lcl, nrd, ncd
integer(psb_ipk_) :: debug_level, debug_unit integer(psb_ipk_) :: debug_level, debug_unit
character(len=20) :: name='coo_getrow' character(len=20) :: name='coo_getrow'
@ -2219,6 +2221,8 @@ contains
nza = a%get_nzeros() nza = a%get_nzeros()
irw = imin irw = imin
lrw = min(imax,a%get_nrows()) lrw = min(imax,a%get_nrows())
icl = jmin
lcl = min(jmax,a%get_ncols())
if (irw<0) then if (irw<0) then
info = psb_err_pivot_too_small_ info = psb_err_pivot_too_small_
return return

8
base/serial/impl/psb_s_csr_impl.f90
Unescape Escape View File

@ -2029,7 +2029,7 @@ contains
logical, intent(in) :: append logical, intent(in) :: append
integer(psb_ipk_) :: info integer(psb_ipk_) :: info
integer(psb_ipk_), optional :: iren(:) integer(psb_ipk_), optional :: iren(:)
integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, nrd,ncd integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, icl,lcl,nrd,ncd
integer(psb_ipk_) :: debug_level, debug_unit integer(psb_ipk_) :: debug_level, debug_unit
character(len=20) :: name='csr_getptn' character(len=20) :: name='csr_getptn'
@ -2039,6 +2039,8 @@ contains
nza = a%get_nzeros() nza = a%get_nzeros()
irw = imin irw = imin
lrw = min(imax,a%get_nrows()) lrw = min(imax,a%get_nrows())
icl = jmin
lcl = min(jmax,a%get_ncols())
if (irw<0) then if (irw<0) then
info = psb_err_pivot_too_small_ info = psb_err_pivot_too_small_
return return
@ -2209,7 +2211,7 @@ contains
logical, intent(in) :: append logical, intent(in) :: append
integer(psb_ipk_) :: info integer(psb_ipk_) :: info
integer(psb_ipk_), optional :: iren(:) integer(psb_ipk_), optional :: iren(:)
integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, nrd, ncd integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, icl,lcl, nrd, ncd
integer(psb_ipk_) :: debug_level, debug_unit integer(psb_ipk_) :: debug_level, debug_unit
character(len=20) :: name='coo_getrow' character(len=20) :: name='coo_getrow'
@ -2219,6 +2221,8 @@ contains
nza = a%get_nzeros() nza = a%get_nzeros()
irw = imin irw = imin
lrw = min(imax,a%get_nrows()) lrw = min(imax,a%get_nrows())
icl = jmin
lcl = min(jmax,a%get_ncols())
if (irw<0) then if (irw<0) then
info = psb_err_pivot_too_small_ info = psb_err_pivot_too_small_
return return

8
base/serial/impl/psb_z_csr_impl.f90
Unescape Escape View File

@ -2029,7 +2029,7 @@ contains
logical, intent(in) :: append logical, intent(in) :: append
integer(psb_ipk_) :: info integer(psb_ipk_) :: info
integer(psb_ipk_), optional :: iren(:) integer(psb_ipk_), optional :: iren(:)
integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, nrd,ncd integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, icl,lcl,nrd,ncd
integer(psb_ipk_) :: debug_level, debug_unit integer(psb_ipk_) :: debug_level, debug_unit
character(len=20) :: name='csr_getptn' character(len=20) :: name='csr_getptn'
@ -2039,6 +2039,8 @@ contains
nza = a%get_nzeros() nza = a%get_nzeros()
irw = imin irw = imin
lrw = min(imax,a%get_nrows()) lrw = min(imax,a%get_nrows())
icl = jmin
lcl = min(jmax,a%get_ncols())
if (irw<0) then if (irw<0) then
info = psb_err_pivot_too_small_ info = psb_err_pivot_too_small_
return return
@ -2209,7 +2211,7 @@ contains
logical, intent(in) :: append logical, intent(in) :: append
integer(psb_ipk_) :: info integer(psb_ipk_) :: info
integer(psb_ipk_), optional :: iren(:) integer(psb_ipk_), optional :: iren(:)
integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, nrd, ncd integer(psb_ipk_) :: nzin_, nza, idx,i,j,k, nzt, irw, lrw, icl,lcl, nrd, ncd
integer(psb_ipk_) :: debug_level, debug_unit integer(psb_ipk_) :: debug_level, debug_unit
character(len=20) :: name='coo_getrow' character(len=20) :: name='coo_getrow'
@ -2219,6 +2221,8 @@ contains
nza = a%get_nzeros() nza = a%get_nzeros()
irw = imin irw = imin
lrw = min(imax,a%get_nrows()) lrw = min(imax,a%get_nrows())
icl = jmin
lcl = min(jmax,a%get_ncols())
if (irw<0) then if (irw<0) then
info = psb_err_pivot_too_small_ info = psb_err_pivot_too_small_
return return

Write Preview
Loading…
Cancel
Save