1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
|
---
:name: zlantp
:md5sum: 5643fadc8657c6b20518835ef43210f4
:category: :function
:type: doublereal
:arguments:
- norm:
:type: char
:intent: input
- uplo:
:type: char
:intent: input
- diag:
:type: char
:intent: input
- n:
:type: integer
:intent: input
- ap:
:type: doublecomplex
:intent: input
:dims:
- n*(n+1)/2
- work:
:type: doublereal
:intent: workspace
:dims:
- MAX(1,lwork)
:substitutions:
lwork: "lsame_(&norm,\"I\") ? n : 0"
:fortran_help: " DOUBLE PRECISION FUNCTION ZLANTP( NORM, UPLO, DIAG, N, AP, WORK )\n\n\
* Purpose\n\
* =======\n\
*\n\
* ZLANTP returns the value of the one norm, or the Frobenius norm, or\n\
* the infinity norm, or the element of largest absolute value of a\n\
* triangular matrix A, supplied in packed form.\n\
*\n\
* Description\n\
* ===========\n\
*\n\
* ZLANTP returns the value\n\
*\n\
* ZLANTP = ( max(abs(A(i,j))), NORM = 'M' or 'm'\n\
* (\n\
* ( norm1(A), NORM = '1', 'O' or 'o'\n\
* (\n\
* ( normI(A), NORM = 'I' or 'i'\n\
* (\n\
* ( normF(A), NORM = 'F', 'f', 'E' or 'e'\n\
*\n\
* where norm1 denotes the one norm of a matrix (maximum column sum),\n\
* normI denotes the infinity norm of a matrix (maximum row sum) and\n\
* normF denotes the Frobenius norm of a matrix (square root of sum of\n\
* squares). Note that max(abs(A(i,j))) is not a consistent matrix norm.\n\
*\n\n\
* Arguments\n\
* =========\n\
*\n\
* NORM (input) CHARACTER*1\n\
* Specifies the value to be returned in ZLANTP as described\n\
* above.\n\
*\n\
* UPLO (input) CHARACTER*1\n\
* Specifies whether the matrix A is upper or lower triangular.\n\
* = 'U': Upper triangular\n\
* = 'L': Lower triangular\n\
*\n\
* DIAG (input) CHARACTER*1\n\
* Specifies whether or not the matrix A is unit triangular.\n\
* = 'N': Non-unit triangular\n\
* = 'U': Unit triangular\n\
*\n\
* N (input) INTEGER\n\
* The order of the matrix A. N >= 0. When N = 0, ZLANTP is\n\
* set to zero.\n\
*\n\
* AP (input) COMPLEX*16 array, dimension (N*(N+1)/2)\n\
* The upper or lower triangular matrix A, packed columnwise in\n\
* a linear array. The j-th column of A is stored in the array\n\
* AP as follows:\n\
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;\n\
* if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.\n\
* Note that when DIAG = 'U', the elements of the array AP\n\
* corresponding to the diagonal elements of the matrix A are\n\
* not referenced, but are assumed to be one.\n\
*\n\
* WORK (workspace) DOUBLE PRECISION array, dimension (MAX(1,LWORK)),\n\
* where LWORK >= N when NORM = 'I'; otherwise, WORK is not\n\
* referenced.\n\
*\n\n\
* =====================================================================\n\
*\n"
|
