From 11c039ea76d03e550c64394658ec0f5ecdd5cc49 Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Fri, 13 Jul 2007 17:07:43 +0000 Subject: [PATCH] Fixed performance problem: need to add fixcoo for ILU when sweeps > 1. --- mld_dbjac_bld.f90 | 1 + mld_zbjac_bld.f90 | 1 + 2 files changed, 2 insertions(+) diff --git a/mld_dbjac_bld.f90 b/mld_dbjac_bld.f90 index a79c3688..095e6085 100644 --- 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. diff --git a/mld_zbjac_bld.f90 b/mld_zbjac_bld.f90 index af6e4637..a271f120 100644 --- 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.