Estimates of sensitivities of selected eigenvalues and eigenvectors of complex upper triangular matrix

はじめる前に
新機能一覧
Maple ワークシートの作成
Mapleワークシートを共有
Maple ウィンドウのカスタマイズ
Maple システムのカスタマイズ
コネクティビティ
Mathematics
数学
物理
プログラミング
グラフィックス
Student パッケージ
Science and Engineering
科学・エンジニアリング
アプリケーションと例題
リファレンス
システム
MapleSim
MapleSim ツールボックス
MapleSim ヘルプ
Tasks
Toolboxes
- BlockImporter
- Global Optimization
- Grid Computing
- Maple-NAG Connector
  - General information
  - Computational Routines
    - a02 (Complex Arithmetic)
    - c02 (Zeros of Polynomials)
    - c05 (Roots of One or More Transcendental Equations)
    - c06 (Fourier Transforms)
    - d01 (Quadrature)
    - d02 (Ordinary Differential Equations)
    - d03 (Partial Differential Equations)
    - d06 (Mesh Generation)
    - e01 (Interpolation)
    - e02 (Curve and Surface Fitting)
    - e04 (Minimizing or Maximizing a Function)
    - f01 (Matrix Factorizations)
    - f02 (Eigenvalues and Eigenvectors)
    - f03 (Determinants)
    - f04 (Simultaneous Linear Equations)
    - f06 (Linear Algebra Support Functions)
    - f07 (Linear Equations (LAPACK))
    - f08 (Least-squares and Eigenvalue Problems (LAPACK))
      - f08 Introduction
      - f08aec (nag_dgeqrf)
      - f08afc (nag_dorgqr)
      - f08agc (nag_dormqr)
      - f08ahc (nag_dgelqf)
      - f08ajc (nag_dorglq)
      - f08akc (nag_dormlq)
      - f08asc (nag_zgeqrf)
      - f08atc (nag_zungqr)
      - f08auc (nag_zunmqr)
      - f08avc (nag_zgelqf)
      - f08awc (nag_zunglq)
      - f08axc (nag_zunmlq)
      - f08bec (nag_dgeqpf)
      - f08bsc (nag_zgeqpf)
      - f08fcc (nag_dsyevd)
      - f08fec (nag_dsytrd)
      - f08ffc (nag_dorgtr)
      - f08fgc (nag_dormtr)
      - f08fqc (nag_zheevd)
      - f08fsc (nag_zhetrd)
      - f08ftc (nag_zungtr)
      - f08fuc (nag_zunmtr)
      - f08gcc (nag_dspevd)
      - f08gec (nag_dsptrd)
      - f08gfc (nag_dopgtr)
      - f08ggc (nag_dopmtr)
      - f08gqc (nag_zhpevd)
      - f08gsc (nag_zhptrd)
      - f08gtc (nag_zupgtr)
      - f08guc (nag_zupmtr)
      - f08hcc (nag_dsbevd)
      - f08hec (nag_dsbtrd)
      - f08hqc (nag_zhbevd)
      - f08hsc (nag_zhbtrd)
      - f08jcc (nag_dstevd)
      - f08jec (nag_dsteqr)
      - f08jfc (nag_dsterf)
      - f08jgc (nag_dpteqr)
      - f08jjc (nag_dstebz)
      - f08jkc (nag_dstein)
      - f08jsc (nag_zsteqr)
      - f08juc (nag_zpteqr)
      - f08jxc (nag_zstein)
      - f08kec (nag_dgebrd)
      - f08kfc (nag_dorgbr)
      - f08kgc (nag_dormbr)
      - f08ksc (nag_zgebrd)
      - f08ktc (nag_zungbr)
      - f08kuc (nag_zunmbr)
      - f08lec (nag_dgbbrd)
      - f08lsc (nag_zgbbrd)
      - f08mec (nag_dbdsqr)
      - f08msc (nag_zbdsqr)
      - f08nec (nag_dgehrd)
      - f08nfc (nag_dorghr)
      - f08ngc (nag_dormhr)
      - f08nhc (nag_dgebal)
      - f08njc (nag_dgebak)
      - f08nsc (nag_zgehrd)
      - f08ntc (nag_zunghr)
      - f08nuc (nag_zunmhr)
      - f08nvc (nag_zgebal)
      - f08nwc (nag_zgebak)
      - f08pec (nag_dhseqr)
      - f08pkc (nag_dhsein)
      - f08psc (nag_zhseqr)
      - f08pxc (nag_zhsein)
      - f08qfc (nag_dtrexc)
      - f08qgc (nag_dtrsen)
      - f08qhc (nag_dtrsyl)
      - f08qkc (nag_dtrevc)
      - f08qlc (nag_dtrsna)
      - f08qtc (nag_ztrexc)
      - f08quc (nag_ztrsen)
      - f08qvc (nag_ztrsyl)
      - f08qxc (nag_ztrevc)
      - f08qyc (nag_ztrsna)
      - f08sec (nag_dsygst)
      - f08ssc (nag_zhegst)
      - f08tec (nag_dspgst)
      - f08tsc (nag_zhpgst)
      - f08uec (nag_dsbgst)
      - f08ufc (nag_dpbstf)
      - f08usc (nag_zhbgst)
      - f08utc (nag_zpbstf)
      - f08wec (nag_dgghrd)
      - f08whc (nag_dggbal)
      - f08wjc (nag_dggbak)
      - f08wsc (nag_zgghrd)
      - f08wvc (nag_zggbal)
      - f08wwc (nag_zggbak)
      - f08xec (nag_dhgeqz)
      - f08xsc (nag_zhgeqz)
      - f08ykc (nag_dtgevc)
      - f08yxc (nag_ztgevc)
    - f11 (Large Scale Linear Systems)
    - f16 (NAG Interface to BLAS)
    - g01 (Simple Calculations on Statistical Data)
    - g02 (Correlation and Regression Analysis)
    - g03 (Multivariate Methods)
    - g04 (Analysis of Variance)
    - g05 (Random Number Generators)
    - g07 (Univariate Estimation)
    - g08 (Nonparametric Statistics)
    - g10 (Smoothing in Statistics)
    - g11 (Contingency Table Analysis)
    - g12 (Survival Analysis)
    - g13 (Time Series Analysis)
    - h (Operations Research)
    - m01 (Sorting)
    - s (Approximations of Special Functions)
    - x01 (Mathematical Constants)
    - x02 (Machine Constants)
    - x04 (Input-Output Utilities)
  - Getting Started
  - Example Worksheet
