1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
---
:name: zlags2
:md5sum: 3e5b9ee0cacaa512ab8a822358e62805
:category: :subroutine
:arguments:
- upper:
:type: logical
:intent: input
- a1:
:type: doublereal
:intent: input
- a2:
:type: doublecomplex
:intent: input
- a3:
:type: doublereal
:intent: input
- b1:
:type: doublereal
:intent: input
- b2:
:type: doublecomplex
:intent: input
- b3:
:type: doublereal
:intent: input
- csu:
:type: doublereal
:intent: output
- snu:
:type: doublecomplex
:intent: output
- csv:
:type: doublereal
:intent: output
- snv:
:type: doublecomplex
:intent: output
- csq:
:type: doublereal
:intent: output
- snq:
:type: doublecomplex
:intent: output
:substitutions: {}
:fortran_help: " SUBROUTINE ZLAGS2( UPPER, A1, A2, A3, B1, B2, B3, CSU, SNU, CSV, SNV, CSQ, SNQ )\n\n\
* Purpose\n\
* =======\n\
*\n\
* ZLAGS2 computes 2-by-2 unitary matrices U, V and Q, such\n\
* that if ( UPPER ) then\n\
*\n\
* U'*A*Q = U'*( A1 A2 )*Q = ( x 0 )\n\
* ( 0 A3 ) ( x x )\n\
* and\n\
* V'*B*Q = V'*( B1 B2 )*Q = ( x 0 )\n\
* ( 0 B3 ) ( x x )\n\
*\n\
* or if ( .NOT.UPPER ) then\n\
*\n\
* U'*A*Q = U'*( A1 0 )*Q = ( x x )\n\
* ( A2 A3 ) ( 0 x )\n\
* and\n\
* V'*B*Q = V'*( B1 0 )*Q = ( x x )\n\
* ( B2 B3 ) ( 0 x )\n\
* where\n\
*\n\
* U = ( CSU SNU ), V = ( CSV SNV ),\n\
* ( -CONJG(SNU) CSU ) ( -CONJG(SNV) CSV )\n\
*\n\
* Q = ( CSQ SNQ )\n\
* ( -CONJG(SNQ) CSQ )\n\
*\n\
* Z' denotes the conjugate transpose of Z.\n\
*\n\
* The rows of the transformed A and B are parallel. Moreover, if the\n\
* input 2-by-2 matrix A is not zero, then the transformed (1,1) entry\n\
* of A is not zero. If the input matrices A and B are both not zero,\n\
* then the transformed (2,2) element of B is not zero, except when the\n\
* first rows of input A and B are parallel and the second rows are\n\
* zero.\n\
*\n\n\
* Arguments\n\
* =========\n\
*\n\
* UPPER (input) LOGICAL\n\
* = .TRUE.: the input matrices A and B are upper triangular.\n\
* = .FALSE.: the input matrices A and B are lower triangular.\n\
*\n\
* A1 (input) DOUBLE PRECISION\n\
* A2 (input) COMPLEX*16\n\
* A3 (input) DOUBLE PRECISION\n\
* On entry, A1, A2 and A3 are elements of the input 2-by-2\n\
* upper (lower) triangular matrix A.\n\
*\n\
* B1 (input) DOUBLE PRECISION\n\
* B2 (input) COMPLEX*16\n\
* B3 (input) DOUBLE PRECISION\n\
* On entry, B1, B2 and B3 are elements of the input 2-by-2\n\
* upper (lower) triangular matrix B.\n\
*\n\
* CSU (output) DOUBLE PRECISION\n\
* SNU (output) COMPLEX*16\n\
* The desired unitary matrix U.\n\
*\n\
* CSV (output) DOUBLE PRECISION\n\
* SNV (output) COMPLEX*16\n\
* The desired unitary matrix V.\n\
*\n\
* CSQ (output) DOUBLE PRECISION\n\
* SNQ (output) COMPLEX*16\n\
* The desired unitary matrix Q.\n\
*\n\n\
* =====================================================================\n\
*\n"
|
