|
2329
|
1 SUBROUTINE DLARFX( SIDE, M, N, V, TAU, C, LDC, WORK ) |
|
|
2 * |
|
3333
|
3 * -- LAPACK auxiliary routine (version 3.0) -- |
|
2329
|
4 * Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd., |
|
|
5 * Courant Institute, Argonne National Lab, and Rice University |
|
|
6 * February 29, 1992 |
|
|
7 * |
|
|
8 * .. Scalar Arguments .. |
|
|
9 CHARACTER SIDE |
|
|
10 INTEGER LDC, M, N |
|
|
11 DOUBLE PRECISION TAU |
|
|
12 * .. |
|
|
13 * .. Array Arguments .. |
|
|
14 DOUBLE PRECISION C( LDC, * ), V( * ), WORK( * ) |
|
|
15 * .. |
|
|
16 * |
|
|
17 * Purpose |
|
|
18 * ======= |
|
|
19 * |
|
|
20 * DLARFX applies a real elementary reflector H to a real m by n |
|
|
21 * matrix C, from either the left or the right. H is represented in the |
|
|
22 * form |
|
|
23 * |
|
|
24 * H = I - tau * v * v' |
|
|
25 * |
|
|
26 * where tau is a real scalar and v is a real vector. |
|
|
27 * |
|
|
28 * If tau = 0, then H is taken to be the unit matrix |
|
|
29 * |
|
|
30 * This version uses inline code if H has order < 11. |
|
|
31 * |
|
|
32 * Arguments |
|
|
33 * ========= |
|
|
34 * |
|
|
35 * SIDE (input) CHARACTER*1 |
|
|
36 * = 'L': form H * C |
|
|
37 * = 'R': form C * H |
|
|
38 * |
|
|
39 * M (input) INTEGER |
|
|
40 * The number of rows of the matrix C. |
|
|
41 * |
|
|
42 * N (input) INTEGER |
|
|
43 * The number of columns of the matrix C. |
|
|
44 * |
|
|
45 * V (input) DOUBLE PRECISION array, dimension (M) if SIDE = 'L' |
|
|
46 * or (N) if SIDE = 'R' |
|
|
47 * The vector v in the representation of H. |
|
|
48 * |
|
|
49 * TAU (input) DOUBLE PRECISION |
|
|
50 * The value tau in the representation of H. |
|
|
51 * |
|
|
52 * C (input/output) DOUBLE PRECISION array, dimension (LDC,N) |
|
|
53 * On entry, the m by n matrix C. |
|
|
54 * On exit, C is overwritten by the matrix H * C if SIDE = 'L', |
|
|
55 * or C * H if SIDE = 'R'. |
|
|
56 * |
|
|
57 * LDC (input) INTEGER |
|
|
58 * The leading dimension of the array C. LDA >= (1,M). |
|
|
59 * |
|
|
60 * WORK (workspace) DOUBLE PRECISION array, dimension |
|
|
61 * (N) if SIDE = 'L' |
|
|
62 * or (M) if SIDE = 'R' |
|
|
63 * WORK is not referenced if H has order < 11. |
|
|
64 * |
|
|
65 * ===================================================================== |
|
|
66 * |
|
|
67 * .. Parameters .. |
|
|
68 DOUBLE PRECISION ZERO, ONE |
|
|
69 PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 ) |
|
|
70 * .. |
|
|
71 * .. Local Scalars .. |
|
|
72 INTEGER J |
|
|
73 DOUBLE PRECISION SUM, T1, T10, T2, T3, T4, T5, T6, T7, T8, T9, |
|
|
74 $ V1, V10, V2, V3, V4, V5, V6, V7, V8, V9 |
|
|
75 * .. |
|
|
76 * .. External Functions .. |
|
|
77 LOGICAL LSAME |
|
|
78 EXTERNAL LSAME |
|
|
79 * .. |
|
|
80 * .. External Subroutines .. |
|
|
81 EXTERNAL DGEMV, DGER |
|
|
82 * .. |
|
|
83 * .. Executable Statements .. |
|
|
84 * |
|
|
85 IF( TAU.EQ.ZERO ) |
|
|
86 $ RETURN |
|
|
87 IF( LSAME( SIDE, 'L' ) ) THEN |
|
|
88 * |
|
|
89 * Form H * C, where H has order m. |
|
|
90 * |
|
|
91 GO TO ( 10, 30, 50, 70, 90, 110, 130, 150, |
|
|
92 $ 170, 190 )M |
|
|
93 * |
|
|
94 * Code for general M |
|
|
95 * |
|
|
96 * w := C'*v |
|
|
97 * |
|
|
98 CALL DGEMV( 'Transpose', M, N, ONE, C, LDC, V, 1, ZERO, WORK, |
|
|
99 $ 1 ) |
|
|
100 * |
|
|
101 * C := C - tau * v * w' |
|
|
102 * |
|
|
103 CALL DGER( M, N, -TAU, V, 1, WORK, 1, C, LDC ) |
|
|
104 GO TO 410 |
|
|
105 10 CONTINUE |
|
|
106 * |
|
|
107 * Special code for 1 x 1 Householder |
|
|
108 * |
|
|
109 T1 = ONE - TAU*V( 1 )*V( 1 ) |
|
|
110 DO 20 J = 1, N |
|
|
111 C( 1, J ) = T1*C( 1, J ) |
|
|
112 20 CONTINUE |
|
|
113 GO TO 410 |
|
|
114 30 CONTINUE |
|
|
115 * |
|
|
116 * Special code for 2 x 2 Householder |
|
|
117 * |
|
|
118 V1 = V( 1 ) |
|
|
119 T1 = TAU*V1 |
|
|
120 V2 = V( 2 ) |
|
|
121 T2 = TAU*V2 |
|
|
122 DO 40 J = 1, N |
|
|
123 SUM = V1*C( 1, J ) + V2*C( 2, J ) |
|
|
124 C( 1, J ) = C( 1, J ) - SUM*T1 |
|
|
125 C( 2, J ) = C( 2, J ) - SUM*T2 |
|
|
126 40 CONTINUE |
|
|
127 GO TO 410 |
|
|
128 50 CONTINUE |
|
|
129 * |
|
|
130 * Special code for 3 x 3 Householder |
|
|
131 * |
|
|
132 V1 = V( 1 ) |
|
|
133 T1 = TAU*V1 |
|
|
134 V2 = V( 2 ) |
|
|
135 T2 = TAU*V2 |
|
|
136 V3 = V( 3 ) |
|
|
137 T3 = TAU*V3 |
|
|
138 DO 60 J = 1, N |
|
|
139 SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) |
|
|
140 C( 1, J ) = C( 1, J ) - SUM*T1 |
|
|
141 C( 2, J ) = C( 2, J ) - SUM*T2 |
|
|
142 C( 3, J ) = C( 3, J ) - SUM*T3 |
|
|
143 60 CONTINUE |
|
|
144 GO TO 410 |
|
|
145 70 CONTINUE |
|
|
146 * |
|
|
147 * Special code for 4 x 4 Householder |
|
|
148 * |
|
|
149 V1 = V( 1 ) |
|
|
150 T1 = TAU*V1 |
|
|
151 V2 = V( 2 ) |
|
|
152 T2 = TAU*V2 |
|
|
153 V3 = V( 3 ) |
|
|
154 T3 = TAU*V3 |
|
|
155 V4 = V( 4 ) |
|
|
156 T4 = TAU*V4 |
|
|
157 DO 80 J = 1, N |
|
|
158 SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + |
|
|
159 $ V4*C( 4, J ) |
|
|
160 C( 1, J ) = C( 1, J ) - SUM*T1 |
|
|
161 C( 2, J ) = C( 2, J ) - SUM*T2 |
|
|
162 C( 3, J ) = C( 3, J ) - SUM*T3 |
|
|
163 C( 4, J ) = C( 4, J ) - SUM*T4 |
|
|
164 80 CONTINUE |
|
|
165 GO TO 410 |
|
|
166 90 CONTINUE |
|
|
167 * |
|
|
168 * Special code for 5 x 5 Householder |
|
|
169 * |
|
|
170 V1 = V( 1 ) |
|
|
171 T1 = TAU*V1 |
|
|
172 V2 = V( 2 ) |
|
|
173 T2 = TAU*V2 |
|
|
174 V3 = V( 3 ) |
|
|
175 T3 = TAU*V3 |
|
|
176 V4 = V( 4 ) |
|
|
177 T4 = TAU*V4 |
|
|
178 V5 = V( 5 ) |
|
|
179 T5 = TAU*V5 |
|
|
180 DO 100 J = 1, N |
|
|
181 SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + |
|
|
182 $ V4*C( 4, J ) + V5*C( 5, J ) |
|
|
183 C( 1, J ) = C( 1, J ) - SUM*T1 |
|
|
184 C( 2, J ) = C( 2, J ) - SUM*T2 |
|
|
185 C( 3, J ) = C( 3, J ) - SUM*T3 |
|
|
186 C( 4, J ) = C( 4, J ) - SUM*T4 |
|
|
187 C( 5, J ) = C( 5, J ) - SUM*T5 |
|
|
188 100 CONTINUE |
|
|
189 GO TO 410 |
|
|
190 110 CONTINUE |
|
|
191 * |
|
|
192 * Special code for 6 x 6 Householder |
|
|
193 * |
|
|
194 V1 = V( 1 ) |
|
|
195 T1 = TAU*V1 |
|
|
196 V2 = V( 2 ) |
|
|
197 T2 = TAU*V2 |
|
|
198 V3 = V( 3 ) |
|
|
199 T3 = TAU*V3 |
|
|
200 V4 = V( 4 ) |
|
|
201 T4 = TAU*V4 |
|
|
202 V5 = V( 5 ) |
|
|
203 T5 = TAU*V5 |
|
|
204 V6 = V( 6 ) |
|
|
205 T6 = TAU*V6 |
|
|
206 DO 120 J = 1, N |
|
|
207 SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + |
|
|
208 $ V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) |
|
|
209 C( 1, J ) = C( 1, J ) - SUM*T1 |
|
|
210 C( 2, J ) = C( 2, J ) - SUM*T2 |
|
|
211 C( 3, J ) = C( 3, J ) - SUM*T3 |
|
|
212 C( 4, J ) = C( 4, J ) - SUM*T4 |
|
|
213 C( 5, J ) = C( 5, J ) - SUM*T5 |
|
|
214 C( 6, J ) = C( 6, J ) - SUM*T6 |
|
|
215 120 CONTINUE |
|
|
216 GO TO 410 |
|
|
217 130 CONTINUE |
|
|
218 * |
|
|
219 * Special code for 7 x 7 Householder |
|
|
220 * |
|
|
221 V1 = V( 1 ) |
|
|
222 T1 = TAU*V1 |
|
|
223 V2 = V( 2 ) |
|
|
224 T2 = TAU*V2 |
|
|
225 V3 = V( 3 ) |
|
|
226 T3 = TAU*V3 |
|
|
227 V4 = V( 4 ) |
|
|
228 T4 = TAU*V4 |
|
|
229 V5 = V( 5 ) |
|
|
230 T5 = TAU*V5 |
|
|
231 V6 = V( 6 ) |
|
|
232 T6 = TAU*V6 |
|
|
233 V7 = V( 7 ) |
|
|
234 T7 = TAU*V7 |
|
|
235 DO 140 J = 1, N |
|
|
236 SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + |
|
|
237 $ V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) + |
|
|
238 $ V7*C( 7, J ) |
|
|
239 C( 1, J ) = C( 1, J ) - SUM*T1 |
|
|
240 C( 2, J ) = C( 2, J ) - SUM*T2 |
|
|
241 C( 3, J ) = C( 3, J ) - SUM*T3 |
|
|
242 C( 4, J ) = C( 4, J ) - SUM*T4 |
|
|
243 C( 5, J ) = C( 5, J ) - SUM*T5 |
|
|
244 C( 6, J ) = C( 6, J ) - SUM*T6 |
|
|
245 C( 7, J ) = C( 7, J ) - SUM*T7 |
|
|
246 140 CONTINUE |
|
|
247 GO TO 410 |
|
|
248 150 CONTINUE |
|
|
249 * |
|
|
250 * Special code for 8 x 8 Householder |
|
|
251 * |
|
|
252 V1 = V( 1 ) |
|
|
253 T1 = TAU*V1 |
|
|
254 V2 = V( 2 ) |
|
|
255 T2 = TAU*V2 |
|
|
256 V3 = V( 3 ) |
|
|
257 T3 = TAU*V3 |
|
|
258 V4 = V( 4 ) |
|
|
259 T4 = TAU*V4 |
|
|
260 V5 = V( 5 ) |
|
|
261 T5 = TAU*V5 |
|
|
262 V6 = V( 6 ) |
|
|
263 T6 = TAU*V6 |
|
|
264 V7 = V( 7 ) |
|
|
265 T7 = TAU*V7 |
|
|
266 V8 = V( 8 ) |
|
|
267 T8 = TAU*V8 |
|
|
268 DO 160 J = 1, N |
|
|
269 SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + |
|
|
270 $ V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) + |
|
|
271 $ V7*C( 7, J ) + V8*C( 8, J ) |
|
|
272 C( 1, J ) = C( 1, J ) - SUM*T1 |
|
|
273 C( 2, J ) = C( 2, J ) - SUM*T2 |
|
|
274 C( 3, J ) = C( 3, J ) - SUM*T3 |
|
|
275 C( 4, J ) = C( 4, J ) - SUM*T4 |
|
|
276 C( 5, J ) = C( 5, J ) - SUM*T5 |
|
|
277 C( 6, J ) = C( 6, J ) - SUM*T6 |
|
|
278 C( 7, J ) = C( 7, J ) - SUM*T7 |
|
|
279 C( 8, J ) = C( 8, J ) - SUM*T8 |
|
|
280 160 CONTINUE |
|
|
281 GO TO 410 |
|
|
282 170 CONTINUE |
|
|
283 * |
|
|
284 * Special code for 9 x 9 Householder |
|
|
285 * |
|
|
286 V1 = V( 1 ) |
|
|
287 T1 = TAU*V1 |
|
|
288 V2 = V( 2 ) |
|
|
289 T2 = TAU*V2 |
|
|
290 V3 = V( 3 ) |
|
|
291 T3 = TAU*V3 |
|
|
292 V4 = V( 4 ) |
|
|
293 T4 = TAU*V4 |
|
|
294 V5 = V( 5 ) |
|
|
295 T5 = TAU*V5 |
|
|
296 V6 = V( 6 ) |
|
|
297 T6 = TAU*V6 |
|
|
298 V7 = V( 7 ) |
|
|
299 T7 = TAU*V7 |
|
|
300 V8 = V( 8 ) |
|
|
301 T8 = TAU*V8 |
|
|
302 V9 = V( 9 ) |
|
|
303 T9 = TAU*V9 |
|
|
304 DO 180 J = 1, N |
|
|
305 SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + |
|
|
306 $ V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) + |
|
|
307 $ V7*C( 7, J ) + V8*C( 8, J ) + V9*C( 9, J ) |
|
|
308 C( 1, J ) = C( 1, J ) - SUM*T1 |
|
|
309 C( 2, J ) = C( 2, J ) - SUM*T2 |
|
|
310 C( 3, J ) = C( 3, J ) - SUM*T3 |
|
|
311 C( 4, J ) = C( 4, J ) - SUM*T4 |
|
|
312 C( 5, J ) = C( 5, J ) - SUM*T5 |
|
|
313 C( 6, J ) = C( 6, J ) - SUM*T6 |
|
|
314 C( 7, J ) = C( 7, J ) - SUM*T7 |
|
|
315 C( 8, J ) = C( 8, J ) - SUM*T8 |
|
|
316 C( 9, J ) = C( 9, J ) - SUM*T9 |
|
|
317 180 CONTINUE |
|
|
318 GO TO 410 |
|
|
319 190 CONTINUE |
|
|
320 * |
|
|
321 * Special code for 10 x 10 Householder |
|
|
322 * |
|
|
323 V1 = V( 1 ) |
|
|
324 T1 = TAU*V1 |
|
|
325 V2 = V( 2 ) |
|
|
326 T2 = TAU*V2 |
|
|
327 V3 = V( 3 ) |
|
|
328 T3 = TAU*V3 |
|
|
329 V4 = V( 4 ) |
|
|
330 T4 = TAU*V4 |
|
|
331 V5 = V( 5 ) |
|
|
332 T5 = TAU*V5 |
|
|
333 V6 = V( 6 ) |
|
|
334 T6 = TAU*V6 |
|
|
335 V7 = V( 7 ) |
|
|
336 T7 = TAU*V7 |
|
|
337 V8 = V( 8 ) |
|
|
338 T8 = TAU*V8 |
|
|
339 V9 = V( 9 ) |
|
|
340 T9 = TAU*V9 |
|
|
341 V10 = V( 10 ) |
|
|
342 T10 = TAU*V10 |
|
|
343 DO 200 J = 1, N |
|
|
344 SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + |
|
|
345 $ V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) + |
|
|
346 $ V7*C( 7, J ) + V8*C( 8, J ) + V9*C( 9, J ) + |
|
|
347 $ V10*C( 10, J ) |
|
|
348 C( 1, J ) = C( 1, J ) - SUM*T1 |
|
|
349 C( 2, J ) = C( 2, J ) - SUM*T2 |
|
|
350 C( 3, J ) = C( 3, J ) - SUM*T3 |
|
|
351 C( 4, J ) = C( 4, J ) - SUM*T4 |
|
|
352 C( 5, J ) = C( 5, J ) - SUM*T5 |
|
|
353 C( 6, J ) = C( 6, J ) - SUM*T6 |
|
|
354 C( 7, J ) = C( 7, J ) - SUM*T7 |
|
|
355 C( 8, J ) = C( 8, J ) - SUM*T8 |
|
|
356 C( 9, J ) = C( 9, J ) - SUM*T9 |
|
|
357 C( 10, J ) = C( 10, J ) - SUM*T10 |
|
|
358 200 CONTINUE |
|
|
359 GO TO 410 |
|
|
360 ELSE |
|
|
361 * |
|
|
362 * Form C * H, where H has order n. |
|
|
363 * |
|
|
364 GO TO ( 210, 230, 250, 270, 290, 310, 330, 350, |
|
|
365 $ 370, 390 )N |
|
|
366 * |
|
|
367 * Code for general N |
|
|
368 * |
|
|
369 * w := C * v |
|
|
370 * |
|
|
371 CALL DGEMV( 'No transpose', M, N, ONE, C, LDC, V, 1, ZERO, |
|
|
372 $ WORK, 1 ) |
|
|
373 * |
|
|
374 * C := C - tau * w * v' |
|
|
375 * |
|
|
376 CALL DGER( M, N, -TAU, WORK, 1, V, 1, C, LDC ) |
|
|
377 GO TO 410 |
|
|
378 210 CONTINUE |
|
|
379 * |
|
|
380 * Special code for 1 x 1 Householder |
|
|
381 * |
|
|
382 T1 = ONE - TAU*V( 1 )*V( 1 ) |
|
|
383 DO 220 J = 1, M |
|
|
384 C( J, 1 ) = T1*C( J, 1 ) |
|
|
385 220 CONTINUE |
|
|
386 GO TO 410 |
|
|
387 230 CONTINUE |
|
|
388 * |
|
|
389 * Special code for 2 x 2 Householder |
|
|
390 * |
|
|
391 V1 = V( 1 ) |
|
|
392 T1 = TAU*V1 |
|
|
393 V2 = V( 2 ) |
|
|
394 T2 = TAU*V2 |
|
|
395 DO 240 J = 1, M |
|
|
396 SUM = V1*C( J, 1 ) + V2*C( J, 2 ) |
|
|
397 C( J, 1 ) = C( J, 1 ) - SUM*T1 |
|
|
398 C( J, 2 ) = C( J, 2 ) - SUM*T2 |
|
|
399 240 CONTINUE |
|
|
400 GO TO 410 |
|
|
401 250 CONTINUE |
|
|
402 * |
|
|
403 * Special code for 3 x 3 Householder |
|
|
404 * |
|
|
405 V1 = V( 1 ) |
|
|
406 T1 = TAU*V1 |
|
|
407 V2 = V( 2 ) |
|
|
408 T2 = TAU*V2 |
|
|
409 V3 = V( 3 ) |
|
|
410 T3 = TAU*V3 |
|
|
411 DO 260 J = 1, M |
|
|
412 SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) |
|
|
413 C( J, 1 ) = C( J, 1 ) - SUM*T1 |
|
|
414 C( J, 2 ) = C( J, 2 ) - SUM*T2 |
|
|
415 C( J, 3 ) = C( J, 3 ) - SUM*T3 |
|
|
416 260 CONTINUE |
|
|
417 GO TO 410 |
|
|
418 270 CONTINUE |
|
|
419 * |
|
|
420 * Special code for 4 x 4 Householder |
|
|
421 * |
|
|
422 V1 = V( 1 ) |
|
|
423 T1 = TAU*V1 |
|
|
424 V2 = V( 2 ) |
|
|
425 T2 = TAU*V2 |
|
|
426 V3 = V( 3 ) |
|
|
427 T3 = TAU*V3 |
|
|
428 V4 = V( 4 ) |
|
|
429 T4 = TAU*V4 |
|
|
430 DO 280 J = 1, M |
|
|
431 SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + |
|
|
432 $ V4*C( J, 4 ) |
|
|
433 C( J, 1 ) = C( J, 1 ) - SUM*T1 |
|
|
434 C( J, 2 ) = C( J, 2 ) - SUM*T2 |
|
|
435 C( J, 3 ) = C( J, 3 ) - SUM*T3 |
|
|
436 C( J, 4 ) = C( J, 4 ) - SUM*T4 |
|
|
437 280 CONTINUE |
|
|
438 GO TO 410 |
|
|
439 290 CONTINUE |
|
|
440 * |
|
|
441 * Special code for 5 x 5 Householder |
|
|
442 * |
|
|
443 V1 = V( 1 ) |
|
|
444 T1 = TAU*V1 |
|
|
445 V2 = V( 2 ) |
|
|
446 T2 = TAU*V2 |
|
|
447 V3 = V( 3 ) |
|
|
448 T3 = TAU*V3 |
|
|
449 V4 = V( 4 ) |
|
|
450 T4 = TAU*V4 |
|
|
451 V5 = V( 5 ) |
|
|
452 T5 = TAU*V5 |
|
|
453 DO 300 J = 1, M |
|
|
454 SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + |
|
|
455 $ V4*C( J, 4 ) + V5*C( J, 5 ) |
|
|
456 C( J, 1 ) = C( J, 1 ) - SUM*T1 |
|
|
457 C( J, 2 ) = C( J, 2 ) - SUM*T2 |
|
|
458 C( J, 3 ) = C( J, 3 ) - SUM*T3 |
|
|
459 C( J, 4 ) = C( J, 4 ) - SUM*T4 |
|
|
460 C( J, 5 ) = C( J, 5 ) - SUM*T5 |
|
|
461 300 CONTINUE |
|
|
462 GO TO 410 |
|
|
463 310 CONTINUE |
|
|
464 * |
|
|
465 * Special code for 6 x 6 Householder |
|
|
466 * |
|
|
467 V1 = V( 1 ) |
|
|
468 T1 = TAU*V1 |
|
|
469 V2 = V( 2 ) |
|
|
470 T2 = TAU*V2 |
|
|
471 V3 = V( 3 ) |
|
|
472 T3 = TAU*V3 |
|
|
473 V4 = V( 4 ) |
|
|
474 T4 = TAU*V4 |
|
|
475 V5 = V( 5 ) |
|
|
476 T5 = TAU*V5 |
|
|
477 V6 = V( 6 ) |
|
|
478 T6 = TAU*V6 |
|
|
479 DO 320 J = 1, M |
|
|
480 SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + |
|
|
481 $ V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) |
|
|
482 C( J, 1 ) = C( J, 1 ) - SUM*T1 |
|
|
483 C( J, 2 ) = C( J, 2 ) - SUM*T2 |
|
|
484 C( J, 3 ) = C( J, 3 ) - SUM*T3 |
|
|
485 C( J, 4 ) = C( J, 4 ) - SUM*T4 |
|
|
486 C( J, 5 ) = C( J, 5 ) - SUM*T5 |
|
|
487 C( J, 6 ) = C( J, 6 ) - SUM*T6 |
|
|
488 320 CONTINUE |
|
|
489 GO TO 410 |
|
|
490 330 CONTINUE |
|
|
491 * |
|
|
492 * Special code for 7 x 7 Householder |
|
|
493 * |
|
|
494 V1 = V( 1 ) |
|
|
495 T1 = TAU*V1 |
|
|
496 V2 = V( 2 ) |
|
|
497 T2 = TAU*V2 |
|
|
498 V3 = V( 3 ) |
|
|
499 T3 = TAU*V3 |
|
|
500 V4 = V( 4 ) |
|
|
501 T4 = TAU*V4 |
|
|
502 V5 = V( 5 ) |
|
|
503 T5 = TAU*V5 |
|
|
504 V6 = V( 6 ) |
|
|
505 T6 = TAU*V6 |
|
|
506 V7 = V( 7 ) |
|
|
507 T7 = TAU*V7 |
|
|
508 DO 340 J = 1, M |
|
|
509 SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + |
|
|
510 $ V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) + |
|
|
511 $ V7*C( J, 7 ) |
|
|
512 C( J, 1 ) = C( J, 1 ) - SUM*T1 |
|
|
513 C( J, 2 ) = C( J, 2 ) - SUM*T2 |
|
|
514 C( J, 3 ) = C( J, 3 ) - SUM*T3 |
|
|
515 C( J, 4 ) = C( J, 4 ) - SUM*T4 |
|
|
516 C( J, 5 ) = C( J, 5 ) - SUM*T5 |
|
|
517 C( J, 6 ) = C( J, 6 ) - SUM*T6 |
|
|
518 C( J, 7 ) = C( J, 7 ) - SUM*T7 |
|
|
519 340 CONTINUE |
|
|
520 GO TO 410 |
|
|
521 350 CONTINUE |
|
|
522 * |
|
|
523 * Special code for 8 x 8 Householder |
|
|
524 * |
|
|
525 V1 = V( 1 ) |
|
|
526 T1 = TAU*V1 |
|
|
527 V2 = V( 2 ) |
|
|
528 T2 = TAU*V2 |
|
|
529 V3 = V( 3 ) |
|
|
530 T3 = TAU*V3 |
|
|
531 V4 = V( 4 ) |
|
|
532 T4 = TAU*V4 |
|
|
533 V5 = V( 5 ) |
|
|
534 T5 = TAU*V5 |
|
|
535 V6 = V( 6 ) |
|
|
536 T6 = TAU*V6 |
|
|
537 V7 = V( 7 ) |
|
|
538 T7 = TAU*V7 |
|
|
539 V8 = V( 8 ) |
|
|
540 T8 = TAU*V8 |
|
|
541 DO 360 J = 1, M |
|
|
542 SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + |
|
|
543 $ V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) + |
|
|
544 $ V7*C( J, 7 ) + V8*C( J, 8 ) |
|
|
545 C( J, 1 ) = C( J, 1 ) - SUM*T1 |
|
|
546 C( J, 2 ) = C( J, 2 ) - SUM*T2 |
|
|
547 C( J, 3 ) = C( J, 3 ) - SUM*T3 |
|
|
548 C( J, 4 ) = C( J, 4 ) - SUM*T4 |
|
|
549 C( J, 5 ) = C( J, 5 ) - SUM*T5 |
|
|
550 C( J, 6 ) = C( J, 6 ) - SUM*T6 |
|
|
551 C( J, 7 ) = C( J, 7 ) - SUM*T7 |
|
|
552 C( J, 8 ) = C( J, 8 ) - SUM*T8 |
|
|
553 360 CONTINUE |
|
|
554 GO TO 410 |
|
|
555 370 CONTINUE |
|
|
556 * |
|
|
557 * Special code for 9 x 9 Householder |
|
|
558 * |
|
|
559 V1 = V( 1 ) |
|
|
560 T1 = TAU*V1 |
|
|
561 V2 = V( 2 ) |
|
|
562 T2 = TAU*V2 |
|
|
563 V3 = V( 3 ) |
|
|
564 T3 = TAU*V3 |
|
|
565 V4 = V( 4 ) |
|
|
566 T4 = TAU*V4 |
|
|
567 V5 = V( 5 ) |
|
|
568 T5 = TAU*V5 |
|
|
569 V6 = V( 6 ) |
|
|
570 T6 = TAU*V6 |
|
|
571 V7 = V( 7 ) |
|
|
572 T7 = TAU*V7 |
|
|
573 V8 = V( 8 ) |
|
|
574 T8 = TAU*V8 |
|
|
575 V9 = V( 9 ) |
|
|
576 T9 = TAU*V9 |
|
|
577 DO 380 J = 1, M |
|
|
578 SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + |
|
|
579 $ V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) + |
|
|
580 $ V7*C( J, 7 ) + V8*C( J, 8 ) + V9*C( J, 9 ) |
|
|
581 C( J, 1 ) = C( J, 1 ) - SUM*T1 |
|
|
582 C( J, 2 ) = C( J, 2 ) - SUM*T2 |
|
|
583 C( J, 3 ) = C( J, 3 ) - SUM*T3 |
|
|
584 C( J, 4 ) = C( J, 4 ) - SUM*T4 |
|
|
585 C( J, 5 ) = C( J, 5 ) - SUM*T5 |
|
|
586 C( J, 6 ) = C( J, 6 ) - SUM*T6 |
|
|
587 C( J, 7 ) = C( J, 7 ) - SUM*T7 |
|
|
588 C( J, 8 ) = C( J, 8 ) - SUM*T8 |
|
|
589 C( J, 9 ) = C( J, 9 ) - SUM*T9 |
|
|
590 380 CONTINUE |
|
|
591 GO TO 410 |
|
|
592 390 CONTINUE |
|
|
593 * |
|
|
594 * Special code for 10 x 10 Householder |
|
|
595 * |
|
|
596 V1 = V( 1 ) |
|
|
597 T1 = TAU*V1 |
|
|
598 V2 = V( 2 ) |
|
|
599 T2 = TAU*V2 |
|
|
600 V3 = V( 3 ) |
|
|
601 T3 = TAU*V3 |
|
|
602 V4 = V( 4 ) |
|
|
603 T4 = TAU*V4 |
|
|
604 V5 = V( 5 ) |
|
|
605 T5 = TAU*V5 |
|
|
606 V6 = V( 6 ) |
|
|
607 T6 = TAU*V6 |
|
|
608 V7 = V( 7 ) |
|
|
609 T7 = TAU*V7 |
|
|
610 V8 = V( 8 ) |
|
|
611 T8 = TAU*V8 |
|
|
612 V9 = V( 9 ) |
|
|
613 T9 = TAU*V9 |
|
|
614 V10 = V( 10 ) |
|
|
615 T10 = TAU*V10 |
|
|
616 DO 400 J = 1, M |
|
|
617 SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + |
|
|
618 $ V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) + |
|
|
619 $ V7*C( J, 7 ) + V8*C( J, 8 ) + V9*C( J, 9 ) + |
|
|
620 $ V10*C( J, 10 ) |
|
|
621 C( J, 1 ) = C( J, 1 ) - SUM*T1 |
|
|
622 C( J, 2 ) = C( J, 2 ) - SUM*T2 |
|
|
623 C( J, 3 ) = C( J, 3 ) - SUM*T3 |
|
|
624 C( J, 4 ) = C( J, 4 ) - SUM*T4 |
|
|
625 C( J, 5 ) = C( J, 5 ) - SUM*T5 |
|
|
626 C( J, 6 ) = C( J, 6 ) - SUM*T6 |
|
|
627 C( J, 7 ) = C( J, 7 ) - SUM*T7 |
|
|
628 C( J, 8 ) = C( J, 8 ) - SUM*T8 |
|
|
629 C( J, 9 ) = C( J, 9 ) - SUM*T9 |
|
|
630 C( J, 10 ) = C( J, 10 ) - SUM*T10 |
|
|
631 400 CONTINUE |
|
|
632 GO TO 410 |
|
|
633 END IF |
|
|
634 410 CONTINUE |
|
|
635 RETURN |
|
|
636 * |
|
|
637 * End of DLARFX |
|
|
638 * |
|
|
639 END |
