Fix bug in heapsort on complex data

1 year ago · 329a295315
parent 3c53066a04
commit 329a295315
2 changed files with 4 additions and 0 deletions
--- a/base/serial/sort/psb_c_hsort_impl.f90
+++ b/base/serial/sort/psb_c_hsort_impl.f90
@ -899,7 +899,9 @@ subroutine psi_c_idx_heap_get_first(key,index,last,heap,idxs,dir,info)
  endif

  key     = heap(1)
+  index   = idxs(1)
  heap(1) = heap(last)
+  idxs(1) = idxs(last)
  last    = last - 1

  select case(dir)
--- a/base/serial/sort/psb_z_hsort_impl.f90
+++ b/base/serial/sort/psb_z_hsort_impl.f90
@ -899,7 +899,9 @@ subroutine psi_z_idx_heap_get_first(key,index,last,heap,idxs,dir,info)
  endif

  key     = heap(1)
+  index   = idxs(1)
  heap(1) = heap(last)
+  idxs(1) = idxs(last)
  last    = last - 1

  select case(dir)