f06waf (dlansf) returns the value of the -norm, the -norm, the Frobenius norm, or the maximum absolute value of the elements of a real symmetric matrix stored in Rectangular Full Packed (RFP) format.
The routine may be called by its
LAPACK
name dlansf.
3
Description
Given a real by symmetric matrix, ,
f06waf (dlansf) calculates one of the values given by
(the -norm of ),
(the -norm of ),
(the Frobenius norm of ), or
(the maximum absolute element value of ).
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. Software37, 2
5
Arguments
1: – Character(1)Input
On entry: specifies the value to be returned.
or
The -norm.
The -norm.
or
The Frobenius (or Euclidean) norm.
The value (not a norm).
Constraint:
, , , , or .
2: – Character(1)Input
On entry: specifies whether the RFP representation of is normal or transposed.
The matrix is stored in normal RFP format.
The matrix is stored in transposed RFP format.
Constraint:
or .
3: – Character(1)Input
On entry: specifies whether the upper or lower triangular part of is stored.
The upper triangular part of is stored.
The lower triangular part of is stored.
Constraint:
or .
4: – IntegerInput
On entry: , the order of the matrix .
When , f06waf (dlansf) returns zero.
Constraint:
.
5: – Real (Kind=nag_wp) arrayInput
On entry: the upper or lower triangular part (as specified by uplo) of the by symmetric matrix , 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: – Real (Kind=nag_wp) arrayWorkspace
Note: the dimension of the array work
must be at least
if , or , and at least otherwise.
f06waf (dlansf) is not threaded in any implementation.
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.