NAG Library Routine Document

F16GHF (BLAS_ZWAXPBY)

Basic Linear Algebra Subprograms Technical (BLAST) Forum (2001) Basic Linear Algebra Subprograms Technical (BLAST) Forum Standard University of Tennessee, Knoxville, Tennessee http://www.netlib.org/blas/blast-forum/blas-report.pdf

5 Parameters

1: N – INTEGERInput: On entry: $n$ , the number of elements in $x$ , $y$ and $w$ .
2: ALPHA – COMPLEX (KIND=nag_wp)Input: On entry: the scalar $α$ .
3: X( $1 + (N - 1) \times |INCX|$ ) – COMPLEX (KIND=nag_wp) arrayInput: On entry: the $n$ -element vector $x$ .
If $INCX > 0$ , $x_{i}$ must be stored in $X (1 + (i - 1) \times INCX)$ , for $i = 1, 2, \dots, N$ .

If $INCX < 0$ , $x_{i}$ must be stored in $X (1 - (N - i) \times INCX)$ , for $i = 1, 2, \dots, N$ .

Intermediate elements of X are not referenced.
4: INCX – INTEGERInput: On entry: the increment in the subscripts of X between successive elements of $x$ .
Constraint: $INCX \neq 0$ .
5: BETA – COMPLEX (KIND=nag_wp)Input: On entry: the scalar $β$ .
6: Y( $1 + (N - 1) \times |INCY|$ ) – COMPLEX (KIND=nag_wp) arrayInput: On entry: the $n$ -element vector $y$ .
If $INCY > 0$ , $y_{i}$ must be stored in $Y (1 + (i - 1) \times INCY)$ , for $i = 1, 2, \dots, N$ .

If $INCY < 0$ , $y_{i}$ must be stored in $Y (1 - (N - i) \times INCY)$ , for $i = 1, 2, \dots, N$ .

Intermediate elements of Y are not referenced.
7: INCY – INTEGERInput: On entry: the increment in the subscripts of Y between successive elements of $y$ .
Constraint: $INCY \neq 0$ .
8: W( $1 + (N - 1) \times |INCW|$ ) – COMPLEX (KIND=nag_wp) arrayOutput: On exit: the $n$ -element vector $w$ .
If $INCW > 0$ , $w_{i}$ is in $W (1 + (i - 1) \times INCW)$ , for $i = 1, 2, \dots, N$ .

If $INCW < 0$ , $w_{i}$ is in $W (1 + (N - i) \times INCW)$ , for $i = 1, 2, \dots, N$ .

Intermediate elements of W are not referenced.
9: INCW – INTEGERInput: On entry: the increment in the subscripts of W between successive elements of $w$ .
Constraint: $INCW \neq 0$ .

6 Error Indicators and Warnings

INCX = 0

INCY = 0

INCW = 0

, an error message is printed and program execution is terminated.

7 Accuracy

The BLAS standard requires accurate implementations which avoid unnecessary over/underflow (see Section 2.7 of Basic Linear Algebra Subprograms Technical (BLAST) Forum (2001)).

8 Further Comments

None.

9 Example

This example computes the result of a scaled vector accumulation for

\begin{matrix} α = 3 + 2 i, & x = {(- 4 + 2.1 i, 3.7 + 4.5 i, - 6 + 1.2 i)}^{T}, \\ β = - i, & y = {(- 3 - 2.4 i, 6.4 - 5 i, - 5.1)}^{T} . \end{matrix}

NAG Library Routine DocumentF16GHF (BLAS_ZWAXPBY)

+− Contents

1 Purpose

2 Specification

3 Description

4 References

5 Parameters

6 Error Indicators and Warnings

7 Accuracy

8 Further Comments

9 Example

9.1 Program Text

9.2 Program Data

9.3 Program Results

NAG Library Routine Document

F16GHF (BLAS_ZWAXPBY)