Merge pull request #5306 from martin-frbg/lapack1131

martin-frbg · web-flow · commit f1097d1cba9d · 2025-06-15T15:00:05.000-07:00
Fix missing initialization leading to bypassing corner cases in C/ZGEQP3RK (Reference-LAPACK PR #1131)
diff --git a/lapack-netlib/SRC/cgeqp3rk.f b/lapack-netlib/SRC/cgeqp3rk.f
@@ -761,6 +761,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/lapack-netlib/SRC/zgeqp3rk.f b/lapack-netlib/SRC/zgeqp3rk.f
@@ -760,6 +760,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
`@@ -761,6 +761,7 @@ SUBROUTINE CGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,`
`761`	`761`	`* for the whole original matrix stored in A(1:M,1:N).`
`762`	`762`	`*`
`763`	`763`	`KP1 = ISAMAX( N, RWORK( 1 ), 1 )`
	`764`	`+ MAXC2NRM = RWORK( KP1 )`
`764`	`765`	`*`
`765`	`766`	`* ==================================================================.`
`766`	`767`	`*`
Original file line number	Diff line number	Diff line change
`@@ -760,6 +760,7 @@ SUBROUTINE ZGEQP3RK( M, N, NRHS, KMAX, ABSTOL, RELTOL, A, LDA,`
`760`	`760`	`* for the whole original matrix stored in A(1:M,1:N).`
`761`	`761`	`*`
`762`	`762`	`KP1 = IDAMAX( N, RWORK( 1 ), 1 )`
	`763`	`+ MAXC2NRM = RWORK( KP1 )`
`763`	`764`	`*`
`764`	`765`	`* ==================================================================.`
`765`	`766`	`*`