| REAL*4 function samaxf(n,dx,incx) |
| c |
| c finds the index of element having max. absolute value. |
| c jack dongarra, linpack, 3/11/78. |
| c modified 3/93 to return if incx .le. 0. |
| c modified 12/3/93, array(1) declarations changed to array(*) |
| c |
| REAL*4 dx(*) |
| integer i,incx,ix,n |
| c |
| samaxf = 0. |
| if( n.lt.1 .or. incx.le.0 ) return |
| |
| if(incx.eq.1)go to 20 |
| c |
| c code for increment not equal to 1 |
| c |
| ix = 1 |
| samaxf = abs(dx(1)) |
| ix = ix + incx |
| do 10 i = 2,n |
| if(abs(dx(ix)).le.samaxf) go to 5 |
| samaxf = abs(dx(ix)) |
| 5 ix = ix + incx |
| 10 continue |
| return |
| c |
| c code for increment equal to 1 |
| c |
| 20 samaxf = abs(dx(1)) |
| do 30 i = 2,n |
| if(abs(dx(i)).le.samaxf) go to 30 |
| samaxf = abs(dx(i)) |
| 30 continue |
| return |
| end |
