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

Added I to L conversion because of GNU PR52162.

Browse Source
ILmat
Salvatore Filippone 8 years ago
parent 7cb828bef1
commit 6faed403ee
4 changed files with 112 additions and 32 deletions
Show Stats Download Patch File Download Diff File

36
base/modules/serial/psb_c_mat_mod.F90
Unescape Escape View File

@ -1314,16 +1314,26 @@ contains
! Local
integer(psb_ipk_), allocatable :: lia(:), lja(:)
info = psb_success_
!
! Note: in principle we could use reallocate on assignment,
! but GCC bug 52162 forces us to take defensive programming.
!
if (allocated(ia)) then
lia = ia
call psb_realloc(size(ia),lia,info)
if (info == psb_success_) lia(:) = ia(:)
end if
if (allocated(ja)) then
lja = ja
call psb_realloc(size(ja),lja,info)
if (info == psb_success_) lja(:) = ja(:)
end if
call a%csget(imin,imax,nz,lia,lja,info,&
& jmin,jmax,iren,append,nzin,rscale,cscale)
ia = lia
ja = lja
call psb_ensure_size(size(lia),ia,info)
if (info == psb_success_) ia(:) = lia(:)
call psb_ensure_size(size(lja),ja,info)
if (info == psb_success_) ja(:) = lja(:)
end subroutine psb_c_lcsgetptn
@ -1343,16 +1353,26 @@ contains
! Local
integer(psb_ipk_), allocatable :: lia(:), lja(:)
!
! Note: in principle we could use reallocate on assignment,
! but GCC bug 52162 forces us to take defensive programming.
!
if (allocated(ia)) then
lia = ia
call psb_realloc(size(ia),lia,info)
if (info == psb_success_) lia(:) = ia(:)
end if
if (allocated(ja)) then
lja = ja
call psb_realloc(size(ja),lja,info)
if (info == psb_success_) lja(:) = ja(:)
end if
call a%csget(imin,imax,nz,lia,lja,val,info,&
& jmin,jmax,iren,append,nzin,rscale,cscale)
ia = lia
ja = lja
call psb_ensure_size(size(lia),ia,info)
if (info == psb_success_) ia(:) = lia(:)
call psb_ensure_size(size(lja),ja,info)
if (info == psb_success_) ja(:) = lja(:)
end subroutine psb_c_lcsgetrow
#endif

36
base/modules/serial/psb_d_mat_mod.F90
Unescape Escape View File

@ -1314,16 +1314,26 @@ contains
! Local
integer(psb_ipk_), allocatable :: lia(:), lja(:)
info = psb_success_
!
! Note: in principle we could use reallocate on assignment,
! but GCC bug 52162 forces us to take defensive programming.
!
if (allocated(ia)) then
lia = ia
call psb_realloc(size(ia),lia,info)
if (info == psb_success_) lia(:) = ia(:)
end if
if (allocated(ja)) then
lja = ja
call psb_realloc(size(ja),lja,info)
if (info == psb_success_) lja(:) = ja(:)
end if
call a%csget(imin,imax,nz,lia,lja,info,&
& jmin,jmax,iren,append,nzin,rscale,cscale)
ia = lia
ja = lja
call psb_ensure_size(size(lia),ia,info)
if (info == psb_success_) ia(:) = lia(:)
call psb_ensure_size(size(lja),ja,info)
if (info == psb_success_) ja(:) = lja(:)
end subroutine psb_d_lcsgetptn
@ -1343,16 +1353,26 @@ contains
! Local
integer(psb_ipk_), allocatable :: lia(:), lja(:)
!
! Note: in principle we could use reallocate on assignment,
! but GCC bug 52162 forces us to take defensive programming.
!
if (allocated(ia)) then
lia = ia
call psb_realloc(size(ia),lia,info)
if (info == psb_success_) lia(:) = ia(:)
end if
if (allocated(ja)) then
lja = ja
call psb_realloc(size(ja),lja,info)
if (info == psb_success_) lja(:) = ja(:)
end if
call a%csget(imin,imax,nz,lia,lja,val,info,&
& jmin,jmax,iren,append,nzin,rscale,cscale)
ia = lia
ja = lja
call psb_ensure_size(size(lia),ia,info)
if (info == psb_success_) ia(:) = lia(:)
call psb_ensure_size(size(lja),ja,info)
if (info == psb_success_) ja(:) = lja(:)
end subroutine psb_d_lcsgetrow
#endif

