From 1e5e86ee89c8a6f6daf6ef43687debca37496075 Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Fri, 10 Aug 2018 15:50:03 +0100 Subject: [PATCH] Use optional argument in GETROW in SOC1_MAP_BLD. --- mlprec/impl/aggregator/mld_c_soc1_map_bld.f90 | 8 ++++---- mlprec/impl/aggregator/mld_d_soc1_map_bld.f90 | 8 ++++---- mlprec/impl/aggregator/mld_s_soc1_map_bld.f90 | 8 ++++---- mlprec/impl/aggregator/mld_z_soc1_map_bld.f90 | 8 ++++---- 4 files changed, 16 insertions(+), 16 deletions(-) diff --git a/mlprec/impl/aggregator/mld_c_soc1_map_bld.f90 b/mlprec/impl/aggregator/mld_c_soc1_map_bld.f90 index 1ec3aa18..89f3703c 100644 --- a/mlprec/impl/aggregator/mld_c_soc1_map_bld.f90 +++ b/mlprec/impl/aggregator/mld_c_soc1_map_bld.f90 @@ -109,7 +109,7 @@ subroutine mld_c_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) ncol = desc_a%get_local_cols() nr = a%get_nrows() - allocate(ilaggr(nr),neigh(nr),ideg(nr),idxs(nr),stat=info) + allocate(ilaggr(nr),neigh(nr),ideg(nr),idxs(nr),irow(nr),icol(nr),val(nr),stat=info) if(info /= psb_success_) then info=psb_err_alloc_request_ call psb_errpush(info,name,i_err=(/2*nr,izero,izero,izero,izero/),& @@ -149,7 +149,7 @@ subroutine mld_c_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) == -(nr+1)) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='csget') @@ -201,7 +201,7 @@ subroutine mld_c_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) == -(nr+1)) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='psb_sp_getrow') @@ -237,7 +237,7 @@ subroutine mld_c_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) < 0) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='psb_sp_getrow') diff --git a/mlprec/impl/aggregator/mld_d_soc1_map_bld.f90 b/mlprec/impl/aggregator/mld_d_soc1_map_bld.f90 index 4c9e935f..8bf12435 100644 --- a/mlprec/impl/aggregator/mld_d_soc1_map_bld.f90 +++ b/mlprec/impl/aggregator/mld_d_soc1_map_bld.f90 @@ -109,7 +109,7 @@ subroutine mld_d_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) ncol = desc_a%get_local_cols() nr = a%get_nrows() - allocate(ilaggr(nr),neigh(nr),ideg(nr),idxs(nr),stat=info) + allocate(ilaggr(nr),neigh(nr),ideg(nr),idxs(nr),irow(nr),icol(nr),val(nr),stat=info) if(info /= psb_success_) then info=psb_err_alloc_request_ call psb_errpush(info,name,i_err=(/2*nr,izero,izero,izero,izero/),& @@ -149,7 +149,7 @@ subroutine mld_d_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) == -(nr+1)) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='csget') @@ -201,7 +201,7 @@ subroutine mld_d_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) == -(nr+1)) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='psb_sp_getrow') @@ -237,7 +237,7 @@ subroutine mld_d_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) < 0) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='psb_sp_getrow') diff --git a/mlprec/impl/aggregator/mld_s_soc1_map_bld.f90 b/mlprec/impl/aggregator/mld_s_soc1_map_bld.f90 index 4729964d..d1ad4e6f 100644 --- a/mlprec/impl/aggregator/mld_s_soc1_map_bld.f90 +++ b/mlprec/impl/aggregator/mld_s_soc1_map_bld.f90 @@ -109,7 +109,7 @@ subroutine mld_s_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) ncol = desc_a%get_local_cols() nr = a%get_nrows() - allocate(ilaggr(nr),neigh(nr),ideg(nr),idxs(nr),stat=info) + allocate(ilaggr(nr),neigh(nr),ideg(nr),idxs(nr),irow(nr),icol(nr),val(nr),stat=info) if(info /= psb_success_) then info=psb_err_alloc_request_ call psb_errpush(info,name,i_err=(/2*nr,izero,izero,izero,izero/),& @@ -149,7 +149,7 @@ subroutine mld_s_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) == -(nr+1)) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='csget') @@ -201,7 +201,7 @@ subroutine mld_s_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) == -(nr+1)) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='psb_sp_getrow') @@ -237,7 +237,7 @@ subroutine mld_s_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) < 0) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='psb_sp_getrow') diff --git a/mlprec/impl/aggregator/mld_z_soc1_map_bld.f90 b/mlprec/impl/aggregator/mld_z_soc1_map_bld.f90 index 9a4bc47b..aa175f9a 100644 --- a/mlprec/impl/aggregator/mld_z_soc1_map_bld.f90 +++ b/mlprec/impl/aggregator/mld_z_soc1_map_bld.f90 @@ -109,7 +109,7 @@ subroutine mld_z_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) ncol = desc_a%get_local_cols() nr = a%get_nrows() - allocate(ilaggr(nr),neigh(nr),ideg(nr),idxs(nr),stat=info) + allocate(ilaggr(nr),neigh(nr),ideg(nr),idxs(nr),irow(nr),icol(nr),val(nr),stat=info) if(info /= psb_success_) then info=psb_err_alloc_request_ call psb_errpush(info,name,i_err=(/2*nr,izero,izero,izero,izero/),& @@ -149,7 +149,7 @@ subroutine mld_z_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) == -(nr+1)) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='csget') @@ -201,7 +201,7 @@ subroutine mld_z_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) == -(nr+1)) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='psb_sp_getrow') @@ -237,7 +237,7 @@ subroutine mld_z_soc1_map_bld(iorder,theta,a,desc_a,nlaggr,ilaggr,info) i = idxs(ii) if (ilaggr(i) < 0) then - call a%csget(i,i,nz,irow,icol,val,info) + call a%csget(i,i,nz,irow,icol,val,info,chksz=.false.) if (info /= psb_success_) then info=psb_err_from_subroutine_ call psb_errpush(info,name,a_err='psb_sp_getrow')