zla_wwaddw()

subroutine zla_wwaddw	(	integer	n,
		complex*16, dimension( * )	x,
		complex*16, dimension( * )	y,
		complex*16, dimension( * )	w
	)

ZLA_WWADDW adds a vector into a doubled-single vector.

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

Purpose:

    ZLA_WWADDW adds a vector W into a doubled-single vector (X, Y).

    This works for all extant IBM's hex and binary floating point
    arithmetic, but not for decimal.

Parameters

[in]	N	N is INTEGER The length of vectors X, Y, and W.
[in,out]	X	X is COMPLEX*16 array, dimension (N) The first part of the doubled-single accumulation vector.
[in,out]	Y	Y is COMPLEX*16 array, dimension (N) The second part of the doubled-single accumulation vector.
[in]	W	W is COMPLEX*16 array, dimension (N) The vector to be added.

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Definition at line 80 of file zla_wwaddw.f.

*
*  -- LAPACK computational routine --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*
*     .. Scalar Arguments ..
      INTEGER            N
*     ..
*     .. Array Arguments ..
      COMPLEX*16         X( * ), Y( * ), W( * )
*     ..
*
*  =====================================================================
*
*     .. Local Scalars ..
      COMPLEX*16         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
*
*     End of ZLA_WWADDW
*

Here is the caller graph for this function: