*** empty log message ***

psblas3-type-indexed
Alfredo Buttari 21 years ago
parent bbf9a4c53a
commit be2da4718d

@ -45,7 +45,7 @@ subroutine psb_dcslu(a,desc_a,p,upd,info)
external mpi_wtime external mpi_wtime
logical, parameter :: debugprt=.false., debug=.false., aggr_dump=.false. logical, parameter :: debugprt=.false., debug=.false., aggr_dump=.false.
integer istpb, istpe, ifctb, ifcte, err_act, irank, icomm, nztota, nztotb,& integer istpb, istpe, ifctb, ifcte, err_act, irank, icomm, nztota, nztotb,&
& nztmp, nzl, nnr, ir, mglob, mtype, n_row, nrow_a,n_col, nhalo,lovr & nztmp, nzl, nnr, ir, mglob, mtype, n_row, nrow_a,n_col, nhalo,lovr, ind, iind, pi,nr,ns
integer ::icontxt,nprow,npcol,me,mycol integer ::icontxt,nprow,npcol,me,mycol
character(len=20) :: name, ch_err character(len=20) :: name, ch_err
@ -176,8 +176,6 @@ subroutine psb_dcslu(a,desc_a,p,upd,info)
call blacs_barrier(icontxt,'All') call blacs_barrier(icontxt,'All')
endif endif
! if (me==0) write(0,*) 'setup time',t2-t1, blck%fida, p%iprcparm(p_type_),blck%m,upd
! write(0,'(i3,1x,a,4(1x,g14.5))') me,' setup time',t2-t1
if (p%iprcparm(iren_) > 0) then if (p%iprcparm(iren_) > 0) then
@ -208,6 +206,7 @@ subroutine psb_dcslu(a,desc_a,p,upd,info)
t3 = mpi_wtime() t3 = mpi_wtime()
if (debugprt) then if (debugprt) then
call blacs_barrier(icontxt,'All')
open(40+me) open(40+me)
call psb_csprt(40+me,atmp,head='% Local matrix') call psb_csprt(40+me,atmp,head='% Local matrix')
close(40+me) close(40+me)
@ -239,12 +238,12 @@ subroutine psb_dcslu(a,desc_a,p,upd,info)
else if (p%iprcparm(iren_) == 0) then else if (p%iprcparm(iren_) == 0) then
t3 = mpi_wtime() t3 = mpi_wtime()
! ierr = MPE_Log_event( ifctb, 0, "st SIMPLE" )
! This is where we have mo renumbering, thus no need ! This is where we have mo renumbering, thus no need
! for ATMP ! for ATMP
if (debugprt) then if (debugprt) then
open(40+me) open(40+me)
call blacs_barrier(icontxt,'All')
call psb_csprt(40+me,a,iv=p%desc_data%loc_to_glob,& call psb_csprt(40+me,a,iv=p%desc_data%loc_to_glob,&
& head='% Local matrix') & head='% Local matrix')
if (p%iprcparm(p_type_)==asm_) then if (p%iprcparm(p_type_)==asm_) then
@ -254,8 +253,6 @@ subroutine psb_dcslu(a,desc_a,p,upd,info)
close(40+me) close(40+me)
endif endif
t5= mpi_wtime() t5= mpi_wtime()
if (debug) write(0,*) me,' Going for dsplu' if (debug) write(0,*) me,' Going for dsplu'
call psb_dsplu(a,p%av(l_pr_),p%av(u_pr_),p%d,info,blck=blck) call psb_dsplu(a,p%av(l_pr_),p%av(u_pr_),p%d,info,blck=blck)

@ -143,10 +143,9 @@ Subroutine psb_dcsrsetup(ptype,novr,a,blk,desc_data,upd,desc_p,info,outfmt)
!!$ !!$
Call blacs_gridinfo(icontxt,nprow,npcol,me,mycol) Call blacs_gridinfo(icontxt,nprow,npcol,me,mycol)
If(debug)Write(0,*)'BEGIN dcsrsetup',me,upd If(debug)Write(0,*)'BEGIN dcsrsetup',me,upd,novr
t1 = mpi_wtime() t1 = mpi_wtime()
If (upd == 'F') Then If (upd == 'F') Then
! !
! Build the auiliary descriptor',desc_p%matrix_data(psb_n_row_) ! Build the auiliary descriptor',desc_p%matrix_data(psb_n_row_)

@ -239,9 +239,6 @@ subroutine psb_dprecbld(a,p,desc_a,info,upd)
call psb_check_def(p%baseprecv(2)%iprcparm(jac_sweeps_),'Jacobi sweeps',& call psb_check_def(p%baseprecv(2)%iprcparm(jac_sweeps_),'Jacobi sweeps',&
& 1,is_legal_jac_sweeps) & 1,is_legal_jac_sweeps)
call blacs_barrier(icontxt,'All') ! to be removed
write(0,'(i2," Calling mlprecbld")')me
call blacs_barrier(icontxt,'All') ! to be removed
call psb_mlprec_bld(a,desc_a,p%baseprecv(2),info) call psb_mlprec_bld(a,desc_a,p%baseprecv(2),info)
if(info /= 0) then if(info /= 0) then
info=4010 info=4010

Loading…
Cancel
Save