From fb35ef9a9eee4e0c91b5130519733a0b592741bf Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Fri, 7 Dec 2012 10:58:29 +0000 Subject: [PATCH] mld2p4-2: mlprec/impl/mld_daggrmat_minnrg_asb.f90 mlprec/impl/mld_zaggrmat_minnrg_asb.f90 mlprec/impl/mld_caggrmat_minnrg_asb.f90 mlprec/impl/mld_saggrmat_minnrg_asb.f90 Finished with impl, now for the subdirs. --- mlprec/impl/mld_caggrmat_minnrg_asb.f90 | 4 ++-- mlprec/impl/mld_daggrmat_minnrg_asb.f90 | 4 ++-- mlprec/impl/mld_saggrmat_minnrg_asb.f90 | 4 ++-- mlprec/impl/mld_zaggrmat_minnrg_asb.f90 | 4 ++-- 4 files changed, 8 insertions(+), 8 deletions(-) diff --git a/mlprec/impl/mld_caggrmat_minnrg_asb.f90 b/mlprec/impl/mld_caggrmat_minnrg_asb.f90 index 335fd75c..7aac81a8 100644 --- a/mlprec/impl/mld_caggrmat_minnrg_asb.f90 +++ b/mlprec/impl/mld_caggrmat_minnrg_asb.f90 @@ -410,7 +410,7 @@ subroutine mld_caggrmat_minnrg_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re call psb_sphalo(atmp,desc_a,am4,info,& & colcnv=.true.,rowscale=.true.) nrt = am4%get_nrows() - call am4%csclip(atmp2,info,1,nrt,1,ncol) + call am4%csclip(atmp2,info,ione,nrt,ione,ncol) call atmp2%cscnv(info,type='CSR') if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2) call am4%free() @@ -419,7 +419,7 @@ subroutine mld_caggrmat_minnrg_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re ! This is to compute the transpose. It ONLY works if the ! original A has a symmetric pattern. call atmp%transc(atmp2) - call atmp2%csclip(dat,info,1,nrow,1,ncol) + call atmp2%csclip(dat,info,ione,nrow,ione,ncol) call dat%cscnv(info,type='csr') call dat%scal(adinv,info) diff --git a/mlprec/impl/mld_daggrmat_minnrg_asb.f90 b/mlprec/impl/mld_daggrmat_minnrg_asb.f90 index 679a09a5..745005d4 100644 --- a/mlprec/impl/mld_daggrmat_minnrg_asb.f90 +++ b/mlprec/impl/mld_daggrmat_minnrg_asb.f90 @@ -410,7 +410,7 @@ subroutine mld_daggrmat_minnrg_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re call psb_sphalo(atmp,desc_a,am4,info,& & colcnv=.true.,rowscale=.true.) nrt = am4%get_nrows() - call am4%csclip(atmp2,info,1,nrt,1,ncol) + call am4%csclip(atmp2,info,ione,nrt,ione,ncol) call atmp2%cscnv(info,type='CSR') if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2) call am4%free() @@ -419,7 +419,7 @@ subroutine mld_daggrmat_minnrg_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re ! This is to compute the transpose. It ONLY works if the ! original A has a symmetric pattern. call atmp%transc(atmp2) - call atmp2%csclip(dat,info,1,nrow,1,ncol) + call atmp2%csclip(dat,info,ione,nrow,ione,ncol) call dat%cscnv(info,type='csr') call dat%scal(adinv,info) diff --git a/mlprec/impl/mld_saggrmat_minnrg_asb.f90 b/mlprec/impl/mld_saggrmat_minnrg_asb.f90 index 6709dac7..fc36b715 100644 --- a/mlprec/impl/mld_saggrmat_minnrg_asb.f90 +++ b/mlprec/impl/mld_saggrmat_minnrg_asb.f90 @@ -410,7 +410,7 @@ subroutine mld_saggrmat_minnrg_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re call psb_sphalo(atmp,desc_a,am4,info,& & colcnv=.true.,rowscale=.true.) nrt = am4%get_nrows() - call am4%csclip(atmp2,info,1,nrt,1,ncol) + call am4%csclip(atmp2,info,ione,nrt,ione,ncol) call atmp2%cscnv(info,type='CSR') if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2) call am4%free() @@ -419,7 +419,7 @@ subroutine mld_saggrmat_minnrg_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re ! This is to compute the transpose. It ONLY works if the ! original A has a symmetric pattern. call atmp%transc(atmp2) - call atmp2%csclip(dat,info,1,nrow,1,ncol) + call atmp2%csclip(dat,info,ione,nrow,ione,ncol) call dat%cscnv(info,type='csr') call dat%scal(adinv,info) diff --git a/mlprec/impl/mld_zaggrmat_minnrg_asb.f90 b/mlprec/impl/mld_zaggrmat_minnrg_asb.f90 index 87bb18e4..1f171ea9 100644 --- a/mlprec/impl/mld_zaggrmat_minnrg_asb.f90 +++ b/mlprec/impl/mld_zaggrmat_minnrg_asb.f90 @@ -410,7 +410,7 @@ subroutine mld_zaggrmat_minnrg_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re call psb_sphalo(atmp,desc_a,am4,info,& & colcnv=.true.,rowscale=.true.) nrt = am4%get_nrows() - call am4%csclip(atmp2,info,1,nrt,1,ncol) + call am4%csclip(atmp2,info,ione,nrt,ione,ncol) call atmp2%cscnv(info,type='CSR') if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2) call am4%free() @@ -419,7 +419,7 @@ subroutine mld_zaggrmat_minnrg_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re ! This is to compute the transpose. It ONLY works if the ! original A has a symmetric pattern. call atmp%transc(atmp2) - call atmp2%csclip(dat,info,1,nrow,1,ncol) + call atmp2%csclip(dat,info,ione,nrow,ione,ncol) call dat%cscnv(info,type='csr') call dat%scal(adinv,info)