Mercurial > octave-nkf
diff libcruft/lapack/clarfx.f @ 7789:82be108cc558
First attempt at single precision tyeps
* * *
corrections to qrupdate single precision routines
* * *
prefer demotion to single over promotion to double
* * *
Add single precision support to log2 function
* * *
Trivial PROJECT file update
* * *
Cache optimized hermitian/transpose methods
* * *
Add tests for tranpose/hermitian and ChangeLog entry for new transpose code
| author | David Bateman <dbateman@free.fr> |
|---|---|
| date | Sun, 27 Apr 2008 22:34:17 +0200 |
| parents | |
| children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/libcruft/lapack/clarfx.f Sun Apr 27 22:34:17 2008 +0200 @@ -0,0 +1,640 @@ + SUBROUTINE CLARFX( SIDE, M, N, V, TAU, C, LDC, WORK ) +* +* -- LAPACK auxiliary routine (version 3.1) -- +* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. +* November 2006 +* +* .. Scalar Arguments .. + CHARACTER SIDE + INTEGER LDC, M, N + COMPLEX TAU +* .. +* .. Array Arguments .. + COMPLEX C( LDC, * ), V( * ), WORK( * ) +* .. +* +* Purpose +* ======= +* +* CLARFX applies a complex elementary reflector H to a complex m by n +* matrix C, from either the left or the right. H is represented in the +* form +* +* H = I - tau * v * v' +* +* where tau is a complex scalar and v is a complex vector. +* +* If tau = 0, then H is taken to be the unit matrix +* +* This version uses inline code if H has order < 11. +* +* Arguments +* ========= +* +* SIDE (input) CHARACTER*1 +* = 'L': form H * C +* = 'R': form C * H +* +* M (input) INTEGER +* The number of rows of the matrix C. +* +* N (input) INTEGER +* The number of columns of the matrix C. +* +* V (input) COMPLEX array, dimension (M) if SIDE = 'L' +* or (N) if SIDE = 'R' +* The vector v in the representation of H. +* +* TAU (input) COMPLEX +* The value tau in the representation of H. +* +* C (input/output) COMPLEX array, dimension (LDC,N) +* On entry, the m by n matrix C. +* On exit, C is overwritten by the matrix H * C if SIDE = 'L', +* or C * H if SIDE = 'R'. +* +* LDC (input) INTEGER +* The leading dimension of the array C. LDA >= max(1,M). +* +* WORK (workspace) COMPLEX array, dimension (N) if SIDE = 'L' +* or (M) if SIDE = 'R' +* WORK is not referenced if H has order < 11. +* +* ===================================================================== +* +* .. Parameters .. + COMPLEX ZERO, ONE + PARAMETER ( ZERO = ( 0.0E+0, 0.0E+0 ), + $ ONE = ( 1.0E+0, 0.0E+0 ) ) +* .. +* .. Local Scalars .. + INTEGER J + COMPLEX SUM, T1, T10, T2, T3, T4, T5, T6, T7, T8, T9, + $ V1, V10, V2, V3, V4, V5, V6, V7, V8, V9 +* .. +* .. External Functions .. + LOGICAL LSAME + EXTERNAL LSAME +* .. +* .. External Subroutines .. + EXTERNAL CGEMV, CGERC +* .. +* .. Intrinsic Functions .. + INTRINSIC CONJG +* .. +* .. Executable Statements .. +* + IF( TAU.EQ.ZERO ) + $ RETURN + IF( LSAME( SIDE, 'L' ) ) THEN +* +* Form H * C, where H has order m. +* + GO TO ( 10, 30, 50, 70, 90, 110, 130, 150, + $ 170, 190 )M +* +* Code for general M +* +* w := C'*v +* + CALL CGEMV( 'Conjugate transpose', M, N, ONE, C, LDC, V, 1, + $ ZERO, WORK, 1 ) +* +* C := C - tau * v * w' +* + CALL CGERC( M, N, -TAU, V, 1, WORK, 1, C, LDC ) + GO TO 410 + 10 CONTINUE +* +* Special code for 1 x 1 Householder +* + T1 = ONE - TAU*V( 1 )*CONJG( V( 1 ) ) + DO 20 J = 1, N + C( 1, J ) = T1*C( 1, J ) + 20 CONTINUE + GO TO 410 + 30 CONTINUE +* +* Special code for 2 x 2 Householder +* + V1 = CONJG( V( 1 ) ) + T1 = TAU*CONJG( V1 ) + V2 = CONJG( V( 2 ) ) + T2 = TAU*CONJG( V2 ) + DO 40 J = 1, N + SUM = V1*C( 1, J ) + V2*C( 2, J ) + C( 1, J ) = C( 1, J ) - SUM*T1 + C( 2, J ) = C( 2, J ) - SUM*T2 + 40 CONTINUE + GO TO 410 + 50 CONTINUE +* +* Special code for 3 x 3 Householder +* + V1 = CONJG( V( 1 ) ) + T1 = TAU*CONJG( V1 ) + V2 = CONJG( V( 2 ) ) + T2 = TAU*CONJG( V2 ) + V3 = CONJG( V( 3 ) ) + T3 = TAU*CONJG( V3 ) + DO 60 J = 1, N + SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + C( 1, J ) = C( 1, J ) - SUM*T1 + C( 2, J ) = C( 2, J ) - SUM*T2 + C( 3, J ) = C( 3, J ) - SUM*T3 + 60 CONTINUE + GO TO 410 + 70 CONTINUE +* +* Special code for 4 x 4 Householder +* + V1 = CONJG( V( 1 ) ) + T1 = TAU*CONJG( V1 ) + V2 = CONJG( V( 2 ) ) + T2 = TAU*CONJG( V2 ) + V3 = CONJG( V( 3 ) ) + T3 = TAU*CONJG( V3 ) + V4 = CONJG( V( 4 ) ) + T4 = TAU*CONJG( V4 ) + DO 80 J = 1, N + SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + + $ V4*C( 4, J ) + C( 1, J ) = C( 1, J ) - SUM*T1 + C( 2, J ) = C( 2, J ) - SUM*T2 + C( 3, J ) = C( 3, J ) - SUM*T3 + C( 4, J ) = C( 4, J ) - SUM*T4 + 80 CONTINUE + GO TO 410 + 90 CONTINUE +* +* Special code for 5 x 5 Householder +* + V1 = CONJG( V( 1 ) ) + T1 = TAU*CONJG( V1 ) + V2 = CONJG( V( 2 ) ) + T2 = TAU*CONJG( V2 ) + V3 = CONJG( V( 3 ) ) + T3 = TAU*CONJG( V3 ) + V4 = CONJG( V( 4 ) ) + T4 = TAU*CONJG( V4 ) + V5 = CONJG( V( 5 ) ) + T5 = TAU*CONJG( V5 ) + DO 100 J = 1, N + SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + + $ V4*C( 4, J ) + V5*C( 5, J ) + C( 1, J ) = C( 1, J ) - SUM*T1 + C( 2, J ) = C( 2, J ) - SUM*T2 + C( 3, J ) = C( 3, J ) - SUM*T3 + C( 4, J ) = C( 4, J ) - SUM*T4 + C( 5, J ) = C( 5, J ) - SUM*T5 + 100 CONTINUE + GO TO 410 + 110 CONTINUE +* +* Special code for 6 x 6 Householder +* + V1 = CONJG( V( 1 ) ) + T1 = TAU*CONJG( V1 ) + V2 = CONJG( V( 2 ) ) + T2 = TAU*CONJG( V2 ) + V3 = CONJG( V( 3 ) ) + T3 = TAU*CONJG( V3 ) + V4 = CONJG( V( 4 ) ) + T4 = TAU*CONJG( V4 ) + V5 = CONJG( V( 5 ) ) + T5 = TAU*CONJG( V5 ) + V6 = CONJG( V( 6 ) ) + T6 = TAU*CONJG( V6 ) + DO 120 J = 1, N + SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + + $ V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) + C( 1, J ) = C( 1, J ) - SUM*T1 + C( 2, J ) = C( 2, J ) - SUM*T2 + C( 3, J ) = C( 3, J ) - SUM*T3 + C( 4, J ) = C( 4, J ) - SUM*T4 + C( 5, J ) = C( 5, J ) - SUM*T5 + C( 6, J ) = C( 6, J ) - SUM*T6 + 120 CONTINUE + GO TO 410 + 130 CONTINUE +* +* Special code for 7 x 7 Householder +* + V1 = CONJG( V( 1 ) ) + T1 = TAU*CONJG( V1 ) + V2 = CONJG( V( 2 ) ) + T2 = TAU*CONJG( V2 ) + V3 = CONJG( V( 3 ) ) + T3 = TAU*CONJG( V3 ) + V4 = CONJG( V( 4 ) ) + T4 = TAU*CONJG( V4 ) + V5 = CONJG( V( 5 ) ) + T5 = TAU*CONJG( V5 ) + V6 = CONJG( V( 6 ) ) + T6 = TAU*CONJG( V6 ) + V7 = CONJG( V( 7 ) ) + T7 = TAU*CONJG( V7 ) + DO 140 J = 1, N + SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + + $ V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) + + $ V7*C( 7, J ) + C( 1, J ) = C( 1, J ) - SUM*T1 + C( 2, J ) = C( 2, J ) - SUM*T2 + C( 3, J ) = C( 3, J ) - SUM*T3 + C( 4, J ) = C( 4, J ) - SUM*T4 + C( 5, J ) = C( 5, J ) - SUM*T5 + C( 6, J ) = C( 6, J ) - SUM*T6 + C( 7, J ) = C( 7, J ) - SUM*T7 + 140 CONTINUE + GO TO 410 + 150 CONTINUE +* +* Special code for 8 x 8 Householder +* + V1 = CONJG( V( 1 ) ) + T1 = TAU*CONJG( V1 ) + V2 = CONJG( V( 2 ) ) + T2 = TAU*CONJG( V2 ) + V3 = CONJG( V( 3 ) ) + T3 = TAU*CONJG( V3 ) + V4 = CONJG( V( 4 ) ) + T4 = TAU*CONJG( V4 ) + V5 = CONJG( V( 5 ) ) + T5 = TAU*CONJG( V5 ) + V6 = CONJG( V( 6 ) ) + T6 = TAU*CONJG( V6 ) + V7 = CONJG( V( 7 ) ) + T7 = TAU*CONJG( V7 ) + V8 = CONJG( V( 8 ) ) + T8 = TAU*CONJG( V8 ) + DO 160 J = 1, N + SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + + $ V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) + + $ V7*C( 7, J ) + V8*C( 8, J ) + C( 1, J ) = C( 1, J ) - SUM*T1 + C( 2, J ) = C( 2, J ) - SUM*T2 + C( 3, J ) = C( 3, J ) - SUM*T3 + C( 4, J ) = C( 4, J ) - SUM*T4 + C( 5, J ) = C( 5, J ) - SUM*T5 + C( 6, J ) = C( 6, J ) - SUM*T6 + C( 7, J ) = C( 7, J ) - SUM*T7 + C( 8, J ) = C( 8, J ) - SUM*T8 + 160 CONTINUE + GO TO 410 + 170 CONTINUE +* +* Special code for 9 x 9 Householder +* + V1 = CONJG( V( 1 ) ) + T1 = TAU*CONJG( V1 ) + V2 = CONJG( V( 2 ) ) + T2 = TAU*CONJG( V2 ) + V3 = CONJG( V( 3 ) ) + T3 = TAU*CONJG( V3 ) + V4 = CONJG( V( 4 ) ) + T4 = TAU*CONJG( V4 ) + V5 = CONJG( V( 5 ) ) + T5 = TAU*CONJG( V5 ) + V6 = CONJG( V( 6 ) ) + T6 = TAU*CONJG( V6 ) + V7 = CONJG( V( 7 ) ) + T7 = TAU*CONJG( V7 ) + V8 = CONJG( V( 8 ) ) + T8 = TAU*CONJG( V8 ) + V9 = CONJG( V( 9 ) ) + T9 = TAU*CONJG( V9 ) + DO 180 J = 1, N + SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + + $ V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) + + $ V7*C( 7, J ) + V8*C( 8, J ) + V9*C( 9, J ) + C( 1, J ) = C( 1, J ) - SUM*T1 + C( 2, J ) = C( 2, J ) - SUM*T2 + C( 3, J ) = C( 3, J ) - SUM*T3 + C( 4, J ) = C( 4, J ) - SUM*T4 + C( 5, J ) = C( 5, J ) - SUM*T5 + C( 6, J ) = C( 6, J ) - SUM*T6 + C( 7, J ) = C( 7, J ) - SUM*T7 + C( 8, J ) = C( 8, J ) - SUM*T8 + C( 9, J ) = C( 9, J ) - SUM*T9 + 180 CONTINUE + GO TO 410 + 190 CONTINUE +* +* Special code for 10 x 10 Householder +* + V1 = CONJG( V( 1 ) ) + T1 = TAU*CONJG( V1 ) + V2 = CONJG( V( 2 ) ) + T2 = TAU*CONJG( V2 ) + V3 = CONJG( V( 3 ) ) + T3 = TAU*CONJG( V3 ) + V4 = CONJG( V( 4 ) ) + T4 = TAU*CONJG( V4 ) + V5 = CONJG( V( 5 ) ) + T5 = TAU*CONJG( V5 ) + V6 = CONJG( V( 6 ) ) + T6 = TAU*CONJG( V6 ) + V7 = CONJG( V( 7 ) ) + T7 = TAU*CONJG( V7 ) + V8 = CONJG( V( 8 ) ) + T8 = TAU*CONJG( V8 ) + V9 = CONJG( V( 9 ) ) + T9 = TAU*CONJG( V9 ) + V10 = CONJG( V( 10 ) ) + T10 = TAU*CONJG( V10 ) + DO 200 J = 1, N + SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) + + $ V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) + + $ V7*C( 7, J ) + V8*C( 8, J ) + V9*C( 9, J ) + + $ V10*C( 10, J ) + C( 1, J ) = C( 1, J ) - SUM*T1 + C( 2, J ) = C( 2, J ) - SUM*T2 + C( 3, J ) = C( 3, J ) - SUM*T3 + C( 4, J ) = C( 4, J ) - SUM*T4 + C( 5, J ) = C( 5, J ) - SUM*T5 + C( 6, J ) = C( 6, J ) - SUM*T6 + C( 7, J ) = C( 7, J ) - SUM*T7 + C( 8, J ) = C( 8, J ) - SUM*T8 + C( 9, J ) = C( 9, J ) - SUM*T9 + C( 10, J ) = C( 10, J ) - SUM*T10 + 200 CONTINUE + GO TO 410 + ELSE +* +* Form C * H, where H has order n. +* + GO TO ( 210, 230, 250, 270, 290, 310, 330, 350, + $ 370, 390 )N +* +* Code for general N +* +* w := C * v +* + CALL CGEMV( 'No transpose', M, N, ONE, C, LDC, V, 1, ZERO, + $ WORK, 1 ) +* +* C := C - tau * w * v' +* + CALL CGERC( M, N, -TAU, WORK, 1, V, 1, C, LDC ) + GO TO 410 + 210 CONTINUE +* +* Special code for 1 x 1 Householder +* + T1 = ONE - TAU*V( 1 )*CONJG( V( 1 ) ) + DO 220 J = 1, M + C( J, 1 ) = T1*C( J, 1 ) + 220 CONTINUE + GO TO 410 + 230 CONTINUE +* +* Special code for 2 x 2 Householder +* + V1 = V( 1 ) + T1 = TAU*CONJG( V1 ) + V2 = V( 2 ) + T2 = TAU*CONJG( V2 ) + DO 240 J = 1, M + SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + C( J, 1 ) = C( J, 1 ) - SUM*T1 + C( J, 2 ) = C( J, 2 ) - SUM*T2 + 240 CONTINUE + GO TO 410 + 250 CONTINUE +* +* Special code for 3 x 3 Householder +* + V1 = V( 1 ) + T1 = TAU*CONJG( V1 ) + V2 = V( 2 ) + T2 = TAU*CONJG( V2 ) + V3 = V( 3 ) + T3 = TAU*CONJG( V3 ) + DO 260 J = 1, M + SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + C( J, 1 ) = C( J, 1 ) - SUM*T1 + C( J, 2 ) = C( J, 2 ) - SUM*T2 + C( J, 3 ) = C( J, 3 ) - SUM*T3 + 260 CONTINUE + GO TO 410 + 270 CONTINUE +* +* Special code for 4 x 4 Householder +* + V1 = V( 1 ) + T1 = TAU*CONJG( V1 ) + V2 = V( 2 ) + T2 = TAU*CONJG( V2 ) + V3 = V( 3 ) + T3 = TAU*CONJG( V3 ) + V4 = V( 4 ) + T4 = TAU*CONJG( V4 ) + DO 280 J = 1, M + SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + + $ V4*C( J, 4 ) + C( J, 1 ) = C( J, 1 ) - SUM*T1 + C( J, 2 ) = C( J, 2 ) - SUM*T2 + C( J, 3 ) = C( J, 3 ) - SUM*T3 + C( J, 4 ) = C( J, 4 ) - SUM*T4 + 280 CONTINUE + GO TO 410 + 290 CONTINUE +* +* Special code for 5 x 5 Householder +* + V1 = V( 1 ) + T1 = TAU*CONJG( V1 ) + V2 = V( 2 ) + T2 = TAU*CONJG( V2 ) + V3 = V( 3 ) + T3 = TAU*CONJG( V3 ) + V4 = V( 4 ) + T4 = TAU*CONJG( V4 ) + V5 = V( 5 ) + T5 = TAU*CONJG( V5 ) + DO 300 J = 1, M + SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + + $ V4*C( J, 4 ) + V5*C( J, 5 ) + C( J, 1 ) = C( J, 1 ) - SUM*T1 + C( J, 2 ) = C( J, 2 ) - SUM*T2 + C( J, 3 ) = C( J, 3 ) - SUM*T3 + C( J, 4 ) = C( J, 4 ) - SUM*T4 + C( J, 5 ) = C( J, 5 ) - SUM*T5 + 300 CONTINUE + GO TO 410 + 310 CONTINUE +* +* Special code for 6 x 6 Householder +* + V1 = V( 1 ) + T1 = TAU*CONJG( V1 ) + V2 = V( 2 ) + T2 = TAU*CONJG( V2 ) + V3 = V( 3 ) + T3 = TAU*CONJG( V3 ) + V4 = V( 4 ) + T4 = TAU*CONJG( V4 ) + V5 = V( 5 ) + T5 = TAU*CONJG( V5 ) + V6 = V( 6 ) + T6 = TAU*CONJG( V6 ) + DO 320 J = 1, M + SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + + $ V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) + C( J, 1 ) = C( J, 1 ) - SUM*T1 + C( J, 2 ) = C( J, 2 ) - SUM*T2 + C( J, 3 ) = C( J, 3 ) - SUM*T3 + C( J, 4 ) = C( J, 4 ) - SUM*T4 + C( J, 5 ) = C( J, 5 ) - SUM*T5 + C( J, 6 ) = C( J, 6 ) - SUM*T6 + 320 CONTINUE + GO TO 410 + 330 CONTINUE +* +* Special code for 7 x 7 Householder +* + V1 = V( 1 ) + T1 = TAU*CONJG( V1 ) + V2 = V( 2 ) + T2 = TAU*CONJG( V2 ) + V3 = V( 3 ) + T3 = TAU*CONJG( V3 ) + V4 = V( 4 ) + T4 = TAU*CONJG( V4 ) + V5 = V( 5 ) + T5 = TAU*CONJG( V5 ) + V6 = V( 6 ) + T6 = TAU*CONJG( V6 ) + V7 = V( 7 ) + T7 = TAU*CONJG( V7 ) + DO 340 J = 1, M + SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + + $ V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) + + $ V7*C( J, 7 ) + C( J, 1 ) = C( J, 1 ) - SUM*T1 + C( J, 2 ) = C( J, 2 ) - SUM*T2 + C( J, 3 ) = C( J, 3 ) - SUM*T3 + C( J, 4 ) = C( J, 4 ) - SUM*T4 + C( J, 5 ) = C( J, 5 ) - SUM*T5 + C( J, 6 ) = C( J, 6 ) - SUM*T6 + C( J, 7 ) = C( J, 7 ) - SUM*T7 + 340 CONTINUE + GO TO 410 + 350 CONTINUE +* +* Special code for 8 x 8 Householder +* + V1 = V( 1 ) + T1 = TAU*CONJG( V1 ) + V2 = V( 2 ) + T2 = TAU*CONJG( V2 ) + V3 = V( 3 ) + T3 = TAU*CONJG( V3 ) + V4 = V( 4 ) + T4 = TAU*CONJG( V4 ) + V5 = V( 5 ) + T5 = TAU*CONJG( V5 ) + V6 = V( 6 ) + T6 = TAU*CONJG( V6 ) + V7 = V( 7 ) + T7 = TAU*CONJG( V7 ) + V8 = V( 8 ) + T8 = TAU*CONJG( V8 ) + DO 360 J = 1, M + SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + + $ V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) + + $ V7*C( J, 7 ) + V8*C( J, 8 ) + C( J, 1 ) = C( J, 1 ) - SUM*T1 + C( J, 2 ) = C( J, 2 ) - SUM*T2 + C( J, 3 ) = C( J, 3 ) - SUM*T3 + C( J, 4 ) = C( J, 4 ) - SUM*T4 + C( J, 5 ) = C( J, 5 ) - SUM*T5 + C( J, 6 ) = C( J, 6 ) - SUM*T6 + C( J, 7 ) = C( J, 7 ) - SUM*T7 + C( J, 8 ) = C( J, 8 ) - SUM*T8 + 360 CONTINUE + GO TO 410 + 370 CONTINUE +* +* Special code for 9 x 9 Householder +* + V1 = V( 1 ) + T1 = TAU*CONJG( V1 ) + V2 = V( 2 ) + T2 = TAU*CONJG( V2 ) + V3 = V( 3 ) + T3 = TAU*CONJG( V3 ) + V4 = V( 4 ) + T4 = TAU*CONJG( V4 ) + V5 = V( 5 ) + T5 = TAU*CONJG( V5 ) + V6 = V( 6 ) + T6 = TAU*CONJG( V6 ) + V7 = V( 7 ) + T7 = TAU*CONJG( V7 ) + V8 = V( 8 ) + T8 = TAU*CONJG( V8 ) + V9 = V( 9 ) + T9 = TAU*CONJG( V9 ) + DO 380 J = 1, M + SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + + $ V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) + + $ V7*C( J, 7 ) + V8*C( J, 8 ) + V9*C( J, 9 ) + C( J, 1 ) = C( J, 1 ) - SUM*T1 + C( J, 2 ) = C( J, 2 ) - SUM*T2 + C( J, 3 ) = C( J, 3 ) - SUM*T3 + C( J, 4 ) = C( J, 4 ) - SUM*T4 + C( J, 5 ) = C( J, 5 ) - SUM*T5 + C( J, 6 ) = C( J, 6 ) - SUM*T6 + C( J, 7 ) = C( J, 7 ) - SUM*T7 + C( J, 8 ) = C( J, 8 ) - SUM*T8 + C( J, 9 ) = C( J, 9 ) - SUM*T9 + 380 CONTINUE + GO TO 410 + 390 CONTINUE +* +* Special code for 10 x 10 Householder +* + V1 = V( 1 ) + T1 = TAU*CONJG( V1 ) + V2 = V( 2 ) + T2 = TAU*CONJG( V2 ) + V3 = V( 3 ) + T3 = TAU*CONJG( V3 ) + V4 = V( 4 ) + T4 = TAU*CONJG( V4 ) + V5 = V( 5 ) + T5 = TAU*CONJG( V5 ) + V6 = V( 6 ) + T6 = TAU*CONJG( V6 ) + V7 = V( 7 ) + T7 = TAU*CONJG( V7 ) + V8 = V( 8 ) + T8 = TAU*CONJG( V8 ) + V9 = V( 9 ) + T9 = TAU*CONJG( V9 ) + V10 = V( 10 ) + T10 = TAU*CONJG( V10 ) + DO 400 J = 1, M + SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) + + $ V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) + + $ V7*C( J, 7 ) + V8*C( J, 8 ) + V9*C( J, 9 ) + + $ V10*C( J, 10 ) + C( J, 1 ) = C( J, 1 ) - SUM*T1 + C( J, 2 ) = C( J, 2 ) - SUM*T2 + C( J, 3 ) = C( J, 3 ) - SUM*T3 + C( J, 4 ) = C( J, 4 ) - SUM*T4 + C( J, 5 ) = C( J, 5 ) - SUM*T5 + C( J, 6 ) = C( J, 6 ) - SUM*T6 + C( J, 7 ) = C( J, 7 ) - SUM*T7 + C( J, 8 ) = C( J, 8 ) - SUM*T8 + C( J, 9 ) = C( J, 9 ) - SUM*T9 + C( J, 10 ) = C( J, 10 ) - SUM*T10 + 400 CONTINUE + GO TO 410 + END IF + 410 RETURN +* +* End of CLARFX +* + END