Fixed performance problem: need to add fixcoo for ILU when sweeps > 1.

18 years ago · 11c039ea76
parent c420d65d9e
commit 11c039ea76
2 changed files with 2 additions and 0 deletions
--- a/mld_dbjac_bld.f90
+++ b/mld_dbjac_bld.f90
@ -282,6 +282,7 @@ subroutine mld_dbjac_bld(a,desc_a,p,upd,info)
        n_row = psb_cd_get_local_rows(p%desc_data)
        n_col = psb_cd_get_local_cols(p%desc_data)
        call psb_rwextd(n_row,atmp,info,b=blck,rowscale=.false.) 
+        call psb_fixcoo(atmp,info)
        !------------------------------------------------------------------
        ! Split AC=M+N  N off-diagonal part
        ! Output in COO format. 
--- a/mld_zbjac_bld.f90
+++ b/mld_zbjac_bld.f90
@ -274,6 +274,7 @@ subroutine mld_zbjac_bld(a,desc_a,p,upd,info)
        n_row = psb_cd_get_local_rows(p%desc_data)
        n_col = psb_cd_get_local_cols(p%desc_data)
        call psb_rwextd(n_row,atmp,info,b=blck,rowscale=.false.) 
+        call psb_fixcoo(atmp,info)
        !------------------------------------------------------------------
        ! Split AC=M+N  N off-diagonal part
        ! Output in COO format.