F06FUF (PDF version)
F06 Chapter Contents
F06 Chapter Introduction
NAG Library Manual

NAG Library Routine Document

F06FUF

Note:  before using this routine, please read the Users' Note for your implementation to check the interpretation of bold italicised terms and other implementation-dependent details.

+ Contents

    1  Purpose
    7  Accuracy
    9  Example

1  Purpose

F06FUF applies a LINPACK (as opposed to NAG) style real elementary reflection to a real vector.

2  Specification

SUBROUTINE F06FUF ( N, Z, INCZ, Z1, ALPHA, X, INCX)
INTEGER  N, INCZ, INCX
REAL (KIND=nag_wp)  Z(*), Z1, ALPHA, X(*)

3  Description

F06FUF applies a real elementary reflection (Householder matrix) P, as generated by F06FSF, to a given real vector:
α x P α x
where x is an n-element real vector and α a real scalar.

4  References

None.

5  Parameters

1:     N – INTEGERInput
On entry: n, the number of elements in x and z.
2:     Z(*) – REAL (KIND=nag_wp) arrayInput
Note: the dimension of the array Z must be at least max1, 1+N-1 ×INCZ .
On entry: the vector z, as returned by F06FSF.
If INCZ>0, zi must be stored in Z1+i-1×INCZ, for i=1,2,,N.
If INCZ<0, zi must be stored in Z1-N-i×INCZ, for i=1,2,,N.
3:     INCZ – INTEGERInput
On entry: the increment in the subscripts of Z between successive elements of z.
4:     Z1 – REAL (KIND=nag_wp)Input
On entry: the scalar ζ, as returned by F06FSF.
If ζ=0, P is assumed to be the unit matrix and the transformation is skipped.
5:     ALPHA – REAL (KIND=nag_wp)Input/Output
On entry: the original scalar α.
On exit: the transformed scalar α.
6:     X(*) – REAL (KIND=nag_wp) arrayInput/Output
Note: the dimension of the array X must be at least max1, 1+N-1 ×INCX .
On entry: the original vector x.
If INCX>0, xi must be stored in X1+i-1×INCX, for i=1,2,,N.
If INCX<0, xi must be stored in X1-N-i×INCX, for i=1,2,,N.
On exit: the transformed vector x stored in the same array elements used to supply the original vector x.
7:     INCX – INTEGERInput
On entry: the increment in the subscripts of X between successive elements of x.

6  Error Indicators and Warnings

None.

7  Accuracy

Not applicable.

8  Further Comments

None.

9  Example

None.

F06FUF (PDF version)
F06 Chapter Contents
F06 Chapter Introduction
NAG Library Manual

© The Numerical Algorithms Group Ltd, Oxford, UK. 2012