1CLANHP(1) LAPACK auxiliary routine (version 3.2) CLANHP(1)
2
3
4
6 CLANHP - returns the value of the one norm, or the Frobenius norm, or
7 the infinity norm, or the element of largest absolute value of a com‐
8 plex hermitian matrix A, supplied in packed form
9
11 REAL FUNCTION CLANHP( NORM, UPLO, N, AP, WORK )
12
13 CHARACTER NORM, UPLO
14
15 INTEGER N
16
17 REAL WORK( * )
18
19 COMPLEX AP( * )
20
22 CLANHP returns the value of the one norm, or the Frobenius norm, or
23 the infinity norm, or the element of largest absolute value of a
24 complex hermitian matrix A, supplied in packed form.
25
27 CLANHP returns the value
28 CLANHP = ( max(abs(A(i,j))), NORM = 'M' or 'm'
29 (
30 ( norm1(A), NORM = '1', 'O' or 'o'
31 (
32 ( normI(A), NORM = 'I' or 'i'
33 (
34 ( normF(A), NORM = 'F', 'f', 'E' or 'e' where
35 norm1 denotes the one norm of a matrix (maximum column sum), normI
36 denotes the infinity norm of a matrix (maximum row sum) and normF
37 denotes the Frobenius norm of a matrix (square root of sum of
38 squares). Note that max(abs(A(i,j))) is not a consistent matrix
39 norm.
40
42 NORM (input) CHARACTER*1
43 Specifies the value to be returned in CLANHP as described
44 above.
45
46 UPLO (input) CHARACTER*1
47 Specifies whether the upper or lower triangular part of the
48 hermitian matrix A is supplied. = 'U': Upper triangular part
49 of A is supplied
50 = 'L': Lower triangular part of A is supplied
51
52 N (input) INTEGER
53 The order of the matrix A. N >= 0. When N = 0, CLANHP is set
54 to zero.
55
56 AP (input) COMPLEX array, dimension (N*(N+1)/2)
57 The upper or lower triangle of the hermitian matrix A, packed
58 columnwise in a linear array. The j-th column of A is stored
59 in the array AP as follows: if UPLO = 'U', AP(i + (j-1)*j/2) =
60 A(i,j) for 1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) =
61 A(i,j) for j<=i<=n. Note that the imaginary parts of the
62 diagonal elements need not be set and are assumed to be zero.
63
64 WORK (workspace) REAL array, dimension (MAX(1,LWORK)),
65 where LWORK >= N when NORM = 'I' or '1' or 'O'; otherwise, WORK
66 is not referenced.
67
68
69
70 LAPACK auxiliary routine (versionNo3v.e2m)ber 2008 CLANHP(1)