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

mld2p4-2:

Browse Source 
mlprec/mld_c_move_alloc_mod.f90
 mlprec/mld_c_prec_type.f90
 mlprec/mld_cmlprec_bld.f90
 mlprec/mld_d_move_alloc_mod.f90
 mlprec/mld_d_prec_type.f90
 mlprec/mld_s_move_alloc_mod.f90
 mlprec/mld_s_prec_type.f90
 mlprec/mld_smlprec_bld.f90
 mlprec/mld_z_move_alloc_mod.f90
 mlprec/mld_z_prec_type.f90
 mlprec/mld_zmlprec_bld.f90

Fixed internal docs, also preprocessed.
stopcriterion
Salvatore Filippone 13 years ago
parent ffe4fe1f0c
commit 089893d3f6
11 changed files with 1507 additions and 1041 deletions
Show Stats Download Patch File Download Diff File

6
mlprec/mld_c_move_alloc_mod.f90
Unescape Escape View File

@ -36,9 +36,9 @@
!!$ POSSIBILITY OF SUCH DAMAGE. !!$ POSSIBILITY OF SUCH DAMAGE.
!!$ !!$
!!$ !!$
! File: mld_move_alloc_mod.f90 ! File: mld_c_move_alloc_mod.f90
! !
! Module: mld_move_alloc_mod ! Module: mld_c_move_alloc_mod
! !
! This module defines move_alloc-like routines, and related interfaces, ! This module defines move_alloc-like routines, and related interfaces,
! for the preconditioner data structures. . ! for the preconditioner data structures. .
@ -61,7 +61,7 @@ contains
type(mld_conelev_type), intent(inout) :: a, b type(mld_conelev_type), intent(inout) :: a, b
integer, intent(out) :: info integer, intent(out) :: info
call mld_precfree(b,info) call b%free(info)
call move_alloc(a%sm,b%sm) call move_alloc(a%sm,b%sm)
if (info == psb_success_) call psb_move_alloc(a%ac,b%ac,info) if (info == psb_success_) call psb_move_alloc(a%ac,b%ac,info)
if (info == psb_success_) call psb_move_alloc(a%desc_ac,b%desc_ac,info) if (info == psb_success_) call psb_move_alloc(a%desc_ac,b%desc_ac,info)

834
mlprec/mld_c_prec_type.f90
View File

File diff suppressed because it is too large Load Diff

12
mlprec/mld_cmlprec_bld.f90
Unescape Escape View File

@ -86,7 +86,7 @@ subroutine mld_cmlprec_bld(a,desc_a,p,info,amold,vmold)
integer :: ipv(mld_ifpsz_), val integer :: ipv(mld_ifpsz_), val
integer :: int_err(5) integer :: int_err(5)
character :: upd_ character :: upd_
type(mld_dml_parms) :: prm type(mld_sml_parms) :: prm
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name, ch_err character(len=20) :: name, ch_err
@ -249,7 +249,7 @@ subroutine mld_cmlprec_bld(a,desc_a,p,info,amold,vmold)
end do end do
call mld_move_alloc(p%precv(iszv),t_prec%precv(newsz),info) call mld_move_alloc(p%precv(iszv),t_prec%precv(newsz),info)
do i=newsz+1, iszv do i=newsz+1, iszv
call mld_precfree(p%precv(i),info) call p%precv(i)%free(info)
end do end do
call mld_move_alloc(t_prec,p,info) call mld_move_alloc(t_prec,p,info)
! Ignore errors from transfer ! Ignore errors from transfer
@ -260,8 +260,8 @@ subroutine mld_cmlprec_bld(a,desc_a,p,info,amold,vmold)
! Fix the pointers, but the level 1 should ! Fix the pointers, but the level 1 should
! be already OK ! be already OK
do i=2, iszv - 1 do i=2, iszv - 1
p%precv(i)%base_a => p%precv(i)%ac p%precv(i)%base_a => p%precv(i)%ac
p%precv(i)%base_desc => p%precv(i)%desc_ac p%precv(i)%base_desc => p%precv(i)%desc_ac
p%precv(i)%map%p_desc_X => p%precv(i-1)%base_desc p%precv(i)%map%p_desc_X => p%precv(i-1)%base_desc
p%precv(i)%map%p_desc_Y => p%precv(i)%base_desc p%precv(i)%map%p_desc_Y => p%precv(i)%base_desc
end do end do
@ -293,13 +293,13 @@ subroutine mld_cmlprec_bld(a,desc_a,p,info,amold,vmold)
& 'Jacobi sweeps',1,is_legal_jac_sweeps) & 'Jacobi sweeps',1,is_legal_jac_sweeps)
if (.not.allocated(p%precv(i)%sm)) then if (.not.allocated(p%precv(i)%sm)) then
!! Error: should have called mld_cprecinit !! Error: should have called mld_dprecinit
info=3111 info=3111
call psb_errpush(info,name) call psb_errpush(info,name)
goto 9999 goto 9999
end if end if
if (.not.allocated(p%precv(i)%sm%sv)) then if (.not.allocated(p%precv(i)%sm%sv)) then
!! Error: should have called mld_cprecinit !! Error: should have called mld_dprecinit
info=3111 info=3111
call psb_errpush(info,name) call psb_errpush(info,name)
goto 9999 goto 9999

