1ZLAEV2(1)           LAPACK auxiliary routine (version 3.2)           ZLAEV2(1)
2
3
4

NAME

6       ZLAEV2  -  computes the eigendecomposition of a 2-by-2 Hermitian matrix
7       [ A B ]  [ CONJG(B) C ]
8

SYNOPSIS

10       SUBROUTINE ZLAEV2( A, B, C, RT1, RT2, CS1, SN1 )
11
12           DOUBLE         PRECISION CS1, RT1, RT2
13
14           COMPLEX*16     A, B, C, SN1
15

PURPOSE

17       ZLAEV2 computes the eigendecomposition of a 2-by-2 Hermitian matrix
18          [  A         B  ]
19          [  CONJG(B)  C  ].  On return, RT1 is the eigenvalue of larger abso‐
20       lute  value,  RT2  is  the  eigenvalue  of  smaller absolute value, and
21       (CS1,SN1) is the unit right eigenvector for RT1, giving the  decomposi‐
22       tion
23       [  CS1   CONJG(SN1)  ] [    A     B ] [ CS1 -CONJG(SN1) ] = [ RT1  0  ]
24       [-SN1     CS1     ] [ CONJG(B) C ] [ SN1     CS1     ]   [  0  RT2 ].
25

ARGUMENTS

27       A      (input) COMPLEX*16
28              The (1,1) element of the 2-by-2 matrix.
29
30       B      (input) COMPLEX*16
31              The (1,2) element and the conjugate of the (2,1) element of  the
32              2-by-2 matrix.
33
34       C      (input) COMPLEX*16
35              The (2,2) element of the 2-by-2 matrix.
36
37       RT1    (output) DOUBLE PRECISION
38              The eigenvalue of larger absolute value.
39
40       RT2    (output) DOUBLE PRECISION
41              The eigenvalue of smaller absolute value.
42
43       CS1    (output) DOUBLE PRECISION
44              SN1    (output) COMPLEX*16 The vector (CS1, SN1) is a unit right
45              eigenvector for RT1.
46

FURTHER DETAILS

48       RT1 is accurate to a few ulps barring over/underflow.
49       RT2 may be inaccurate if there is massive cancellation in the  determi‐
50       nant  A*C-B*B; higher precision or correctly rounded or correctly trun‐
51       cated arithmetic would be needed  to  compute  RT2  accurately  in  all
52       cases.
53       CS1  and  SN1 are accurate to a few ulps barring over/underflow.  Over‐
54       flow is possible only if RT1 is within  a  factor  of  5  of  overflow.
55       Underflow is harmless if the input data is 0 or exceeds
56          underflow_threshold / macheps.
57
58
59
60 LAPACK auxiliary routine (versionNo3v.e2m)ber 2008                       ZLAEV2(1)
Impressum