diff --git a/Changelog b/Changelog index 61bb7994..e0affe25 100644 --- a/Changelog +++ b/Changelog @@ -1,8 +1,15 @@ Changelog. A lot less detailed than usual, at least for past history. +2006/11/28: Merged the mods for descriptors of large index spaces to + avoid having the GLOB_TO_LOC array. Took the chance to + reorganize the descriptor build routines and define some + access functions for decsriptor features and entries, so + as not to use the descriptor components directly. Tested + with AS, 2- and 3- level Post smoothers. + 2006/11/09: The allocatable version works, but under gcc42 there is a - strange problem when using -fbounds-check. Be careful! + compiler bug when using -fbounds-check. 2006/11/08: Merged the allocatable version; hope everything works! diff --git a/src/modules/psb_spmat_type.f90 b/src/modules/psb_spmat_type.f90 index 5389b0bc..d1927f26 100644 --- a/src/modules/psb_spmat_type.f90 +++ b/src/modules/psb_spmat_type.f90 @@ -1295,6 +1295,11 @@ contains return endif irw = iaux + if (irw > a%m) then + write(0,*) 'SPINFO: Accessing out of bounds? ',irw,a%m + ires = 0 + return + endif if (toupper(a%fida) == 'CSR') then ires = a%ia2(irw+1)-a%ia2(irw) else if ((toupper(a%fida) == 'COO').or.(toupper(a%fida) == 'COI')) then diff --git a/src/tools/psb_dcdovrbld.f90 b/src/tools/psb_dcdovrbld.f90 index 1494ee96..d9b6d827 100644 --- a/src/tools/psb_dcdovrbld.f90 +++ b/src/tools/psb_dcdovrbld.f90 @@ -217,8 +217,7 @@ Subroutine psb_dcdovrbld(n_ovr,desc_p,desc_a,a,& idxs = 0 counter = 1 counter_t = 1 - - + t1 = mpi_wtime() Do While (halo(counter) /= -1) tot_elem=0 @@ -435,6 +434,13 @@ Subroutine psb_dcdovrbld(n_ovr,desc_p,desc_a,a,& call psb_errpush(info,name,a_err=ch_err) goto 9999 end if + if (.false.) then + open(70+me,position='append') + write(70+me,*) ' Workr on iteration ',i_ovr + write(70+me,'(8(i7,1x))') workr(1:iszr) + call flush(70+me) + close(70+me) + end if if (debug) write(0,*) 'ISZR :',iszr diff --git a/src/tools/psb_dsphalo.f90 b/src/tools/psb_dsphalo.f90 index eb77ff0f..1ea3b455 100644 --- a/src/tools/psb_dsphalo.f90 +++ b/src/tools/psb_dsphalo.f90 @@ -295,6 +295,7 @@ Subroutine psb_dsphalo(a,desc_a,blk,info,rwcnv,clcnv,outfmt) Enddo blk%fida='COO' blk%infoa(psb_nnz_)=l1 + if (debugprt) then open(50+me) call psb_csprt(50+me,blk,head='% SPHALO border .')