miller's simplex subrtn

From: Yiorgos Lykidis (abacus_at_freemail.gr)
Date: 05/17/04 

Date: Mon, 17 May 2004 20:58:55 +0300

We are having problems in using your simplex fortran subroutine smplx.f90
written by A. Miller. A DO loop results in an "array bounds exceeded
message", since the A vector has a dimension equal to ka (stored in variable
m of the previous sbrtn smplx1). Do you have any idea of what could be going
wrong?

SUBROUTINE smplx1
bi is defined as bi(m,m)
....
CALL crout1 (bi(1:,1), m, m, iend, indx, y, ierr)
....

SUBROUTINE crout1(a, ka, n, iend, indx, temp, ierr)
....
maxdim = ka*n
....
  DO kj = kk, maxdim, ka
    c = a(kj)
! when kj > ka then a(kj) is not defined!!!
    a(kj) = a(lj)
    a(lj) = c
    lj = lj + ka
  END DO

We are using CVF66C.

Yiorgos