/*
 * Author:  G. Jungman
 * RCS:     $Id: source_iamax_r.h,v 1.3 1999/10/18 00:45:12 jungman Exp $
 */

  BASE_TYPE max = 0.0;
  size_t i = 0;
  CBLAS_INDEX n;
  CBLAS_INDEX result;
  for(n=0; n<N; n++) {
    if(fabs(X[i]) > max) {
      max = fabs(X[i]);
      result = n;
    }
    i += incX;
  }
  return result;