fdurastante/amg4psblas
1
0
Fork 0
mirror of https://github.com/sfilippone/amg4psblas.git synced 2026-10-06 14:44:55 +00:00
Code Issues Packages Projects Releases Wiki Activity

Merge branch 'unify_aggr_bld' into development

Browse Source
This commit is contained in:
Salvatore Filippone
2020-02-20 12:30:24 +00:00
parent 497984f9e3 7f3f0cb2e5
commit 343a0e3661
8 changed files with 44 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files
mlprec/impl/aggregator/mld_c_dec_aggregator_mat_asb.f90
+7
View File
@@ -191,7 +191,14 @@ subroutine mld_c_dec_aggregator_mat_asb(ag,parms,a,desc_a,ilaggr,nlaggr,&
!
! If we are here, it means we assume that an IPK version of the
! coarse matrix can hold all indices. User beware!
!
!
! op_prol/op_restr come from par_spmm_bld with local sizes
! suitable for DIST option, fix relevant sizes
!
call op_prol%set_ncols(ntaggr)
call op_restr%set_nrows(ntaggr)
call psb_cdall(ictxt,desc_ac,info,mg=ntaggr,repl=.true.)
if (info == psb_success_) call psb_cdasb(desc_ac,info)
if (info == psb_success_) call ac%mv_to(acoo)
mlprec/impl/aggregator/mld_caggrmat_nosmth_bld.f90
+4
View File
@@ -146,6 +146,9 @@ subroutine mld_caggrmat_nosmth_bld(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
naggrp1 = sum(nlaggr(1:me+1))
if (.false.) then
!
! If we ever want to switch to explicit product when building this..
!
call a%cp_to(acsr)
call op_prol%mv_to(coo_prol)
inaggr = naggr
@@ -161,6 +164,7 @@ subroutine mld_caggrmat_nosmth_bld(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
else
call a%cp_to(ac_coo)
nzt = ac_coo%get_nzeros()
k = 0
mlprec/impl/aggregator/mld_d_dec_aggregator_mat_asb.f90
+7
View File
@@ -191,7 +191,14 @@ subroutine mld_d_dec_aggregator_mat_asb(ag,parms,a,desc_a,ilaggr,nlaggr,&
!
! If we are here, it means we assume that an IPK version of the
! coarse matrix can hold all indices. User beware!
!
!
! op_prol/op_restr come from par_spmm_bld with local sizes
! suitable for DIST option, fix relevant sizes
!
call op_prol%set_ncols(ntaggr)
call op_restr%set_nrows(ntaggr)
call psb_cdall(ictxt,desc_ac,info,mg=ntaggr,repl=.true.)
if (info == psb_success_) call psb_cdasb(desc_ac,info)
if (info == psb_success_) call ac%mv_to(acoo)
mlprec/impl/aggregator/mld_daggrmat_nosmth_bld.f90
+4
View File
@@ -146,6 +146,9 @@ subroutine mld_daggrmat_nosmth_bld(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
naggrp1 = sum(nlaggr(1:me+1))
if (.false.) then
!
! If we ever want to switch to explicit product when building this..
!
call a%cp_to(acsr)
call op_prol%mv_to(coo_prol)
inaggr = naggr
@@ -161,6 +164,7 @@ subroutine mld_daggrmat_nosmth_bld(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
else
call a%cp_to(ac_coo)
nzt = ac_coo%get_nzeros()
k = 0
mlprec/impl/aggregator/mld_s_dec_aggregator_mat_asb.f90
+7
View File
@@ -191,7 +191,14 @@ subroutine mld_s_dec_aggregator_mat_asb(ag,parms,a,desc_a,ilaggr,nlaggr,&
!
! If we are here, it means we assume that an IPK version of the
! coarse matrix can hold all indices. User beware!
!
!
! op_prol/op_restr come from par_spmm_bld with local sizes
! suitable for DIST option, fix relevant sizes
!
call op_prol%set_ncols(ntaggr)
call op_restr%set_nrows(ntaggr)
call psb_cdall(ictxt,desc_ac,info,mg=ntaggr,repl=.true.)
if (info == psb_success_) call psb_cdasb(desc_ac,info)
if (info == psb_success_) call ac%mv_to(acoo)
mlprec/impl/aggregator/mld_saggrmat_nosmth_bld.f90
+4
View File
@@ -146,6 +146,9 @@ subroutine mld_saggrmat_nosmth_bld(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
naggrp1 = sum(nlaggr(1:me+1))
if (.false.) then
!
! If we ever want to switch to explicit product when building this..
!
call a%cp_to(acsr)
call op_prol%mv_to(coo_prol)
inaggr = naggr
@@ -161,6 +164,7 @@ subroutine mld_saggrmat_nosmth_bld(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
else
call a%cp_to(ac_coo)
nzt = ac_coo%get_nzeros()
k = 0
mlprec/impl/aggregator/mld_z_dec_aggregator_mat_asb.f90
+7
View File
@@ -191,7 +191,14 @@ subroutine mld_z_dec_aggregator_mat_asb(ag,parms,a,desc_a,ilaggr,nlaggr,&
!
! If we are here, it means we assume that an IPK version of the
! coarse matrix can hold all indices. User beware!
!
!
! op_prol/op_restr come from par_spmm_bld with local sizes
! suitable for DIST option, fix relevant sizes
!
call op_prol%set_ncols(ntaggr)
call op_restr%set_nrows(ntaggr)
call psb_cdall(ictxt,desc_ac,info,mg=ntaggr,repl=.true.)
if (info == psb_success_) call psb_cdasb(desc_ac,info)
if (info == psb_success_) call ac%mv_to(acoo)
mlprec/impl/aggregator/mld_zaggrmat_nosmth_bld.f90
+4
View File
@@ -146,6 +146,9 @@ subroutine mld_zaggrmat_nosmth_bld(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
naggrp1 = sum(nlaggr(1:me+1))
if (.false.) then
!
! If we ever want to switch to explicit product when building this..
!
call a%cp_to(acsr)
call op_prol%mv_to(coo_prol)
inaggr = naggr
@@ -161,6 +164,7 @@ subroutine mld_zaggrmat_nosmth_bld(a,desc_a,ilaggr,nlaggr,parms,ac,op_prol,op_re
else
call a%cp_to(ac_coo)
nzt = ac_coo%get_nzeros()
k = 0