Fixed array shape problem for the double list version of the sweep

MintonGroup · Oct 25, 2021 · abd655c · abd655c
1 parent a3b8b0d
commit abd655c
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/src/encounter/encounter_check.f90 b/src/encounter/encounter_check.f90
@@ -938,9 +938,11 @@ module subroutine encounter_check_sweep_aabb_double_list(self, n1, n2, nenc, ind
       call util_index_array(ind_arr, ntot)
       allocate(ibeg(SWEEPDIM * ntot))
       allocate(iend(SWEEPDIM * ntot))
-      do dim = 1, SWEEPDIM
-         ibeg((dim - 1) * ntot + 1:dim * ntot) = self%aabb(dim)%ibeg(:)
-         iend((dim - 1) * ntot + 1:dim * ntot) = self%aabb(dim)%iend(:)
+      do i = 1, ntot
+         do dim = 1, SWEEPDIM
+            ibeg((i - 1) * SWEEPDIM + dim) = self%aabb(dim)%ibeg(i)
+            iend((i - 1) * SWEEPDIM + dim) = self%aabb(dim)%iend(i)
+         end do
       end do
       ! Sweep the intervals for each of the massive bodies along one dimension
       ! This will build a ragged pair of index lists inside of the lenc data structure