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

Dump of Krylov base

Browse Source
lambdaI
Cirdans-Home 4 years ago
parent 0bda91069e
commit f3412de169
4 changed files with 108 additions and 96 deletions
Show Stats Download Patch File Download Diff File

51
krylov/psb_ckrylovsubspace_mod.F90
Unescape Escape View File

@ -284,6 +284,7 @@ contains
integer(psb_lpk_) :: mglob
character(len=20) :: name, ch_err
complex(psb_spk_), allocatable :: vglobal(:,:)
complex(psb_spk_), allocatable :: w(:)
info = psb_success_
@ -302,30 +303,32 @@ contains
end if
end if
! mglob = desc_a%get_global_rows()
! allocate(vglobal(mglob,kryl%k), stat=info)
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
! do i=1,kryl%k
! call psb_gather(vglobal(:,i),kryl%v(i),desc_a,info)
! end do
! if (iam == psb_root_) then
! call mm_array_write(vglobal,"Krylov basis",info,filename=filename//"_v")
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
! end if
! if (allocated(vglobal)) deallocate(vglobal, stat=info)
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
mglob = desc_a%get_global_rows()
allocate(vglobal(mglob,kryl%k), stat=info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
do i=1,kryl%k
call psb_gather(w,kryl%v(i),desc_a,info)
vglobal(:,i) = w
end do
if (iam == psb_root_) then
call mm_array_write(vglobal,"Krylov basis",info,filename=filename//"_v")
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
end if
if (allocated(vglobal)) deallocate(vglobal, stat=info)
if (allocated(w)) deallocate(w, stat=info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
call psb_erractionrestore(err_act)

51
krylov/psb_dkrylovsubspace_mod.F90
Unescape Escape View File

@ -284,6 +284,7 @@ contains
integer(psb_lpk_) :: mglob
character(len=20) :: name, ch_err
real(psb_dpk_), allocatable :: vglobal(:,:)
real(psb_dpk_), allocatable :: w(:)
info = psb_success_
@ -302,30 +303,32 @@ contains
end if
end if
! mglob = desc_a%get_global_rows()
! allocate(vglobal(mglob,kryl%k), stat=info)
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
! do i=1,kryl%k
! call psb_gather(vglobal(:,i),kryl%v(i),desc_a,info)
! end do
! if (iam == psb_root_) then
! call mm_array_write(vglobal,"Krylov basis",info,filename=filename//"_v")
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
! end if
! if (allocated(vglobal)) deallocate(vglobal, stat=info)
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
mglob = desc_a%get_global_rows()
allocate(vglobal(mglob,kryl%k), stat=info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
do i=1,kryl%k
call psb_gather(w,kryl%v(i),desc_a,info)
vglobal(:,i) = w
end do
if (iam == psb_root_) then
call mm_array_write(vglobal,"Krylov basis",info,filename=filename//"_v")
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
end if
if (allocated(vglobal)) deallocate(vglobal, stat=info)
if (allocated(w)) deallocate(w, stat=info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
call psb_erractionrestore(err_act)

51
krylov/psb_skrylovsubspace_mod.F90
Unescape Escape View File

@ -284,6 +284,7 @@ contains
integer(psb_lpk_) :: mglob
character(len=20) :: name, ch_err
real(psb_spk_), allocatable :: vglobal(:,:)
real(psb_spk_), allocatable :: w(:)
info = psb_success_
@ -302,30 +303,32 @@ contains
end if
end if
! mglob = desc_a%get_global_rows()
! allocate(vglobal(mglob,kryl%k), stat=info)
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
! do i=1,kryl%k
! call psb_gather(vglobal(:,i),kryl%v(i),desc_a,info)
! end do
! if (iam == psb_root_) then
! call mm_array_write(vglobal,"Krylov basis",info,filename=filename//"_v")
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
! end if
! if (allocated(vglobal)) deallocate(vglobal, stat=info)
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
mglob = desc_a%get_global_rows()
allocate(vglobal(mglob,kryl%k), stat=info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
do i=1,kryl%k
call psb_gather(w,kryl%v(i),desc_a,info)
vglobal(:,i) = w
end do
if (iam == psb_root_) then
call mm_array_write(vglobal,"Krylov basis",info,filename=filename//"_v")
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
end if
if (allocated(vglobal)) deallocate(vglobal, stat=info)
if (allocated(w)) deallocate(w, stat=info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
call psb_erractionrestore(err_act)

51
krylov/psb_zkrylovsubspace_mod.F90
Unescape Escape View File

@ -284,6 +284,7 @@ contains
integer(psb_lpk_) :: mglob
character(len=20) :: name, ch_err
complex(psb_dpk_), allocatable :: vglobal(:,:)
complex(psb_dpk_), allocatable :: w(:)
info = psb_success_
@ -302,30 +303,32 @@ contains
end if
end if
! mglob = desc_a%get_global_rows()
! allocate(vglobal(mglob,kryl%k), stat=info)
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
! do i=1,kryl%k
! call psb_gather(vglobal(:,i),kryl%v(i),desc_a,info)
! end do
! if (iam == psb_root_) then
! call mm_array_write(vglobal,"Krylov basis",info,filename=filename//"_v")
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
! end if
! if (allocated(vglobal)) deallocate(vglobal, stat=info)
! if (info /= psb_success_) then
! info=psb_err_from_subroutine_non_
! call psb_errpush(info,name)
! goto 9999
! end if
mglob = desc_a%get_global_rows()
allocate(vglobal(mglob,kryl%k), stat=info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
do i=1,kryl%k
call psb_gather(w,kryl%v(i),desc_a,info)
vglobal(:,i) = w
end do
if (iam == psb_root_) then
call mm_array_write(vglobal,"Krylov basis",info,filename=filename//"_v")
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
end if
if (allocated(vglobal)) deallocate(vglobal, stat=info)
if (allocated(w)) deallocate(w, stat=info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_non_
call psb_errpush(info,name)
goto 9999
end if
call psb_erractionrestore(err_act)

Write Preview
Loading…
Cancel
Save