Skip to content

Commit 326bb27

Browse files
perazzperazz
committed
bound-check calls to stdlib_*laset from *gesdd 
1 parent 627d72c commit 326bb27

6 files changed

+24
-24
lines changed

‎src/stdlib_linalg_lapack_c.fypp‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -65732,7 +65732,7 @@ module stdlib_linalg_lapack_c
6573265732
call stdlib_cgeqrf( m, n, a, lda, work( itau ), work( nwork ),lwork-nwork+1, &
6573365733
ierr )
6573465734
! zero out below r
65735-
call stdlib_claset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
65735+
if (n>1) call stdlib_claset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
6573665736
ie = 1
6573765737
itauq = 1
6573865738
itaup = itauq + n
@@ -65918,7 +65918,7 @@ module stdlib_linalg_lapack_c
6591865918
call stdlib_cungqr( m, m, n, u, ldu, work( itau ),work( nwork ), lwork-nwork+&
6591965919
1, ierr )
6592065920
! produce r in a, zeroing out below it
65921-
call stdlib_claset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
65921+
if (n>1) call stdlib_claset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
6592265922
ie = 1
6592365923
itauq = itau
6592465924
itaup = itauq + n
@@ -66294,7 +66294,7 @@ module stdlib_linalg_lapack_c
6629466294
call stdlib_cgelqf( m, n, a, lda, work( itau ), work( nwork ),lwork-nwork+1, &
6629566295
ierr )
6629666296
! zero out above l
66297-
call stdlib_claset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
66297+
if (m>1) call stdlib_claset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
6629866298
ie = 1
6629966299
itauq = 1
6630066300
itaup = itauq + m
@@ -66485,7 +66485,7 @@ module stdlib_linalg_lapack_c
6648566485
call stdlib_cunglq( n, n, m, vt, ldvt, work( itau ),work( nwork ), lwork-&
6648666486
nwork+1, ierr )
6648766487
! produce l in a, zeroing out above it
66488-
call stdlib_claset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
66488+
if (m>1) call stdlib_claset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
6648966489
ie = 1
6649066490
itauq = itau
6649166491
itaup = itauq + m

‎src/stdlib_linalg_lapack_d.fypp‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -75410,7 +75410,7 @@ module stdlib_linalg_lapack_d
7541075410
call stdlib_dgeqrf( m, n, a, lda, work( itau ), work( nwork ),lwork - nwork + &
7541175411
1, ierr )
7541275412
! zero out below r
75413-
call stdlib_dlaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
75413+
if (n>1) call stdlib_dlaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
7541475414
ie = 1
7541575415
itauq = ie + n
7541675416
itaup = itauq + n
@@ -75560,7 +75560,7 @@ module stdlib_linalg_lapack_d
7556075560
call stdlib_dorgqr( m, m, n, u, ldu, work( itau ),work( nwork ), lwork - &
7556175561
nwork + 1, ierr )
7556275562
! produce r in a, zeroing out other entries
75563-
call stdlib_dlaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
75563+
if (n>1) call stdlib_dlaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
7556475564
ie = itau
7556575565
itauq = ie + n
7556675566
itaup = itauq + n
@@ -75726,7 +75726,7 @@ module stdlib_linalg_lapack_d
7572675726
call stdlib_dgelqf( m, n, a, lda, work( itau ), work( nwork ),lwork - nwork + &
7572775727
1, ierr )
7572875728
! zero out above l
75729-
call stdlib_dlaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
75729+
if (m>1) call stdlib_dlaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
7573075730
ie = 1
7573175731
itauq = ie + m
7573275732
itaup = itauq + m
@@ -75880,7 +75880,7 @@ module stdlib_linalg_lapack_d
7588075880
call stdlib_dorglq( n, n, m, vt, ldvt, work( itau ),work( nwork ), lwork - &
7588175881
nwork + 1, ierr )
7588275882
! produce l in a, zeroing out other entries
75883-
call stdlib_dlaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
75883+
if (m>1) call stdlib_dlaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
7588475884
ie = itau
7588575885
itauq = ie + m
7588675886
itaup = itauq + m

