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

NAME

6       CLAESY  -  computes the eigendecomposition of a 2-by-2 symmetric matrix
7       ( ( A, B );( B, C ) ) provided the norm of the matrix  of  eigenvectors
8       is larger than some threshold value
9

SYNOPSIS

11       SUBROUTINE CLAESY( A, B, C, RT1, RT2, EVSCAL, CS1, SN1 )
12
13           COMPLEX        A, B, C, CS1, EVSCAL, RT1, RT2, SN1
14

PURPOSE

16       CLAESY computes the eigendecomposition of a 2-by-2 symmetric matrix
17          (  (  A, B );( B, C ) ) provided the norm of the matrix of eigenvec‐
18       tors is larger than some threshold value.  RT1  is  the  eigenvalue  of
19       larger  absolute  value,  and  RT2  of  smaller absolute value.  If the
20       eigenvectors are computed, then on return ( CS1,  SN1  )  is  the  unit
21       eigenvector  for  RT1,  hence  [   CS1      SN1    ] . [ A  B ] . [ CS1
22       -SN1   ] = [ RT1  0  ] [ -SN1     CS1   ]   [ B  C ]    [  SN1      CS1
23       ]   [  0  RT2 ]
24

ARGUMENTS

26       A       (input) COMPLEX
27               The ( 1, 1 ) element of input matrix.
28
29       B       (input) COMPLEX
30               The  ( 1, 2 ) element of input matrix.  The ( 2, 1 ) element is
31               also given by B, since the 2-by-2 matrix is symmetric.
32
33       C       (input) COMPLEX
34               The ( 2, 2 ) element of input matrix.
35
36       RT1     (output) COMPLEX
37               The eigenvalue of larger modulus.
38
39       RT2     (output) COMPLEX
40               The eigenvalue of smaller modulus.
41
42       EVSCAL  (output) COMPLEX
43               The complex value by which the eigenvector matrix was scaled to
44               make  it orthonormal.  If EVSCAL is zero, the eigenvectors were
45               not computed.  This means one of two things:  the 2-by-2 matrix
46               could  not be diagonalized, or the norm of the matrix of eigen‐
47               vectors before scaling was  larger  than  the  threshold  value
48               THRESH (set below).
49
50       CS1     (output) COMPLEX
51               SN1     (output) COMPLEX If EVSCAL .NE. 0,  ( CS1, SN1 ) is the
52               unit right eigenvector for RT1.
53
54
55
56 LAPACK auxiliary routine (versionNo3v.e2m)ber 2008                       CLAESY(1)
Impressum