エラーメッセージガイド
マニュアル
数学アプリ

NAG[f08qyc] NAG[nag_ztrsna] - Estimates of sensitivities of selected eigenvalues and eigenvectors of complex upper triangular matrix

Calling Sequence

f08qyc(job, how_many, select, t, vl, vr, s, sep, m, 'n'=n, 'mm'=mm, 'fail'=fail)

nag_ztrsna(. . .)

Parameters

job - String;

On entry: indicates whether condition numbers are required for eigenvalues and/or eigenvectors.

Condition numbers for eigenvalues only are computed.

Condition numbers for eigenvectors only are computed.

Condition numbers for both eigenvalues and eigenvectors are computed.

Constraint: "Nag_EigVals", "Nag_EigVecs" or "Nag_DoBoth". .

how_many - String;

On entry: indicates how many condition numbers are to be computed.

Condition numbers for all eigenpairs are computed.

Condition numbers for selected eigenpairs (as specified by select) are computed.

Constraint: "Nag_ComputeAll" or "Nag_ComputeSelected". .

select - Vector(1..dim, datatype=boolean);

Note: the dimension, dim, of the array select must be at least

when ;

otherwise.

On entry: specifies the eigenpairs for which condition numbers are to be computed if . To select condition numbers for the eigenpair corresponding to the eigenvalue , must be set to true.

If , select is not referenced.

t - Matrix(1..dim1, 1..dim2, datatype=complex[8], order=order);

Note: this array may be supplied in Fortran_order or C_order , as specified by order. All array parameters must use a consistent order.

On entry: the by upper triangular matrix , as returned by f08psc (nag_zhseqr).

vl - Matrix(1..dim1, 1..dim2, datatype=complex[8], order=order);

Note: this array may be supplied in Fortran_order or C_order , as specified by order. All array parameters must use a consistent order.

On entry: if "Nag_EigVals" or "Nag_DoBoth", vl must contain the left eigenvectors of (or of any matrix with unitary) corresponding to the eigenpairs specified by how_many and select. The eigenvectors must be stored in consecutive rows or columns (depending on the value of storage order) of vl, as returned by f08qxc (nag_ztrevc) or f08pxc (nag_zhsein).

If , vl is not referenced.

vr - Matrix(1..dim1, 1..dim2, datatype=complex[8], order=order);

Note: this array may be supplied in Fortran_order or C_order , as specified by order. All array parameters must use a consistent order.

On entry: if "Nag_EigVals" or "Nag_DoBoth", vr must contain the right eigenvectors of (or of any matrix with unitary) corresponding to the eigenpairs specified by how_many and select. The eigenvectors must be stored in consecutive rows or columns (depending on the value of storage order) of vr, as returned by f08qxc (nag_ztrevc) or f08pxc (nag_zhsein).

If , vr is not referenced.

s - Vector(1..dim, datatype=float[8]);

Note: the dimension, dim, of the array s must be at least

when "Nag_EigVals" or "Nag_DoBoth";

when .

On exit: the reciprocal condition numbers of the selected eigenvalues if "Nag_EigVals" or "Nag_DoBoth", stored in consecutive elements of the array. Thus , and the th rows or columns of vl and vr all correspond to the same eigenpair (but not in general the th eigenpair unless all eigenpairs have been selected).

s is not referenced if .

sep - Vector(1..dim, datatype=float[8]);

Note: the dimension, dim, of the array sep must be at least

when "Nag_EigVecs" or "Nag_DoBoth";

when .

