diff --git a/psb_dbjac_bld.f90 b/psb_dbjac_bld.f90 index 214081a7..f1e52379 100644 --- a/psb_dbjac_bld.f90 +++ b/psb_dbjac_bld.f90 @@ -111,6 +111,15 @@ subroutine psb_dbjac_bld(a,desc_a,p,upd,info) call psb_asmatbld(p%iprcparm(p_type_),p%iprcparm(n_ovr_),a,& & blck,desc_a,upd,p%desc_data,info,outfmt=coofmt) + if (debugprt) then + open(60+me) + call psb_csprt(60+me,a,head='% A') + close(60+me) + open(70+me) + call psb_csprt(70+me,blck,head='% BLCK') + close(70+me) + endif + if(info/=0) then call psb_errpush(4010,name,a_err='psb_asmatbld') goto 9999 @@ -145,7 +154,7 @@ subroutine psb_dbjac_bld(a,desc_a,p,upd,info) call psb_ipcoo2csr(p%av(ap_nd_),info) if(info /= 0) then - call psb_errpush(4010,name,a_err='psb_ipcoo2csr') + call psb_errpush(4010,name,a_err='psb_ipcoo2csr 1') goto 9999 end if @@ -177,7 +186,7 @@ subroutine psb_dbjac_bld(a,desc_a,p,upd,info) call psb_ipcoo2csr(atmp,info) if (info /= 0) then - call psb_errpush(4010,name,a_err='psb_ipcoo2csr') + call psb_errpush(4010,name,a_err='psb_ipcoo2csr 2') goto 9999 end if @@ -208,7 +217,7 @@ subroutine psb_dbjac_bld(a,desc_a,p,upd,info) call psb_ipcoo2csr(atmp,info) if (info /= 0) then - call psb_errpush(4010,name,a_err='psb_ipcoo2csr') + call psb_errpush(4010,name,a_err='psb_ipcoo2csr 3') goto 9999 end if @@ -278,10 +287,9 @@ subroutine psb_dbjac_bld(a,desc_a,p,upd,info) ! Output in COO format. call psb_sp_clip(atmp,p%av(ap_nd_),info,& & jmin=atmp%m+1,rscale=.false.,cscale=.false.) - call psb_ipcoo2csr(p%av(ap_nd_),info) if(info /= 0) then - call psb_errpush(4010,name,a_err='psb_ipcoo2csr') + call psb_errpush(4010,name,a_err='psb_ipcoo2csr 4') goto 9999 end if @@ -300,7 +308,7 @@ subroutine psb_dbjac_bld(a,desc_a,p,upd,info) call psb_ipcoo2csr(blck,info,rwshr=.true.) if(info/=0) then - call psb_errpush(4010,name,a_err='psb_ipcoo2csr') + call psb_errpush(4010,name,a_err='psb_ipcoo2csr 5') goto 9999 end if @@ -349,7 +357,7 @@ subroutine psb_dbjac_bld(a,desc_a,p,upd,info) call psb_ipcoo2csr(p%av(ap_nd_),info) if(info /= 0) then - call psb_errpush(4010,name,a_err='psb_ipcoo2csr') + call psb_errpush(4010,name,a_err='psb_ipcoo2csr 6') goto 9999 end if @@ -400,7 +408,7 @@ subroutine psb_dbjac_bld(a,desc_a,p,upd,info) call psb_ipcoo2csr(p%av(ap_nd_),info) if(info /= 0) then - call psb_errpush(4010,name,a_err='psb_ipcoo2csr') + call psb_errpush(4010,name,a_err='psb_ipcoo2csr 7') goto 9999 end if @@ -455,7 +463,7 @@ subroutine psb_dbjac_bld(a,desc_a,p,upd,info) call psb_ipcoo2csr(p%av(ap_nd_),info) if(info /= 0) then - call psb_errpush(4010,name,a_err='psb_ipcoo2csr') + call psb_errpush(4010,name,a_err='psb_ipcoo2csr 8') goto 9999 end if diff --git a/psb_dprecbld.f90 b/psb_dprecbld.f90 index 295f5ebc..cfa4be5f 100644 --- a/psb_dprecbld.f90 +++ b/psb_dprecbld.f90 @@ -52,7 +52,7 @@ subroutine psb_dprecbld(a,desc_a,p,info,upd) integer :: int_err(5) character :: iupd - logical, parameter :: debug=.false. + logical, parameter :: debug=.false., dump_aggr=.false. integer,parameter :: iroot=0,iout=60,ilout=40 character(len=20) :: name, ch_err @@ -126,11 +126,16 @@ subroutine psb_dprecbld(a,desc_a,p,info,upd) call psb_mlprc_bld(p%baseprecv(i-1)%base_a,p%baseprecv(i-1)%base_desc,& & p%baseprecv(i),info) + if (dump_aggr) then + call psb_csprt(90+me,p%baseprecv(i)%base_a,head='% Smoothed aggregate.') + end if + if (info /= 0) then info=4010 call psb_errpush(info,name) goto 9999 endif + if (debug) then write(0,*) 'Return from ',i-1,' call to mlprcbld ',info endif