Merge pull request #1541 from ExtremeFLOW/martin/fix/aabb_sort

Add option to create aabb tree directly from aabb and fix sorting bug
ExtremeFLOW · Nov 11, 2024 · bbf0a12 · bbf0a12
2 parents e33f49e + 22ef147
commit bbf0a12
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 2 deletions.
diff --git a/src/mesh/aabb_tree.f90 b/src/mesh/aabb_tree.f90
@@ -439,8 +439,9 @@ function sort(array) result(indices)
        minidx = -1
        do imin = 1, size(array)
           if (.not. visited(imin) .and. minidx .eq. -1) minidx = imin
-
-          if (visited(imin) .and. array(imin) .lt. array(minidx)) minidx = imin
+          if (minidx .gt. -1) then
+             if (visited(imin) .and. array(imin) .lt. array(minidx)) minidx = imin
+          end if
        end do
 
        indices(i) = minidx
@@ -641,6 +642,7 @@ subroutine aabb_tree_query_overlaps(this, object, object_index, overlaps)
           end if
        end if
     end do
+    call simple_stack%free()
   end subroutine aabb_tree_query_overlaps
 
   ! -------------------------------------------------------------------------- !

diff --git a/src/mesh/search_tree/aabb.f90 b/src/mesh/search_tree/aabb.f90
@@ -187,6 +187,8 @@ function get_aabb(object, padding) result(box)
 
     select type(object)
 
+      type is (aabb_t)
+       call box%init(object%get_min(), object%get_max())
       type is (tri_t)
        box = get_aabb_element(object)
       type is (hex_t)