From c29d3aae1d59c67fb66c490672d10e32cd4a1ac7 Mon Sep 17 00:00:00 2001
From: Salvatore Filippone <salvatore.filippone@uniroma2.it>
Date: Fri, 4 Dec 2020 14:45:44 +0100
Subject: [PATCH] Additional info on data exchanges on output cdprt

---
 base/tools/psb_cdprt.f90 | 38 ++++++++++++++++++++++++++++++++++++--
 1 file changed, 36 insertions(+), 2 deletions(-)

diff --git a/base/tools/psb_cdprt.f90 b/base/tools/psb_cdprt.f90
index 0ee64f2a..1e14ed64 100644
--- a/base/tools/psb_cdprt.f90
+++ b/base/tools/psb_cdprt.f90
@@ -56,7 +56,8 @@ subroutine psb_cdprt(iout,desc_p,glob,short, verbosity)
   integer(psb_ipk_)   :: me, np
   integer(psb_ipk_)   :: total_snd, total_rcv, total_xhcg, global_halo, global_points
   integer(psb_ipk_)   :: local_snd, local_rcv, local_xhcg, local_halo, local_points
-  real(psb_dpk_)      :: av2s, v2s
+  integer(psb_ipk_)   :: max_xchg
+  real(psb_dpk_)      :: av2s, v2s, anxchg
   
   if (present(glob)) then 
     glob_ = glob
@@ -94,15 +95,24 @@ subroutine psb_cdprt(iout,desc_p,glob,short, verbosity)
   av2s          = v2s
   call psb_sum(ctxt, global_halo)
   call psb_sum(ctxt, av2s)
-  av2s = av2s / np 
+  av2s = av2s / np
+  max_xchg = get_nxchg(desc_p)
+  anxchg   = max_xchg
+  call psb_max(ctxt, max_xchg)
+  call psb_sum(ctxt, anxchg)
+  anxchg = anxchg/np
+  
   if (me == psb_root_) then
     write(iout,*) ' Communication descriptor details '
     write(iout,*) '               Descriptor format:       ',desc_p%get_fmt()    
     write(iout,*) ' Global descriptor data: points:',global_points,' halo:',global_halo
+    write(iout,*) ' Maximum number of exchanges     :',max_xchg
+    write(iout,*) ' Average number of exchanges     :',anxchg
     write(iout,*)
     write(iout,*) ' Average volume to surface ratio :',av2s
     write(iout,*)
   end if
+  flush(iout)
   call psb_barrier(ctxt)
   
   if (verb_ <= 1) return
@@ -123,12 +133,16 @@ subroutine psb_cdprt(iout,desc_p,glob,short, verbosity)
     call psb_barrier(ctxt)
   end do
 
+  if (verb_ <= 2) return
+
   if (me==psb_root_) write(iout,*) 'Communication data for : comm_halo'
   do i=0, np-1
     if (me == i)  &
        &  call  print_my_xchg(iout,desc_p,verbosity=verb_,data=psb_comm_halo_,glob=glob_)
     call psb_barrier(ctxt)
   end do
+
+  if (verb_ <= 3) return
   
   if (me==psb_root_) write(iout,*) 'Communication data for : comm_ext'
   do i=0, np-1
@@ -140,6 +154,26 @@ subroutine psb_cdprt(iout,desc_p,glob,short, verbosity)
   return
 
 contains
+
+  function get_nxchg(desc,data) result(res)
+    implicit none
+    type(psb_desc_type), intent(in), target :: desc
+    integer(psb_ipk_), intent(in), optional :: data
+    integer(psb_ipk_) :: res
+    
+    class(psb_i_base_vect_type), pointer :: vpnt
+    integer(psb_ipk_) :: ip, nerv, nesd, data_, totxch, idxr, idxs, info
+    
+    if (present(data)) then
+      data_ = data
+    else
+      data_ = psb_comm_halo_
+    end if
+    call psb_cd_v_get_list(data_,desc_p,vpnt,totxch,idxr,idxs,info)
+    res = totxch
+  end function get_nxchg
+  
+    
   subroutine print_my_xchg(iout,desc_p,data,glob,short, verbosity)
     implicit none 
     type(psb_desc_type), intent(in), target   :: desc_p