Changed psb_const_mod to include psb_const.fh, to have constants

defined in just one place.

src/serial/csr/dcocrupd.f

@@ -56,7 +56,6 @@ C     .. Local scalars ..
       INTEGER           J, NNZ, IP1, NNZI
 C     .. Local arrays ..
       IERROR = 0
-c$$$      write(0,*) 'dcocrupd ',infoa(upd_),ibits(infoa(upd_),2,1)
       IF (IBITS(INFOA(PSB_UPD_),2,1).EQ.1) THEN 
 C
 C     Smart update capability

src/serial/csr/dcrcrupd.f

@@ -57,7 +57,6 @@ C     .. Local scalars ..
      +  NRC, IPH, JPH, JPA, LPA, IRET, LNK, NNZ, IP1
 C     .. Local arrays ..
       IERROR = 0
-c$$$      write(0,*) 'dcrcrupd ',infoa(upd_),ibits(infoa(upd_),2,1)
       IF (IBITS(INFOA(PSB_UPD_),2,1).EQ.1) THEN 
 C
 C     Smart update capability

src/serial/psb_dcoins.f90

@@ -108,7 +108,7 @@ subroutine psb_dcoins(nz,ia,ja,val,a,gtl,imin,imax,jmin,jmax,info)
      endif
 
      if ((nza+nz)>isza) then 
-        call psb_sp_reall(a,nza+nz,info)
+        call psb_sp_reall(a,max(nza+nz,int(1.5*isza)),info)
         if(info.ne.izero) then
            info=4010
            ch_err='psb_sp_reall'