From d3feb81c17112768d648a33caed7afc5599e256d Mon Sep 17 00:00:00 2001 From: Christoph Conrads Date: Fri, 24 Nov 2023 16:33:51 +0100 Subject: [PATCH] xBBCSD: fix uninitialized read fixes #943 --- SRC/cbbcsd.f | 2 +- SRC/dbbcsd.f | 2 +- SRC/sbbcsd.f | 2 +- SRC/zbbcsd.f | 2 +- 4 files changed, 4 insertions(+), 4 deletions(-) diff --git a/SRC/cbbcsd.f b/SRC/cbbcsd.f index e67be4ae36..8be43f5554 100644 --- a/SRC/cbbcsd.f +++ b/SRC/cbbcsd.f @@ -805,7 +805,7 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q, CALL SLARTGP( B22BULGE, B22E(I-1), RWORK(IU2SN+I-1), $ RWORK(IU2CS+I-1), R ) ELSE IF( NU .LT. MU ) THEN - CALL SLARTGS( B21E(I), B21E(I+1), NU, RWORK(IU2CS+I-1), + CALL SLARTGS( B21E(I), B21D(I+1), NU, RWORK(IU2CS+I-1), $ RWORK(IU2SN+I-1) ) ELSE CALL SLARTGS( B22D(I), B22E(I), MU, RWORK(IU2CS+I-1), diff --git a/SRC/dbbcsd.f b/SRC/dbbcsd.f index 0c94743822..d382ae792d 100644 --- a/SRC/dbbcsd.f +++ b/SRC/dbbcsd.f @@ -805,7 +805,7 @@ SUBROUTINE DBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q, CALL DLARTGP( B22BULGE, B22E(I-1), WORK(IU2SN+I-1), $ WORK(IU2CS+I-1), R ) ELSE IF( NU .LT. MU ) THEN - CALL DLARTGS( B21E(I), B21E(I+1), NU, WORK(IU2CS+I-1), + CALL DLARTGS( B21E(I), B21D(I+1), NU, WORK(IU2CS+I-1), $ WORK(IU2SN+I-1) ) ELSE CALL DLARTGS( B22D(I), B22E(I), MU, WORK(IU2CS+I-1), diff --git a/SRC/sbbcsd.f b/SRC/sbbcsd.f index 0e84bff8be..2ca5242bae 100644 --- a/SRC/sbbcsd.f +++ b/SRC/sbbcsd.f @@ -805,7 +805,7 @@ SUBROUTINE SBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q, CALL SLARTGP( B22BULGE, B22E(I-1), WORK(IU2SN+I-1), $ WORK(IU2CS+I-1), R ) ELSE IF( NU .LT. MU ) THEN - CALL SLARTGS( B21E(I), B21E(I+1), NU, WORK(IU2CS+I-1), + CALL SLARTGS( B21E(I), B21D(I+1), NU, WORK(IU2CS+I-1), $ WORK(IU2SN+I-1) ) ELSE CALL SLARTGS( B22D(I), B22E(I), MU, WORK(IU2CS+I-1), diff --git a/SRC/zbbcsd.f b/SRC/zbbcsd.f index 4083b2738f..02b7ee7b98 100644 --- a/SRC/zbbcsd.f +++ b/SRC/zbbcsd.f @@ -804,7 +804,7 @@ SUBROUTINE ZBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q, CALL DLARTGP( B22BULGE, B22E(I-1), RWORK(IU2SN+I-1), $ RWORK(IU2CS+I-1), R ) ELSE IF( NU .LT. MU ) THEN - CALL DLARTGS( B21E(I), B21E(I+1), NU, RWORK(IU2CS+I-1), + CALL DLARTGS( B21E(I), B21D(I+1), NU, RWORK(IU2CS+I-1), $ RWORK(IU2SN+I-1) ) ELSE CALL DLARTGS( B22D(I), B22E(I), MU, RWORK(IU2CS+I-1),