Fix documentation for sprn and usage of reinit in sphalo.

docs/pdf/toolsrout.tex

@@ -507,7 +507,7 @@ Specified as: an integer variable.
 %
 \subroutine{psb\_sprn}{Reinit sparse matrix structure for psblas routines.}
 
-\syntax{call psb\_sprn}{a, decsc\_a, info}
+\syntax{call psb\_sprn}{a, decsc\_a, info, clear}
 
 \begin{description}
 \item[\bf On Entry]
@@ -519,6 +519,10 @@ Specified as: a structured data of type \spdata.
 Scope:{\bf local}.\\
 Type:{\bf required}.\\
 Specified as: a structured data of type \descdata.
+\item[clear] Choose whether to zero out matrix coefficients\\
+Scope:{\bf local}.\\
+Type:{\bf optional}.\\
+Default: true.
 \end{description}
 
 \begin{description}

src/modules/psb_spmat_type.f90

@@ -173,7 +173,7 @@ contains
       ! right after allocate, with spins doing the right thing.
       ! hopefully :-)
 
-    case( psb_spmat_upd_) 
+    case(psb_spmat_upd_) 
 
     case default
       info=591     

src/tools/psb_dsphalo.f90

@@ -194,17 +194,18 @@ Subroutine psb_dsphalo(a,desc_a,blk,info,rwcnv,clcnv,outfmt)
   end if
   mat_recv = iszr
   iszs=sum(sdsz)
+  call psb_nullify_sp(tmp)
   call psb_sp_all(0,0,tmp,max(iszs,1),info)
   tmp%fida='COO'
-  
+  call psb_sp_setifld(psb_spmat_asb_,psb_state_,tmp,info)
+
   t2 = mpi_wtime()
 
   l1  = 0
   ipx = 1
   counter=1
   idx = 0
-  call psb_sp_reinit(tmp,info)
-  tmp%infoa(psb_nnz_) = 0 
+
   Do 
     proc=desc_a%halo_index(counter)
     if (proc == -1) exit 

src/tools/psb_zsphalo.f90

@@ -194,17 +194,18 @@ Subroutine psb_zsphalo(a,desc_a,blk,info,rwcnv,clcnv,outfmt)
   end if
   mat_recv = iszr
   iszs=sum(sdsz)
+  call psb_nullify_sp(tmp)
   call psb_sp_all(0,0,tmp,max(iszs,1),info)
   tmp%fida='COO'
-  
+  call psb_sp_setifld(psb_spmat_asb_,psb_state_,tmp,info)
+
   t2 = mpi_wtime()
 
   l1  = 0
   ipx = 1
   counter=1
   idx = 0
-  call psb_sp_reinit(tmp,info)
-  tmp%infoa(psb_nnz_) = 0 
+
   Do 
     proc=desc_a%halo_index(counter)
     if (proc == -1) exit