fdurastante
/
amg4psblas
mirror of https://github.com/sfilippone/amg4psblas.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

mld2p4:

Browse Source 
mlprec/impl/level/mld_c_base_onelev_check.f90
 mlprec/impl/level/mld_d_base_onelev_check.f90
 mlprec/impl/level/mld_s_base_onelev_check.f90
 mlprec/impl/level/mld_z_base_onelev_check.f90
 mlprec/impl/mld_c_hierarchy_bld.f90
 mlprec/impl/mld_d_hierarchy_bld.f90
 mlprec/impl/mld_s_hierarchy_bld.f90
 mlprec/impl/mld_z_hierarchy_bld.f90

Number of sweeps may be 0 to represent PRE/POST.
Correction factor to scale target aggregate size.
stopcriterion
Salvatore Filippone 9 years ago
parent 086ff71d48
commit 1295fdcc60
8 changed files with 16 additions and 12 deletions
Show Stats Download Patch File Download Diff File

6
mlprec/impl/level/mld_c_base_onelev_check.f90
Unescape Escape View File

@ -53,11 +53,11 @@ subroutine mld_c_base_onelev_check(lv,info)
info = psb_success_
call mld_check_def(lv%parms%sweeps,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
call mld_check_def(lv%parms%sweeps_pre,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
call mld_check_def(lv%parms%sweeps_post,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
if (allocated(lv%sm)) then

6
mlprec/impl/level/mld_d_base_onelev_check.f90
Unescape Escape View File

@ -53,11 +53,11 @@ subroutine mld_d_base_onelev_check(lv,info)
info = psb_success_
call mld_check_def(lv%parms%sweeps,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
call mld_check_def(lv%parms%sweeps_pre,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
call mld_check_def(lv%parms%sweeps_post,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
if (allocated(lv%sm)) then

6
mlprec/impl/level/mld_s_base_onelev_check.f90
Unescape Escape View File

@ -53,11 +53,11 @@ subroutine mld_s_base_onelev_check(lv,info)
info = psb_success_
call mld_check_def(lv%parms%sweeps,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
call mld_check_def(lv%parms%sweeps_pre,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
call mld_check_def(lv%parms%sweeps_post,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
if (allocated(lv%sm)) then

6
mlprec/impl/level/mld_z_base_onelev_check.f90
Unescape Escape View File

@ -53,11 +53,11 @@ subroutine mld_z_base_onelev_check(lv,info)
info = psb_success_
call mld_check_def(lv%parms%sweeps,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
call mld_check_def(lv%parms%sweeps_pre,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
call mld_check_def(lv%parms%sweeps_post,&
& 'Jacobi sweeps',ione,is_int_positive)
& 'Jacobi sweeps',ione,is_int_non_negative)
if (allocated(lv%sm)) then

1
mlprec/impl/mld_c_hierarchy_bld.f90
Unescape Escape View File

@ -202,6 +202,7 @@ subroutine mld_c_hierarchy_bld(a,desc_a,p,info,amold,vmold,imold)
casize = desc_a%get_global_rows()
casize = int((sone*casize)**(sone/(sone*3)),psb_ipk_)
casize = max(casize,ione)
casize = casize*40_psb_ipk_
end if
call mld_bld_mlhier_aggsize(casize,mxplevs,mnaggratio,a,desc_a,p%precv,info)
else

1
mlprec/impl/mld_d_hierarchy_bld.f90
Unescape Escape View File

@ -202,6 +202,7 @@ subroutine mld_d_hierarchy_bld(a,desc_a,p,info,amold,vmold,imold)
casize = desc_a%get_global_rows()
casize = int((done*casize)**(done/(done*3)),psb_ipk_)
casize = max(casize,ione)
casize = casize*40_psb_ipk_
end if
call mld_bld_mlhier_aggsize(casize,mxplevs,mnaggratio,a,desc_a,p%precv,info)
else

1
mlprec/impl/mld_s_hierarchy_bld.f90
Unescape Escape View File

@ -202,6 +202,7 @@ subroutine mld_s_hierarchy_bld(a,desc_a,p,info,amold,vmold,imold)
casize = desc_a%get_global_rows()
casize = int((sone*casize)**(sone/(sone*3)),psb_ipk_)
casize = max(casize,ione)
casize = casize*40_psb_ipk_
end if
call mld_bld_mlhier_aggsize(casize,mxplevs,mnaggratio,a,desc_a,p%precv,info)
else

1
mlprec/impl/mld_z_hierarchy_bld.f90
Unescape Escape View File

@ -202,6 +202,7 @@ subroutine mld_z_hierarchy_bld(a,desc_a,p,info,amold,vmold,imold)
casize = desc_a%get_global_rows()
casize = int((done*casize)**(done/(done*3)),psb_ipk_)
casize = max(casize,ione)
casize = casize*40_psb_ipk_
end if
call mld_bld_mlhier_aggsize(casize,mxplevs,mnaggratio,a,desc_a,p%precv,info)
else

Write Preview
Loading…
Cancel
Save