Sponsor:

Your company here, and a link to your site. Click to find out more.

cpst01.f - Man Page

TESTING/LIN/cpst01.f

Synopsis

Functions/Subroutines

subroutine cpst01 (uplo, n, a, lda, afac, ldafac, perm, ldperm, piv, rwork, resid, rank)
CPST01

Function/Subroutine Documentation

subroutine cpst01 (character uplo, integer n, complex, dimension( lda, * ) a, integer lda, complex, dimension( ldafac, * ) afac, integer ldafac, complex, dimension( ldperm, * ) perm, integer ldperm, integer, dimension( * ) piv, real, dimension( * ) rwork, real resid, integer rank)

CPST01

Purpose:

 CPST01 reconstructs an Hermitian positive semidefinite matrix A
 from its L or U factors and the permutation matrix P and computes
 the residual
    norm( P*L*L'*P' - A ) / ( N * norm(A) * EPS ) or
    norm( P*U'*U*P' - A ) / ( N * norm(A) * EPS ),
 where EPS is the machine epsilon, L' is the conjugate transpose of L,
 and U' is the conjugate transpose of U.
Parameters

UPLO

          UPLO is CHARACTER*1
          Specifies whether the upper or lower triangular part of the
          Hermitian matrix A is stored:
          = 'U':  Upper triangular
          = 'L':  Lower triangular

N

          N is INTEGER
          The number of rows and columns of the matrix A.  N >= 0.

A

          A is COMPLEX array, dimension (LDA,N)
          The original Hermitian matrix A.

LDA

          LDA is INTEGER
          The leading dimension of the array A.  LDA >= max(1,N)

AFAC

          AFAC is COMPLEX array, dimension (LDAFAC,N)
          The factor L or U from the L*L' or U'*U
          factorization of A.

LDAFAC

          LDAFAC is INTEGER
          The leading dimension of the array AFAC.  LDAFAC >= max(1,N).

PERM

          PERM is COMPLEX array, dimension (LDPERM,N)
          Overwritten with the reconstructed matrix, and then with the
          difference P*L*L'*P' - A (or P*U'*U*P' - A)

LDPERM

          LDPERM is INTEGER
          The leading dimension of the array PERM.
          LDAPERM >= max(1,N).

PIV

          PIV is INTEGER array, dimension (N)
          PIV is such that the nonzero entries are
          P( PIV( K ), K ) = 1.

RWORK

          RWORK is REAL array, dimension (N)

RESID

          RESID is REAL
          If UPLO = 'L', norm(L*L' - A) / ( N * norm(A) * EPS )
          If UPLO = 'U', norm(U'*U - A) / ( N * norm(A) * EPS )

RANK

          RANK is INTEGER
          number of nonzero singular values of A.
Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Definition at line 134 of file cpst01.f.

Author

Generated automatically by Doxygen for LAPACK from the source code.

Referenced By

The man page cpst01(3) is an alias of cpst01.f(3).

Tue Nov 28 2023 12:08:42 Version 3.12.0 LAPACK