LAPACK  3.8.0
LAPACK: Linear Algebra PACKage

◆ dla_wwaddw()

subroutine dla_wwaddw ( integer  N,
double precision, dimension( * )  X,
double precision, dimension( * )  Y,
double precision, dimension( * )  W 
)

DLA_WWADDW adds a vector into a doubled-single vector.

Download DLA_WWADDW + dependencies [TGZ] [ZIP] [TXT]

Purpose:
    DLA_WWADDW adds a vector W into a doubled-single vector (X, Y).

    This works for all extant IBM's hex and binary floating point
    arithmetics, but not for decimal.
Parameters
[in]N
          N is INTEGER
            The length of vectors X, Y, and W.
[in,out]X
          X is DOUBLE PRECISION array, dimension (N)
            The first part of the doubled-single accumulation vector.
[in,out]Y
          Y is DOUBLE PRECISION array, dimension (N)
            The second part of the doubled-single accumulation vector.
[in]W
          W is DOUBLE PRECISION array, dimension (N)
            The vector to be added.
Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
Date
December 2016

Definition at line 83 of file dla_wwaddw.f.

83 *
84 * -- LAPACK computational routine (version 3.7.0) --
85 * -- LAPACK is a software package provided by Univ. of Tennessee, --
86 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
87 * December 2016
88 *
89 * .. Scalar Arguments ..
90  INTEGER n
91 * ..
92 * .. Array Arguments ..
93  DOUBLE PRECISION x( * ), y( * ), w( * )
94 * ..
95 *
96 * =====================================================================
97 *
98 * .. Local Scalars ..
99  DOUBLE PRECISION s
100  INTEGER i
101 * ..
102 * .. Executable Statements ..
103 *
104  DO 10 i = 1, n
105  s = x(i) + w(i)
106  s = (s + s) - s
107  y(i) = ((x(i) - s) + w(i)) + y(i)
108  x(i) = s
109  10 CONTINUE
110  RETURN
Here is the caller graph for this function: