From cf4af74883744ccf8130830519a5dc7a118cea34 Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Thu, 12 Apr 2018 13:52:53 +0100 Subject: [PATCH] Fixed compilation in krylov. --- krylov/psb_cbicg.f90 | 7 ++++--- krylov/psb_ccg.F90 | 7 ++++--- krylov/psb_ccgs.f90 | 7 ++++--- krylov/psb_ccgstab.f90 | 7 ++++--- krylov/psb_ccgstabl.f90 | 7 ++++--- krylov/psb_cfcg.F90 | 7 ++++--- krylov/psb_cgcr.f90 | 7 ++++--- krylov/psb_crgmres.f90 | 15 ++++++++------- krylov/psb_dbicg.f90 | 7 ++++--- krylov/psb_dcg.F90 | 7 ++++--- krylov/psb_dcgs.f90 | 7 ++++--- krylov/psb_dcgstab.f90 | 7 ++++--- krylov/psb_dcgstabl.f90 | 7 ++++--- krylov/psb_dfcg.F90 | 7 ++++--- krylov/psb_dgcr.f90 | 7 ++++--- krylov/psb_drgmres.f90 | 15 ++++++++------- krylov/psb_sbicg.f90 | 7 ++++--- krylov/psb_scg.F90 | 7 ++++--- krylov/psb_scgs.f90 | 7 ++++--- krylov/psb_scgstab.f90 | 7 ++++--- krylov/psb_scgstabl.f90 | 7 ++++--- krylov/psb_sfcg.F90 | 7 ++++--- krylov/psb_sgcr.f90 | 7 ++++--- krylov/psb_srgmres.f90 | 15 ++++++++------- krylov/psb_zbicg.f90 | 7 ++++--- krylov/psb_zcg.F90 | 7 ++++--- krylov/psb_zcgs.f90 | 7 ++++--- krylov/psb_zcgstab.f90 | 7 ++++--- krylov/psb_zcgstabl.f90 | 7 ++++--- krylov/psb_zfcg.F90 | 7 ++++--- krylov/psb_zgcr.f90 | 7 ++++--- krylov/psb_zrgmres.f90 | 15 ++++++++------- 32 files changed, 144 insertions(+), 112 deletions(-) diff --git a/krylov/psb_cbicg.f90 b/krylov/psb_cbicg.f90 index bb957ef11..0a4fd1185 100644 --- a/krylov/psb_cbicg.f90 +++ b/krylov/psb_cbicg.f90 @@ -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') diff --git a/krylov/psb_ccg.F90 b/krylov/psb_ccg.F90 index 35434e163..af970f560 100644 --- a/krylov/psb_ccg.F90 +++ b/krylov/psb_ccg.F90 @@ -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') diff --git a/krylov/psb_ccgs.f90 b/krylov/psb_ccgs.f90 index b6542aaee..3043c4c3f 100644 --- a/krylov/psb_ccgs.f90 +++ b/krylov/psb_ccgs.f90 @@ -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') diff --git a/krylov/psb_ccgstab.f90 b/krylov/psb_ccgstab.f90 index 9cbc3f90b..c8c1709c1 100644 --- a/krylov/psb_ccgstab.f90 +++ b/krylov/psb_ccgstab.f90 @@ -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') diff --git a/krylov/psb_ccgstabl.f90 b/krylov/psb_ccgstabl.f90 index 8ab0453ac..093785b2a 100644 --- a/krylov/psb_ccgstabl.f90 +++ b/krylov/psb_ccgstabl.f90 @@ -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') diff --git a/krylov/psb_cfcg.F90 b/krylov/psb_cfcg.F90 index 0e4c0271d..d96ab5417 100644 --- a/krylov/psb_cfcg.F90 +++ b/krylov/psb_cfcg.F90 @@ -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') diff --git a/krylov/psb_cgcr.f90 b/krylov/psb_cgcr.f90 index 1ce78fc38..5eb715938 100644 --- a/krylov/psb_cgcr.f90 +++ b/krylov/psb_cgcr.f90 @@ -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') diff --git a/krylov/psb_crgmres.f90 b/krylov/psb_crgmres.f90 index d7e752513..f32f3d13a 100644 --- a/krylov/psb_crgmres.f90 +++ b/krylov/psb_crgmres.f90 @@ -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') diff --git a/krylov/psb_dbicg.f90 b/krylov/psb_dbicg.f90 index be2ed48f4..9a2c109c6 100644 --- a/krylov/psb_dbicg.f90 +++ b/krylov/psb_dbicg.f90 @@ -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') diff --git a/krylov/psb_dcg.F90 b/krylov/psb_dcg.F90 index f99877b67..efdb1334d 100644 --- a/krylov/psb_dcg.F90 +++ b/krylov/psb_dcg.F90 @@ -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') diff --git a/krylov/psb_dcgs.f90 b/krylov/psb_dcgs.f90 index dca1748b8..22bd66710 100644 --- a/krylov/psb_dcgs.f90 +++ b/krylov/psb_dcgs.f90 @@ -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') diff --git a/krylov/psb_dcgstab.f90 b/krylov/psb_dcgstab.f90 index 73b0662d1..449e5511f 100644 --- a/krylov/psb_dcgstab.f90 +++ b/krylov/psb_dcgstab.f90 @@ -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') diff --git a/krylov/psb_dcgstabl.f90 b/krylov/psb_dcgstabl.f90 index c85f562b2..1c57cfbeb 100644 --- a/krylov/psb_dcgstabl.f90 +++ b/krylov/psb_dcgstabl.f90 @@ -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') diff --git a/krylov/psb_dfcg.F90 b/krylov/psb_dfcg.F90 index 2821829a0..ba5d25a81 100644 --- a/krylov/psb_dfcg.F90 +++ b/krylov/psb_dfcg.F90 @@ -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') diff --git a/krylov/psb_dgcr.f90 b/krylov/psb_dgcr.f90 index acbaa5a73..cef088d75 100644 --- a/krylov/psb_dgcr.f90 +++ b/krylov/psb_dgcr.f90 @@ -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') diff --git a/krylov/psb_drgmres.f90 b/krylov/psb_drgmres.f90 index 368a64822..16f58d873 100644 --- a/krylov/psb_drgmres.f90 +++ b/krylov/psb_drgmres.f90 @@ -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') diff --git a/krylov/psb_sbicg.f90 b/krylov/psb_sbicg.f90 index 2afaa7160..5560cfc53 100644 --- a/krylov/psb_sbicg.f90 +++ b/krylov/psb_sbicg.f90 @@ -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') diff --git a/krylov/psb_scg.F90 b/krylov/psb_scg.F90 index 5eade9a19..0e6c88a73 100644 --- a/krylov/psb_scg.F90 +++ b/krylov/psb_scg.F90 @@ -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') diff --git a/krylov/psb_scgs.f90 b/krylov/psb_scgs.f90 index d053444a4..f2e0cfe97 100644 --- a/krylov/psb_scgs.f90 +++ b/krylov/psb_scgs.f90 @@ -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') diff --git a/krylov/psb_scgstab.f90 b/krylov/psb_scgstab.f90 index c0ca345cc..62eb965ad 100644 --- a/krylov/psb_scgstab.f90 +++ b/krylov/psb_scgstab.f90 @@ -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') diff --git a/krylov/psb_scgstabl.f90 b/krylov/psb_scgstabl.f90 index 0f888980c..bae82d5ee 100644 --- a/krylov/psb_scgstabl.f90 +++ b/krylov/psb_scgstabl.f90 @@ -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') diff --git a/krylov/psb_sfcg.F90 b/krylov/psb_sfcg.F90 index 2fae24b42..3e8c39841 100644 --- a/krylov/psb_sfcg.F90 +++ b/krylov/psb_sfcg.F90 @@ -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') diff --git a/krylov/psb_sgcr.f90 b/krylov/psb_sgcr.f90 index 08333af28..2da4f49a1 100644 --- a/krylov/psb_sgcr.f90 +++ b/krylov/psb_sgcr.f90 @@ -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') diff --git a/krylov/psb_srgmres.f90 b/krylov/psb_srgmres.f90 index 311bc0fb3..261636f23 100644 --- a/krylov/psb_srgmres.f90 +++ b/krylov/psb_srgmres.f90 @@ -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') diff --git a/krylov/psb_zbicg.f90 b/krylov/psb_zbicg.f90 index 969280d39..520c9c5c4 100644 --- a/krylov/psb_zbicg.f90 +++ b/krylov/psb_zbicg.f90 @@ -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') diff --git a/krylov/psb_zcg.F90 b/krylov/psb_zcg.F90 index 858b7822a..379ca0f36 100644 --- a/krylov/psb_zcg.F90 +++ b/krylov/psb_zcg.F90 @@ -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') diff --git a/krylov/psb_zcgs.f90 b/krylov/psb_zcgs.f90 index 3130427ad..1e83df369 100644 --- a/krylov/psb_zcgs.f90 +++ b/krylov/psb_zcgs.f90 @@ -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') diff --git a/krylov/psb_zcgstab.f90 b/krylov/psb_zcgstab.f90 index 68381c482..4fca5c038 100644 --- a/krylov/psb_zcgstab.f90 +++ b/krylov/psb_zcgstab.f90 @@ -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') diff --git a/krylov/psb_zcgstabl.f90 b/krylov/psb_zcgstabl.f90 index 5ddbf3092..5efb11395 100644 --- a/krylov/psb_zcgstabl.f90 +++ b/krylov/psb_zcgstabl.f90 @@ -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') diff --git a/krylov/psb_zfcg.F90 b/krylov/psb_zfcg.F90 index 0e2ad2dab..c6aa6807f 100644 --- a/krylov/psb_zfcg.F90 +++ b/krylov/psb_zfcg.F90 @@ -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') diff --git a/krylov/psb_zgcr.f90 b/krylov/psb_zgcr.f90 index 40b4cf7fe..55ba874c9 100644 --- a/krylov/psb_zgcr.f90 +++ b/krylov/psb_zgcr.f90 @@ -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') diff --git a/krylov/psb_zrgmres.f90 b/krylov/psb_zrgmres.f90 index deb4111f9..38f02bd91 100644 --- a/krylov/psb_zrgmres.f90 +++ b/krylov/psb_zrgmres.f90 @@ -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')