stopcriterion
			
			
		
Salvatore Filippone 16 years ago
parent 4de41865e8
commit 1ccb5158ed

@ -420,7 +420,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(cone,mlprec_wrk(ilev-1)%x2l,&
call psb_map_X2Y(cone,mlprec_wrk(ilev-1)%x2l,&
& czero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -451,7 +451,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(cone,mlprec_wrk(ilev)%y2l,&
call psb_map_Y2X(cone,mlprec_wrk(ilev)%y2l,&
& cone,mlprec_wrk(ilev-1)%y2l,&
& precv(ilev)%map_desc,info,work=work)
@ -681,7 +681,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(cone,mlprec_wrk(ilev-1)%tx,&
call psb_map_X2Y(cone,mlprec_wrk(ilev-1)%tx,&
& czero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -720,7 +720,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(cone,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(cone,mlprec_wrk(ilev+1)%y2l,&
& cone,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)
@ -927,7 +927,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(cone,mlprec_wrk(ilev-1)%x2l,&
call psb_map_X2Y(cone,mlprec_wrk(ilev-1)%x2l,&
& czero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -981,7 +981,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(cone,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(cone,mlprec_wrk(ilev+1)%y2l,&
& czero,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)
@ -1238,7 +1238,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(cone,mlprec_wrk(ilev-1)%ty,&
call psb_map_X2Y(cone,mlprec_wrk(ilev-1)%ty,&
& czero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -1281,7 +1281,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(cone,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(cone,mlprec_wrk(ilev+1)%y2l,&
& cone,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)

@ -268,8 +268,8 @@ subroutine mld_cprecbld(a,desc_a,p,info)
do i=2, iszv - 1
p%precv(i)%base_a => p%precv(i)%ac
p%precv(i)%base_desc => p%precv(i)%desc_ac
p%precv(i)%map_desc%p_desc_fw => p%precv(i-1)%base_desc
p%precv(i)%map_desc%p_desc_bk => p%precv(i)%base_desc
p%precv(i)%map_desc%p_desc_X => p%precv(i-1)%base_desc
p%precv(i)%map_desc%p_desc_Y => p%precv(i)%base_desc
end do

@ -419,7 +419,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(done,mlprec_wrk(ilev-1)%x2l,&
call psb_map_X2Y(done,mlprec_wrk(ilev-1)%x2l,&
& dzero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -450,7 +450,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(done,mlprec_wrk(ilev)%y2l,&
call psb_map_Y2X(done,mlprec_wrk(ilev)%y2l,&
& done,mlprec_wrk(ilev-1)%y2l,&
& precv(ilev)%map_desc,info,work=work)
@ -679,7 +679,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(done,mlprec_wrk(ilev-1)%tx,&
call psb_map_X2Y(done,mlprec_wrk(ilev-1)%tx,&
& dzero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -718,7 +718,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(done,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(done,mlprec_wrk(ilev+1)%y2l,&
& done,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)
@ -924,7 +924,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(done,mlprec_wrk(ilev-1)%x2l,&
call psb_map_X2Y(done,mlprec_wrk(ilev-1)%x2l,&
& dzero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -978,7 +978,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(done,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(done,mlprec_wrk(ilev+1)%y2l,&
& dzero,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)
@ -1234,7 +1234,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(done,mlprec_wrk(ilev-1)%ty,&
call psb_map_X2Y(done,mlprec_wrk(ilev-1)%ty,&
& dzero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -1277,7 +1277,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(done,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(done,mlprec_wrk(ilev+1)%y2l,&
& done,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)

@ -268,8 +268,8 @@ subroutine mld_dprecbld(a,desc_a,p,info)
do i=2, iszv - 1
p%precv(i)%base_a => p%precv(i)%ac
p%precv(i)%base_desc => p%precv(i)%desc_ac
p%precv(i)%map_desc%p_desc_fw => p%precv(i-1)%base_desc
p%precv(i)%map_desc%p_desc_bk => p%precv(i)%base_desc
p%precv(i)%map_desc%p_desc_X => p%precv(i-1)%base_desc
p%precv(i)%map_desc%p_desc_Y => p%precv(i)%base_desc
end do

@ -419,7 +419,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(sone,mlprec_wrk(ilev-1)%x2l,&
call psb_map_X2Y(sone,mlprec_wrk(ilev-1)%x2l,&
& szero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -450,7 +450,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(sone,mlprec_wrk(ilev)%y2l,&
call psb_map_Y2X(sone,mlprec_wrk(ilev)%y2l,&
& sone,mlprec_wrk(ilev-1)%y2l,&
& precv(ilev)%map_desc,info,work=work)
@ -679,7 +679,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(sone,mlprec_wrk(ilev-1)%tx,&
call psb_map_X2Y(sone,mlprec_wrk(ilev-1)%tx,&
& szero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -718,7 +718,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(sone,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(sone,mlprec_wrk(ilev+1)%y2l,&
& sone,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)
@ -924,7 +924,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(sone,mlprec_wrk(ilev-1)%x2l,&
call psb_map_X2Y(sone,mlprec_wrk(ilev-1)%x2l,&
& szero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -978,7 +978,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(sone,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(sone,mlprec_wrk(ilev+1)%y2l,&
& szero,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)
@ -1234,7 +1234,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(sone,mlprec_wrk(ilev-1)%ty,&
call psb_map_X2Y(sone,mlprec_wrk(ilev-1)%ty,&
& szero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -1277,7 +1277,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(sone,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(sone,mlprec_wrk(ilev+1)%y2l,&
& sone,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)

@ -268,8 +268,8 @@ subroutine mld_sprecbld(a,desc_a,p,info)
do i=2, iszv - 1
p%precv(i)%base_a => p%precv(i)%ac
p%precv(i)%base_desc => p%precv(i)%desc_ac
p%precv(i)%map_desc%p_desc_fw => p%precv(i-1)%base_desc
p%precv(i)%map_desc%p_desc_bk => p%precv(i)%base_desc
p%precv(i)%map_desc%p_desc_X => p%precv(i-1)%base_desc
p%precv(i)%map_desc%p_desc_Y => p%precv(i)%base_desc
end do

@ -162,8 +162,8 @@ contains
do i=2, isz
b%precv(i)%base_a => b%precv(i)%ac
b%precv(i)%base_desc => b%precv(i)%desc_ac
b%precv(i)%map_desc%p_desc_fw => b%precv(i-1)%base_desc
b%precv(i)%map_desc%p_desc_bk => b%precv(i)%base_desc
b%precv(i)%map_desc%p_desc_X => b%precv(i-1)%base_desc
b%precv(i)%map_desc%p_desc_Y => b%precv(i)%base_desc
end do
end subroutine mld_sprec_transfer
@ -256,8 +256,8 @@ contains
do i=2, isz
b%precv(i)%base_a => b%precv(i)%ac
b%precv(i)%base_desc => b%precv(i)%desc_ac
b%precv(i)%map_desc%p_desc_fw => b%precv(i-1)%base_desc
b%precv(i)%map_desc%p_desc_bk => b%precv(i)%base_desc
b%precv(i)%map_desc%p_desc_X => b%precv(i-1)%base_desc
b%precv(i)%map_desc%p_desc_Y => b%precv(i)%base_desc
end do
end subroutine mld_dprec_transfer
@ -350,8 +350,8 @@ contains
do i=2, isz
b%precv(i)%base_a => b%precv(i)%ac
b%precv(i)%base_desc => b%precv(i)%desc_ac
b%precv(i)%map_desc%p_desc_fw => b%precv(i-1)%base_desc
b%precv(i)%map_desc%p_desc_bk => b%precv(i)%base_desc
b%precv(i)%map_desc%p_desc_X => b%precv(i-1)%base_desc
b%precv(i)%map_desc%p_desc_Y => b%precv(i)%base_desc
end do
end subroutine mld_cprec_transfer
@ -444,8 +444,8 @@ contains
do i=2, isz
b%precv(i)%base_a => b%precv(i)%ac
b%precv(i)%base_desc => b%precv(i)%desc_ac
b%precv(i)%map_desc%p_desc_fw => b%precv(i-1)%base_desc
b%precv(i)%map_desc%p_desc_bk => b%precv(i)%base_desc
b%precv(i)%map_desc%p_desc_X => b%precv(i-1)%base_desc
b%precv(i)%map_desc%p_desc_Y => b%precv(i)%base_desc
end do
end subroutine mld_zprec_transfer

@ -420,7 +420,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(zone,mlprec_wrk(ilev-1)%x2l,&
call psb_map_X2Y(zone,mlprec_wrk(ilev-1)%x2l,&
& zzero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -451,7 +451,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(zone,mlprec_wrk(ilev)%y2l,&
call psb_map_Y2X(zone,mlprec_wrk(ilev)%y2l,&
& zone,mlprec_wrk(ilev-1)%y2l,&
& precv(ilev)%map_desc,info,work=work)
@ -681,7 +681,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(zone,mlprec_wrk(ilev-1)%tx,&
call psb_map_X2Y(zone,mlprec_wrk(ilev-1)%tx,&
& zzero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -720,7 +720,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(zone,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(zone,mlprec_wrk(ilev+1)%y2l,&
& zone,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)
@ -927,7 +927,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(zone,mlprec_wrk(ilev-1)%x2l,&
call psb_map_X2Y(zone,mlprec_wrk(ilev-1)%x2l,&
& zzero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -981,7 +981,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(zone,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(zone,mlprec_wrk(ilev+1)%y2l,&
& zzero,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)
@ -1238,7 +1238,7 @@ contains
end if
! Apply prolongator transpose, i.e. restriction
call psb_forward_map(zone,mlprec_wrk(ilev-1)%ty,&
call psb_map_X2Y(zone,mlprec_wrk(ilev-1)%ty,&
& zzero,mlprec_wrk(ilev)%x2l,&
& precv(ilev)%map_desc,info,work=work)
@ -1281,7 +1281,7 @@ contains
!
! Apply prolongator
!
call psb_backward_map(zone,mlprec_wrk(ilev+1)%y2l,&
call psb_map_Y2X(zone,mlprec_wrk(ilev+1)%y2l,&
& zone,mlprec_wrk(ilev)%y2l,&
& precv(ilev+1)%map_desc,info,work=work)

@ -268,8 +268,8 @@ subroutine mld_zprecbld(a,desc_a,p,info)
do i=2, iszv - 1
p%precv(i)%base_a => p%precv(i)%ac
p%precv(i)%base_desc => p%precv(i)%desc_ac
p%precv(i)%map_desc%p_desc_fw => p%precv(i-1)%base_desc
p%precv(i)%map_desc%p_desc_bk => p%precv(i)%base_desc
p%precv(i)%map_desc%p_desc_X => p%precv(i-1)%base_desc
p%precv(i)%map_desc%p_desc_Y => p%precv(i)%base_desc
end do

Loading…
Cancel
Save