4
mlprec/mld_d_move_alloc_mod.f90
Unescape Escape View File

@ -36,9 +36,9 @@
!!$ POSSIBILITY OF SUCH DAMAGE. !!$ POSSIBILITY OF SUCH DAMAGE.
!!$ !!$
!!$ !!$
! File: mld_move_alloc_mod.f90 ! File: mld_d_move_alloc_mod.f90
! !
! Module: mld_move_alloc_mod ! Module: mld_d_move_alloc_mod
! !
! This module defines move_alloc-like routines, and related interfaces, ! This module defines move_alloc-like routines, and related interfaces,
! for the preconditioner data structures. . ! for the preconditioner data structures. .

11
mlprec/mld_d_prec_type.f90
Unescape Escape View File

@ -36,17 +36,16 @@
!!$ POSSIBILITY OF SUCH DAMAGE. !!$ POSSIBILITY OF SUCH DAMAGE.
!!$ !!$
!!$ !!$
! File: mld_prec_type.f90 ! File: mld_d_prec_type.f90
! !
! Module: mld_prec_type ! Module: mld_d_prec_type
! !
! This module defines: ! This module defines:
! - the mld_prec_type data structure containing the preconditioner and related ! - the mld_d_prec_type data structure containing the preconditioner and related
! data structures; ! data structures;
! !
! It contains routines for ! It contains routines for
! - converting character constants defining the preconditioner into integer ! - Building and applying;
! constants;
! - checking if the preconditioner is correctly defined; ! - checking if the preconditioner is correctly defined;
! - printing a description of the preconditioner; ! - printing a description of the preconditioner;
! - deallocating the preconditioner data structure. ! - deallocating the preconditioner data structure.
@ -644,7 +643,7 @@ contains
! !
! Subroutines: mld_Tbase_precfree, mld_T_onelev_precfree, mld_Tprec_free ! Subroutines: mld_Tbase_precfree, mld_T_onelev_precfree, mld_Tprec_free
! Version: real/complex ! Version: real
! !
! These routines deallocate the mld_Tbaseprec_type, mld_Tonelev_type and ! These routines deallocate the mld_Tbaseprec_type, mld_Tonelev_type and
! mld_Tprec_type data structures. ! mld_Tprec_type data structures.

6
mlprec/mld_s_move_alloc_mod.f90
Unescape Escape View File

