Added fix for condition when the narrow phase returns no true encounters

MintonGroup · Oct 25, 2021 · 7f95561 · 7f95561
1 parent 862357f
commit 7f95561
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/encounter/encounter_check.f90 b/src/encounter/encounter_check.f90
@@ -289,6 +289,12 @@ subroutine encounter_check_reduce_broadphase(n, nenc, index1, index2, lencounter
       logical, dimension(:), allocatable :: ltmp
 
       nenc = count(lencounter(:)) ! Count the true number of encounters
+      if (nenc == 0) then
+         if (allocated(index1)) deallocate(index1)
+         if (allocated(index2)) deallocate(index2)
+         if (allocated(lvdotr)) deallocate(lvdotr)
+         return
+      end if
 
       allocate(itmp(nenc))
       itmp(:) = pack(index1(:), lencounter(:))