Do not do the swap if the index is the same -- fixing ESCOMP#397

ekluzek · Jun 4, 2023 · 4addd32 · 4addd32
1 parent bdf4d46
commit 4addd32
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/route/build/src/nr_utils.f90 b/route/build/src/nr_utils.f90
@@ -108,7 +108,7 @@ SUBROUTINE indexx(arr,index)
  jstack=jstack-2
  else
  k=(l+r)/2
- call swap(index(k),index(l+1))
+ if ( k /= l+1 ) call swap(index(k),index(l+1))
  call icomp_xchg(index(l),index(r))
  call icomp_xchg(index(l+1),index(r))
  call icomp_xchg(index(l),index(l+1))
@@ -126,7 +126,7 @@ SUBROUTINE indexx(arr,index)
  if (arr(index(j)) <= a) exit
  end do
  if (j < i) exit
- call swap(index(i),index(j))
+ if ( i /= j ) call swap(index(i),index(j))
  end do
  index(l+1)=index(j)
  index(j)=indext