Merge (minimal) changes from SLUDist branch.

Make.inc

@@ -35,7 +35,9 @@ BLACS=-lmpiblacs-gcc42 -L$(HOME)/LIB
 ##########################################################
 SLUDIR=/usr/local/SuperLU_3.0
 SLU=-lslu_lx_gcc42 -L$(SLUDIR)
-SLUDEF=-DHave_SLU_ -I$(SLUDIR)
+SLUDistDIR=/usr/local/SLUDist_2.0
+SLUDist=-lslud-bck-gcc42 -L$(SLUDistDIR) -L $(HOME)/LIB
+SLUDEF=-DHave_SLU_ -I$(SLUDIR) -DHave_SLUDist_ -I$(SLUDistDIR)
 
 UMFDIR=$(HOME)/LIB/Umfpack_gcc41
 UMF=-lumfpack -lamd -L$(UMFDIR)
@@ -45,7 +47,7 @@ UMFDEF=-DHave_UMF_ -I$(UMFDIR)
 #  check URL http://www.cs.umn.edu:~karypis
 #  
 METIS_LIB = -L$(HOME)/NUMERICAL/metis-4.0 -lmetis 
-LDLIBS=$(BLACS) $(SLU) $(UMF) $(BLAS) $(METIS_LIB)
+LDLIBS=$(BLACS) $(SLU) $(SLUDist) $(UMF) $(BLAS) $(METIS_LIB)
 
 #  Add -DLargeFptr for 64-bit addresses
 CDEFINES=-DAdd_ $(SLUDEF) $(UMFDEF)

Make.inc.gcc42

@@ -35,7 +35,9 @@ BLACS=-lmpiblacs-gcc42 -L$(HOME)/LIB
 ##########################################################
 SLUDIR=/usr/local/SuperLU_3.0
 SLU=-lslu_lx_gcc42 -L$(SLUDIR)
-SLUDEF=-DHave_SLU_ -I$(SLUDIR)
+SLUDistDIR=/usr/local/SLUDist_2.0
+SLUDist=-lslud-bck-gcc42 -L$(SLUDistDIR) -L $(HOME)/LIB
+SLUDEF=-DHave_SLU_ -I$(SLUDIR) -DHave_SLUDist_ -I$(SLUDistDIR)
 
 UMFDIR=$(HOME)/LIB/Umfpack_gcc41
 UMF=-lumfpack -lamd -L$(UMFDIR)
@@ -45,7 +47,7 @@ UMFDEF=-DHave_UMF_ -I$(UMFDIR)
 #  check URL http://www.cs.umn.edu:~karypis
 #  
 METIS_LIB = -L$(HOME)/NUMERICAL/metis-4.0 -lmetis 
-LDLIBS=$(BLACS) $(SLU) $(UMF) $(BLAS) $(METIS_LIB)
+LDLIBS=$(BLACS) $(SLU) $(SLUDist) $(UMF) $(BLAS) $(METIS_LIB)
 
 #  Add -DLargeFptr for 64-bit addresses
 CDEFINES=-DAdd_ $(SLUDEF) $(UMFDEF)

Makefile

@@ -1,6 +1,6 @@
 include Make.inc
-#PREC=../mld2p4-dev
-PREC=prec
+PREC=../mld2p4-dev
+#PREC=prec
 
 library:
 	( [ -d lib ] || mkdir lib)

base/internals/psi_zswapdata.f90

@@ -165,7 +165,8 @@ subroutine psi_zswapdatam(flag,n,beta,y,desc_a,work,info,data)
     end if
   end if
 
-
+  idxr = max(idxr,1)
+  idxs = max(idxs,1)
   if((idxr+idxs) < size(work)) then
     sndbuf => work(1:idxs)
     rcvbuf => work(idxs+1:idxs+idxr)
@@ -581,7 +582,8 @@ subroutine psi_zswapdatav(flag,beta,y,desc_a,work,info,data)
   end if
 
   
-
+  idxr = max(idxr,1)
+  idxs = max(idxs,1)
   if((idxr+idxs) < size(work)) then
     sndbuf => work(1:idxs)
     rcvbuf => work(idxs+1:idxs+idxr)

base/modules/psb_realloc_mod.F90

@@ -1213,7 +1213,7 @@ Contains
     if (allocated(vin)) then 
       call move_alloc(vin,vout)
     else if (allocated(vout)) then 
-      write(0,*) 'transfer: Clearing output'
+!!$      write(0,*) 'transfer: Clearing output'
       deallocate(vout)
     end if
 
@@ -1361,7 +1361,7 @@ Contains
     if (allocated(vin)) then 
       call move_alloc(vin,vout)
     else if (allocated(vout)) then 
-      write(0,*) 'transfer: Clearing output'
+!!$      write(0,*) 'transfer: Clearing output'
       deallocate(vout)
     end if
 #else