diff --git a/mlprec/impl/mld_caggrmat_nosmth_asb.f90 b/mlprec/impl/mld_caggrmat_nosmth_asb.f90 index cd8f26b7..91b01614 100644 --- a/mlprec/impl/mld_caggrmat_nosmth_asb.f90 +++ b/mlprec/impl/mld_caggrmat_nosmth_asb.f90 @@ -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) diff --git a/mlprec/impl/mld_daggrmat_nosmth_asb.f90 b/mlprec/impl/mld_daggrmat_nosmth_asb.f90 index 0cb2f9a4..9fbec777 100644 --- a/mlprec/impl/mld_daggrmat_nosmth_asb.f90 +++ b/mlprec/impl/mld_daggrmat_nosmth_asb.f90 @@ -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) diff --git a/mlprec/impl/mld_saggrmat_nosmth_asb.f90 b/mlprec/impl/mld_saggrmat_nosmth_asb.f90 index 0ecef6ff..78608ea0 100644 --- a/mlprec/impl/mld_saggrmat_nosmth_asb.f90 +++ b/mlprec/impl/mld_saggrmat_nosmth_asb.f90 @@ -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) diff --git a/mlprec/impl/mld_zaggrmat_nosmth_asb.f90 b/mlprec/impl/mld_zaggrmat_nosmth_asb.f90 index 70b627b7..580b4e73 100644 --- a/mlprec/impl/mld_zaggrmat_nosmth_asb.f90 +++ b/mlprec/impl/mld_zaggrmat_nosmth_asb.f90 @@ -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)