diff options
author | lawrence.mulholland <lawrence.mulholland@8a072113-8704-0410-8d35-dd094bca7971> | 2016-01-27 17:24:52 +0000 |
---|---|---|
committer | lawrence.mulholland <lawrence.mulholland@8a072113-8704-0410-8d35-dd094bca7971> | 2016-01-27 17:24:52 +0000 |
commit | 6e3304a458b4c9e33b9a30a8e5a4345c0b2339ab (patch) | |
tree | 1544a43a18a78567c7f67dea750444a668a0eb06 /SRC/dgesvd.f | |
parent | 5f9dd6e3f5fc477d420dba82be1227b46856ecab (diff) | |
download | lapack-6e3304a458b4c9e33b9a30a8e5a4345c0b2339ab.tar.gz lapack-6e3304a458b4c9e33b9a30a8e5a4345c0b2339ab.tar.bz2 lapack-6e3304a458b4c9e33b9a30a8e5a4345c0b2339ab.zip |
protect calls with args n-1, a(2,1) against n=1 case
Diffstat (limited to 'SRC/dgesvd.f')
-rw-r--r-- | SRC/dgesvd.f | 29 |
1 files changed, 20 insertions, 9 deletions
diff --git a/SRC/dgesvd.f b/SRC/dgesvd.f index 898570b6..f3034042 100644 --- a/SRC/dgesvd.f +++ b/SRC/dgesvd.f @@ -692,7 +692,10 @@ * * Zero out below R * - CALL DLASET( 'L', N-1, N-1, ZERO, ZERO, A( 2, 1 ), LDA ) + IF( N .GT. 1 ) THEN + CALL DLASET( 'L', N-1, N-1, ZERO, ZERO, A( 2, 1 ), + $ LDA ) + END IF IE = 1 ITAUQ = IE + N ITAUP = ITAUQ + N @@ -1121,8 +1124,10 @@ * * Zero out below R in A * - CALL DLASET( 'L', N-1, N-1, ZERO, ZERO, A( 2, 1 ), - $ LDA ) + IF( N .GT. 1 ) THEN + CALL DLASET( 'L', N-1, N-1, ZERO, ZERO, + $ A( 2, 1 ), LDA ) + END IF * * Bidiagonalize R in A * (Workspace: need 4*N, prefer 3*N+2*N*NB) @@ -1284,8 +1289,10 @@ * * Zero out below R in A * - CALL DLASET( 'L', N-1, N-1, ZERO, ZERO, A( 2, 1 ), - $ LDA ) + IF( N .GT. 1 ) THEN + CALL DLASET( 'L', N-1, N-1, ZERO, ZERO, + $ A( 2, 1 ), LDA ) + END IF * * Bidiagonalize R in A * (Workspace: need 4*N, prefer 3*N+2*N*NB) @@ -1587,8 +1594,10 @@ * * Zero out below R in A * - CALL DLASET( 'L', N-1, N-1, ZERO, ZERO, A( 2, 1 ), - $ LDA ) + IF( N .GT. 1 ) THEN + CALL DLASET( 'L', N-1, N-1, ZERO, ZERO, + $ A( 2, 1 ), LDA ) + END IF * * Bidiagonalize R in A * (Workspace: need 4*N, prefer 3*N+2*N*NB) @@ -1755,8 +1764,10 @@ * * Zero out below R in A * - CALL DLASET( 'L', N-1, N-1, ZERO, ZERO, A( 2, 1 ), - $ LDA ) + IF( N .GT. 1 ) THEN + CALL DLASET( 'L', N-1, N-1, ZERO, ZERO, + $ A( 2, 1 ), LDA ) + END IF * * Bidiagonalize R in A * (Workspace: need 4*N, prefer 3*N+2*N*NB) |