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