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

mld2p4:

Browse Source 
mlprec/impl/mld_caggrmat_nosmth_asb.f90
 mlprec/impl/mld_daggrmat_nosmth_asb.f90
 mlprec/impl/mld_saggrmat_nosmth_asb.f90
 mlprec/impl/mld_zaggrmat_nosmth_asb.f90

Do not clip elements off the diagonal block.
stopcriterion
Salvatore Filippone 9 years ago
parent b0d7272f11
commit a128383525
4 changed files with 16 additions and 24 deletions
Show Stats Download Patch File Download Diff File

10
mlprec/impl/mld_caggrmat_nosmth_asb.f90
Unescape Escape View File

@ -157,12 +157,10 @@ subroutine mld_caggrmat_nosmth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
nzt = ac_coo%get_nzeros()
k = 0
do i=1, nzt
if (ac_coo%ja(i) <= nrow) then
k = k + 1
ac_coo%ia(k) = ilaggr(ac_coo%ia(i))
ac_coo%ja(k) = ilaggr(ac_coo%ja(i))
ac_coo%val(k) = ac_coo%val(i)
end if
k = k + 1
ac_coo%ia(k) = ilaggr(ac_coo%ia(i))
ac_coo%ja(k) = ilaggr(ac_coo%ja(i))
ac_coo%val(k) = ac_coo%val(i)
enddo
call ac_coo%set_nrows(naggr)
call ac_coo%set_ncols(naggr)

10
mlprec/impl/mld_daggrmat_nosmth_asb.f90
Unescape Escape View File

@ -157,12 +157,10 @@ subroutine mld_daggrmat_nosmth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
nzt = ac_coo%get_nzeros()
k = 0
do i=1, nzt
if (ac_coo%ja(i) <= nrow) then
k = k + 1
ac_coo%ia(k) = ilaggr(ac_coo%ia(i))
ac_coo%ja(k) = ilaggr(ac_coo%ja(i))
ac_coo%val(k) = ac_coo%val(i)
end if
k = k + 1
ac_coo%ia(k) = ilaggr(ac_coo%ia(i))
ac_coo%ja(k) = ilaggr(ac_coo%ja(i))
ac_coo%val(k) = ac_coo%val(i)
enddo
call ac_coo%set_nrows(naggr)
call ac_coo%set_ncols(naggr)

10
mlprec/impl/mld_saggrmat_nosmth_asb.f90
Unescape Escape View File

@ -157,12 +157,10 @@ subroutine mld_saggrmat_nosmth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
nzt = ac_coo%get_nzeros()
k = 0
do i=1, nzt
if (ac_coo%ja(i) <= nrow) then
k = k + 1
ac_coo%ia(k) = ilaggr(ac_coo%ia(i))
ac_coo%ja(k) = ilaggr(ac_coo%ja(i))
ac_coo%val(k) = ac_coo%val(i)
end if
k = k + 1
ac_coo%ia(k) = ilaggr(ac_coo%ia(i))
ac_coo%ja(k) = ilaggr(ac_coo%ja(i))
ac_coo%val(k) = ac_coo%val(i)
enddo
call ac_coo%set_nrows(naggr)
call ac_coo%set_ncols(naggr)

10
mlprec/impl/mld_zaggrmat_nosmth_asb.f90
Unescape Escape View File

@ -157,12 +157,10 @@ subroutine mld_zaggrmat_nosmth_asb(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
nzt = ac_coo%get_nzeros()
k = 0
do i=1, nzt
if (ac_coo%ja(i) <= nrow) then
k = k + 1
ac_coo%ia(k) = ilaggr(ac_coo%ia(i))
ac_coo%ja(k) = ilaggr(ac_coo%ja(i))
ac_coo%val(k) = ac_coo%val(i)
end if
k = k + 1
ac_coo%ia(k) = ilaggr(ac_coo%ia(i))
ac_coo%ja(k) = ilaggr(ac_coo%ja(i))
ac_coo%val(k) = ac_coo%val(i)
enddo
call ac_coo%set_nrows(naggr)
call ac_coo%set_ncols(naggr)

Write Preview
Loading…
Cancel
Save