36
base/modules/serial/psb_s_mat_mod.F90
Unescape Escape View File

@ -1314,16 +1314,26 @@ contains
! Local
integer(psb_ipk_), allocatable :: lia(:), lja(:)
info = psb_success_
!
! Note: in principle we could use reallocate on assignment,
! but GCC bug 52162 forces us to take defensive programming.
!
if (allocated(ia)) then
lia = ia
call psb_realloc(size(ia),lia,info)
if (info == psb_success_) lia(:) = ia(:)
end if
if (allocated(ja)) then
lja = ja
call psb_realloc(size(ja),lja,info)
if (info == psb_success_) lja(:) = ja(:)
end if
call a%csget(imin,imax,nz,lia,lja,info,&
& jmin,jmax,iren,append,nzin,rscale,cscale)
ia = lia
ja = lja
call psb_ensure_size(size(lia),ia,info)
if (info == psb_success_) ia(:) = lia(:)
call psb_ensure_size(size(lja),ja,info)
if (info == psb_success_) ja(:) = lja(:)
end subroutine psb_s_lcsgetptn
@ -1343,16 +1353,26 @@ contains
! Local
integer(psb_ipk_), allocatable :: lia(:), lja(:)
!
! Note: in principle we could use reallocate on assignment,
! but GCC bug 52162 forces us to take defensive programming.
!
if (allocated(ia)) then
lia = ia
call psb_realloc(size(ia),lia,info)
if (info == psb_success_) lia(:) = ia(:)
end if
if (allocated(ja)) then
lja = ja
call psb_realloc(size(ja),lja,info)
if (info == psb_success_) lja(:) = ja(:)
end if
call a%csget(imin,imax,nz,lia,lja,val,info,&
& jmin,jmax,iren,append,nzin,rscale,cscale)
ia = lia
ja = lja
call psb_ensure_size(size(lia),ia,info)
if (info == psb_success_) ia(:) = lia(:)
call psb_ensure_size(size(lja),ja,info)
if (info == psb_success_) ja(:) = lja(:)
end subroutine psb_s_lcsgetrow
#endif

36
base/modules/serial/psb_z_mat_mod.F90
Unescape Escape View File

@ -1314,16 +1314,26 @@ contains
! Local
integer(psb_ipk_), allocatable :: lia(:), lja(:)
info = psb_success_
!
! Note: in principle we could use reallocate on assignment,
! but GCC bug 52162 forces us to take defensive programming.
!
if (allocated(ia)) then
lia = ia
call psb_realloc(size(ia),lia,info)
if (info == psb_success_) lia(:) = ia(:)
end if
if (allocated(ja)) then
lja = ja
call psb_realloc(size(ja),lja,info)
if (info == psb_success_) lja(:) = ja(:)
end if
call a%csget(imin,imax,nz,lia,lja,info,&
& jmin,jmax,iren,append,nzin,rscale,cscale)
ia = lia
ja = lja
call psb_ensure_size(size(lia),ia,info)
if (info == psb_success_) ia(:) = lia(:)
call psb_ensure_size(size(lja),ja,info)
if (info == psb_success_) ja(:) = lja(:)
end subroutine psb_z_lcsgetptn
@ -1343,16 +1353,26 @@ contains
! Local
integer(psb_ipk_), allocatable :: lia(:), lja(:)
!
! Note: in principle we could use reallocate on assignment,
! but GCC bug 52162 forces us to take defensive programming.
!
if (allocated(ia)) then
lia = ia
call psb_realloc(size(ia),lia,info)
if (info == psb_success_) lia(:) = ia(:)
end if
if (allocated(ja)) then
lja = ja
call psb_realloc(size(ja),lja,info)
if (info == psb_success_) lja(:) = ja(:)
end if
call a%csget(imin,imax,nz,lia,lja,val,info,&
& jmin,jmax,iren,append,nzin,rscale,cscale)
ia = lia
ja = lja
call psb_ensure_size(size(lia),ia,info)
if (info == psb_success_) ia(:) = lia(:)
call psb_ensure_size(size(lja),ja,info)
if (info == psb_success_) ja(:) = lja(:)
end subroutine psb_z_lcsgetrow
#endif

Write Preview
Loading…
Cancel
Save