Fixed optional itmax usage.

psblas3-type-indexed
Salvatore Filippone 19 years ago
parent d283f6ffa8
commit 316d3859a2

@ -186,7 +186,7 @@ subroutine psb_dbicg(a,prec,b,x,eps,desc_a,info,&
!!$
!!$ r0 = b-ax0
!!$
if (itx.ge.itmax) exit restart
if (itx.ge.litmax) exit restart
it = 0
call psb_axpby(one,b,zero,r,desc_a,info)
call psb_spmm(-one,a,x,one,r,desc_a,info,work=aux)
@ -301,7 +301,7 @@ subroutine psb_dbicg(a,prec,b,x,eps,desc_a,info,&
if (rerr<=eps) then
exit restart
end if
if (itx.ge.itmax) exit restart
if (itx.ge.litmax) exit restart
end do iteration
end do restart

@ -182,7 +182,7 @@ Subroutine psb_dcgs(a,prec,b,x,eps,desc_a,info,&
!!$
!!$ r0 = b-ax0
!!$
If (itx.Ge.itmax) Exit restart
If (itx.Ge.litmax) Exit restart
it = 0
Call psb_axpby(one,b,zero,r,desc_a,info)
Call psb_spmm(-one,a,x,one,r,desc_a,info,work=aux)
@ -294,7 +294,7 @@ Subroutine psb_dcgs(a,prec,b,x,eps,desc_a,info,&
If (rerr<=eps) Then
Exit restart
End If
If (itx.Ge.itmax) Exit restart
If (itx.Ge.litmax) Exit restart
End Do iteration
End Do restart

@ -184,7 +184,7 @@ Subroutine psb_dcgstab(a,prec,b,x,eps,desc_a,info,&
!!$
!!$ r0 = b-Ax0
!!$
If (itx >= itmax) Exit restart
If (itx >= litmax) Exit restart
it = 0
Call psb_axpby(one,b,zero,r,desc_a,info)
Call psb_spmm(-one,a,x,one,r,desc_a,info,work=aux)
@ -347,7 +347,7 @@ Subroutine psb_dcgstab(a,prec,b,x,eps,desc_a,info,&
Exit restart
End If
If (itx.Ge.itmax) Exit restart
If (itx.Ge.litmax) Exit restart
End Do iteration
End Do restart

@ -202,7 +202,7 @@ Subroutine psb_dcgstabl(a,prec,b,x,eps,desc_a,info,&
!!$ r0 = b-ax0
!!$
If (debug) Write(0,*) 'restart: ',itx,it
If (itx.Ge.itmax) Exit restart
If (itx.Ge.litmax) Exit restart
it = 0
Call psb_axpby(one,b,zero,r,desc_a,info)
Call psb_spmm(-one,a,x,one,r,desc_a,info,work=aux)
@ -356,7 +356,7 @@ Subroutine psb_dcgstabl(a,prec,b,x,eps,desc_a,info,&
Exit restart
End If
If (itx.Ge.itmax) Exit restart
If (itx.Ge.litmax) Exit restart
End Do iteration
End Do restart

@ -235,7 +235,7 @@ Subroutine psb_dgmresr(a,prec,b,x,eps,desc_a,info,&
Exit restart
End If
If (itx.Ge.itmax) Exit restart
If (itx.Ge.litmax) Exit restart
v(:,1) = v(:,1) * scal

Loading…
Cancel
Save