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

Fixed compilation in krylov.

Browse Source
ILmat
Salvatore Filippone 8 years ago
parent 1bb802328b
commit cf4af74883
32 changed files with 144 additions and 112 deletions
Show Stats Download Patch File Download Diff File

7
krylov/psb_cbicg.f90
Unescape Escape View File

@ -117,8 +117,9 @@ subroutine psb_cbicg_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_c_vect_type), pointer :: ww, q, r, p,&
& zt, pt, z, rt, qt
integer(psb_ipk_) :: int_err(5)
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,&
& n_row, n_col, istop_, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
logical, parameter :: exchange=.true., noexchange=.false.
integer(psb_ipk_), parameter :: irmax = 8
@ -174,13 +175,13 @@ subroutine psb_cbicg_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_ccg.F90
Unescape Escape View File

@ -118,8 +118,9 @@ subroutine psb_ccg_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_c_vect_type), allocatable, target :: wwrk(:)
type(psb_c_vect_type), pointer :: q, p, r, z, w
complex(psb_spk_) :: alpha, beta, rho, rho_old, sigma,alpha_old,beta_old
integer(psb_ipk_) :: itmax_, istop_, naux, mglob, it, itx, itrace_,&
integer(psb_ipk_) :: itmax_, istop_, naux, it, itx, itrace_,&
& n_col, n_row,err_act, int_err(5), ieg,nspl, istebz
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
real(psb_dpk_) :: derr
@ -160,9 +161,9 @@ subroutine psb_ccg_vect(a,prec,b,x,eps,desc_a,info,&
istop_ = 2
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_ccgs.f90
Unescape Escape View File

@ -114,8 +114,9 @@ Subroutine psb_ccgs_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_c_vect_type), allocatable, target :: wwrk(:)
type(psb_c_vect_type), pointer :: ww, q, r, p, v,&
& s, z, f, rt, qt, uv
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,int_err(5),&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,int_err(5),&
& n_row, n_col,istop_, itx, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, ictxt
integer(psb_ipk_) :: debug_level, debug_unit
complex(psb_spk_) :: alpha, beta, rho, rho_old, sigma
@ -155,8 +156,8 @@ Subroutine psb_ccgs_vect(a,prec,b,x,eps,desc_a,info,&
istop_ = 2
Endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_ccgstab.f90
Unescape Escape View File

@ -113,8 +113,9 @@ Subroutine psb_ccgstab_vect(a,prec,b,x,eps,desc_a,info,itmax,iter,err,itrace,ist
complex(psb_spk_), allocatable, target :: aux(:),wwrk(:,:)
type(psb_c_vect_type) :: q, r, p, v, s, t, z, f
integer(psb_ipk_) :: itmax_, naux, mglob, it,itrace_,&
integer(psb_ipk_) :: itmax_, naux, it,itrace_,&
& n_row, n_col
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
Logical, Parameter :: exchange=.True., noexchange=.False., debug1 = .False.
integer(psb_ipk_), Parameter :: irmax = 8
@ -165,13 +166,13 @@ Subroutine psb_ccgstab_vect(a,prec,b,x,eps,desc_a,info,itmax,iter,err,itrace,ist
! = write(0,*) 'Warning: different dynamic types for X and B '
! = end if
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (psb_errstatus_fatal()) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if (psb_errstatus_fatal()) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_ccgstabl.f90
Unescape Escape View File

@ -127,8 +127,9 @@ Subroutine psb_ccgstabl_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_c_vect_type), Pointer :: ww, q, r, rt0, p, v, &
& s, t, z, f
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,&
& n_row, n_col, nl, err_act
integer(psb_lpk_) :: mglob
Logical, Parameter :: exchange=.True., noexchange=.False.
integer(psb_ipk_), Parameter :: irmax = 8
integer(psb_ipk_) :: itx, i, istop_,j, k, int_err(5)
@ -204,8 +205,8 @@ Subroutine psb_ccgstabl_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_cfcg.F90
Unescape Escape View File

@ -126,7 +126,8 @@ subroutine psb_cfcg_vect(a,prec,b,x,eps,desc_a,info,&
complex(psb_spk_) :: alpha, tau, tau1, beta, delta
real(psb_dpk_) :: derr
integer(psb_ipk_) :: i, idx, nc2l, it, itx, istop_, itmax_, itrace_
integer(psb_ipk_) :: n_col, mglob, naux, err_act
integer(psb_ipk_) :: n_col, naux, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
complex(psb_spk_), allocatable, target :: aux(:)
@ -165,9 +166,9 @@ subroutine psb_cfcg_vect(a,prec,b,x,eps,desc_a,info,&
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_cgcr.f90
Unescape Escape View File

@ -130,7 +130,8 @@ subroutine psb_cgcr_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_c_vect_type) :: r
real(psb_dpk_) :: r_norm, b_norm, a_norm, derr
integer(psb_ipk_) :: n_col, mglob, naux, err_act
integer(psb_ipk_) :: n_col, naux, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
integer(psb_ipk_) :: i, j, it, itx, istop_, itmax_, itrace_, nl, m, nrst
@ -183,9 +184,9 @@ subroutine psb_cgcr_vect(a,prec,b,x,eps,desc_a,info,&
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

15
krylov/psb_crgmres.f90
Unescape Escape View File

@ -29,10 +29,10 @@
! SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
! INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
! CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
!!$ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
!!$ POSSIBILITY OF SUCH DAMAGE.
!!$
!!$
! ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
! POSSIBILITY OF SUCH DAMAGE.
!
!
! CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
! C C
! C References: C
@ -130,8 +130,9 @@ subroutine psb_crgmres_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_c_vect_type) :: w, w1, xt
real(psb_spk_) :: tmp
complex(psb_spk_) :: scal, gm, rti, rti1
integer(psb_ipk_) ::litmax, naux, mglob, it,k, itrace_,&
integer(psb_ipk_) ::litmax, naux, it, k, itrace_,&
& n_row, n_col, nl, int_err(5)
integer(psb_lpk_) :: mglob
Logical, Parameter :: exchange=.True., noexchange=.False., use_srot=.true.
integer(psb_ipk_), Parameter :: irmax = 8
integer(psb_ipk_) :: itx, i, istop_, err_act
@ -216,13 +217,13 @@ subroutine psb_crgmres_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_dbicg.f90
Unescape Escape View File

@ -117,8 +117,9 @@ subroutine psb_dbicg_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_d_vect_type), pointer :: ww, q, r, p,&
& zt, pt, z, rt, qt
integer(psb_ipk_) :: int_err(5)
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,&
& n_row, n_col, istop_, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
logical, parameter :: exchange=.true., noexchange=.false.
integer(psb_ipk_), parameter :: irmax = 8
@ -174,13 +175,13 @@ subroutine psb_dbicg_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_dcg.F90
Unescape Escape View File

@ -118,8 +118,9 @@ subroutine psb_dcg_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_d_vect_type), allocatable, target :: wwrk(:)
type(psb_d_vect_type), pointer :: q, p, r, z, w
real(psb_dpk_) :: alpha, beta, rho, rho_old, sigma,alpha_old,beta_old
integer(psb_ipk_) :: itmax_, istop_, naux, mglob, it, itx, itrace_,&
integer(psb_ipk_) :: itmax_, istop_, naux, it, itx, itrace_,&
& n_col, n_row,err_act, int_err(5), ieg,nspl, istebz
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
real(psb_dpk_) :: derr
@ -160,9 +161,9 @@ subroutine psb_dcg_vect(a,prec,b,x,eps,desc_a,info,&
istop_ = 2
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_dcgs.f90
Unescape Escape View File

@ -114,8 +114,9 @@ Subroutine psb_dcgs_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_d_vect_type), allocatable, target :: wwrk(:)
type(psb_d_vect_type), pointer :: ww, q, r, p, v,&
& s, z, f, rt, qt, uv
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,int_err(5),&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,int_err(5),&
& n_row, n_col,istop_, itx, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, ictxt
integer(psb_ipk_) :: debug_level, debug_unit
real(psb_dpk_) :: alpha, beta, rho, rho_old, sigma
@ -155,8 +156,8 @@ Subroutine psb_dcgs_vect(a,prec,b,x,eps,desc_a,info,&
istop_ = 2
Endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_dcgstab.f90
Unescape Escape View File

@ -113,8 +113,9 @@ Subroutine psb_dcgstab_vect(a,prec,b,x,eps,desc_a,info,itmax,iter,err,itrace,ist
real(psb_dpk_), allocatable, target :: aux(:),wwrk(:,:)
type(psb_d_vect_type) :: q, r, p, v, s, t, z, f
integer(psb_ipk_) :: itmax_, naux, mglob, it,itrace_,&
integer(psb_ipk_) :: itmax_, naux, it,itrace_,&
& n_row, n_col
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
Logical, Parameter :: exchange=.True., noexchange=.False., debug1 = .False.
integer(psb_ipk_), Parameter :: irmax = 8
@ -165,13 +166,13 @@ Subroutine psb_dcgstab_vect(a,prec,b,x,eps,desc_a,info,itmax,iter,err,itrace,ist
! = write(0,*) 'Warning: different dynamic types for X and B '
! = end if
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (psb_errstatus_fatal()) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if (psb_errstatus_fatal()) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_dcgstabl.f90
Unescape Escape View File

@ -127,8 +127,9 @@ Subroutine psb_dcgstabl_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_d_vect_type), Pointer :: ww, q, r, rt0, p, v, &
& s, t, z, f
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,&
& n_row, n_col, nl, err_act
integer(psb_lpk_) :: mglob
Logical, Parameter :: exchange=.True., noexchange=.False.
integer(psb_ipk_), Parameter :: irmax = 8
integer(psb_ipk_) :: itx, i, istop_,j, k, int_err(5)
@ -204,8 +205,8 @@ Subroutine psb_dcgstabl_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_dfcg.F90
Unescape Escape View File

@ -126,7 +126,8 @@ subroutine psb_dfcg_vect(a,prec,b,x,eps,desc_a,info,&
real(psb_dpk_) :: alpha, tau, tau1, beta, delta
real(psb_dpk_) :: derr
integer(psb_ipk_) :: i, idx, nc2l, it, itx, istop_, itmax_, itrace_
integer(psb_ipk_) :: n_col, mglob, naux, err_act
integer(psb_ipk_) :: n_col, naux, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
real(psb_dpk_), allocatable, target :: aux(:)
@ -165,9 +166,9 @@ subroutine psb_dfcg_vect(a,prec,b,x,eps,desc_a,info,&
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_dgcr.f90
Unescape Escape View File

@ -130,7 +130,8 @@ subroutine psb_dgcr_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_d_vect_type) :: r
real(psb_dpk_) :: r_norm, b_norm, a_norm, derr
integer(psb_ipk_) :: n_col, mglob, naux, err_act
integer(psb_ipk_) :: n_col, naux, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
integer(psb_ipk_) :: i, j, it, itx, istop_, itmax_, itrace_, nl, m, nrst
@ -183,9 +184,9 @@ subroutine psb_dgcr_vect(a,prec,b,x,eps,desc_a,info,&
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

15
krylov/psb_drgmres.f90
Unescape Escape View File

@ -29,10 +29,10 @@
! SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
! INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
! CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
!!$ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
!!$ POSSIBILITY OF SUCH DAMAGE.
!!$
!!$
! ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
! POSSIBILITY OF SUCH DAMAGE.
!
!
! CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
! C C
! C References: C
@ -130,8 +130,9 @@ subroutine psb_drgmres_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_d_vect_type) :: w, w1, xt
real(psb_dpk_) :: tmp
real(psb_dpk_) :: scal, gm, rti, rti1
integer(psb_ipk_) ::litmax, naux, mglob, it,k, itrace_,&
integer(psb_ipk_) ::litmax, naux, it, k, itrace_,&
& n_row, n_col, nl, int_err(5)
integer(psb_lpk_) :: mglob
Logical, Parameter :: exchange=.True., noexchange=.False., use_srot=.true.
integer(psb_ipk_), Parameter :: irmax = 8
integer(psb_ipk_) :: itx, i, istop_, err_act
@ -216,13 +217,13 @@ subroutine psb_drgmres_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_sbicg.f90
Unescape Escape View File

@ -117,8 +117,9 @@ subroutine psb_sbicg_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_s_vect_type), pointer :: ww, q, r, p,&
& zt, pt, z, rt, qt
integer(psb_ipk_) :: int_err(5)
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,&
& n_row, n_col, istop_, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
logical, parameter :: exchange=.true., noexchange=.false.
integer(psb_ipk_), parameter :: irmax = 8
@ -174,13 +175,13 @@ subroutine psb_sbicg_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_scg.F90
Unescape Escape View File

@ -118,8 +118,9 @@ subroutine psb_scg_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_s_vect_type), allocatable, target :: wwrk(:)
type(psb_s_vect_type), pointer :: q, p, r, z, w
real(psb_spk_) :: alpha, beta, rho, rho_old, sigma,alpha_old,beta_old
integer(psb_ipk_) :: itmax_, istop_, naux, mglob, it, itx, itrace_,&
integer(psb_ipk_) :: itmax_, istop_, naux, it, itx, itrace_,&
& n_col, n_row,err_act, int_err(5), ieg,nspl, istebz
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
real(psb_dpk_) :: derr
@ -160,9 +161,9 @@ subroutine psb_scg_vect(a,prec,b,x,eps,desc_a,info,&
istop_ = 2
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_scgs.f90
Unescape Escape View File

@ -114,8 +114,9 @@ Subroutine psb_scgs_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_s_vect_type), allocatable, target :: wwrk(:)
type(psb_s_vect_type), pointer :: ww, q, r, p, v,&
& s, z, f, rt, qt, uv
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,int_err(5),&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,int_err(5),&
& n_row, n_col,istop_, itx, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, ictxt
integer(psb_ipk_) :: debug_level, debug_unit
real(psb_spk_) :: alpha, beta, rho, rho_old, sigma
@ -155,8 +156,8 @@ Subroutine psb_scgs_vect(a,prec,b,x,eps,desc_a,info,&
istop_ = 2
Endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_scgstab.f90
Unescape Escape View File

@ -113,8 +113,9 @@ Subroutine psb_scgstab_vect(a,prec,b,x,eps,desc_a,info,itmax,iter,err,itrace,ist
real(psb_spk_), allocatable, target :: aux(:),wwrk(:,:)
type(psb_s_vect_type) :: q, r, p, v, s, t, z, f
integer(psb_ipk_) :: itmax_, naux, mglob, it,itrace_,&
integer(psb_ipk_) :: itmax_, naux, it,itrace_,&
& n_row, n_col
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
Logical, Parameter :: exchange=.True., noexchange=.False., debug1 = .False.
integer(psb_ipk_), Parameter :: irmax = 8
@ -165,13 +166,13 @@ Subroutine psb_scgstab_vect(a,prec,b,x,eps,desc_a,info,itmax,iter,err,itrace,ist
! = write(0,*) 'Warning: different dynamic types for X and B '
! = end if
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (psb_errstatus_fatal()) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if (psb_errstatus_fatal()) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_scgstabl.f90
Unescape Escape View File

@ -127,8 +127,9 @@ Subroutine psb_scgstabl_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_s_vect_type), Pointer :: ww, q, r, rt0, p, v, &
& s, t, z, f
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,&
& n_row, n_col, nl, err_act
integer(psb_lpk_) :: mglob
Logical, Parameter :: exchange=.True., noexchange=.False.
integer(psb_ipk_), Parameter :: irmax = 8
integer(psb_ipk_) :: itx, i, istop_,j, k, int_err(5)
@ -204,8 +205,8 @@ Subroutine psb_scgstabl_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_sfcg.F90
Unescape Escape View File

@ -126,7 +126,8 @@ subroutine psb_sfcg_vect(a,prec,b,x,eps,desc_a,info,&
real(psb_spk_) :: alpha, tau, tau1, beta, delta
real(psb_dpk_) :: derr
integer(psb_ipk_) :: i, idx, nc2l, it, itx, istop_, itmax_, itrace_
integer(psb_ipk_) :: n_col, mglob, naux, err_act
integer(psb_ipk_) :: n_col, naux, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
real(psb_spk_), allocatable, target :: aux(:)
@ -165,9 +166,9 @@ subroutine psb_sfcg_vect(a,prec,b,x,eps,desc_a,info,&
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_sgcr.f90
Unescape Escape View File

@ -130,7 +130,8 @@ subroutine psb_sgcr_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_s_vect_type) :: r
real(psb_dpk_) :: r_norm, b_norm, a_norm, derr
integer(psb_ipk_) :: n_col, mglob, naux, err_act
integer(psb_ipk_) :: n_col, naux, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
integer(psb_ipk_) :: i, j, it, itx, istop_, itmax_, itrace_, nl, m, nrst
@ -183,9 +184,9 @@ subroutine psb_sgcr_vect(a,prec,b,x,eps,desc_a,info,&
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

15
krylov/psb_srgmres.f90
Unescape Escape View File

@ -29,10 +29,10 @@
! SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
! INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
! CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
!!$ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
!!$ POSSIBILITY OF SUCH DAMAGE.
!!$
!!$
! ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
! POSSIBILITY OF SUCH DAMAGE.
!
!
! CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
! C C
! C References: C
@ -130,8 +130,9 @@ subroutine psb_srgmres_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_s_vect_type) :: w, w1, xt
real(psb_spk_) :: tmp
real(psb_spk_) :: scal, gm, rti, rti1
integer(psb_ipk_) ::litmax, naux, mglob, it,k, itrace_,&
integer(psb_ipk_) ::litmax, naux, it, k, itrace_,&
& n_row, n_col, nl, int_err(5)
integer(psb_lpk_) :: mglob
Logical, Parameter :: exchange=.True., noexchange=.False., use_srot=.true.
integer(psb_ipk_), Parameter :: irmax = 8
integer(psb_ipk_) :: itx, i, istop_, err_act
@ -216,13 +217,13 @@ subroutine psb_srgmres_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_zbicg.f90
Unescape Escape View File

@ -117,8 +117,9 @@ subroutine psb_zbicg_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_z_vect_type), pointer :: ww, q, r, p,&
& zt, pt, z, rt, qt
integer(psb_ipk_) :: int_err(5)
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,&
& n_row, n_col, istop_, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
logical, parameter :: exchange=.true., noexchange=.false.
integer(psb_ipk_), parameter :: irmax = 8
@ -174,13 +175,13 @@ subroutine psb_zbicg_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_zcg.F90
Unescape Escape View File

@ -118,8 +118,9 @@ subroutine psb_zcg_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_z_vect_type), allocatable, target :: wwrk(:)
type(psb_z_vect_type), pointer :: q, p, r, z, w
complex(psb_dpk_) :: alpha, beta, rho, rho_old, sigma,alpha_old,beta_old
integer(psb_ipk_) :: itmax_, istop_, naux, mglob, it, itx, itrace_,&
integer(psb_ipk_) :: itmax_, istop_, naux, it, itx, itrace_,&
& n_col, n_row,err_act, int_err(5), ieg,nspl, istebz
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
real(psb_dpk_) :: derr
@ -160,9 +161,9 @@ subroutine psb_zcg_vect(a,prec,b,x,eps,desc_a,info,&
istop_ = 2
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_zcgs.f90
Unescape Escape View File

@ -114,8 +114,9 @@ Subroutine psb_zcgs_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_z_vect_type), allocatable, target :: wwrk(:)
type(psb_z_vect_type), pointer :: ww, q, r, p, v,&
& s, z, f, rt, qt, uv
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,int_err(5),&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,int_err(5),&
& n_row, n_col,istop_, itx, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, ictxt
integer(psb_ipk_) :: debug_level, debug_unit
complex(psb_dpk_) :: alpha, beta, rho, rho_old, sigma
@ -155,8 +156,8 @@ Subroutine psb_zcgs_vect(a,prec,b,x,eps,desc_a,info,&
istop_ = 2
Endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_zcgstab.f90
Unescape Escape View File

@ -113,8 +113,9 @@ Subroutine psb_zcgstab_vect(a,prec,b,x,eps,desc_a,info,itmax,iter,err,itrace,ist
complex(psb_dpk_), allocatable, target :: aux(:),wwrk(:,:)
type(psb_z_vect_type) :: q, r, p, v, s, t, z, f
integer(psb_ipk_) :: itmax_, naux, mglob, it,itrace_,&
integer(psb_ipk_) :: itmax_, naux, it,itrace_,&
& n_row, n_col
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
Logical, Parameter :: exchange=.True., noexchange=.False., debug1 = .False.
integer(psb_ipk_), Parameter :: irmax = 8
@ -165,13 +166,13 @@ Subroutine psb_zcgstab_vect(a,prec,b,x,eps,desc_a,info,itmax,iter,err,itrace,ist
! = write(0,*) 'Warning: different dynamic types for X and B '
! = end if
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (psb_errstatus_fatal()) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if (psb_errstatus_fatal()) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

7
krylov/psb_zcgstabl.f90
Unescape Escape View File

@ -127,8 +127,9 @@ Subroutine psb_zcgstabl_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_z_vect_type), Pointer :: ww, q, r, rt0, p, v, &
& s, t, z, f
integer(psb_ipk_) :: itmax_, naux, mglob, it, itrace_,&
integer(psb_ipk_) :: itmax_, naux, it, itrace_,&
& n_row, n_col, nl, err_act
integer(psb_lpk_) :: mglob
Logical, Parameter :: exchange=.True., noexchange=.False.
integer(psb_ipk_), Parameter :: irmax = 8
integer(psb_ipk_) :: itx, i, istop_,j, k, int_err(5)
@ -204,8 +205,8 @@ Subroutine psb_zcgstabl_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_) call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_zfcg.F90
Unescape Escape View File

@ -126,7 +126,8 @@ subroutine psb_zfcg_vect(a,prec,b,x,eps,desc_a,info,&
complex(psb_dpk_) :: alpha, tau, tau1, beta, delta
real(psb_dpk_) :: derr
integer(psb_ipk_) :: i, idx, nc2l, it, itx, istop_, itmax_, itrace_
integer(psb_ipk_) :: n_col, mglob, naux, err_act
integer(psb_ipk_) :: n_col, naux, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
complex(psb_dpk_), allocatable, target :: aux(:)
@ -165,9 +166,9 @@ subroutine psb_zfcg_vect(a,prec,b,x,eps,desc_a,info,&
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

7
krylov/psb_zgcr.f90
Unescape Escape View File

@ -130,7 +130,8 @@ subroutine psb_zgcr_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_z_vect_type) :: r
real(psb_dpk_) :: r_norm, b_norm, a_norm, derr
integer(psb_ipk_) :: n_col, mglob, naux, err_act
integer(psb_ipk_) :: n_col, naux, err_act
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_ipk_) :: np, me, ictxt
integer(psb_ipk_) :: i, j, it, itx, istop_, itmax_, itrace_, nl, m, nrst
@ -183,9 +184,9 @@ subroutine psb_zgcr_vect(a,prec,b,x,eps,desc_a,info,&
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if (info == psb_success_)&
& call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
& call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X/B')

15
krylov/psb_zrgmres.f90
Unescape Escape View File

@ -29,10 +29,10 @@
! SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
! INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
! CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
!!$ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
!!$ POSSIBILITY OF SUCH DAMAGE.
!!$
!!$
! ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
! POSSIBILITY OF SUCH DAMAGE.
!
!
! CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
! C C
! C References: C
@ -130,8 +130,9 @@ subroutine psb_zrgmres_vect(a,prec,b,x,eps,desc_a,info,&
type(psb_z_vect_type) :: w, w1, xt
real(psb_dpk_) :: tmp
complex(psb_dpk_) :: scal, gm, rti, rti1
integer(psb_ipk_) ::litmax, naux, mglob, it,k, itrace_,&
integer(psb_ipk_) ::litmax, naux, it, k, itrace_,&
& n_row, n_col, nl, int_err(5)
integer(psb_lpk_) :: mglob
Logical, Parameter :: exchange=.True., noexchange=.False., use_srot=.true.
integer(psb_ipk_), Parameter :: irmax = 8
integer(psb_ipk_) :: itx, i, istop_, err_act
@ -216,13 +217,13 @@ subroutine psb_zrgmres_vect(a,prec,b,x,eps,desc_a,info,&
goto 9999
endif
call psb_chkvect(mglob,ione,x%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,x%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on X')
goto 9999
end if
call psb_chkvect(mglob,ione,b%get_nrows(),ione,ione,desc_a,info)
call psb_chkvect(mglob,lone,b%get_nrows(),lone,lone,desc_a,info)
if(info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='psb_chkvect on B')

Write Preview
Loading…
Cancel
Save