2
/* -- translated by f2c (version 19940927).
3
You must link the resulting object file with the libraries:
4
-lf2c -lm (in that order)
9
integer idamax_(integer *n, doublereal *dx, integer *incx)
13
/* System generated locals */
14
integer ret_val, i__1;
18
static doublereal dmax__;
22
/* finds the index of element having max. absolute value.
23
jack dongarra, linpack, 3/11/78.
24
modified 3/93 to return if incx .le. 0.
25
modified 12/3/93, array(1) declarations changed to array(*)
31
#define DX(I) dx[(I)-1]
35
if (*n < 1 || *incx <= 0) {
46
/* code for increment not equal to 1 */
52
for (i = 2; i <= *n; ++i) {
53
if ((d__1 = DX(ix), abs(d__1)) <= dmax__) {
57
dmax__ = (d__1 = DX(ix), abs(d__1));
64
/* code for increment equal to 1 */
69
for (i = 2; i <= *n; ++i) {
70
if ((d__1 = DX(i), abs(d__1)) <= dmax__) {
74
dmax__ = (d__1 = DX(i), abs(d__1));