Added call to psb_ovrl at the end of the preconditioner application,

to handle data distributions with a "real" overlap.
psblas3-type-indexed
Salvatore Filippone 17 years ago
parent 83213e1fc5
commit c30ec19bc4

@ -73,7 +73,7 @@ subroutine psb_dprc_aply(prec,x,y,desc_data,info,trans, work)
end if end if
call psb_gprec_aply(done,prec,x,dzero,y,desc_data,trans_,work_,info) call psb_gprec_aply(done,prec,x,dzero,y,desc_data,trans_,work_,info)
call psb_ovrl(y,desc_data,info,update=psb_avg_)
if (present(work)) then if (present(work)) then
else else
deallocate(work_) deallocate(work_)

@ -74,6 +74,7 @@ subroutine psb_zprc_aply(prec,x,y,desc_data,info,trans, work)
end if end if
call psb_gprec_aply(zone,prec,x,zzero,y,desc_data,trans_, work_,info) call psb_gprec_aply(zone,prec,x,zzero,y,desc_data,trans_, work_,info)
call psb_ovrl(y,desc_data,info,update=psb_avg_)
if (present(work)) then if (present(work)) then
else else

Loading…
Cancel
Save