Merge pull request #1131 from angsch/czgeqp3rk

langou · web-flow · commit 4758c863f59e · 2025-06-09T19:17:45.000-06:00
Fix bypassing corner cases due to uninitialized variable in [c,z]geqp3rk
diff --git a/SRC/cgeqp3rk.f b/SRC/cgeqp3rk.f
@@ -751,6 +751,7 @@ SUBROUTINE CGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
 *     for the whole original matrix stored in A(1:M,1:N).
 *
       KP1 = ISAMAX( N, RWORK( 1 ), 1 )
+      MAXC2NRM = RWORK( KP1 )
 *
 *     ==================================================================.
 *
diff --git a/SRC/zgeqp3rk.f b/SRC/zgeqp3rk.f
@@ -750,6 +750,7 @@ SUBROUTINE ZGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,
 *     for the whole original matrix stored in A(1:M,1:N).
 *
       KP1 = IDAMAX( N, RWORK( 1 ), 1 )
+      MAXC2NRM = RWORK( KP1 )
 *
 *     ==================================================================.
 *

Original file line number	Diff line number	Diff line change
`@@ -751,6 +751,7 @@ SUBROUTINE CGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,`
`751`	`751`	`* for the whole original matrix stored in A(1:M,1:N).`
`752`	`752`	`*`
`753`	`753`	`KP1 = ISAMAX( N, RWORK( 1 ), 1 )`
	`754`	`+ MAXC2NRM = RWORK( KP1 )`
`754`	`755`	`*`
`755`	`756`	`* ==================================================================.`
`756`	`757`	`*`
Original file line number	Diff line number	Diff line change
`@@ -750,6 +750,7 @@ SUBROUTINE ZGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,`
`750`	`750`	`* for the whole original matrix stored in A(1:M,1:N).`
`751`	`751`	`*`
`752`	`752`	`KP1 = IDAMAX( N, RWORK( 1 ), 1 )`
	`753`	`+ MAXC2NRM = RWORK( KP1 )`
`753`	`754`	`*`
`754`	`755`	`* ==================================================================.`
`755`	`756`	`*`