F06WNF (ZLANHF) (PDF version)
F06 Chapter Contents
F06 Chapter Introduction
NAG Library Manual

NAG Library Routine Document

F06WNF (ZLANHF)

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

1  Purpose

F06WNF (ZLANHF) returns the value of the 1-norm, the -norm, the Frobenius norm, or the maximum absolute value of the elements of a complex Hermitian matrix A stored in Rectangular Full Packed (RFP) format.

2  Specification

FUNCTION F06WNF ( NORM, TRANSR, UPLO, N, AR, WORK)
REAL (KIND=nag_wp) F06WNF
INTEGER  N
REAL (KIND=nag_wp)  WORK(*)
COMPLEX (KIND=nag_wp)  AR(N*(N+1)/2)
CHARACTER(1)  NORM, TRANSR, UPLO
The routine may be called by its LAPACK name zlanhf.

3  Description

Given a complex n by n symmetric matrix, A, F06WNF (ZLANHF) calculates one of the values given by
A1=maxji=1naij (the 1-norm of A),
A=maxij= 1naij (the -norm of A),
AF=i=1nj=1naij21/2 (the Frobenius norm of A),   or
maxi,jaij (the maximum absolute element value of A).
A is stored in compact form using the RFP format. The RFP storage format is described in Section 3.3.3 in the F07 Chapter Introduction.

4  References

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
Gustavson F G, Waśniewski J, Dongarra J J and Langou J (2010) Rectangular full packed format for Cholesky's algorithm: factorization, solution, and inversion ACM Trans. Math. Software 37, 2

5  Parameters

1:     NORM – CHARACTER(1)Input
On entry: specifies the value to be returned.
NORM='1' or 'O'
The 1-norm.
NORM='I'
The -norm.
NORM='F' or 'E'
The Frobenius (or Euclidean) norm.
NORM='M'
The value maxi,jaij (not a norm).
Constraint: NORM='1', 'O', 'I', 'F', 'E' or 'M'.
2:     TRANSR – CHARACTER(1)Input
On entry: specifies whether the normal RFP representation of A or its conjugate transpose is stored.
TRANSR='N'
The matrix A is stored in normal RFP format.
TRANSR='C'
The conjugate transpose of the RFP representation of the matrix A is stored.
Constraint: TRANSR='N' or 'C'.
3:     UPLO – CHARACTER(1)Input
On entry: specifies whether the upper or lower triangular part of A is stored.
UPLO='U'
The upper triangular part of A is stored.
UPLO='L'
The lower triangular part of A is stored.
Constraint: UPLO='U' or 'L'.
4:     N – INTEGERInput
On entry: n, the order of the matrix A.
When N=0, F06WNF (ZLANHF) returns zero.
Constraint: N0.
5:     ARN×N+1/2 – COMPLEX (KIND=nag_wp) arrayInput
On entry: the upper or lower triangular part (as specified by UPLO) of the n by n Hermitian matrix A, in either normal or transposed RFP format (as specified by TRANSR). The storage format is described in detail in Section 3.3.3 in the F07 Chapter Introduction.
6:     WORK* – REAL (KIND=nag_wp) arrayWorkspace
Note: the dimension of the array WORK must be at least max1,N  if NORM='1', 'O' or 'I', and at least 1 otherwise.

6  Error Indicators and Warnings

None.

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  Parallelism and Performance

Not applicable.

9  Further Comments

None.

10  Example

This example reads in the lower triangular part of a symmetric matrix, converts this to RFP format, then calculates the norm of the matrix for each of the available norm types.

10.1  Program Text

Program Text (f06wnfe.f90)

10.2  Program Data

Program Data (f06wnfe.d)

10.3  Program Results

Program Results (f06wnfe.r)


F06WNF (ZLANHF) (PDF version)
F06 Chapter Contents
F06 Chapter Introduction
NAG Library Manual

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