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

September 2012

Definition at line 83 of file dla_wwaddw.f.

*
*  -- LAPACK computational routine (version 3.4.2) --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*     September 2012
*
*     .. Scalar Arguments ..
      INTEGER            n
*     ..
*     .. Array Arguments ..
      DOUBLE PRECISION   x( * ), y( * ), w( * )
*     ..
*
*  =====================================================================
*
*     .. Local Scalars ..
      DOUBLE PRECISION   s
      INTEGER            i
*     ..
*     .. Executable Statements ..
*
      DO 10 i = 1, n
        s = x(i) + w(i)
        s = (s + s) - s
        y(i) = ((x(i) - s) + w(i)) + y(i)
        x(i) = s
 10   CONTINUE
      RETURN

Here is the caller graph for this function: