From 257fd0d00adc2dd6ed84b65b5a115d941bdbbb2c Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Thu, 9 Jul 2020 11:23:53 +0200 Subject: [PATCH] Fix matrix assembly for REPL case. --- mlprec/impl/aggregator/mld_c_dec_aggregator_mat_asb.f90 | 4 ++-- mlprec/impl/aggregator/mld_d_dec_aggregator_mat_asb.f90 | 4 ++-- mlprec/impl/aggregator/mld_s_dec_aggregator_mat_asb.f90 | 4 ++-- mlprec/impl/aggregator/mld_z_dec_aggregator_mat_asb.f90 | 4 ++-- 4 files changed, 8 insertions(+), 8 deletions(-) diff --git a/mlprec/impl/aggregator/mld_c_dec_aggregator_mat_asb.f90 b/mlprec/impl/aggregator/mld_c_dec_aggregator_mat_asb.f90 index 01d6e6e2..a98771e3 100644 --- a/mlprec/impl/aggregator/mld_c_dec_aggregator_mat_asb.f90 +++ b/mlprec/impl/aggregator/mld_c_dec_aggregator_mat_asb.f90 @@ -143,15 +143,15 @@ subroutine mld_c_dec_aggregator_mat_asb(ag,parms,a,desc_a,& call op_prol%mv_to(tmpcoo) nzl = tmpcoo%get_nzeros() call psb_loc_to_glob(tmpcoo%ja(1:nzl),desc_ac,info,'I') + call tmpcoo%set_ncols(i_nr) call op_prol%mv_from(tmpcoo) call op_restr%mv_to(tmpcoo) nzl = tmpcoo%get_nzeros() call psb_loc_to_glob(tmpcoo%ia(1:nzl),desc_ac,info,'I') + call tmpcoo%set_nrows(i_nr) call op_restr%mv_from(tmpcoo) - call op_prol%set_ncols(i_nr) - call op_restr%set_nrows(i_nr) call psb_gather(tmp_ac,ac,desc_ac,info,root=-ione,& & dupl=psb_dupl_add_,keeploc=.false.) diff --git a/mlprec/impl/aggregator/mld_d_dec_aggregator_mat_asb.f90 b/mlprec/impl/aggregator/mld_d_dec_aggregator_mat_asb.f90 index 00ae22b0..61c2a1c3 100644 --- a/mlprec/impl/aggregator/mld_d_dec_aggregator_mat_asb.f90 +++ b/mlprec/impl/aggregator/mld_d_dec_aggregator_mat_asb.f90 @@ -143,15 +143,15 @@ subroutine mld_d_dec_aggregator_mat_asb(ag,parms,a,desc_a,& call op_prol%mv_to(tmpcoo) nzl = tmpcoo%get_nzeros() call psb_loc_to_glob(tmpcoo%ja(1:nzl),desc_ac,info,'I') + call tmpcoo%set_ncols(i_nr) call op_prol%mv_from(tmpcoo) call op_restr%mv_to(tmpcoo) nzl = tmpcoo%get_nzeros() call psb_loc_to_glob(tmpcoo%ia(1:nzl),desc_ac,info,'I') + call tmpcoo%set_nrows(i_nr) call op_restr%mv_from(tmpcoo) - call op_prol%set_ncols(i_nr) - call op_restr%set_nrows(i_nr) call psb_gather(tmp_ac,ac,desc_ac,info,root=-ione,& & dupl=psb_dupl_add_,keeploc=.false.) diff --git a/mlprec/impl/aggregator/mld_s_dec_aggregator_mat_asb.f90 b/mlprec/impl/aggregator/mld_s_dec_aggregator_mat_asb.f90 index 74bba297..8a7dab3e 100644 --- a/mlprec/impl/aggregator/mld_s_dec_aggregator_mat_asb.f90 +++ b/mlprec/impl/aggregator/mld_s_dec_aggregator_mat_asb.f90 @@ -143,15 +143,15 @@ subroutine mld_s_dec_aggregator_mat_asb(ag,parms,a,desc_a,& call op_prol%mv_to(tmpcoo) nzl = tmpcoo%get_nzeros() call psb_loc_to_glob(tmpcoo%ja(1:nzl),desc_ac,info,'I') + call tmpcoo%set_ncols(i_nr) call op_prol%mv_from(tmpcoo) call op_restr%mv_to(tmpcoo) nzl = tmpcoo%get_nzeros() call psb_loc_to_glob(tmpcoo%ia(1:nzl),desc_ac,info,'I') + call tmpcoo%set_nrows(i_nr) call op_restr%mv_from(tmpcoo) - call op_prol%set_ncols(i_nr) - call op_restr%set_nrows(i_nr) call psb_gather(tmp_ac,ac,desc_ac,info,root=-ione,& & dupl=psb_dupl_add_,keeploc=.false.) diff --git a/mlprec/impl/aggregator/mld_z_dec_aggregator_mat_asb.f90 b/mlprec/impl/aggregator/mld_z_dec_aggregator_mat_asb.f90 index d0c50884..133fe8a0 100644 --- a/mlprec/impl/aggregator/mld_z_dec_aggregator_mat_asb.f90 +++ b/mlprec/impl/aggregator/mld_z_dec_aggregator_mat_asb.f90 @@ -143,15 +143,15 @@ subroutine mld_z_dec_aggregator_mat_asb(ag,parms,a,desc_a,& call op_prol%mv_to(tmpcoo) nzl = tmpcoo%get_nzeros() call psb_loc_to_glob(tmpcoo%ja(1:nzl),desc_ac,info,'I') + call tmpcoo%set_ncols(i_nr) call op_prol%mv_from(tmpcoo) call op_restr%mv_to(tmpcoo) nzl = tmpcoo%get_nzeros() call psb_loc_to_glob(tmpcoo%ia(1:nzl),desc_ac,info,'I') + call tmpcoo%set_nrows(i_nr) call op_restr%mv_from(tmpcoo) - call op_prol%set_ncols(i_nr) - call op_restr%set_nrows(i_nr) call psb_gather(tmp_ac,ac,desc_ac,info,root=-ione,& & dupl=psb_dupl_add_,keeploc=.false.)