real function sasum	(	integer	N,
		real, dimension(*)	SX,
		integer	INCX
	)

SASUM

Purpose:

    SASUM takes the sum of the absolute values.
    uses unrolled loops for increment equal to one.

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Date

December 2016

Further Details:

     jack dongarra, linpack, 3/11/78.
     modified 3/93 to return if incx .le. 0.
     modified 12/3/93, array(1) declarations changed to array(*)

Definition at line 54 of file sasum.f.

*
*  -- Reference BLAS level1 routine (version 3.7.0) --
*  -- Reference BLAS is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*     December 2016
*
*     .. Scalar Arguments ..
      INTEGER incx,n
*     ..
*     .. Array Arguments ..
      REAL sx(*)
*     ..
*
*  =====================================================================
*
*     .. Local Scalars ..
      REAL stemp
      INTEGER i,m,mp1,nincx
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC abs,mod
*     ..
      sasum = 0.0e0
      stemp = 0.0e0
      IF (n.LE.0 .OR. incx.LE.0) RETURN
      IF (incx.EQ.1) THEN
*        code for increment equal to 1
*
*
*        clean-up loop
*
         m = mod(n,6)
         IF (m.NE.0) THEN
            DO i = 1,m
               stemp = stemp + abs(sx(i))
            END DO
            IF (n.LT.6) THEN
               sasum = stemp
               RETURN
            END IF
         END IF
         mp1 = m + 1
         DO i = mp1,n,6
            stemp = stemp + abs(sx(i)) + abs(sx(i+1)) +
     $              abs(sx(i+2)) + abs(sx(i+3)) +
     $              abs(sx(i+4)) + abs(sx(i+5))
         END DO
      ELSE
*
*        code for increment not equal to 1
*
         nincx = n*incx
         DO i = 1,nincx,incx
            stemp = stemp + abs(sx(i))
         END DO
      END IF
      sasum = stemp
      RETURN

Here is the caller graph for this function: