LAPACK 3.12.0
LAPACK: Linear Algebra PACKage
Loading...
Searching...
No Matches
zla_wwaddw.f
Go to the documentation of this file.
1*> \brief \b ZLA_WWADDW adds a vector into a doubled-single vector.
2*
3* =========== DOCUMENTATION ===========
4*
5* Online html documentation available at
6* http://www.netlib.org/lapack/explore-html/
7*
8*> \htmlonly
9*> Download ZLA_WWADDW + dependencies
10*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zla_wwaddw.f">
11*> [TGZ]</a>
12*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zla_wwaddw.f">
13*> [ZIP]</a>
14*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zla_wwaddw.f">
15*> [TXT]</a>
16*> \endhtmlonly
17*
18* Definition:
19* ===========
20*
21* SUBROUTINE ZLA_WWADDW( N, X, Y, W )
22*
23* .. Scalar Arguments ..
24* INTEGER N
25* ..
26* .. Array Arguments ..
27* COMPLEX*16 X( * ), Y( * ), W( * )
28* ..
29*
30*
31*> \par Purpose:
32* =============
33*>
34*> \verbatim
35*>
36*> ZLA_WWADDW adds a vector W into a doubled-single vector (X, Y).
37*>
38*> This works for all extant IBM's hex and binary floating point
39*> arithmetic, but not for decimal.
40*> \endverbatim
41*
42* Arguments:
43* ==========
44*
45*> \param[in] N
46*> \verbatim
47*> N is INTEGER
48*> The length of vectors X, Y, and W.
49*> \endverbatim
50*>
51*> \param[in,out] X
52*> \verbatim
53*> X is COMPLEX*16 array, dimension (N)
54*> The first part of the doubled-single accumulation vector.
55*> \endverbatim
56*>
57*> \param[in,out] Y
58*> \verbatim
59*> Y is COMPLEX*16 array, dimension (N)
60*> The second part of the doubled-single accumulation vector.
61*> \endverbatim
62*>
63*> \param[in] W
64*> \verbatim
65*> W is COMPLEX*16 array, dimension (N)
66*> The vector to be added.
67*> \endverbatim
68*
69* Authors:
70* ========
71*
72*> \author Univ. of Tennessee
73*> \author Univ. of California Berkeley
74*> \author Univ. of Colorado Denver
75*> \author NAG Ltd.
76*
77*> \ingroup la_wwaddw
78*
79* =====================================================================
80 SUBROUTINE zla_wwaddw( N, X, Y, W )
81*
82* -- LAPACK computational routine --
83* -- LAPACK is a software package provided by Univ. of Tennessee, --
84* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
85*
86* .. Scalar Arguments ..
87 INTEGER N
88* ..
89* .. Array Arguments ..
90 COMPLEX*16 X( * ), Y( * ), W( * )
91* ..
92*
93* =====================================================================
94*
95* .. Local Scalars ..
96 COMPLEX*16 S
97 INTEGER I
98* ..
99* .. Executable Statements ..
100 DO 10 i = 1, n
101 s = x(i) + w(i)
102 s = (s + s) - s
103 y(i) = ((x(i) - s) + w(i)) + y(i)
104 x(i) = s
105 10 CONTINUE
106 RETURN
107*
108* End of ZLA_WWADDW
109*
110 END
subroutine zla_wwaddw(n, x, y, w)
ZLA_WWADDW adds a vector into a doubled-single vector.
Definition zla_wwaddw.f:81