ismin.f90 Source File

!
! $Header$
!
FUNCTION ismin(n,sx,incx)
  !
  IMPLICIT NONE
  !
  integer :: n,i,incx,ismin,ix
  real :: sx((n-1)*incx+1),sxmin
  !
  ix=1
  ismin=1
  sxmin=sx(1)
  DO i=1,n-1
     ix=ix+incx
     if(sx(ix).lt.sxmin) then
         sxmin=sx(ix)
         ismin=i+1
     endif
  ENDDO
  !
  return
end function ismin
!