@@ -382,7 +382,8 @@ SUBROUTINE CGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
382
382
30 CONTINUE
383
383
40 CONTINUE
384
384
IF ( P.GT. L )
385
- $ CALL CLASET( ' Full' , P- L, N, CZERO, CZERO, B( L+1 , 1 ), LDB )
385
+ $ CALL CLASET( ' Full' , P- L, N, CZERO, CZERO, B( L+1 , 1 ),
386
+ $ LDB )
386
387
*
387
388
IF ( WANTQ ) THEN
388
389
*
@@ -444,17 +445,18 @@ SUBROUTINE CGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
444
445
*
445
446
* Update A12 := U**H*A12, where A12 = A( 1:M, N-L+1:N )
446
447
*
447
- CALL CUNM2R( ' Left' , ' Conjugate transpose' , M, L, MIN ( M, N- L ),
448
- $ A, LDA, TAU, A( 1 , N- L+1 ), LDA, WORK, INFO )
448
+ CALL CUNM2R( ' Left' , ' Conjugate transpose' , M, L,
449
+ $ MIN ( M, N- L ), A, LDA, TAU, A( 1 , N- L+1 ), LDA, WORK,
450
+ $ INFO )
449
451
*
450
452
IF ( WANTU ) THEN
451
453
*
452
454
* Copy the details of U, and form U
453
455
*
454
456
CALL CLASET( ' Full' , M, M, CZERO, CZERO, U, LDU )
455
457
IF ( M.GT. 1 )
456
- $ CALL CLACPY( ' Lower' , M-1 , N- L, A( 2 , 1 ), LDA, U( 2 , 1 ),
457
- $ LDU )
458
+ $ CALL CLACPY( ' Lower' , M-1 , N- L, A( 2 , 1 ), LDA,
459
+ $ U( 2 , 1 ), LDU )
458
460
CALL CUNG2R( M, M, MIN ( M, N- L ), U, LDU, TAU, WORK, INFO )
459
461
END IF
460
462
*
@@ -474,7 +476,8 @@ SUBROUTINE CGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
474
476
90 CONTINUE
475
477
100 CONTINUE
476
478
IF ( M.GT. K )
477
- $ CALL CLASET( ' Full' , M- K, N- L, CZERO, CZERO, A( K+1 , 1 ), LDA )
479
+ $ CALL CLASET( ' Full' , M- K, N- L, CZERO, CZERO, A( K+1 , 1 ),
480
+ $ LDA )
478
481
*
479
482
IF ( N- L.GT. K ) THEN
480
483
*
@@ -511,9 +514,9 @@ SUBROUTINE CGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
511
514
*
512
515
* Update U(:,K+1:M) := U(:,K+1:M)*U1
513
516
*
514
- CALL CUNM2R( ' Right' , ' No transpose' , M, M- K, MIN ( M - K, L ),
515
- $ A( K +1 , N - L +1 ), LDA, TAU, U( 1 , K +1 ), LDU ,
516
- $ WORK, INFO )
517
+ CALL CUNM2R( ' Right' , ' No transpose' , M, M- K,
518
+ $ MIN ( M - K, L ), A( K + 1 , N - L +1 ), LDA, TAU ,
519
+ $ U( 1 , K +1 ), LDU, WORK, INFO )
517
520
END IF
518
521
*
519
522
* Clean up
0 commit comments