LAPACK  3.10.0
LAPACK: Linear Algebra PACKage

◆ sscal()

subroutine sscal ( integer  N,
real  SA,
real, dimension(*)  SX,
integer  INCX 
)

SSCAL

Purpose:
    SSCAL scales a vector by a constant.
    uses unrolled loops for increment equal to 1.
Parameters
[in]N
          N is INTEGER
         number of elements in input vector(s)
[in]SA
          SA is REAL
           On entry, SA specifies the scalar alpha.
[in,out]SX
          SX is REAL array, dimension ( 1 + ( N - 1 )*abs( INCX ) )
[in]INCX
          INCX is INTEGER
         storage spacing between elements of SX
Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
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 78 of file sscal.f.

79 *
80 * -- Reference BLAS level1 routine --
81 * -- Reference BLAS is a software package provided by Univ. of Tennessee, --
82 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
83 *
84 * .. Scalar Arguments ..
85  REAL SA
86  INTEGER INCX,N
87 * ..
88 * .. Array Arguments ..
89  REAL SX(*)
90 * ..
91 *
92 * =====================================================================
93 *
94 * .. Local Scalars ..
95  INTEGER I,M,MP1,NINCX
96 * ..
97 * .. Intrinsic Functions ..
98  INTRINSIC mod
99 * ..
100  IF (n.LE.0 .OR. incx.LE.0) RETURN
101  IF (incx.EQ.1) THEN
102 *
103 * code for increment equal to 1
104 *
105 *
106 * clean-up loop
107 *
108  m = mod(n,5)
109  IF (m.NE.0) THEN
110  DO i = 1,m
111  sx(i) = sa*sx(i)
112  END DO
113  IF (n.LT.5) RETURN
114  END IF
115  mp1 = m + 1
116  DO i = mp1,n,5
117  sx(i) = sa*sx(i)
118  sx(i+1) = sa*sx(i+1)
119  sx(i+2) = sa*sx(i+2)
120  sx(i+3) = sa*sx(i+3)
121  sx(i+4) = sa*sx(i+4)
122  END DO
123  ELSE
124 *
125 * code for increment not equal to 1
126 *
127  nincx = n*incx
128  DO i = 1,nincx,incx
129  sx(i) = sa*sx(i)
130  END DO
131  END IF
132  RETURN
133 *
134 * End of SSCAL
135 *