@ -36,9 +36,9 @@
!!$ POSSIBILITY OF SUCH DAMAGE. !!$ POSSIBILITY OF SUCH DAMAGE.
!!$ !!$
!!$ !!$
! File: mld_move_alloc_mod.f90 ! File: mld_s_move_alloc_mod.f90
! !
! Module: mld_move_alloc_mod ! Module: mld_s_move_alloc_mod
! !
! This module defines move_alloc-like routines, and related interfaces, ! This module defines move_alloc-like routines, and related interfaces,
! for the preconditioner data structures. . ! for the preconditioner data structures. .
@ -61,7 +61,7 @@ contains
type(mld_sonelev_type), intent(inout) :: a, b type(mld_sonelev_type), intent(inout) :: a, b
integer, intent(out) :: info integer, intent(out) :: info
call mld_precfree(b,info) call b%free(info)
call move_alloc(a%sm,b%sm) call move_alloc(a%sm,b%sm)
if (info == psb_success_) call psb_move_alloc(a%ac,b%ac,info) if (info == psb_success_) call psb_move_alloc(a%ac,b%ac,info)
if (info == psb_success_) call psb_move_alloc(a%desc_ac,b%desc_ac,info) if (info == psb_success_) call psb_move_alloc(a%desc_ac,b%desc_ac,info)

823
mlprec/mld_s_prec_type.f90
View File

File diff suppressed because it is too large Load Diff

10
mlprec/mld_smlprec_bld.f90
Unescape Escape View File

@ -249,7 +249,7 @@ subroutine mld_smlprec_bld(a,desc_a,p,info,amold,vmold)
end do end do
call mld_move_alloc(p%precv(iszv),t_prec%precv(newsz),info) call mld_move_alloc(p%precv(iszv),t_prec%precv(newsz),info)
do i=newsz+1, iszv do i=newsz+1, iszv
call mld_precfree(p%precv(i),info) call p%precv(i)%free(info)
end do end do
call mld_move_alloc(t_prec,p,info) call mld_move_alloc(t_prec,p,info)
! Ignore errors from transfer ! Ignore errors from transfer
@ -260,8 +260,8 @@ subroutine mld_smlprec_bld(a,desc_a,p,info,amold,vmold)
! Fix the pointers, but the level 1 should ! Fix the pointers, but the level 1 should
! be already OK ! be already OK
do i=2, iszv - 1 do i=2, iszv - 1
p%precv(i)%base_a => p%precv(i)%ac p%precv(i)%base_a => p%precv(i)%ac
p%precv(i)%base_desc => p%precv(i)%desc_ac p%precv(i)%base_desc => p%precv(i)%desc_ac
p%precv(i)%map%p_desc_X => p%precv(i-1)%base_desc p%precv(i)%map%p_desc_X => p%precv(i-1)%base_desc
p%precv(i)%map%p_desc_Y => p%precv(i)%base_desc p%precv(i)%map%p_desc_Y => p%precv(i)%base_desc
end do end do
@ -293,13 +293,13 @@ subroutine mld_smlprec_bld(a,desc_a,p,info,amold,vmold)
& 'Jacobi sweeps',1,is_legal_jac_sweeps) & 'Jacobi sweeps',1,is_legal_jac_sweeps)
if (.not.allocated(p%precv(i)%sm)) then if (.not.allocated(p%precv(i)%sm)) then
!! Error: should have called mld_sprecinit !! Error: should have called mld_dprecinit
info=3111 info=3111
call psb_errpush(info,name) call psb_errpush(info,name)
goto 9999 goto 9999
end if end if
if (.not.allocated(p%precv(i)%sm%sv)) then if (.not.allocated(p%precv(i)%sm%sv)) then
!! Error: should have called mld_sprecinit !! Error: should have called mld_dprecinit
info=3111 info=3111
call psb_errpush(info,name) call psb_errpush(info,name)
goto 9999 goto 9999

6
mlprec/mld_z_move_alloc_mod.f90
Unescape Escape View File

