1
integer function isamax(n,sx,incx)
3
* $Id: isamax.f 19695 2010-10-29 16:51:02Z d3y133 $
6
c finds the index of element having max. absolute value.
7
c jack dongarra, linpack, 3/11/78.
8
c modified 3/93 to return if incx .le. 0.
9
c modified 12/3/93, array(1) declarations changed to array(*)
15
if( n.lt.1 .or. incx.le.0 ) return
20
c code for increment not equal to 1
26
if(abs(sx(ix)).le.smax) go to 5
33
c code for increment equal to 1
37
if(abs(sx(i)).le.smax) go to 30
3
* =========== DOCUMENTATION ===========
5
* Online html documentation available at
6
* http://www.netlib.org/lapack/explore-html/
11
* INTEGER FUNCTION ISAMAX(N,SX,INCX)
13
* .. Scalar Arguments ..
16
* .. Array Arguments ..
26
*> ISAMAX finds the index of element having max. absolute value.
32
*> \author Univ. of Tennessee
33
*> \author Univ. of California Berkeley
34
*> \author Univ. of Colorado Denver
37
*> \date November 2011
41
*> \par Further Details:
42
* =====================
46
*> jack dongarra, linpack, 3/11/78.
47
*> modified 3/93 to return if incx .le. 0.
48
*> modified 12/3/93, array(1) declarations changed to array(*)
51
* =====================================================================
52
INTEGER FUNCTION ISAMAX(N,SX,INCX)
54
* -- Reference BLAS level1 routine (version 3.4.0) --
55
* -- Reference BLAS is a software package provided by Univ. of Tennessee, --
56
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
59
* .. Scalar Arguments ..
62
* .. Array Arguments ..
66
* =====================================================================
72
* .. Intrinsic Functions ..
76
IF (N.LT.1 .OR. INCX.LE.0) RETURN
81
* code for increment equal to 1
85
IF (ABS(SX(I)).GT.SMAX) THEN
92
* code for increment not equal to 1
98
IF (ABS(SX(IX)).GT.SMAX) THEN