‎src/stdlib_linalg_lapack_q.fypp‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -10442,7 +10442,7 @@ module stdlib_linalg_lapack_q
1044210442
call stdlib_qgeqrf( m, n, a, lda, work( itau ), work( nwork ),lwork - nwork + &
1044310443
1, ierr )
1044410444
! zero out below r
10445-
call stdlib_qlaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
10445+
if (n>1) call stdlib_qlaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
1044610446
ie = 1
1044710447
itauq = ie + n
1044810448
itaup = itauq + n
@@ -10592,7 +10592,7 @@ module stdlib_linalg_lapack_q
1059210592
call stdlib_qorgqr( m, m, n, u, ldu, work( itau ),work( nwork ), lwork - &
1059310593
nwork + 1, ierr )
1059410594
! produce r in a, zeroing out other entries
10595-
call stdlib_qlaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
10595+
if (n>1) call stdlib_qlaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
1059610596
ie = itau
1059710597
itauq = ie + n
1059810598
itaup = itauq + n
@@ -10758,7 +10758,7 @@ module stdlib_linalg_lapack_q
1075810758
call stdlib_qgelqf( m, n, a, lda, work( itau ), work( nwork ),lwork - nwork + &
1075910759
1, ierr )
1076010760
! zero out above l
10761-
call stdlib_qlaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
10761+
if (m>1) call stdlib_qlaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
1076210762
ie = 1
1076310763
itauq = ie + m
1076410764
itaup = itauq + m
@@ -10912,7 +10912,7 @@ module stdlib_linalg_lapack_q
1091210912
call stdlib_qorglq( n, n, m, vt, ldvt, work( itau ),work( nwork ), lwork - &
1091310913
nwork + 1, ierr )
1091410914
! produce l in a, zeroing out other entries
10915-
call stdlib_qlaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
10915+
if (m>1) call stdlib_qlaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
1091610916
ie = itau
1091710917
itauq = ie + m
1091810918
itaup = itauq + m

‎src/stdlib_linalg_lapack_s.fypp‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -72903,7 +72903,7 @@ module stdlib_linalg_lapack_s
7290372903
call stdlib_sgeqrf( m, n, a, lda, work( itau ), work( nwork ),lwork - nwork + &
7290472904
1, ierr )
7290572905
! zero out below r
72906-
call stdlib_slaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
72906+
if (n>1) call stdlib_slaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
7290772907
ie = 1
7290872908
itauq = ie + n
7290972909
itaup = itauq + n
@@ -73053,7 +73053,7 @@ module stdlib_linalg_lapack_s
7305373053
call stdlib_sorgqr( m, m, n, u, ldu, work( itau ),work( nwork ), lwork - &
7305473054
nwork + 1, ierr )
7305573055
! produce r in a, zeroing out other entries
73056-
call stdlib_slaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
73056+
if (n>1) call stdlib_slaset( 'L', n-1, n-1, zero, zero, a( 2, 1 ), lda )
7305773057
ie = itau
7305873058
itauq = ie + n
7305973059
itaup = itauq + n
@@ -73219,7 +73219,7 @@ module stdlib_linalg_lapack_s
7321973219
call stdlib_sgelqf( m, n, a, lda, work( itau ), work( nwork ),lwork - nwork + &
7322073220
1, ierr )
7322173221
! zero out above l
73222-
call stdlib_slaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
73222+
if (m>1) call stdlib_slaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
7322373223
ie = 1
7322473224
itauq = ie + m
7322573225
itaup = itauq + m
@@ -73373,7 +73373,7 @@ module stdlib_linalg_lapack_s
7337373373
call stdlib_sorglq( n, n, m, vt, ldvt, work( itau ),work( nwork ), lwork - &
7337473374
nwork + 1, ierr )
7337573375
! produce l in a, zeroing out other entries
73376-
call stdlib_slaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
73376+
if (m>1) call stdlib_slaset( 'U', m-1, m-1, zero, zero, a( 1, 2 ), lda )
7337773377
ie = itau
7337873378
itauq = ie + m
7337973379
itaup = itauq + m