On exit: the estimated reciprocal condition numbers of the selected right eigenvectors if "Nag_EigVecs" or "Nag_DoBoth", stored in consecutive elements of the array.

If , sep is not referenced i.

m - assignable;

Note: On exit the variable m will have a value of type integer.

On exit: , the number of selected eigenpairs. If , m is set to .

'n'=n - integer; (optional)

Default value: the dimension of the array t.

On entry: , the order of the matrix .

Constraint: . .

'mm'=mm - integer; (optional)

Default value: the second dimension of the arrays vl, vr.

On entry: the number of elements in the arrays s and sep, and the number of rows or columns (depending on the value of storage order) in the arrays vl and vr (if used). The precise number required, , is if ; if , is the number of selected eigenpairs (see select), in which case .

Constraint: . .

'fail'=fail - table; (optional)

The NAG error argument, see the documentation for NagError.

Description

	Purpose
	nag_ztrsna (f08qyc) estimates condition numbers for specified eigenvalues and/or right eigenvectors of a complex upper triangular matrix.

Description

nag_ztrsna (f08qyc) estimates condition numbers for specified eigenvalues and/or right eigenvectors of a complex upper triangular matrix . These are the same as the condition numbers of the eigenvalues and right eigenvectors of an original matrix (with unitary ), from which may have been derived.

nag_ztrsna (f08qyc) computes the reciprocal of the condition number of an eigenvalue as

s[i]=(|v^H,u|)/(||u||[E]||v||[E])

where and are the right and left eigenvectors of , respectively, corresponding to . This reciprocal condition number always lies between zero (i.e., ill-conditioned) and one (i.e., well-conditioned).

An approximate error estimate for a computed eigenvalue is then given by

epsilon*LinearAlgebra[Norm](T)/s[i]

where is the machine precision.

To estimate the reciprocal of the condition number of the right eigenvector corresponding to , the function first calls f08qtc (nag_ztrexc) to reorder the eigenvalues so that is in the leading position:

The reciprocal condition number of the eigenvector is then estimated as , the smallest singular value of the matrix . This number ranges from zero (i.e., ill-conditioned) to very large (i.e., well-conditioned).

An approximate error estimate for a computed right eigenvector corresponding to is then given by

epsilon*LinearAlgebra[Norm](T)/sep[i]

	Error Indicators and Warnings
	"NE_ALLOC_FAIL" Dynamic memory allocation failed. "NE_BAD_PARAM" On entry, argument had an illegal value. "NE_INT" On entry, . Constraint: , where is the number of selected eigenpairs. On entry, . Constraint: . "NE_INTERNAL_ERROR" An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please consult NAG for assistance.

	Accuracy
	The computed values may over estimate the true value, but seldom by a factor of more than 3.

	Further Comments
	The real analogue of this function is f08qlc (nag_dtrsna).

Examples

job := "Nag_DoBoth":
how_many := "Nag_ComputeAll":
n := 4:
mm := 4:
selector := Vector([false], datatype=boolean):
t := Matrix([[-6.0004 -6.9999*I , 0.3637 -0.3656*I , -0.188 +0.4787*I , 0.8784999999999999 -0.2539*I ], [0 +0*I , -5 +2.006*I , -0.0307 -0.7217*I , -0.229 +0.1313*I ], [0 +0*I , 0 +0*I , 7.9982 -0.9964*I , 0.9357 +0.5359*I ], [0 +0*I , 0 +0*I , 0 +0*I , 3.0023 -3.9998*I ]], datatype=complex[8]):
vl := Matrix([[1 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [0.03566943515948521 -0.04434689513913646*I , 1 +0*I , 0 +0*I , 0 +0*I ], [-0.002207376074206878 +0.03131341253433394*I , -0.009933197113414047 -0.05322864465746686*I , 1 +0*I , 0 +0*I ], [-0.07824388063236407 -0.05827078558456029*I , 0.03187495330451339 -0.001955906687244087*I , 0.1849408889864735 +0.003913502268254905*I , 1 +0*I ]], datatype=complex[8]):
vr := Matrix([[1 +0*I , -0.03566943515948521 -0.04434689513913646*I , -0.0005074605791794693 +0.03277154177278578*I , 0.0792597025312869 -0.06285024830308554*I ], [0 +0*I , 1 +0*I , 0.009933197113414047 -0.05322864465746686*I , -0.03350369718754292 +0.007927119764687306*I ], [0 +0*I , 0 +0*I , 1 +0*I , -0.1849408889864735 +0.003913502268254905*I ], [0 +0*I , 0 +0*I , 0 +0*I , 1 +0*I ]], datatype=complex[8]):
s := Vector(4, datatype=float[8]):
sep := Vector(4, datatype=float[8]):
NAG:-f08qyc(job, how_many, selector, t, vl, vr, s, sep, m, 'n' = n, 'mm' = mm):

sep;

	See Also
	Golub G H and Van Loan C F (1996) Matrix Computations (3rd Edition) Johns Hopkins University Press, Baltimore f08 Chapter Introduction. NAG Toolbox Overview. NAG Web Site.