home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libblas / src_original / sasum.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.5 KB  |  63 lines
  1.       REAL             FUNCTION SASUM( N, SX, INCX )
  2. *
  3. *     takes the sum of the absolute values.
  4. *     uses unrolled loops for increment equal to one.
  5. *     jack dongarra, linpack, 3/11/78.
  6. *
  7. *     .. Scalar Arguments ..
  8.       INTEGER                          INCX, N
  9. *     ..
  10. *     .. Array Arguments ..
  11.       REAL                             SX( 1 )
  12. *     ..
  13. *     .. Local Scalars ..
  14.       INTEGER                          I, IX, M, MP1, NINCX
  15.       REAL                             STEMP
  16. *     ..
  17. *     .. Intrinsic Functions ..
  18.       INTRINSIC                        ABS, MOD
  19. *     ..
  20. *     .. Executable Statements ..
  21. *
  22.       SASUM = 0.0E0
  23.       STEMP = 0.0E0
  24.       IF( N.LE.0 )
  25.      $   RETURN
  26.       IF( INCX.EQ.1 )
  27.      $   GO TO 20
  28. *
  29. *        code for increment not equal to 1
  30. *
  31.       IX = 1
  32.       IF( INCX.LT.0 )
  33.      $   IX = 1 - ( N-1 )*INCX
  34.       NINCX = IX + ( N-1 )*INCX
  35.       DO 10 I = IX, NINCX, INCX
  36.          STEMP = STEMP + ABS( SX( I ) )
  37.    10 CONTINUE
  38.       SASUM = STEMP
  39.       RETURN
  40. *
  41. *        code for increment equal to 1
  42. *
  43. *
  44. *        clean-up loop
  45. *
  46.    20 M = MOD( N, 6 )
  47.       IF( M.EQ.0 )
  48.      $   GO TO 40
  49.       DO 30 I = 1, M
  50.          STEMP = STEMP + ABS( SX( I ) )
  51.    30 CONTINUE
  52.       IF( N.LT.6 )
  53.      $   GO TO 60
  54.    40 MP1 = M + 1
  55.       DO 50 I = MP1, N, 6
  56.          STEMP = STEMP + ABS( SX( I ) ) + ABS( SX( I+1 ) ) +
  57.      $           ABS( SX( I+2 ) ) + ABS( SX( I+3 ) ) +
  58.      $           ABS( SX( I+4 ) ) + ABS( SX( I+5 ) )
  59.    50 CONTINUE
  60.    60 SASUM = STEMP
  61.       RETURN
  62.       END
  63.