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

mld2p4-2:

Browse Source 
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.
stopcriterion
Salvatore Filippone 12 years ago
parent 0b5ebcbfb5
commit fb35ef9a9e
4 changed files with 8 additions and 8 deletions
Show Stats Download Patch File Download Diff File

4
mlprec/impl/mld_caggrmat_minnrg_asb.f90
Unescape Escape View File

@ -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,& call psb_sphalo(atmp,desc_a,am4,info,&
& colcnv=.true.,rowscale=.true.) & colcnv=.true.,rowscale=.true.)
nrt = am4%get_nrows() 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') call atmp2%cscnv(info,type='CSR')
if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2) if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2)
call am4%free() 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 ! This is to compute the transpose. It ONLY works if the
! original A has a symmetric pattern. ! original A has a symmetric pattern.
call atmp%transc(atmp2) 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%cscnv(info,type='csr')
call dat%scal(adinv,info) call dat%scal(adinv,info)

4
mlprec/impl/mld_daggrmat_minnrg_asb.f90
Unescape Escape View File

@ -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,& call psb_sphalo(atmp,desc_a,am4,info,&
& colcnv=.true.,rowscale=.true.) & colcnv=.true.,rowscale=.true.)
nrt = am4%get_nrows() 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') call atmp2%cscnv(info,type='CSR')
if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2) if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2)
call am4%free() 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 ! This is to compute the transpose. It ONLY works if the
! original A has a symmetric pattern. ! original A has a symmetric pattern.
call atmp%transc(atmp2) 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%cscnv(info,type='csr')
call dat%scal(adinv,info) call dat%scal(adinv,info)

4
mlprec/impl/mld_saggrmat_minnrg_asb.f90
Unescape Escape View File

@ -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,& call psb_sphalo(atmp,desc_a,am4,info,&
& colcnv=.true.,rowscale=.true.) & colcnv=.true.,rowscale=.true.)
nrt = am4%get_nrows() 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') call atmp2%cscnv(info,type='CSR')
if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2) if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2)
call am4%free() 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 ! This is to compute the transpose. It ONLY works if the
! original A has a symmetric pattern. ! original A has a symmetric pattern.
call atmp%transc(atmp2) 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%cscnv(info,type='csr')
call dat%scal(adinv,info) call dat%scal(adinv,info)

4
mlprec/impl/mld_zaggrmat_minnrg_asb.f90
Unescape Escape View File

@ -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,& call psb_sphalo(atmp,desc_a,am4,info,&
& colcnv=.true.,rowscale=.true.) & colcnv=.true.,rowscale=.true.)
nrt = am4%get_nrows() 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') call atmp2%cscnv(info,type='CSR')
if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2) if (info == psb_success_) call psb_rwextd(ncol,atmp,info,b=atmp2)
call am4%free() 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 ! This is to compute the transpose. It ONLY works if the
! original A has a symmetric pattern. ! original A has a symmetric pattern.
call atmp%transc(atmp2) 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%cscnv(info,type='csr')
call dat%scal(adinv,info) call dat%scal(adinv,info)

Write Preview
Loading…
Cancel
Save