Avoid divide by zero for desc printout.

7 years ago · ee6992e938
parent 64fdf563fc
commit ee6992e938
1 changed files with 5 additions and 1 deletions
--- a/base/tools/psb_cdprt.f90
+++ b/base/tools/psb_cdprt.f90
@ -93,7 +93,11 @@ subroutine psb_cdprt(iout,desc_p,glob,short, verbosity)
    if (me == i)  then
      write(iout,*) me,': Local descriptor data: points:',local_points,&
           & ' halo:',local_halo
-      write(iout,*) me,': Volume to surface ratio:',real(local_points,psb_dpk_)/real(local_halo,psb_dpk_)
+      if (local_halo>0) then 
+        write(iout,*) me,': Volume to surface ratio:',real(local_points,psb_dpk_)/real(local_halo,psb_dpk_)
+      else
+        write(iout,*) me,': Volume to surface ratio:',0.0_psb_dpk_
+      end if
    end if
    call psb_barrier(ictxt)
  end do