@ -36,9 +36,9 @@
!!$ POSSIBILITY OF SUCH DAMAGE. !!$ POSSIBILITY OF SUCH DAMAGE.
!!$ !!$
!!$ !!$
! File: mld_move_alloc_mod.f90 ! File: mld_z_move_alloc_mod.f90
! !
! Module: mld_move_alloc_mod ! Module: mld_z_move_alloc_mod
! !
! This module defines move_alloc-like routines, and related interfaces, ! This module defines move_alloc-like routines, and related interfaces,
! for the preconditioner data structures. . ! for the preconditioner data structures. .
@ -61,7 +61,7 @@ contains
type(mld_zonelev_type), intent(inout) :: a, b type(mld_zonelev_type), intent(inout) :: a, b
integer, intent(out) :: info integer, intent(out) :: info
call mld_precfree(b,info) call b%free(info)
call move_alloc(a%sm,b%sm) call move_alloc(a%sm,b%sm)
if (info == psb_success_) call psb_move_alloc(a%ac,b%ac,info) if (info == psb_success_) call psb_move_alloc(a%ac,b%ac,info)
if (info == psb_success_) call psb_move_alloc(a%desc_ac,b%desc_ac,info) if (info == psb_success_) call psb_move_alloc(a%desc_ac,b%desc_ac,info)

832
mlprec/mld_z_prec_type.f90
View File

File diff suppressed because it is too large Load Diff

10
mlprec/mld_zmlprec_bld.f90
Unescape Escape View File

@ -249,7 +249,7 @@ subroutine mld_zmlprec_bld(a,desc_a,p,info,amold,vmold)
end do end do
call mld_move_alloc(p%precv(iszv),t_prec%precv(newsz),info) call mld_move_alloc(p%precv(iszv),t_prec%precv(newsz),info)
do i=newsz+1, iszv do i=newsz+1, iszv
call mld_precfree(p%precv(i),info) call p%precv(i)%free(info)
end do end do
call mld_move_alloc(t_prec,p,info) call mld_move_alloc(t_prec,p,info)
! Ignore errors from transfer ! Ignore errors from transfer
@ -260,8 +260,8 @@ subroutine mld_zmlprec_bld(a,desc_a,p,info,amold,vmold)
! Fix the pointers, but the level 1 should ! Fix the pointers, but the level 1 should
! be already OK ! be already OK
do i=2, iszv - 1 do i=2, iszv - 1
p%precv(i)%base_a => p%precv(i)%ac p%precv(i)%base_a => p%precv(i)%ac
p%precv(i)%base_desc => p%precv(i)%desc_ac p%precv(i)%base_desc => p%precv(i)%desc_ac
p%precv(i)%map%p_desc_X => p%precv(i-1)%base_desc p%precv(i)%map%p_desc_X => p%precv(i-1)%base_desc
p%precv(i)%map%p_desc_Y => p%precv(i)%base_desc p%precv(i)%map%p_desc_Y => p%precv(i)%base_desc
end do end do
@ -293,13 +293,13 @@ subroutine mld_zmlprec_bld(a,desc_a,p,info,amold,vmold)
& 'Jacobi sweeps',1,is_legal_jac_sweeps) & 'Jacobi sweeps',1,is_legal_jac_sweeps)
if (.not.allocated(p%precv(i)%sm)) then if (.not.allocated(p%precv(i)%sm)) then
!! Error: should have called mld_zprecinit !! Error: should have called mld_dprecinit
info=3111 info=3111
call psb_errpush(info,name) call psb_errpush(info,name)
goto 9999 goto 9999
end if end if
if (.not.allocated(p%precv(i)%sm%sv)) then if (.not.allocated(p%precv(i)%sm%sv)) then
!! Error: should have called mld_zprecinit !! Error: should have called mld_dprecinit
info=3111 info=3111
call psb_errpush(info,name) call psb_errpush(info,name)
goto 9999 goto 9999

Write Preview
Loading…
Cancel
Save