KTH framework for Nek5000 toolboxes; testing version  0.0.1
dlarfx.f
Go to the documentation of this file.
1  SUBROUTINE dlarfx( SIDE, M, N, V, TAU, C, LDC, WORK )
2 *
3 * -- LAPACK auxiliary routine (version 3.0) --
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
dgemv
subroutine dgemv(TRANS, M, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
Definition: dgemv.f:3
dger
subroutine dger(M, N, ALPHA, X, INCX, Y, INCY, A, LDA)
Definition: dger.f:2
dlarfx
subroutine dlarfx(SIDE, M, N, V, TAU, C, LDC, WORK)
Definition: dlarfx.f:2