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

Merge branch 'development' of https://github.com/sfilippone/mld2p4-2 into development

Browse Source
stopcriterion
Salvatore Filippone 7 years ago
parent 46355253a5 a7bd1fb466
commit cbe0a5754f
6 changed files with 40 additions and 36 deletions
Show Stats Download Patch File Download Diff File

18
mlprec/impl/mld_cmlprec_aply.f90
Unescape Escape View File

@ -520,7 +520,7 @@ contains
call psb_map_X2Y(cone,vx2l,& call psb_map_X2Y(cone,vx2l,&
& czero,p%precv(level+1)%wrk%vx2l,& & czero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -540,7 +540,7 @@ contains
call psb_map_Y2X(cone,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(cone,p%precv(level+1)%wrk%vy2l,&
& cone,vy2l,& & cone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during prolongation') & a_err='Error during prolongation')
@ -653,7 +653,7 @@ contains
call psb_map_X2Y(cone,vty,& call psb_map_X2Y(cone,vty,&
& czero,p%precv(level+1)%wrk%vx2l,& & czero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -664,7 +664,7 @@ contains
call psb_map_X2Y(cone,vx2l,& call psb_map_X2Y(cone,vx2l,&
& czero,p%precv(level+1)%wrk%vx2l,& & czero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -680,7 +680,7 @@ contains
call psb_map_Y2X(cone,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(cone,p%precv(level+1)%wrk%vy2l,&
& cone,vy2l,& & cone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during prolongation') & a_err='Error during prolongation')
@ -697,7 +697,7 @@ contains
if (info == psb_success_) call psb_map_X2Y(cone,vty,& if (info == psb_success_) call psb_map_X2Y(cone,vty,&
& czero,p%precv(level+1)%wrk%vx2l,& & czero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during W-cycle restriction') & a_err='Error during W-cycle restriction')
@ -709,7 +709,7 @@ contains
if (info == psb_success_) call psb_map_Y2X(cone,p%precv(level+1)%wrk%vy2l,& if (info == psb_success_) call psb_map_Y2X(cone,p%precv(level+1)%wrk%vy2l,&
& cone,vy2l,& & cone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
@ -889,7 +889,7 @@ contains
call psb_map_X2Y(cone,vty,& call psb_map_X2Y(cone,vty,&
& czero,p%precv(level + 1)%wrk%vx2l,& & czero,p%precv(level + 1)%wrk%vx2l,&
& p%precv(level + 1)%map,info,work=work,& & p%precv(level + 1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
@ -925,7 +925,7 @@ contains
call psb_map_Y2X(cone,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(cone,p%precv(level+1)%wrk%vy2l,&
& cone,vy2l,& & cone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&

18
mlprec/impl/mld_dmlprec_aply.f90
Unescape Escape View File

@ -520,7 +520,7 @@ contains
call psb_map_X2Y(done,vx2l,& call psb_map_X2Y(done,vx2l,&
& dzero,p%precv(level+1)%wrk%vx2l,& & dzero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -540,7 +540,7 @@ contains
call psb_map_Y2X(done,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(done,p%precv(level+1)%wrk%vy2l,&
& done,vy2l,& & done,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during prolongation') & a_err='Error during prolongation')
@ -653,7 +653,7 @@ contains
call psb_map_X2Y(done,vty,& call psb_map_X2Y(done,vty,&
& dzero,p%precv(level+1)%wrk%vx2l,& & dzero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -664,7 +664,7 @@ contains
call psb_map_X2Y(done,vx2l,& call psb_map_X2Y(done,vx2l,&
& dzero,p%precv(level+1)%wrk%vx2l,& & dzero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -680,7 +680,7 @@ contains
call psb_map_Y2X(done,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(done,p%precv(level+1)%wrk%vy2l,&
& done,vy2l,& & done,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during prolongation') & a_err='Error during prolongation')
@ -697,7 +697,7 @@ contains
if (info == psb_success_) call psb_map_X2Y(done,vty,& if (info == psb_success_) call psb_map_X2Y(done,vty,&
& dzero,p%precv(level+1)%wrk%vx2l,& & dzero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during W-cycle restriction') & a_err='Error during W-cycle restriction')
@ -709,7 +709,7 @@ contains
if (info == psb_success_) call psb_map_Y2X(done,p%precv(level+1)%wrk%vy2l,& if (info == psb_success_) call psb_map_Y2X(done,p%precv(level+1)%wrk%vy2l,&
& done,vy2l,& & done,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
@ -889,7 +889,7 @@ contains
call psb_map_X2Y(done,vty,& call psb_map_X2Y(done,vty,&
& dzero,p%precv(level + 1)%wrk%vx2l,& & dzero,p%precv(level + 1)%wrk%vx2l,&
& p%precv(level + 1)%map,info,work=work,& & p%precv(level + 1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
@ -925,7 +925,7 @@ contains
call psb_map_Y2X(done,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(done,p%precv(level+1)%wrk%vy2l,&
& done,vy2l,& & done,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&

2
mlprec/impl/mld_dumf_interface.c
Unescape Escape View File

@ -106,6 +106,7 @@ int mld_dumf_fact(int n, int nnz,
info = 0; info = 0;
} else { } else {
printf("umfpack_di_symbolic() error returns INFO= %d\n", info); printf("umfpack_di_symbolic() error returns INFO= %d\n", info);
umfpack_di_report_status(Control,info);
*symptr = (void *) NULL; *symptr = (void *) NULL;
*numptr = (void *) NULL; *numptr = (void *) NULL;
return -11; return -11;
@ -127,6 +128,7 @@ int mld_dumf_fact(int n, int nnz,
} else { } else {
printf("umfpack_di_numeric() error returns INFO= %d\n", info); printf("umfpack_di_numeric() error returns INFO= %d\n", info);
umfpack_di_report_status(Control,info);
info = -12; info = -12;
*numptr = NULL; *numptr = NULL;
} }

18
mlprec/impl/mld_smlprec_aply.f90
Unescape Escape View File

@ -520,7 +520,7 @@ contains
call psb_map_X2Y(sone,vx2l,& call psb_map_X2Y(sone,vx2l,&
& szero,p%precv(level+1)%wrk%vx2l,& & szero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -540,7 +540,7 @@ contains
call psb_map_Y2X(sone,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(sone,p%precv(level+1)%wrk%vy2l,&
& sone,vy2l,& & sone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during prolongation') & a_err='Error during prolongation')
@ -653,7 +653,7 @@ contains
call psb_map_X2Y(sone,vty,& call psb_map_X2Y(sone,vty,&
& szero,p%precv(level+1)%wrk%vx2l,& & szero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -664,7 +664,7 @@ contains
call psb_map_X2Y(sone,vx2l,& call psb_map_X2Y(sone,vx2l,&
& szero,p%precv(level+1)%wrk%vx2l,& & szero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -680,7 +680,7 @@ contains
call psb_map_Y2X(sone,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(sone,p%precv(level+1)%wrk%vy2l,&
& sone,vy2l,& & sone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during prolongation') & a_err='Error during prolongation')
@ -697,7 +697,7 @@ contains
if (info == psb_success_) call psb_map_X2Y(sone,vty,& if (info == psb_success_) call psb_map_X2Y(sone,vty,&
& szero,p%precv(level+1)%wrk%vx2l,& & szero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during W-cycle restriction') & a_err='Error during W-cycle restriction')
@ -709,7 +709,7 @@ contains
if (info == psb_success_) call psb_map_Y2X(sone,p%precv(level+1)%wrk%vy2l,& if (info == psb_success_) call psb_map_Y2X(sone,p%precv(level+1)%wrk%vy2l,&
& sone,vy2l,& & sone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
@ -889,7 +889,7 @@ contains
call psb_map_X2Y(sone,vty,& call psb_map_X2Y(sone,vty,&
& szero,p%precv(level + 1)%wrk%vx2l,& & szero,p%precv(level + 1)%wrk%vx2l,&
& p%precv(level + 1)%map,info,work=work,& & p%precv(level + 1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
@ -925,7 +925,7 @@ contains
call psb_map_Y2X(sone,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(sone,p%precv(level+1)%wrk%vy2l,&
& sone,vy2l,& & sone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&

18
mlprec/impl/mld_zmlprec_aply.f90
Unescape Escape View File

@ -520,7 +520,7 @@ contains
call psb_map_X2Y(zone,vx2l,& call psb_map_X2Y(zone,vx2l,&
& zzero,p%precv(level+1)%wrk%vx2l,& & zzero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -540,7 +540,7 @@ contains
call psb_map_Y2X(zone,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(zone,p%precv(level+1)%wrk%vy2l,&
& zone,vy2l,& & zone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during prolongation') & a_err='Error during prolongation')
@ -653,7 +653,7 @@ contains
call psb_map_X2Y(zone,vty,& call psb_map_X2Y(zone,vty,&
& zzero,p%precv(level+1)%wrk%vx2l,& & zzero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -664,7 +664,7 @@ contains
call psb_map_X2Y(zone,vx2l,& call psb_map_X2Y(zone,vx2l,&
& zzero,p%precv(level+1)%wrk%vx2l,& & zzero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during restriction') & a_err='Error during restriction')
@ -680,7 +680,7 @@ contains
call psb_map_Y2X(zone,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(zone,p%precv(level+1)%wrk%vy2l,&
& zone,vy2l,& & zone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during prolongation') & a_err='Error during prolongation')
@ -697,7 +697,7 @@ contains
if (info == psb_success_) call psb_map_X2Y(zone,vty,& if (info == psb_success_) call psb_map_X2Y(zone,vty,&
& zzero,p%precv(level+1)%wrk%vx2l,& & zzero,p%precv(level+1)%wrk%vx2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
& a_err='Error during W-cycle restriction') & a_err='Error during W-cycle restriction')
@ -709,7 +709,7 @@ contains
if (info == psb_success_) call psb_map_Y2X(zone,p%precv(level+1)%wrk%vy2l,& if (info == psb_success_) call psb_map_Y2X(zone,p%precv(level+1)%wrk%vy2l,&
& zone,vy2l,& & zone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
@ -889,7 +889,7 @@ contains
call psb_map_X2Y(zone,vty,& call psb_map_X2Y(zone,vty,&
& zzero,p%precv(level + 1)%wrk%vx2l,& & zzero,p%precv(level + 1)%wrk%vx2l,&
& p%precv(level + 1)%map,info,work=work,& & p%precv(level + 1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=wv(1),vty=p%precv(level+1)%wrk%wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&
@ -925,7 +925,7 @@ contains
call psb_map_Y2X(zone,p%precv(level+1)%wrk%vy2l,& call psb_map_Y2X(zone,p%precv(level+1)%wrk%vy2l,&
& zone,vy2l,& & zone,vy2l,&
& p%precv(level+1)%map,info,work=work,& & p%precv(level+1)%map,info,work=work,&
& vtx=wv(1),vty=wv(2)) & vtx=p%precv(level+1)%wrk%wv(1),vty=wv(1))
if (info /= psb_success_) then if (info /= psb_success_) then
call psb_errpush(psb_err_internal_error_,name,& call psb_errpush(psb_err_internal_error_,name,&

2
mlprec/impl/mld_zumf_interface.c
Unescape Escape View File

@ -107,6 +107,7 @@ int mld_zumf_fact(int n, int nnz,
info = 0; info = 0;
} else { } else {
printf("umfpack_zi_symbolic() error returns INFO= %d\n", info); printf("umfpack_zi_symbolic() error returns INFO= %d\n", info);
umfpack_zi_report_status(Control,info);
*symptr = (void *) NULL; *symptr = (void *) NULL;
*numptr = (void *) NULL; *numptr = (void *) NULL;
return -11; return -11;
@ -128,6 +129,7 @@ int mld_zumf_fact(int n, int nnz,
} else { } else {
printf("umfpack_zi_numeric() error returns INFO= %d\n", info); printf("umfpack_zi_numeric() error returns INFO= %d\n", info);
umfpack_zi_report_status(Control,info);
info = -12; info = -12;
*numptr = NULL; *numptr = NULL;
} }

Write Preview
Loading…
Cancel
Save