phynexis-deps/lapack
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Merge pull request #1202 from nakatamaho/xLAQP2RK

Browse Source 
fix(xlaqp2rk): do not modify [in] argument KMAX; introduce KBOUND
This commit is contained in:
langou
2026-03-23 08:30:31 -06:00
committed by GitHub
parent 6fe88c7ac8 4f13bb0d6c
commit 5dfd515e90
4 changed files with 20 additions and 20 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SRC/claqp2rk.f
+5 -5
View File
@@ -358,8 +358,8 @@
PARAMETER ( CZERO = ( 0.0E+0, 0.0E+0 ) )
* ..
* .. Local Scalars ..
INTEGER I, ITEMP, J, JMAXC2NRM, KK, KP, MINMNFACT,
$ MINMNUPDT
INTEGER I, ITEMP, J, JMAXC2NRM, KK, KP,
$ KBOUND, MINMNFACT, MINMNUPDT
REAL HUGEVAL, TAUNAN, TEMP, TEMP2, TOL3Z
* ..
* .. External Subroutines ..
@@ -390,13 +390,13 @@
*
MINMNFACT = MIN( M-IOFFSET, N )
MINMNUPDT = MIN( M-IOFFSET, N+NRHS )
KMAX = MIN( KMAX, MINMNFACT )
KBOUND = MIN( KMAX, MINMNFACT )
TOL3Z = SQRT( SLAMCH( 'Epsilon' ) )
HUGEVAL = SLAMCH( 'Overflow' )
*
* Compute the factorization, KK is the lomn loop index.
*
DO KK = 1, KMAX
DO KK = 1, KBOUND
*
I = IOFFSET + KK
*
@@ -674,7 +674,7 @@
* i.e. no condition was triggered to exit the routine.
* Set the number of factorized columns.
*
K = KMAX
K = KBOUND
*
* We reached the end of the loop, i.e. all KMAX columns were
* factorized, we need to set MAXC2NRMK and RELMAXC2NRMK before
SRC/dlaqp2rk.f
+5 -5
View File
@@ -355,8 +355,8 @@
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
* ..
* .. Local Scalars ..
INTEGER I, ITEMP, J, JMAXC2NRM, KK, KP, MINMNFACT,
$ MINMNUPDT
INTEGER I, ITEMP, J, JMAXC2NRM, KK, KP,
$ KBOUND, MINMNFACT, MINMNUPDT
DOUBLE PRECISION HUGEVAL, TEMP, TEMP2, TOL3Z
* ..
* .. External Subroutines ..
@@ -387,13 +387,13 @@
*
MINMNFACT = MIN( M-IOFFSET, N )
MINMNUPDT = MIN( M-IOFFSET, N+NRHS )
KMAX = MIN( KMAX, MINMNFACT )
KBOUND = MIN( KMAX, MINMNFACT )
TOL3Z = SQRT( DLAMCH( 'Epsilon' ) )
HUGEVAL = DLAMCH( 'Overflow' )
*
* Compute the factorization, KK is the lomn loop index.
*
DO KK = 1, KMAX
DO KK = 1, KBOUND
*
I = IOFFSET + KK
*
@@ -663,7 +663,7 @@
* i.e. no condition was triggered to exit the routine.
* Set the number of factorized columns.
*
K = KMAX
K = KBOUND
*
* We reached the end of the loop, i.e. all KMAX columns were
* factorized, we need to set MAXC2NRMK and RELMAXC2NRMK before
SRC/slaqp2rk.f
+5 -5
View File
@@ -355,8 +355,8 @@
PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )
* ..
* .. Local Scalars ..
INTEGER I, ITEMP, J, JMAXC2NRM, KK, KP, MINMNFACT,
$ MINMNUPDT
INTEGER I, ITEMP, J, JMAXC2NRM, KK, KP,
$ KBOUND, MINMNFACT, MINMNUPDT
REAL HUGEVAL, TEMP, TEMP2, TOL3Z
* ..
* .. External Subroutines ..
@@ -387,13 +387,13 @@
*
MINMNFACT = MIN( M-IOFFSET, N )
MINMNUPDT = MIN( M-IOFFSET, N+NRHS )
KMAX = MIN( KMAX, MINMNFACT )
KBOUND = MIN( KMAX, MINMNFACT )
TOL3Z = SQRT( SLAMCH( 'Epsilon' ) )
HUGEVAL = SLAMCH( 'Overflow' )
*
* Compute the factorization, KK is the lomn loop index.
*
DO KK = 1, KMAX
DO KK = 1, KBOUND
*
I = IOFFSET + KK
*
@@ -663,7 +663,7 @@
* i.e. no condition was triggered to exit the routine.
* Set the number of factorized columns.
*
K = KMAX
K = KBOUND
*
* We reached the end of the loop, i.e. all KMAX columns were
* factorized, we need to set MAXC2NRMK and RELMAXC2NRMK before
SRC/zlaqp2rk.f
+5 -5
View File
@@ -359,8 +359,8 @@
$ CONE = ( 1.0D+0, 0.0D+0 ) )
* ..
* .. Local Scalars ..
INTEGER I, ITEMP, J, JMAXC2NRM, KK, KP, MINMNFACT,
$ MINMNUPDT
INTEGER I, ITEMP, J, JMAXC2NRM, KK, KP,
$ KBOUND, MINMNFACT, MINMNUPDT
DOUBLE PRECISION HUGEVAL, TAUNAN, TEMP, TEMP2, TOL3Z
* ..
* .. External Subroutines ..
@@ -391,13 +391,13 @@
*
MINMNFACT = MIN( M-IOFFSET, N )
MINMNUPDT = MIN( M-IOFFSET, N+NRHS )
KMAX = MIN( KMAX, MINMNFACT )
KBOUND = MIN( KMAX, MINMNFACT )
TOL3Z = SQRT( DLAMCH( 'Epsilon' ) )
HUGEVAL = DLAMCH( 'Overflow' )
*
* Compute the factorization, KK is the lomn loop index.
*
DO KK = 1, KMAX
DO KK = 1, KBOUND
*
I = IOFFSET + KK
*
@@ -675,7 +675,7 @@
* i.e. no condition was triggered to exit the routine.
* Set the number of factorized columns.
*
K = KMAX
K = KBOUND
*
* We reached the end of the loop, i.e. all KMAX columns were
* factorized, we need to set MAXC2NRMK and RELMAXC2NRMK before