‎src/stdlib_linalg_lapack_w.fypp‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -10802,7 +10802,7 @@ module stdlib_linalg_lapack_w
1080210802
call stdlib_wgeqrf( m, n, a, lda, work( itau ), work( nwork ),lwork-nwork+1, &
1080310803
ierr )
1080410804
! zero out below r
10805-
call stdlib_wlaset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
10805+
if (n>1) call stdlib_wlaset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
1080610806
ie = 1
1080710807
itauq = 1
1080810808
itaup = itauq + n
@@ -10988,7 +10988,7 @@ module stdlib_linalg_lapack_w
1098810988
call stdlib_wungqr( m, m, n, u, ldu, work( itau ),work( nwork ), lwork-nwork+&
1098910989
1, ierr )
1099010990
! produce r in a, zeroing out below it
10991-
call stdlib_wlaset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
10991+
if (n>1) call stdlib_wlaset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
1099210992
ie = 1
1099310993
itauq = itau
1099410994
itaup = itauq + n
@@ -11364,7 +11364,7 @@ module stdlib_linalg_lapack_w
1136411364
call stdlib_wgelqf( m, n, a, lda, work( itau ), work( nwork ),lwork-nwork+1, &
1136511365
ierr )
1136611366
! zero out above l
11367-
call stdlib_wlaset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
11367+
if (m>1) call stdlib_wlaset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
1136811368
ie = 1
1136911369
itauq = 1
1137011370
itaup = itauq + m
@@ -11555,7 +11555,7 @@ module stdlib_linalg_lapack_w
1155511555
call stdlib_wunglq( n, n, m, vt, ldvt, work( itau ),work( nwork ), lwork-&
1155611556
nwork+1, ierr )
1155711557
! produce l in a, zeroing out above it
11558-
call stdlib_wlaset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
11558+
if (m>1) call stdlib_wlaset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
1155911559
ie = 1
1156011560
itauq = itau
1156111561
itaup = itauq + m

‎src/stdlib_linalg_lapack_z.fypp‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -66211,7 +66211,7 @@ module stdlib_linalg_lapack_z
6621166211
call stdlib_zgeqrf( m, n, a, lda, work( itau ), work( nwork ),lwork-nwork+1, &
6621266212
ierr )
6621366213
! zero out below r
66214-
call stdlib_zlaset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
66214+
if (n>1) call stdlib_zlaset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
6621566215
ie = 1
6621666216
itauq = 1
6621766217
itaup = itauq + n
@@ -66397,7 +66397,7 @@ module stdlib_linalg_lapack_z
6639766397
call stdlib_zungqr( m, m, n, u, ldu, work( itau ),work( nwork ), lwork-nwork+&
6639866398
1, ierr )
6639966399
! produce r in a, zeroing out below it
66400-
call stdlib_zlaset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
66400+
if (n>1) call stdlib_zlaset( 'L', n-1, n-1, czero, czero, a( 2, 1 ),lda )
6640166401
ie = 1
6640266402
itauq = itau
6640366403
itaup = itauq + n
@@ -66773,7 +66773,7 @@ module stdlib_linalg_lapack_z
6677366773
call stdlib_zgelqf( m, n, a, lda, work( itau ), work( nwork ),lwork-nwork+1, &
6677466774
ierr )
6677566775
! zero out above l
66776-
call stdlib_zlaset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
66776+
if (m>1) call stdlib_zlaset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
6677766777
ie = 1
6677866778
itauq = 1
6677966779
itaup = itauq + m
@@ -66964,7 +66964,7 @@ module stdlib_linalg_lapack_z
6696466964
call stdlib_zunglq( n, n, m, vt, ldvt, work( itau ),work( nwork ), lwork-&
6696566965
nwork+1, ierr )
6696666966
! produce l in a, zeroing out above it
66967-
call stdlib_zlaset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
66967+
if (m>1) call stdlib_zlaset( 'U', m-1, m-1, czero, czero, a( 1, 2 ),lda )
6696866968
ie = 1
6696966969
itauq = itau
6697066970
itaup = itauq + m

0 commit comments

Comments
 (0)