LAPACK 3.3.0

clarfx.f

Go to the documentation of this file.
00001       SUBROUTINE CLARFX( SIDE, M, N, V, TAU, C, LDC, WORK )
00002       IMPLICIT NONE
00003 *
00004 *  -- LAPACK auxiliary routine (version 3.2) --
00005 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
00006 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
00007 *     November 2006
00008 *
00009 *     .. Scalar Arguments ..
00010       CHARACTER          SIDE
00011       INTEGER            LDC, M, N
00012       COMPLEX            TAU
00013 *     ..
00014 *     .. Array Arguments ..
00015       COMPLEX            C( LDC, * ), V( * ), WORK( * )
00016 *     ..
00017 *
00018 *  Purpose
00019 *  =======
00020 *
00021 *  CLARFX applies a complex elementary reflector H to a complex m by n
00022 *  matrix C, from either the left or the right. H is represented in the
00023 *  form
00024 *
00025 *        H = I - tau * v * v'
00026 *
00027 *  where tau is a complex scalar and v is a complex vector.
00028 *
00029 *  If tau = 0, then H is taken to be the unit matrix
00030 *
00031 *  This version uses inline code if H has order < 11.
00032 *
00033 *  Arguments
00034 *  =========
00035 *
00036 *  SIDE    (input) CHARACTER*1
00037 *          = 'L': form  H * C
00038 *          = 'R': form  C * H
00039 *
00040 *  M       (input) INTEGER
00041 *          The number of rows of the matrix C.
00042 *
00043 *  N       (input) INTEGER
00044 *          The number of columns of the matrix C.
00045 *
00046 *  V       (input) COMPLEX array, dimension (M) if SIDE = 'L'
00047 *                                        or (N) if SIDE = 'R'
00048 *          The vector v in the representation of H.
00049 *
00050 *  TAU     (input) COMPLEX
00051 *          The value tau in the representation of H.
00052 *
00053 *  C       (input/output) COMPLEX array, dimension (LDC,N)
00054 *          On entry, the m by n matrix C.
00055 *          On exit, C is overwritten by the matrix H * C if SIDE = 'L',
00056 *          or C * H if SIDE = 'R'.
00057 *
00058 *  LDC     (input) INTEGER
00059 *          The leading dimension of the array C. LDA >= max(1,M).
00060 *
00061 *  WORK    (workspace) COMPLEX array, dimension (N) if SIDE = 'L'
00062 *                                            or (M) if SIDE = 'R'
00063 *          WORK is not referenced if H has order < 11.
00064 *
00065 *  =====================================================================
00066 *
00067 *     .. Parameters ..
00068       COMPLEX            ZERO, ONE
00069       PARAMETER          ( ZERO = ( 0.0E+0, 0.0E+0 ),
00070      $                   ONE = ( 1.0E+0, 0.0E+0 ) )
00071 *     ..
00072 *     .. Local Scalars ..
00073       INTEGER            J
00074       COMPLEX            SUM, T1, T10, T2, T3, T4, T5, T6, T7, T8, T9,
00075      $                   V1, V10, V2, V3, V4, V5, V6, V7, V8, V9
00076 *     ..
00077 *     .. External Functions ..
00078       LOGICAL            LSAME
00079       EXTERNAL           LSAME
00080 *     ..
00081 *     .. External Subroutines ..
00082       EXTERNAL           CLARF
00083 *     ..
00084 *     .. Intrinsic Functions ..
00085       INTRINSIC          CONJG
00086 *     ..
00087 *     .. Executable Statements ..
00088 *
00089       IF( TAU.EQ.ZERO )
00090      $   RETURN
00091       IF( LSAME( SIDE, 'L' ) ) THEN
00092 *
00093 *        Form  H * C, where H has order m.
00094 *
00095          GO TO ( 10, 30, 50, 70, 90, 110, 130, 150,
00096      $           170, 190 )M
00097 *
00098 *        Code for general M
00099 *
00100          CALL CLARF( SIDE, M, N, V, 1, TAU, C, LDC, WORK )
00101          GO TO 410
00102    10    CONTINUE
00103 *
00104 *        Special code for 1 x 1 Householder
00105 *
00106          T1 = ONE - TAU*V( 1 )*CONJG( V( 1 ) )
00107          DO 20 J = 1, N
00108             C( 1, J ) = T1*C( 1, J )
00109    20    CONTINUE
00110          GO TO 410
00111    30    CONTINUE
00112 *
00113 *        Special code for 2 x 2 Householder
00114 *
00115          V1 = CONJG( V( 1 ) )
00116          T1 = TAU*CONJG( V1 )
00117          V2 = CONJG( V( 2 ) )
00118          T2 = TAU*CONJG( V2 )
00119          DO 40 J = 1, N
00120             SUM = V1*C( 1, J ) + V2*C( 2, J )
00121             C( 1, J ) = C( 1, J ) - SUM*T1
00122             C( 2, J ) = C( 2, J ) - SUM*T2
00123    40    CONTINUE
00124          GO TO 410
00125    50    CONTINUE
00126 *
00127 *        Special code for 3 x 3 Householder
00128 *
00129          V1 = CONJG( V( 1 ) )
00130          T1 = TAU*CONJG( V1 )
00131          V2 = CONJG( V( 2 ) )
00132          T2 = TAU*CONJG( V2 )
00133          V3 = CONJG( V( 3 ) )
00134          T3 = TAU*CONJG( V3 )
00135          DO 60 J = 1, N
00136             SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J )
00137             C( 1, J ) = C( 1, J ) - SUM*T1
00138             C( 2, J ) = C( 2, J ) - SUM*T2
00139             C( 3, J ) = C( 3, J ) - SUM*T3
00140    60    CONTINUE
00141          GO TO 410
00142    70    CONTINUE
00143 *
00144 *        Special code for 4 x 4 Householder
00145 *
00146          V1 = CONJG( V( 1 ) )
00147          T1 = TAU*CONJG( V1 )
00148          V2 = CONJG( V( 2 ) )
00149          T2 = TAU*CONJG( V2 )
00150          V3 = CONJG( V( 3 ) )
00151          T3 = TAU*CONJG( V3 )
00152          V4 = CONJG( V( 4 ) )
00153          T4 = TAU*CONJG( V4 )
00154          DO 80 J = 1, N
00155             SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) +
00156      $            V4*C( 4, J )
00157             C( 1, J ) = C( 1, J ) - SUM*T1
00158             C( 2, J ) = C( 2, J ) - SUM*T2
00159             C( 3, J ) = C( 3, J ) - SUM*T3
00160             C( 4, J ) = C( 4, J ) - SUM*T4
00161    80    CONTINUE
00162          GO TO 410
00163    90    CONTINUE
00164 *
00165 *        Special code for 5 x 5 Householder
00166 *
00167          V1 = CONJG( V( 1 ) )
00168          T1 = TAU*CONJG( V1 )
00169          V2 = CONJG( V( 2 ) )
00170          T2 = TAU*CONJG( V2 )
00171          V3 = CONJG( V( 3 ) )
00172          T3 = TAU*CONJG( V3 )
00173          V4 = CONJG( V( 4 ) )
00174          T4 = TAU*CONJG( V4 )
00175          V5 = CONJG( V( 5 ) )
00176          T5 = TAU*CONJG( V5 )
00177          DO 100 J = 1, N
00178             SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) +
00179      $            V4*C( 4, J ) + V5*C( 5, J )
00180             C( 1, J ) = C( 1, J ) - SUM*T1
00181             C( 2, J ) = C( 2, J ) - SUM*T2
00182             C( 3, J ) = C( 3, J ) - SUM*T3
00183             C( 4, J ) = C( 4, J ) - SUM*T4
00184             C( 5, J ) = C( 5, J ) - SUM*T5
00185   100    CONTINUE
00186          GO TO 410
00187   110    CONTINUE
00188 *
00189 *        Special code for 6 x 6 Householder
00190 *
00191          V1 = CONJG( V( 1 ) )
00192          T1 = TAU*CONJG( V1 )
00193          V2 = CONJG( V( 2 ) )
00194          T2 = TAU*CONJG( V2 )
00195          V3 = CONJG( V( 3 ) )
00196          T3 = TAU*CONJG( V3 )
00197          V4 = CONJG( V( 4 ) )
00198          T4 = TAU*CONJG( V4 )
00199          V5 = CONJG( V( 5 ) )
00200          T5 = TAU*CONJG( V5 )
00201          V6 = CONJG( V( 6 ) )
00202          T6 = TAU*CONJG( V6 )
00203          DO 120 J = 1, N
00204             SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) +
00205      $            V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J )
00206             C( 1, J ) = C( 1, J ) - SUM*T1
00207             C( 2, J ) = C( 2, J ) - SUM*T2
00208             C( 3, J ) = C( 3, J ) - SUM*T3
00209             C( 4, J ) = C( 4, J ) - SUM*T4
00210             C( 5, J ) = C( 5, J ) - SUM*T5
00211             C( 6, J ) = C( 6, J ) - SUM*T6
00212   120    CONTINUE
00213          GO TO 410
00214   130    CONTINUE
00215 *
00216 *        Special code for 7 x 7 Householder
00217 *
00218          V1 = CONJG( V( 1 ) )
00219          T1 = TAU*CONJG( V1 )
00220          V2 = CONJG( V( 2 ) )
00221          T2 = TAU*CONJG( V2 )
00222          V3 = CONJG( V( 3 ) )
00223          T3 = TAU*CONJG( V3 )
00224          V4 = CONJG( V( 4 ) )
00225          T4 = TAU*CONJG( V4 )
00226          V5 = CONJG( V( 5 ) )
00227          T5 = TAU*CONJG( V5 )
00228          V6 = CONJG( V( 6 ) )
00229          T6 = TAU*CONJG( V6 )
00230          V7 = CONJG( V( 7 ) )
00231          T7 = TAU*CONJG( V7 )
00232          DO 140 J = 1, N
00233             SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) +
00234      $            V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) +
00235      $            V7*C( 7, J )
00236             C( 1, J ) = C( 1, J ) - SUM*T1
00237             C( 2, J ) = C( 2, J ) - SUM*T2
00238             C( 3, J ) = C( 3, J ) - SUM*T3
00239             C( 4, J ) = C( 4, J ) - SUM*T4
00240             C( 5, J ) = C( 5, J ) - SUM*T5
00241             C( 6, J ) = C( 6, J ) - SUM*T6
00242             C( 7, J ) = C( 7, J ) - SUM*T7
00243   140    CONTINUE
00244          GO TO 410
00245   150    CONTINUE
00246 *
00247 *        Special code for 8 x 8 Householder
00248 *
00249          V1 = CONJG( V( 1 ) )
00250          T1 = TAU*CONJG( V1 )
00251          V2 = CONJG( V( 2 ) )
00252          T2 = TAU*CONJG( V2 )
00253          V3 = CONJG( V( 3 ) )
00254          T3 = TAU*CONJG( V3 )
00255          V4 = CONJG( V( 4 ) )
00256          T4 = TAU*CONJG( V4 )
00257          V5 = CONJG( V( 5 ) )
00258          T5 = TAU*CONJG( V5 )
00259          V6 = CONJG( V( 6 ) )
00260          T6 = TAU*CONJG( V6 )
00261          V7 = CONJG( V( 7 ) )
00262          T7 = TAU*CONJG( V7 )
00263          V8 = CONJG( V( 8 ) )
00264          T8 = TAU*CONJG( V8 )
00265          DO 160 J = 1, N
00266             SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) +
00267      $            V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) +
00268      $            V7*C( 7, J ) + V8*C( 8, J )
00269             C( 1, J ) = C( 1, J ) - SUM*T1
00270             C( 2, J ) = C( 2, J ) - SUM*T2
00271             C( 3, J ) = C( 3, J ) - SUM*T3
00272             C( 4, J ) = C( 4, J ) - SUM*T4
00273             C( 5, J ) = C( 5, J ) - SUM*T5
00274             C( 6, J ) = C( 6, J ) - SUM*T6
00275             C( 7, J ) = C( 7, J ) - SUM*T7
00276             C( 8, J ) = C( 8, J ) - SUM*T8
00277   160    CONTINUE
00278          GO TO 410
00279   170    CONTINUE
00280 *
00281 *        Special code for 9 x 9 Householder
00282 *
00283          V1 = CONJG( V( 1 ) )
00284          T1 = TAU*CONJG( V1 )
00285          V2 = CONJG( V( 2 ) )
00286          T2 = TAU*CONJG( V2 )
00287          V3 = CONJG( V( 3 ) )
00288          T3 = TAU*CONJG( V3 )
00289          V4 = CONJG( V( 4 ) )
00290          T4 = TAU*CONJG( V4 )
00291          V5 = CONJG( V( 5 ) )
00292          T5 = TAU*CONJG( V5 )
00293          V6 = CONJG( V( 6 ) )
00294          T6 = TAU*CONJG( V6 )
00295          V7 = CONJG( V( 7 ) )
00296          T7 = TAU*CONJG( V7 )
00297          V8 = CONJG( V( 8 ) )
00298          T8 = TAU*CONJG( V8 )
00299          V9 = CONJG( V( 9 ) )
00300          T9 = TAU*CONJG( V9 )
00301          DO 180 J = 1, N
00302             SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) +
00303      $            V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) +
00304      $            V7*C( 7, J ) + V8*C( 8, J ) + V9*C( 9, J )
00305             C( 1, J ) = C( 1, J ) - SUM*T1
00306             C( 2, J ) = C( 2, J ) - SUM*T2
00307             C( 3, J ) = C( 3, J ) - SUM*T3
00308             C( 4, J ) = C( 4, J ) - SUM*T4
00309             C( 5, J ) = C( 5, J ) - SUM*T5
00310             C( 6, J ) = C( 6, J ) - SUM*T6
00311             C( 7, J ) = C( 7, J ) - SUM*T7
00312             C( 8, J ) = C( 8, J ) - SUM*T8
00313             C( 9, J ) = C( 9, J ) - SUM*T9
00314   180    CONTINUE
00315          GO TO 410
00316   190    CONTINUE
00317 *
00318 *        Special code for 10 x 10 Householder
00319 *
00320          V1 = CONJG( V( 1 ) )
00321          T1 = TAU*CONJG( V1 )
00322          V2 = CONJG( V( 2 ) )
00323          T2 = TAU*CONJG( V2 )
00324          V3 = CONJG( V( 3 ) )
00325          T3 = TAU*CONJG( V3 )
00326          V4 = CONJG( V( 4 ) )
00327          T4 = TAU*CONJG( V4 )
00328          V5 = CONJG( V( 5 ) )
00329          T5 = TAU*CONJG( V5 )
00330          V6 = CONJG( V( 6 ) )
00331          T6 = TAU*CONJG( V6 )
00332          V7 = CONJG( V( 7 ) )
00333          T7 = TAU*CONJG( V7 )
00334          V8 = CONJG( V( 8 ) )
00335          T8 = TAU*CONJG( V8 )
00336          V9 = CONJG( V( 9 ) )
00337          T9 = TAU*CONJG( V9 )
00338          V10 = CONJG( V( 10 ) )
00339          T10 = TAU*CONJG( V10 )
00340          DO 200 J = 1, N
00341             SUM = V1*C( 1, J ) + V2*C( 2, J ) + V3*C( 3, J ) +
00342      $            V4*C( 4, J ) + V5*C( 5, J ) + V6*C( 6, J ) +
00343      $            V7*C( 7, J ) + V8*C( 8, J ) + V9*C( 9, J ) +
00344      $            V10*C( 10, J )
00345             C( 1, J ) = C( 1, J ) - SUM*T1
00346             C( 2, J ) = C( 2, J ) - SUM*T2
00347             C( 3, J ) = C( 3, J ) - SUM*T3
00348             C( 4, J ) = C( 4, J ) - SUM*T4
00349             C( 5, J ) = C( 5, J ) - SUM*T5
00350             C( 6, J ) = C( 6, J ) - SUM*T6
00351             C( 7, J ) = C( 7, J ) - SUM*T7
00352             C( 8, J ) = C( 8, J ) - SUM*T8
00353             C( 9, J ) = C( 9, J ) - SUM*T9
00354             C( 10, J ) = C( 10, J ) - SUM*T10
00355   200    CONTINUE
00356          GO TO 410
00357       ELSE
00358 *
00359 *        Form  C * H, where H has order n.
00360 *
00361          GO TO ( 210, 230, 250, 270, 290, 310, 330, 350,
00362      $           370, 390 )N
00363 *
00364 *        Code for general N
00365 *
00366          CALL CLARF( SIDE, M, N, V, 1, TAU, C, LDC, WORK )
00367          GO TO 410
00368   210    CONTINUE
00369 *
00370 *        Special code for 1 x 1 Householder
00371 *
00372          T1 = ONE - TAU*V( 1 )*CONJG( V( 1 ) )
00373          DO 220 J = 1, M
00374             C( J, 1 ) = T1*C( J, 1 )
00375   220    CONTINUE
00376          GO TO 410
00377   230    CONTINUE
00378 *
00379 *        Special code for 2 x 2 Householder
00380 *
00381          V1 = V( 1 )
00382          T1 = TAU*CONJG( V1 )
00383          V2 = V( 2 )
00384          T2 = TAU*CONJG( V2 )
00385          DO 240 J = 1, M
00386             SUM = V1*C( J, 1 ) + V2*C( J, 2 )
00387             C( J, 1 ) = C( J, 1 ) - SUM*T1
00388             C( J, 2 ) = C( J, 2 ) - SUM*T2
00389   240    CONTINUE
00390          GO TO 410
00391   250    CONTINUE
00392 *
00393 *        Special code for 3 x 3 Householder
00394 *
00395          V1 = V( 1 )
00396          T1 = TAU*CONJG( V1 )
00397          V2 = V( 2 )
00398          T2 = TAU*CONJG( V2 )
00399          V3 = V( 3 )
00400          T3 = TAU*CONJG( V3 )
00401          DO 260 J = 1, M
00402             SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 )
00403             C( J, 1 ) = C( J, 1 ) - SUM*T1
00404             C( J, 2 ) = C( J, 2 ) - SUM*T2
00405             C( J, 3 ) = C( J, 3 ) - SUM*T3
00406   260    CONTINUE
00407          GO TO 410
00408   270    CONTINUE
00409 *
00410 *        Special code for 4 x 4 Householder
00411 *
00412          V1 = V( 1 )
00413          T1 = TAU*CONJG( V1 )
00414          V2 = V( 2 )
00415          T2 = TAU*CONJG( V2 )
00416          V3 = V( 3 )
00417          T3 = TAU*CONJG( V3 )
00418          V4 = V( 4 )
00419          T4 = TAU*CONJG( V4 )
00420          DO 280 J = 1, M
00421             SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) +
00422      $            V4*C( J, 4 )
00423             C( J, 1 ) = C( J, 1 ) - SUM*T1
00424             C( J, 2 ) = C( J, 2 ) - SUM*T2
00425             C( J, 3 ) = C( J, 3 ) - SUM*T3
00426             C( J, 4 ) = C( J, 4 ) - SUM*T4
00427   280    CONTINUE
00428          GO TO 410
00429   290    CONTINUE
00430 *
00431 *        Special code for 5 x 5 Householder
00432 *
00433          V1 = V( 1 )
00434          T1 = TAU*CONJG( V1 )
00435          V2 = V( 2 )
00436          T2 = TAU*CONJG( V2 )
00437          V3 = V( 3 )
00438          T3 = TAU*CONJG( V3 )
00439          V4 = V( 4 )
00440          T4 = TAU*CONJG( V4 )
00441          V5 = V( 5 )
00442          T5 = TAU*CONJG( V5 )
00443          DO 300 J = 1, M
00444             SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) +
00445      $            V4*C( J, 4 ) + V5*C( J, 5 )
00446             C( J, 1 ) = C( J, 1 ) - SUM*T1
00447             C( J, 2 ) = C( J, 2 ) - SUM*T2
00448             C( J, 3 ) = C( J, 3 ) - SUM*T3
00449             C( J, 4 ) = C( J, 4 ) - SUM*T4
00450             C( J, 5 ) = C( J, 5 ) - SUM*T5
00451   300    CONTINUE
00452          GO TO 410
00453   310    CONTINUE
00454 *
00455 *        Special code for 6 x 6 Householder
00456 *
00457          V1 = V( 1 )
00458          T1 = TAU*CONJG( V1 )
00459          V2 = V( 2 )
00460          T2 = TAU*CONJG( V2 )
00461          V3 = V( 3 )
00462          T3 = TAU*CONJG( V3 )
00463          V4 = V( 4 )
00464          T4 = TAU*CONJG( V4 )
00465          V5 = V( 5 )
00466          T5 = TAU*CONJG( V5 )
00467          V6 = V( 6 )
00468          T6 = TAU*CONJG( V6 )
00469          DO 320 J = 1, M
00470             SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) +
00471      $            V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 )
00472             C( J, 1 ) = C( J, 1 ) - SUM*T1
00473             C( J, 2 ) = C( J, 2 ) - SUM*T2
00474             C( J, 3 ) = C( J, 3 ) - SUM*T3
00475             C( J, 4 ) = C( J, 4 ) - SUM*T4
00476             C( J, 5 ) = C( J, 5 ) - SUM*T5
00477             C( J, 6 ) = C( J, 6 ) - SUM*T6
00478   320    CONTINUE
00479          GO TO 410
00480   330    CONTINUE
00481 *
00482 *        Special code for 7 x 7 Householder
00483 *
00484          V1 = V( 1 )
00485          T1 = TAU*CONJG( V1 )
00486          V2 = V( 2 )
00487          T2 = TAU*CONJG( V2 )
00488          V3 = V( 3 )
00489          T3 = TAU*CONJG( V3 )
00490          V4 = V( 4 )
00491          T4 = TAU*CONJG( V4 )
00492          V5 = V( 5 )
00493          T5 = TAU*CONJG( V5 )
00494          V6 = V( 6 )
00495          T6 = TAU*CONJG( V6 )
00496          V7 = V( 7 )
00497          T7 = TAU*CONJG( V7 )
00498          DO 340 J = 1, M
00499             SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) +
00500      $            V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) +
00501      $            V7*C( J, 7 )
00502             C( J, 1 ) = C( J, 1 ) - SUM*T1
00503             C( J, 2 ) = C( J, 2 ) - SUM*T2
00504             C( J, 3 ) = C( J, 3 ) - SUM*T3
00505             C( J, 4 ) = C( J, 4 ) - SUM*T4
00506             C( J, 5 ) = C( J, 5 ) - SUM*T5
00507             C( J, 6 ) = C( J, 6 ) - SUM*T6
00508             C( J, 7 ) = C( J, 7 ) - SUM*T7
00509   340    CONTINUE
00510          GO TO 410
00511   350    CONTINUE
00512 *
00513 *        Special code for 8 x 8 Householder
00514 *
00515          V1 = V( 1 )
00516          T1 = TAU*CONJG( V1 )
00517          V2 = V( 2 )
00518          T2 = TAU*CONJG( V2 )
00519          V3 = V( 3 )
00520          T3 = TAU*CONJG( V3 )
00521          V4 = V( 4 )
00522          T4 = TAU*CONJG( V4 )
00523          V5 = V( 5 )
00524          T5 = TAU*CONJG( V5 )
00525          V6 = V( 6 )
00526          T6 = TAU*CONJG( V6 )
00527          V7 = V( 7 )
00528          T7 = TAU*CONJG( V7 )
00529          V8 = V( 8 )
00530          T8 = TAU*CONJG( V8 )
00531          DO 360 J = 1, M
00532             SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) +
00533      $            V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) +
00534      $            V7*C( J, 7 ) + V8*C( J, 8 )
00535             C( J, 1 ) = C( J, 1 ) - SUM*T1
00536             C( J, 2 ) = C( J, 2 ) - SUM*T2
00537             C( J, 3 ) = C( J, 3 ) - SUM*T3
00538             C( J, 4 ) = C( J, 4 ) - SUM*T4
00539             C( J, 5 ) = C( J, 5 ) - SUM*T5
00540             C( J, 6 ) = C( J, 6 ) - SUM*T6
00541             C( J, 7 ) = C( J, 7 ) - SUM*T7
00542             C( J, 8 ) = C( J, 8 ) - SUM*T8
00543   360    CONTINUE
00544          GO TO 410
00545   370    CONTINUE
00546 *
00547 *        Special code for 9 x 9 Householder
00548 *
00549          V1 = V( 1 )
00550          T1 = TAU*CONJG( V1 )
00551          V2 = V( 2 )
00552          T2 = TAU*CONJG( V2 )
00553          V3 = V( 3 )
00554          T3 = TAU*CONJG( V3 )
00555          V4 = V( 4 )
00556          T4 = TAU*CONJG( V4 )
00557          V5 = V( 5 )
00558          T5 = TAU*CONJG( V5 )
00559          V6 = V( 6 )
00560          T6 = TAU*CONJG( V6 )
00561          V7 = V( 7 )
00562          T7 = TAU*CONJG( V7 )
00563          V8 = V( 8 )
00564          T8 = TAU*CONJG( V8 )
00565          V9 = V( 9 )
00566          T9 = TAU*CONJG( V9 )
00567          DO 380 J = 1, M
00568             SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) +
00569      $            V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) +
00570      $            V7*C( J, 7 ) + V8*C( J, 8 ) + V9*C( J, 9 )
00571             C( J, 1 ) = C( J, 1 ) - SUM*T1
00572             C( J, 2 ) = C( J, 2 ) - SUM*T2
00573             C( J, 3 ) = C( J, 3 ) - SUM*T3
00574             C( J, 4 ) = C( J, 4 ) - SUM*T4
00575             C( J, 5 ) = C( J, 5 ) - SUM*T5
00576             C( J, 6 ) = C( J, 6 ) - SUM*T6
00577             C( J, 7 ) = C( J, 7 ) - SUM*T7
00578             C( J, 8 ) = C( J, 8 ) - SUM*T8
00579             C( J, 9 ) = C( J, 9 ) - SUM*T9
00580   380    CONTINUE
00581          GO TO 410
00582   390    CONTINUE
00583 *
00584 *        Special code for 10 x 10 Householder
00585 *
00586          V1 = V( 1 )
00587          T1 = TAU*CONJG( V1 )
00588          V2 = V( 2 )
00589          T2 = TAU*CONJG( V2 )
00590          V3 = V( 3 )
00591          T3 = TAU*CONJG( V3 )
00592          V4 = V( 4 )
00593          T4 = TAU*CONJG( V4 )
00594          V5 = V( 5 )
00595          T5 = TAU*CONJG( V5 )
00596          V6 = V( 6 )
00597          T6 = TAU*CONJG( V6 )
00598          V7 = V( 7 )
00599          T7 = TAU*CONJG( V7 )
00600          V8 = V( 8 )
00601          T8 = TAU*CONJG( V8 )
00602          V9 = V( 9 )
00603          T9 = TAU*CONJG( V9 )
00604          V10 = V( 10 )
00605          T10 = TAU*CONJG( V10 )
00606          DO 400 J = 1, M
00607             SUM = V1*C( J, 1 ) + V2*C( J, 2 ) + V3*C( J, 3 ) +
00608      $            V4*C( J, 4 ) + V5*C( J, 5 ) + V6*C( J, 6 ) +
00609      $            V7*C( J, 7 ) + V8*C( J, 8 ) + V9*C( J, 9 ) +
00610      $            V10*C( J, 10 )
00611             C( J, 1 ) = C( J, 1 ) - SUM*T1
00612             C( J, 2 ) = C( J, 2 ) - SUM*T2
00613             C( J, 3 ) = C( J, 3 ) - SUM*T3
00614             C( J, 4 ) = C( J, 4 ) - SUM*T4
00615             C( J, 5 ) = C( J, 5 ) - SUM*T5
00616             C( J, 6 ) = C( J, 6 ) - SUM*T6
00617             C( J, 7 ) = C( J, 7 ) - SUM*T7
00618             C( J, 8 ) = C( J, 8 ) - SUM*T8
00619             C( J, 9 ) = C( J, 9 ) - SUM*T9
00620             C( J, 10 ) = C( J, 10 ) - SUM*T10
00621   400    CONTINUE
00622          GO TO 410
00623       END IF
00624   410 RETURN
00625 *
00626 *     End of CLARFX
00627 *
00628       END
 All Files Functions