psblas3:

base/comm/psb_cgather.f90
 base/comm/psb_dgather.f90
 base/comm/psb_sgather.f90
 base/comm/psb_zgather.f90
 base/modules/psb_c_comm_mod.f90
 base/modules/psb_d_comm_mod.f90
 base/modules/psb_s_comm_mod.f90
 base/modules/psb_z_comm_mod.f90


Fixed ALLOCATABLE attribute of globx.
Also updated docs.

base/comm/psb_cgather.f90

@@ -50,7 +50,7 @@ subroutine  psb_cgatherm(globx, locx, desc_a, info, iroot)
   implicit none
 
   complex(psb_spk_), intent(in)    :: locx(:,:)
-  complex(psb_spk_), intent(out)   :: globx(:,:)
+  complex(psb_spk_), intent(out), allocatable :: globx(:,:)
   type(psb_desc_type), intent(in) :: desc_a
   integer(psb_ipk_), intent(out)            :: info
   integer(psb_ipk_), intent(in), optional   :: iroot
@@ -100,16 +100,13 @@ subroutine  psb_cgatherm(globx, locx, desc_a, info, iroot)
   ilocx = 1
   jlocx = 1
 
-  lda_globx = size(globx,1)
-  lda_locx  = size(locx, 1)
-
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
-
-  lock=size(locx,2)-jlocx+1
-  globk=size(globx,2)-jglobx+1
-  maxk=min(lock,globk)
-  k = maxk
+  lda_globx = m
+  lda_locx  = size(locx, 1)
+  lock      = size(locx,2)
+  maxk      = lock
+  k         = maxk
 
   call psb_bcast(ictxt,k,root=iiroot)
 
@@ -130,13 +127,20 @@ subroutine  psb_cgatherm(globx, locx, desc_a, info, iroot)
     call psb_errpush(info,name)
     goto 9999
   end if
