Fixed old usage of igebs2d.

base/comm/psb_dscatter.f90

@@ -121,11 +121,7 @@ subroutine  psb_dscatterm(globx, locx, desc_a, info, iroot)
   m = psb_cd_get_global_rows(desc_a)
   n = psb_cd_get_global_cols(desc_a)
   
-  if (me == iiroot) then
-     call igebs2d(ictxt, 'all', ' ', 1, 1, k, 1)
-  else
-     call igebr2d(ictxt, 'all', ' ', 1, 1, k, 1, iiroot, 0)
-  end if
+  call psb_bcast(ictxt,k,root=iiroot)
 
   !  there should be a global check on k here!!!
 
@@ -331,13 +327,6 @@ subroutine  psb_dscatterv(globx, locx, desc_a, info, iroot)
   n = psb_cd_get_global_cols(desc_a)
   
   k = 1
-
-  if (me == iiroot) then
-     call igebs2d(ictxt, 'all', ' ', 1, 1, k, 1)
-  else
-     call igebr2d(ictxt, 'all', ' ', 1, 1, k, 1, iiroot, 0)
-  end if
-
   !  there should be a global check on k here!!!
 
   call psb_chkglobvect(m,n,size(globx),iglobx,jglobx,desc_a,info)

base/comm/psb_iscatter.f90

@@ -120,11 +120,7 @@ subroutine  psb_iscatterm(globx, locx, desc_a, info, iroot)
   m = psb_cd_get_global_rows(desc_a)
   n = psb_cd_get_global_cols(desc_a)
   
-  if (me == iiroot) then
-     call igebs2d(ictxt, 'all', ' ', 1, 1, k, 1)
-  else
-     call igebr2d(ictxt, 'all', ' ', 1, 1, k, 1, iiroot, 0)
-  end if
+  call psb_bcast(ictxt,k,root=iiroot)
 
   !  there should be a global check on k here!!!
 
@@ -331,12 +327,6 @@ subroutine  psb_iscatterv(globx, locx, desc_a, info, iroot)
   
   k = 1
 
-  if (me == iiroot) then
-     call igebs2d(ictxt, 'all', ' ', 1, 1, k, 1)
-  else
-     call igebr2d(ictxt, 'all', ' ', 1, 1, k, 1, iiroot, 0)
-  end if
-
   !  there should be a global check on k here!!!
 
   call psb_chkglobvect(m,n,size(globx),iglobx,jglobx,desc_a,info)

base/comm/psb_zscatter.f90

@@ -122,11 +122,7 @@ subroutine  psb_zscatterm(globx, locx, desc_a, info, iroot)
   m = psb_cd_get_global_rows(desc_a)
   n = psb_cd_get_global_cols(desc_a)
 
-  if (me == iiroot) then
-    call igebs2d(ictxt, 'all', ' ', 1, 1, k, 1)
-  else
-    call igebr2d(ictxt, 'all', ' ', 1, 1, k, 1, iiroot, 0)
-  end if
+  call psb_bcast(ictxt,k,root=iiroot)
 
   !  there should be a global check on k here!!!
 
@@ -334,12 +330,6 @@ subroutine  psb_zscatterv(globx, locx, desc_a, info, iroot)
 
   k = 1
 
-  if (me == iiroot) then
-    call igebs2d(ictxt, 'all', ' ', 1, 1, k, 1)
-  else
-    call igebr2d(ictxt, 'all', ' ', 1, 1, k, 1, iiroot, 0)
-  end if
-
   !  there should be a global check on k here!!!
 
   call psb_chkglobvect(m,n,size(globx),iglobx,jglobx,desc_a,info)