Fixed problem that was causing bodies to overlap

MintonGroup · Feb 20, 2023 · 13f6a8a · 13f6a8a
1 parent 8df230d
commit 13f6a8a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/fraggle/fraggle_generate.f90 b/src/fraggle/fraggle_generate.f90
@@ -420,7 +420,7 @@ module subroutine fraggle_generate_pos_vec(collider, nbody_system, param, lfailu
                if (.not.loverlap(j)) cycle
                loverlap(j) = .false.
                ! Check for overlaps between fragments
-               do i = j + 1, nfrag
+               do concurrent(i = 1:nfrag, i/=j)
                   dis = .mag.(fragments%rc(:,j) - fragments%rc(:,i))
                   loverlap(j) = loverlap(j) .or. (dis <= rbuffer * (fragments%radius(i) + fragments%radius(j))) 
                end do