+  
+  call psb_realloc(m,k,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
 
   globx(:,:)=czero
 
   do j=1,k
     do i=1,desc_a%get_local_rows()
       call psb_loc_to_glob(i,idx,desc_a,info)
-      globx(idx,jglobx+j-1) = locx(i,jlx+j-1)
+      globx(idx,j) = locx(i,jlx+j-1)
     end do
   end do
 
@@ -146,12 +150,12 @@ subroutine  psb_cgatherm(globx, locx, desc_a, info, iroot)
       if (me /= desc_a%ovrlap_elem(i,3)) then 
         idx = desc_a%ovrlap_elem(i,1)
         call psb_loc_to_glob(idx,desc_a,info)
-        globx(idx,jglobx+j-1) = czero
+        globx(idx,j) = czero
       end if
     end do
   end do
 
-  call psb_sum(ictxt,globx(1:m,jglobx:jglobx+k-1),root=root)
+  call psb_sum(ictxt,globx(1:m,1:k),root=root)
 
   call psb_erractionrestore(err_act)
   return  
@@ -223,7 +227,7 @@ subroutine  psb_cgatherv(globx, locx, desc_a, info, iroot)
   implicit none
 
   complex(psb_spk_), intent(in)    :: locx(:)
-  complex(psb_spk_), intent(out)   :: globx(:)
+  complex(psb_spk_), intent(out), allocatable :: globx(:)
   type(psb_desc_type), intent(in) :: desc_a
   integer(psb_ipk_), intent(out)            :: info
   integer(psb_ipk_), intent(in), optional   :: iroot
@@ -268,12 +272,12 @@ subroutine  psb_cgatherv(globx, locx, desc_a, info, iroot)
   jlocx=1
   ilocx = 1
 
-  lda_globx = size(globx)
-  lda_locx  = size(locx)
-
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
 
+  lda_globx = m
+  lda_locx  = size(locx)
+
   k = 1
 
 
@@ -295,6 +299,13 @@ subroutine  psb_cgatherv(globx, locx, desc_a, info, iroot)
     goto 9999
   end if
   
+  call psb_realloc(m,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+  
   globx(:)=czero
 
   do i=1,desc_a%get_local_rows()
@@ -334,7 +345,7 @@ subroutine  psb_cgather_vect(globx, locx, desc_a, info, iroot)
   implicit none
 
   type(psb_c_vect_type), intent(inout) :: locx
-  complex(psb_spk_), intent(out)    :: globx(:)
+  complex(psb_spk_), intent(out), allocatable   :: globx(:)
   type(psb_desc_type), intent(in)   :: desc_a
   integer(psb_ipk_), intent(out)              :: info
   integer(psb_ipk_), intent(in), optional     :: iroot
@@ -379,11 +390,11 @@ subroutine  psb_cgather_vect(globx, locx, desc_a, info, iroot)
   jlocx=1
   ilocx = 1
 
-  lda_globx = size(globx)
-  lda_locx  = locx%get_nrows()
 
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
+  lda_globx = m
+  lda_locx  = locx%get_nrows()
   
   k = 1
 
@@ -406,6 +417,13 @@ subroutine  psb_cgather_vect(globx, locx, desc_a, info, iroot)
     goto 9999
   end if
   
+  call psb_realloc(m,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+  
   globx(:) = czero
   llocx    = locx%get_vect()
 

base/comm/psb_dgather.f90

@@ -50,7 +50,7 @@ subroutine  psb_dgatherm(globx, locx, desc_a, info, iroot)
   implicit none
 
   real(psb_dpk_), intent(in)    :: locx(:,:)
-  real(psb_dpk_), intent(out)   :: globx(:,:)
+  real(psb_dpk_), intent(out), allocatable :: globx(:,:)
   type(psb_desc_type), intent(in) :: desc_a
   integer(psb_ipk_), intent(out)            :: info
   integer(psb_ipk_), intent(in), optional   :: iroot
@@ -100,16 +100,13 @@ subroutine  psb_dgatherm(globx, locx, desc_a, info, iroot)
   ilocx = 1
   jlocx = 1
 
-  lda_globx = size(globx,1)
-  lda_locx  = size(locx, 1)
-
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
-
-  lock=size(locx,2)-jlocx+1
-  globk=size(globx,2)-jglobx+1
-  maxk=min(lock,globk)
-  k = maxk
+  lda_globx = m
+  lda_locx  = size(locx, 1)
+  lock      = size(locx,2)
+  maxk      = lock
+  k         = maxk
 
   call psb_bcast(ictxt,k,root=iiroot)
 
@@ -130,13 +127,20 @@ subroutine  psb_dgatherm(globx, locx, desc_a, info, iroot)
     call psb_errpush(info,name)
     goto 9999
   end if
+  
+  call psb_realloc(m,k,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
 
   globx(:,:)=dzero
 
   do j=1,k
     do i=1,desc_a%get_local_rows()
       call psb_loc_to_glob(i,idx,desc_a,info)
-      globx(idx,jglobx+j-1) = locx(i,jlx+j-1)
+      globx(idx,j) = locx(i,jlx+j-1)
     end do
   end do
 
@@ -146,12 +150,12 @@ subroutine  psb_dgatherm(globx, locx, desc_a, info, iroot)
       if (me /= desc_a%ovrlap_elem(i,3)) then 
         idx = desc_a%ovrlap_elem(i,1)
         call psb_loc_to_glob(idx,desc_a,info)
-        globx(idx,jglobx+j-1) = dzero
+        globx(idx,j) = dzero
       end if
     end do
   end do
 
-  call psb_sum(ictxt,globx(1:m,jglobx:jglobx+k-1),root=root)
+  call psb_sum(ictxt,globx(1:m,1:k),root=root)
 
   call psb_erractionrestore(err_act)
   return  
@@ -223,7 +227,7 @@ subroutine  psb_dgatherv(globx, locx, desc_a, info, iroot)
   implicit none
 
   real(psb_dpk_), intent(in)    :: locx(:)
-  real(psb_dpk_), intent(out)   :: globx(:)
+  real(psb_dpk_), intent(out), allocatable :: globx(:)
   type(psb_desc_type), intent(in) :: desc_a
   integer(psb_ipk_), intent(out)            :: info
   integer(psb_ipk_), intent(in), optional   :: iroot
@@ -268,12 +272,12 @@ subroutine  psb_dgatherv(globx, locx, desc_a, info, iroot)
   jlocx=1
   ilocx = 1
 
-  lda_globx = size(globx)
-  lda_locx  = size(locx)
-
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
 
+  lda_globx = m
+  lda_locx  = size(locx)
+
   k = 1
 
 
@@ -295,6 +299,13 @@ subroutine  psb_dgatherv(globx, locx, desc_a, info, iroot)
     goto 9999
   end if
   
+  call psb_realloc(m,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+  
   globx(:)=dzero
 
   do i=1,desc_a%get_local_rows()
@@ -334,7 +345,7 @@ subroutine  psb_dgather_vect(globx, locx, desc_a, info, iroot)
   implicit none
 
   type(psb_d_vect_type), intent(inout) :: locx
-  real(psb_dpk_), intent(out)    :: globx(:)
+  real(psb_dpk_), intent(out), allocatable   :: globx(:)
   type(psb_desc_type), intent(in)   :: desc_a
   integer(psb_ipk_), intent(out)              :: info
   integer(psb_ipk_), intent(in), optional     :: iroot
@@ -379,11 +390,11 @@ subroutine  psb_dgather_vect(globx, locx, desc_a, info, iroot)
   jlocx=1
   ilocx = 1
 
-  lda_globx = size(globx)
-  lda_locx  = locx%get_nrows()
 
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
+  lda_globx = m
+  lda_locx  = locx%get_nrows()
   
   k = 1
 
@@ -406,6 +417,13 @@ subroutine  psb_dgather_vect(globx, locx, desc_a, info, iroot)
     goto 9999
   end if
   
+  call psb_realloc(m,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+  
   globx(:) = dzero
   llocx    = locx%get_vect()
 

base/comm/psb_sgather.f90

@@ -50,7 +50,7 @@ subroutine  psb_sgatherm(globx, locx, desc_a, info, iroot)
   implicit none
 
   real(psb_spk_), intent(in)    :: locx(:,:)
-  real(psb_spk_), intent(out)   :: globx(:,:)
+  real(psb_spk_), intent(out), allocatable :: globx(:,:)
   type(psb_desc_type), intent(in) :: desc_a
   integer(psb_ipk_), intent(out)            :: info
   integer(psb_ipk_), intent(in), optional   :: iroot
@@ -100,16 +100,13 @@ subroutine  psb_sgatherm(globx, locx, desc_a, info, iroot)
   ilocx = 1
   jlocx = 1
 
-  lda_globx = size(globx,1)
-  lda_locx  = size(locx, 1)
-
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
-
-  lock=size(locx,2)-jlocx+1
-  globk=size(globx,2)-jglobx+1
-  maxk=min(lock,globk)
-  k = maxk
+  lda_globx = m
+  lda_locx  = size(locx, 1)
+  lock      = size(locx,2)
+  maxk      = lock
+  k         = maxk
 
   call psb_bcast(ictxt,k,root=iiroot)
 
@@ -130,13 +127,20 @@ subroutine  psb_sgatherm(globx, locx, desc_a, info, iroot)
     call psb_errpush(info,name)
     goto 9999
   end if
+  
+  call psb_realloc(m,k,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
 
   globx(:,:)=szero
 
   do j=1,k
     do i=1,desc_a%get_local_rows()
       call psb_loc_to_glob(i,idx,desc_a,info)
-      globx(idx,jglobx+j-1) = locx(i,jlx+j-1)
+      globx(idx,j) = locx(i,jlx+j-1)
     end do
   end do
 
@@ -146,12 +150,12 @@ subroutine  psb_sgatherm(globx, locx, desc_a, info, iroot)
       if (me /= desc_a%ovrlap_elem(i,3)) then 
         idx = desc_a%ovrlap_elem(i,1)
         call psb_loc_to_glob(idx,desc_a,info)
-        globx(idx,jglobx+j-1) = szero
+        globx(idx,j) = szero
       end if
     end do
   end do
 
-  call psb_sum(ictxt,globx(1:m,jglobx:jglobx+k-1),root=root)
+  call psb_sum(ictxt,globx(1:m,1:k),root=root)
 
   call psb_erractionrestore(err_act)
   return  
@@ -223,7 +227,7 @@ subroutine  psb_sgatherv(globx, locx, desc_a, info, iroot)
   implicit none
 
   real(psb_spk_), intent(in)    :: locx(:)
-  real(psb_spk_), intent(out)   :: globx(:)
+  real(psb_spk_), intent(out), allocatable :: globx(:)
   type(psb_desc_type), intent(in) :: desc_a
   integer(psb_ipk_), intent(out)            :: info
   integer(psb_ipk_), intent(in), optional   :: iroot
@@ -268,12 +272,12 @@ subroutine  psb_sgatherv(globx, locx, desc_a, info, iroot)
   jlocx=1
   ilocx = 1
 
-  lda_globx = size(globx)
-  lda_locx  = size(locx)
-
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
 
+  lda_globx = m
+  lda_locx  = size(locx)
+
   k = 1
 
 
@@ -295,6 +299,13 @@ subroutine  psb_sgatherv(globx, locx, desc_a, info, iroot)
     goto 9999
   end if
   
+  call psb_realloc(m,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+  
   globx(:)=szero
 
   do i=1,desc_a%get_local_rows()
@@ -334,7 +345,7 @@ subroutine  psb_sgather_vect(globx, locx, desc_a, info, iroot)
   implicit none
 
   type(psb_s_vect_type), intent(inout) :: locx
-  real(psb_spk_), intent(out)    :: globx(:)
+  real(psb_spk_), intent(out), allocatable   :: globx(:)
   type(psb_desc_type), intent(in)   :: desc_a
   integer(psb_ipk_), intent(out)              :: info
   integer(psb_ipk_), intent(in), optional     :: iroot
@@ -379,11 +390,11 @@ subroutine  psb_sgather_vect(globx, locx, desc_a, info, iroot)
   jlocx=1
   ilocx = 1
 
-  lda_globx = size(globx)
-  lda_locx  = locx%get_nrows()
 
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
+  lda_globx = m
+  lda_locx  = locx%get_nrows()
   
   k = 1
 
@@ -406,6 +417,13 @@ subroutine  psb_sgather_vect(globx, locx, desc_a, info, iroot)
     goto 9999
   end if
   
+  call psb_realloc(m,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+  
   globx(:) = szero
   llocx    = locx%get_vect()
 

base/comm/psb_zgather.f90

@@ -50,7 +50,7 @@ subroutine  psb_zgatherm(globx, locx, desc_a, info, iroot)
   implicit none
 
   complex(psb_dpk_), intent(in)    :: locx(:,:)
-  complex(psb_dpk_), intent(out)   :: globx(:,:)
+  complex(psb_dpk_), intent(out), allocatable :: globx(:,:)
   type(psb_desc_type), intent(in) :: desc_a
   integer(psb_ipk_), intent(out)            :: info
   integer(psb_ipk_), intent(in), optional   :: iroot
@@ -100,16 +100,13 @@ subroutine  psb_zgatherm(globx, locx, desc_a, info, iroot)
   ilocx = 1
   jlocx = 1
 
-  lda_globx = size(globx,1)
-  lda_locx  = size(locx, 1)
-
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
-
-  lock=size(locx,2)-jlocx+1
-  globk=size(globx,2)-jglobx+1
-  maxk=min(lock,globk)
-  k = maxk
+  lda_globx = m
+  lda_locx  = size(locx, 1)
+  lock      = size(locx,2)
+  maxk      = lock
+  k         = maxk
 
   call psb_bcast(ictxt,k,root=iiroot)
 
@@ -130,13 +127,20 @@ subroutine  psb_zgatherm(globx, locx, desc_a, info, iroot)
     call psb_errpush(info,name)
     goto 9999
   end if
+  
+  call psb_realloc(m,k,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
 
   globx(:,:)=zzero
 
   do j=1,k
     do i=1,desc_a%get_local_rows()
       call psb_loc_to_glob(i,idx,desc_a,info)
-      globx(idx,jglobx+j-1) = locx(i,jlx+j-1)
+      globx(idx,j) = locx(i,jlx+j-1)
     end do
   end do
 
@@ -146,12 +150,12 @@ subroutine  psb_zgatherm(globx, locx, desc_a, info, iroot)
       if (me /= desc_a%ovrlap_elem(i,3)) then 
         idx = desc_a%ovrlap_elem(i,1)
         call psb_loc_to_glob(idx,desc_a,info)
-        globx(idx,jglobx+j-1) = zzero
+        globx(idx,j) = zzero
       end if
     end do
   end do
 
-  call psb_sum(ictxt,globx(1:m,jglobx:jglobx+k-1),root=root)
+  call psb_sum(ictxt,globx(1:m,1:k),root=root)
 
   call psb_erractionrestore(err_act)
   return  
@@ -223,7 +227,7 @@ subroutine  psb_zgatherv(globx, locx, desc_a, info, iroot)
   implicit none
 
   complex(psb_dpk_), intent(in)    :: locx(:)
-  complex(psb_dpk_), intent(out)   :: globx(:)
+  complex(psb_dpk_), intent(out), allocatable :: globx(:)
   type(psb_desc_type), intent(in) :: desc_a
   integer(psb_ipk_), intent(out)            :: info
   integer(psb_ipk_), intent(in), optional   :: iroot
@@ -268,12 +272,12 @@ subroutine  psb_zgatherv(globx, locx, desc_a, info, iroot)
   jlocx=1
   ilocx = 1
 
-  lda_globx = size(globx)
-  lda_locx  = size(locx)
-
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
 
+  lda_globx = m
+  lda_locx  = size(locx)
+
   k = 1
 
 
@@ -295,6 +299,13 @@ subroutine  psb_zgatherv(globx, locx, desc_a, info, iroot)
     goto 9999
   end if
   
+  call psb_realloc(m,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+  
   globx(:)=zzero
 
   do i=1,desc_a%get_local_rows()
@@ -334,7 +345,7 @@ subroutine  psb_zgather_vect(globx, locx, desc_a, info, iroot)
   implicit none
 
   type(psb_z_vect_type), intent(inout) :: locx
-  complex(psb_dpk_), intent(out)    :: globx(:)
+  complex(psb_dpk_), intent(out), allocatable   :: globx(:)
   type(psb_desc_type), intent(in)   :: desc_a
   integer(psb_ipk_), intent(out)              :: info
   integer(psb_ipk_), intent(in), optional     :: iroot
@@ -379,11 +390,11 @@ subroutine  psb_zgather_vect(globx, locx, desc_a, info, iroot)
   jlocx=1
   ilocx = 1
 
-  lda_globx = size(globx)
-  lda_locx  = locx%get_nrows()
 
   m = desc_a%get_global_rows()
   n = desc_a%get_global_cols()
+  lda_globx = m
+  lda_locx  = locx%get_nrows()
   
   k = 1
 
@@ -406,6 +417,13 @@ subroutine  psb_zgather_vect(globx, locx, desc_a, info, iroot)
     goto 9999
   end if
   
+  call psb_realloc(m,globx,info)
+  if (info /= psb_success_) then 
+    info=psb_err_alloc_dealloc_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+  
   globx(:) = zzero
   llocx    = locx%get_vect()
 

base/modules/psb_c_comm_mod.f90

@@ -128,7 +128,7 @@ module psb_c_comm_mod
     subroutine psb_cgatherm(globx, locx, desc_a, info, root)
       use psb_descriptor_type
       complex(psb_spk_), intent(in)  :: locx(:,:)
-      complex(psb_spk_), intent(out) :: globx(:,:)
+      complex(psb_spk_), intent(out), allocatable  :: globx(:,:)
       type(psb_desc_type), intent(in)  :: desc_a
       integer(psb_ipk_), intent(out)             :: info
       integer(psb_ipk_), intent(in), optional    :: root
@@ -136,7 +136,7 @@ module psb_c_comm_mod
     subroutine  psb_cgatherv(globx, locx, desc_a, info, root)
       use psb_descriptor_type
       complex(psb_spk_), intent(in)  :: locx(:)
-      complex(psb_spk_), intent(out) :: globx(:)
+      complex(psb_spk_), intent(out), allocatable  :: globx(:)
       type(psb_desc_type), intent(in)  :: desc_a
       integer(psb_ipk_), intent(out)             :: info
       integer(psb_ipk_), intent(in), optional    :: root
@@ -145,7 +145,7 @@ module psb_c_comm_mod
       use psb_descriptor_type
       use psb_c_vect_mod
       type(psb_c_vect_type), intent(inout) :: locx
-      complex(psb_spk_), intent(out)     :: globx(:)
+      complex(psb_spk_), intent(out), allocatable :: globx(:)
       type(psb_desc_type), intent(in) :: desc_a
       integer(psb_ipk_), intent(out)            :: info
       integer(psb_ipk_), intent(in), optional   :: root

base/modules/psb_d_comm_mod.f90

@@ -35,47 +35,47 @@ module psb_d_comm_mod
     subroutine  psb_dovrlm(x,desc_a,info,jx,ik,work,update,mode)
       use psb_descriptor_type
       real(psb_dpk_), intent(inout), target   :: x(:,:)
-      type(psb_desc_type), intent(in)         :: desc_a
-      integer(psb_ipk_), intent(out)                    :: info
+      type(psb_desc_type), intent(in)            :: desc_a
+      integer(psb_ipk_), intent(out)                       :: info
       real(psb_dpk_), intent(inout), optional, target :: work(:)
-      integer(psb_ipk_), intent(in), optional           :: update,jx,ik,mode
+      integer(psb_ipk_), intent(in), optional              :: update,jx,ik,mode
     end subroutine psb_dovrlm
     subroutine  psb_dovrlv(x,desc_a,info,work,update,mode)
       use psb_descriptor_type
       real(psb_dpk_), intent(inout), target   :: x(:)
-      type(psb_desc_type), intent(in)         :: desc_a
-      integer(psb_ipk_), intent(out)                    :: info
+      type(psb_desc_type), intent(in)            :: desc_a
+      integer(psb_ipk_), intent(out)                       :: info
       real(psb_dpk_), intent(inout), optional, target :: work(:)
-      integer(psb_ipk_), intent(in), optional           :: update,mode
+      integer(psb_ipk_), intent(in), optional              :: update,mode
     end subroutine psb_dovrlv
     subroutine  psb_dovrl_vect(x,desc_a,info,work,update,mode)
       use psb_descriptor_type
       use psb_d_vect_mod
-      type(psb_d_vect_type), intent(inout)   :: x
+      type(psb_d_vect_type), intent(inout)    :: x
       type(psb_desc_type), intent(in)         :: desc_a
       integer(psb_ipk_), intent(out)                    :: info
       real(psb_dpk_), intent(inout), optional, target :: work(:)
       integer(psb_ipk_), intent(in), optional           :: update,mode
     end subroutine psb_dovrl_vect
-  end interface
+  end interface psb_ovrl
 
   interface psb_halo
     subroutine  psb_dhalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
       use psb_descriptor_type
-      real(psb_dpk_), intent(inout),target    :: x(:,:)
-      type(psb_desc_type), intent(in)         :: desc_a
-      integer(psb_ipk_), intent(out)                    :: info
-      real(psb_dpk_), intent(in), optional    :: alpha
+      real(psb_dpk_), intent(inout), target :: x(:,:)
+      type(psb_desc_type), intent(in)          :: desc_a
+      integer(psb_ipk_), intent(out)                     :: info
+      real(psb_dpk_), intent(in), optional  :: alpha
       real(psb_dpk_), target, optional, intent(inout) :: work(:)
       integer(psb_ipk_), intent(in), optional           :: mode,jx,ik,data
       character, intent(in), optional         :: tran
     end subroutine psb_dhalom
     subroutine  psb_dhalov(x,desc_a,info,alpha,work,tran,mode,data)
       use psb_descriptor_type
-      real(psb_dpk_), intent(inout)           :: x(:)
+      real(psb_dpk_), intent(inout)        :: x(:)
       type(psb_desc_type), intent(in)         :: desc_a
       integer(psb_ipk_), intent(out)                    :: info
-      real(psb_dpk_), intent(in), optional    :: alpha
+      real(psb_dpk_), intent(in), optional :: alpha
       real(psb_dpk_), target, optional, intent(inout) :: work(:)
       integer(psb_ipk_), intent(in), optional           :: mode,data
       character, intent(in), optional         :: tran
@@ -91,27 +91,27 @@ module psb_d_comm_mod
       integer(psb_ipk_), intent(in), optional           :: mode,data
       character, intent(in), optional         :: tran
     end subroutine psb_dhalo_vect
-  end interface
+  end interface psb_halo
 
 
   interface psb_scatter
     subroutine  psb_dscatterm(globx, locx, desc_a, info, root)
       use psb_descriptor_type
-      real(psb_dpk_), intent(out)    :: locx(:,:)
-      real(psb_dpk_), intent(in)     :: globx(:,:)
+      real(psb_dpk_), intent(out) :: locx(:,:)
+      real(psb_dpk_), intent(in)  :: globx(:,:)
       type(psb_desc_type), intent(in)  :: desc_a
       integer(psb_ipk_), intent(out)             :: info
       integer(psb_ipk_), intent(in), optional    :: root
     end subroutine psb_dscatterm
     subroutine  psb_dscatterv(globx, locx, desc_a, info, root)
       use psb_descriptor_type
-      real(psb_dpk_), intent(out)    :: locx(:)
-      real(psb_dpk_), intent(in)     :: globx(:)
+      real(psb_dpk_), intent(out) :: locx(:)
+      real(psb_dpk_), intent(in)  :: globx(:)
       type(psb_desc_type), intent(in)  :: desc_a
       integer(psb_ipk_), intent(out)             :: info
       integer(psb_ipk_), intent(in), optional    :: root
     end subroutine psb_dscatterv
-  end interface
+  end interface psb_scatter
 
   interface psb_gather
     subroutine  psb_dsp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keeploc)
@@ -125,31 +125,31 @@ module psb_d_comm_mod
       integer(psb_ipk_), intent(in), optional   :: root,dupl
       logical, intent(in), optional   :: keepnum,keeploc
     end subroutine psb_dsp_allgather
-    subroutine  psb_dgatherm(globx, locx, desc_a, info, root)
+    subroutine psb_dgatherm(globx, locx, desc_a, info, root)
       use psb_descriptor_type
-      real(psb_dpk_), intent(in)    :: locx(:,:)
-      real(psb_dpk_), intent(out)   :: globx(:,:)
-      type(psb_desc_type), intent(in) :: desc_a
-      integer(psb_ipk_), intent(out)            :: info
-      integer(psb_ipk_), intent(in), optional   :: root
+      real(psb_dpk_), intent(in)  :: locx(:,:)
+      real(psb_dpk_), intent(out), allocatable  :: globx(:,:)
+      type(psb_desc_type), intent(in)  :: desc_a
+      integer(psb_ipk_), intent(out)             :: info
+      integer(psb_ipk_), intent(in), optional    :: root
     end subroutine psb_dgatherm
     subroutine  psb_dgatherv(globx, locx, desc_a, info, root)
       use psb_descriptor_type
-      real(psb_dpk_), intent(in)    :: locx(:)
-      real(psb_dpk_), intent(out)   :: globx(:)
-      type(psb_desc_type), intent(in) :: desc_a
-      integer(psb_ipk_), intent(out)            :: info
-      integer(psb_ipk_), intent(in), optional   :: root
+      real(psb_dpk_), intent(in)  :: locx(:)
+      real(psb_dpk_), intent(out), allocatable  :: globx(:)
+      type(psb_desc_type), intent(in)  :: desc_a
+      integer(psb_ipk_), intent(out)             :: info
+      integer(psb_ipk_), intent(in), optional    :: root
     end subroutine psb_dgatherv
     subroutine  psb_dgather_vect(globx, locx, desc_a, info, root)
       use psb_descriptor_type
       use psb_d_vect_mod
       type(psb_d_vect_type), intent(inout) :: locx
-      real(psb_dpk_), intent(out)     :: globx(:)
+      real(psb_dpk_), intent(out), allocatable :: globx(:)
       type(psb_desc_type), intent(in) :: desc_a
       integer(psb_ipk_), intent(out)            :: info
       integer(psb_ipk_), intent(in), optional   :: root
     end subroutine psb_dgather_vect
-  end interface
-  
+  end interface psb_gather
+
 end module psb_d_comm_mod

base/modules/psb_s_comm_mod.f90

@@ -35,47 +35,47 @@ module psb_s_comm_mod
     subroutine  psb_sovrlm(x,desc_a,info,jx,ik,work,update,mode)
       use psb_descriptor_type
       real(psb_spk_), intent(inout), target   :: x(:,:)
-      type(psb_desc_type), intent(in)         :: desc_a
-      integer(psb_ipk_), intent(out)                    :: info
+      type(psb_desc_type), intent(in)            :: desc_a
+      integer(psb_ipk_), intent(out)                       :: info
       real(psb_spk_), intent(inout), optional, target :: work(:)
-      integer(psb_ipk_), intent(in), optional           :: update,jx,ik,mode
+      integer(psb_ipk_), intent(in), optional              :: update,jx,ik,mode
     end subroutine psb_sovrlm
     subroutine  psb_sovrlv(x,desc_a,info,work,update,mode)
       use psb_descriptor_type
       real(psb_spk_), intent(inout), target   :: x(:)
-      type(psb_desc_type), intent(in)         :: desc_a
-      integer(psb_ipk_), intent(out)                    :: info
+      type(psb_desc_type), intent(in)            :: desc_a
+      integer(psb_ipk_), intent(out)                       :: info
       real(psb_spk_), intent(inout), optional, target :: work(:)
-      integer(psb_ipk_), intent(in), optional           :: update,mode
+      integer(psb_ipk_), intent(in), optional              :: update,mode
     end subroutine psb_sovrlv
     subroutine  psb_sovrl_vect(x,desc_a,info,work,update,mode)
       use psb_descriptor_type
       use psb_s_vect_mod
-      type(psb_s_vect_type), intent(inout)   :: x
+      type(psb_s_vect_type), intent(inout)    :: x
       type(psb_desc_type), intent(in)         :: desc_a
       integer(psb_ipk_), intent(out)                    :: info
       real(psb_spk_), intent(inout), optional, target :: work(:)
       integer(psb_ipk_), intent(in), optional           :: update,mode
     end subroutine psb_sovrl_vect
-  end interface
+  end interface psb_ovrl
 
   interface psb_halo
     subroutine  psb_shalom(x,desc_a,info,alpha,jx,ik,work,tran,mode,data)
       use psb_descriptor_type
-      real(psb_spk_), intent(inout),target    :: x(:,:)
-      type(psb_desc_type), intent(in)         :: desc_a
-      integer(psb_ipk_), intent(out)                    :: info
-      real(psb_spk_), intent(in), optional    :: alpha
+      real(psb_spk_), intent(inout), target :: x(:,:)
+      type(psb_desc_type), intent(in)          :: desc_a
+      integer(psb_ipk_), intent(out)                     :: info
+      real(psb_spk_), intent(in), optional  :: alpha
       real(psb_spk_), target, optional, intent(inout) :: work(:)
       integer(psb_ipk_), intent(in), optional           :: mode,jx,ik,data
       character, intent(in), optional         :: tran
     end subroutine psb_shalom
     subroutine  psb_shalov(x,desc_a,info,alpha,work,tran,mode,data)
       use psb_descriptor_type
-      real(psb_spk_), intent(inout)           :: x(:)
+      real(psb_spk_), intent(inout)        :: x(:)
       type(psb_desc_type), intent(in)         :: desc_a
       integer(psb_ipk_), intent(out)                    :: info
-      real(psb_spk_), intent(in), optional    :: alpha
+      real(psb_spk_), intent(in), optional :: alpha
       real(psb_spk_), target, optional, intent(inout) :: work(:)
       integer(psb_ipk_), intent(in), optional           :: mode,data
       character, intent(in), optional         :: tran
@@ -91,27 +91,27 @@ module psb_s_comm_mod
       integer(psb_ipk_), intent(in), optional           :: mode,data
       character, intent(in), optional         :: tran
     end subroutine psb_shalo_vect
-  end interface
+  end interface psb_halo
 
 
   interface psb_scatter
     subroutine  psb_sscatterm(globx, locx, desc_a, info, root)
       use psb_descriptor_type
-      real(psb_spk_), intent(out)    :: locx(:,:)
-      real(psb_spk_), intent(in)     :: globx(:,:)
+      real(psb_spk_), intent(out) :: locx(:,:)
+      real(psb_spk_), intent(in)  :: globx(:,:)
       type(psb_desc_type), intent(in)  :: desc_a
       integer(psb_ipk_), intent(out)             :: info
       integer(psb_ipk_), intent(in), optional    :: root
     end subroutine psb_sscatterm
     subroutine  psb_sscatterv(globx, locx, desc_a, info, root)
       use psb_descriptor_type
-      real(psb_spk_), intent(out)    :: locx(:)
-      real(psb_spk_), intent(in)     :: globx(:)
+      real(psb_spk_), intent(out) :: locx(:)
+      real(psb_spk_), intent(in)  :: globx(:)
       type(psb_desc_type), intent(in)  :: desc_a
       integer(psb_ipk_), intent(out)             :: info
       integer(psb_ipk_), intent(in), optional    :: root
     end subroutine psb_sscatterv
-  end interface
+  end interface psb_scatter
 
   interface psb_gather
     subroutine  psb_ssp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keeploc)
@@ -125,31 +125,31 @@ module psb_s_comm_mod
       integer(psb_ipk_), intent(in), optional   :: root,dupl
       logical, intent(in), optional   :: keepnum,keeploc
     end subroutine psb_ssp_allgather
-    subroutine  psb_sgatherm(globx, locx, desc_a, info, root)
+    subroutine psb_sgatherm(globx, locx, desc_a, info, root)
       use psb_descriptor_type
-      real(psb_spk_), intent(in)    :: locx(:,:)
-      real(psb_spk_), intent(out)   :: globx(:,:)
-      type(psb_desc_type), intent(in) :: desc_a
-      integer(psb_ipk_), intent(out)            :: info
-      integer(psb_ipk_), intent(in), optional   :: root
+      real(psb_spk_), intent(in)  :: locx(:,:)
+      real(psb_spk_), intent(out), allocatable  :: globx(:,:)
+      type(psb_desc_type), intent(in)  :: desc_a
+      integer(psb_ipk_), intent(out)             :: info
+      integer(psb_ipk_), intent(in), optional    :: root
     end subroutine psb_sgatherm
     subroutine  psb_sgatherv(globx, locx, desc_a, info, root)
       use psb_descriptor_type
-      real(psb_spk_), intent(in)    :: locx(:)
-      real(psb_spk_), intent(out)   :: globx(:)
-      type(psb_desc_type), intent(in) :: desc_a
-      integer(psb_ipk_), intent(out)            :: info
-      integer(psb_ipk_), intent(in), optional   :: root
+      real(psb_spk_), intent(in)  :: locx(:)
+      real(psb_spk_), intent(out), allocatable  :: globx(:)
+      type(psb_desc_type), intent(in)  :: desc_a
+      integer(psb_ipk_), intent(out)             :: info
+      integer(psb_ipk_), intent(in), optional    :: root
     end subroutine psb_sgatherv
     subroutine  psb_sgather_vect(globx, locx, desc_a, info, root)
       use psb_descriptor_type
       use psb_s_vect_mod
       type(psb_s_vect_type), intent(inout) :: locx
-      real(psb_spk_), intent(out)     :: globx(:)
+      real(psb_spk_), intent(out), allocatable :: globx(:)
       type(psb_desc_type), intent(in) :: desc_a
       integer(psb_ipk_), intent(out)            :: info
       integer(psb_ipk_), intent(in), optional   :: root
     end subroutine psb_sgather_vect
-  end interface
+  end interface psb_gather
 
 end module psb_s_comm_mod

base/modules/psb_z_comm_mod.f90

@@ -111,7 +111,7 @@ module psb_z_comm_mod
       integer(psb_ipk_), intent(out)             :: info
       integer(psb_ipk_), intent(in), optional    :: root
     end subroutine psb_zscatterv
-  end interface
+  end interface psb_scatter
 
   interface psb_gather
     subroutine  psb_zsp_allgather(globa, loca, desc_a, info, root, dupl,keepnum,keeploc)
@@ -125,10 +125,10 @@ module psb_z_comm_mod
       integer(psb_ipk_), intent(in), optional   :: root,dupl
       logical, intent(in), optional   :: keepnum,keeploc
     end subroutine psb_zsp_allgather
-    subroutine  psb_zgatherm(globx, locx, desc_a, info, root)
+    subroutine psb_zgatherm(globx, locx, desc_a, info, root)
       use psb_descriptor_type
       complex(psb_dpk_), intent(in)  :: locx(:,:)
-      complex(psb_dpk_), intent(out) :: globx(:,:)
+      complex(psb_dpk_), intent(out), allocatable  :: globx(:,:)
       type(psb_desc_type), intent(in)  :: desc_a
       integer(psb_ipk_), intent(out)             :: info
       integer(psb_ipk_), intent(in), optional    :: root
@@ -136,7 +136,7 @@ module psb_z_comm_mod
     subroutine  psb_zgatherv(globx, locx, desc_a, info, root)
       use psb_descriptor_type
       complex(psb_dpk_), intent(in)  :: locx(:)
-      complex(psb_dpk_), intent(out) :: globx(:)
+      complex(psb_dpk_), intent(out), allocatable  :: globx(:)
       type(psb_desc_type), intent(in)  :: desc_a
       integer(psb_ipk_), intent(out)             :: info
       integer(psb_ipk_), intent(in), optional    :: root
@@ -145,7 +145,7 @@ module psb_z_comm_mod
       use psb_descriptor_type
       use psb_z_vect_mod
       type(psb_z_vect_type), intent(inout) :: locx
-      complex(psb_dpk_), intent(out)     :: globx(:)
+      complex(psb_dpk_), intent(out), allocatable :: globx(:)
       type(psb_desc_type), intent(in) :: desc_a
       integer(psb_ipk_), intent(out)            :: info
       integer(psb_ipk_), intent(in), optional   :: root