mld2p4-2:

mlprec/impl/mld_cmlprec_bld.f90
 mlprec/impl/mld_dmlprec_bld.f90
 mlprec/impl/mld_smlprec_bld.f90
 mlprec/impl/mld_zmlprec_bld.f90

Comment update.
stopcriterion
Salvatore Filippone 12 years ago
parent 43bb186dcd
commit 534adf7142

@ -289,8 +289,11 @@ subroutine mld_cmlprec_bld(a,desc_a,p,info,amold,vmold)
newnode => head newnode => head
do i=1, newsz do i=1, newsz
current => newnode current => newnode
! First do a move_alloc.
! This handles the AC, DESC_AC and MAP fields
if (info == psb_success_) & if (info == psb_success_) &
& call mld_move_alloc(current%item,p%precv(i),info) & call mld_move_alloc(current%item,p%precv(i),info)
! Now set the smoother/solver parts.
if (info == psb_success_) then if (info == psb_success_) then
if (i ==1) then if (i ==1) then
! This is a workaround for a bug in gfortran 4.7.2 ! This is a workaround for a bug in gfortran 4.7.2

@ -289,8 +289,11 @@ subroutine mld_dmlprec_bld(a,desc_a,p,info,amold,vmold)
newnode => head newnode => head
do i=1, newsz do i=1, newsz
current => newnode current => newnode
! First do a move_alloc.
! This handles the AC, DESC_AC and MAP fields
if (info == psb_success_) & if (info == psb_success_) &
& call mld_move_alloc(current%item,p%precv(i),info) & call mld_move_alloc(current%item,p%precv(i),info)
! Now set the smoother/solver parts.
if (info == psb_success_) then if (info == psb_success_) then
if (i ==1) then if (i ==1) then
! This is a workaround for a bug in gfortran 4.7.2 ! This is a workaround for a bug in gfortran 4.7.2

@ -289,8 +289,11 @@ subroutine mld_smlprec_bld(a,desc_a,p,info,amold,vmold)
newnode => head newnode => head
do i=1, newsz do i=1, newsz
current => newnode current => newnode
! First do a move_alloc.
! This handles the AC, DESC_AC and MAP fields
if (info == psb_success_) & if (info == psb_success_) &
& call mld_move_alloc(current%item,p%precv(i),info) & call mld_move_alloc(current%item,p%precv(i),info)
! Now set the smoother/solver parts.
if (info == psb_success_) then if (info == psb_success_) then
if (i ==1) then if (i ==1) then
! This is a workaround for a bug in gfortran 4.7.2 ! This is a workaround for a bug in gfortran 4.7.2

@ -289,8 +289,11 @@ subroutine mld_zmlprec_bld(a,desc_a,p,info,amold,vmold)
newnode => head newnode => head
do i=1, newsz do i=1, newsz
current => newnode current => newnode
! First do a move_alloc.
! This handles the AC, DESC_AC and MAP fields
if (info == psb_success_) & if (info == psb_success_) &
& call mld_move_alloc(current%item,p%precv(i),info) & call mld_move_alloc(current%item,p%precv(i),info)
! Now set the smoother/solver parts.
if (info == psb_success_) then if (info == psb_success_) then
if (i ==1) then if (i ==1) then
! This is a workaround for a bug in gfortran 4.7.2 ! This is a workaround for a bug in gfortran 4.7.2

Loading…
Cancel
Save