Updated array index calculations in cchkqp3rk.f, dchkqp3rk.f,

ACSimon33 · ACSimon33 · commit f6355dc7697a · 2023-12-15T17:54:40.000Z
schkqp3rk.f, and zchkqp3rk.f to use the leading dimension (LDA) instead
of the fixed size (M)
diff --git a/TESTING/LIN/cchkqp3rk.f b/TESTING/LIN/cchkqp3rk.f
@@ -720,8 +720,8 @@ SUBROUTINE CCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
 *
                      DO J = 1, KFACT-1, 1
 *
-                        DTEMP = (( ABS( A( (J-1)*M+J ) ) -
-     $                          ABS( A( (J)*M+J+1 ) ) ) /
+                        DTEMP = (( ABS( A( (J-1)*LDA+J ) ) -
+     $                          ABS( A( (J)*LDA+J+1 ) ) ) /
      $                          ABS( A(1) ) )
 *
                         IF( DTEMP.LT.ZERO ) THEN
diff --git a/TESTING/LIN/dchkqp3rk.f b/TESTING/LIN/dchkqp3rk.f
@@ -716,8 +716,8 @@ SUBROUTINE DCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
 *
                      DO J = 1, KFACT-1, 1
 
-                        DTEMP = (( ABS( A( (J-1)*M+J ) ) -
-     $                          ABS( A( (J)*M+J+1 ) ) ) /
+                        DTEMP = (( ABS( A( (J-1)*LDA+J ) ) -
+     $                          ABS( A( (J)*LDA+J+1 ) ) ) /
      $                          ABS( A(1) ) )
 *
                         IF( DTEMP.LT.ZERO ) THEN
diff --git a/TESTING/LIN/schkqp3rk.f b/TESTING/LIN/schkqp3rk.f
@@ -715,8 +715,8 @@ SUBROUTINE SCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
 *
                      DO J = 1, KFACT-1, 1
 
-                        DTEMP = (( ABS( A( (J-1)*M+J ) ) -
-     $                          ABS( A( (J)*M+J+1 ) ) ) /
+                        DTEMP = (( ABS( A( (J-1)*LDA+J ) ) -
+     $                          ABS( A( (J)*LDA+J+1 ) ) ) /
      $                          ABS( A(1) ) )
 *
                         IF( DTEMP.LT.ZERO ) THEN
diff --git a/TESTING/LIN/zchkqp3rk.f b/TESTING/LIN/zchkqp3rk.f
@@ -720,8 +720,8 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,
 *
                      DO J = 1, KFACT-1, 1
 *
-                        DTEMP = (( ABS( A( (J-1)*M+J ) ) -
-     $                          ABS( A( (J)*M+J+1 ) ) ) /
+                        DTEMP = (( ABS( A( (J-1)*LDA+J ) ) -
+     $                          ABS( A( (J)*LDA+J+1 ) ) ) /
      $                          ABS( A(1) ) )
 *
                         IF( DTEMP.LT.ZERO ) THEN

Original file line number	Diff line number	Diff line change
`@@ -720,8 +720,8 @@ SUBROUTINE CCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,`
`720`	`720`	`*`
`721`	`721`	`DO J = 1, KFACT-1, 1`
`722`	`722`	`*`
`723`		`- DTEMP = (( ABS( A( (J-1)*M+J ) ) -`
`724`		`- $ ABS( A( (J)*M+J+1 ) ) ) /`
	`723`	`+ DTEMP = (( ABS( A( (J-1)*LDA+J ) ) -`
	`724`	`+ $ ABS( A( (J)*LDA+J+1 ) ) ) /`
`725`	`725`	`$ ABS( A(1) ) )`
`726`	`726`	`*`
`727`	`727`	`IF( DTEMP.LT.ZERO ) THEN`
Original file line number	Diff line number	Diff line change
`@@ -716,8 +716,8 @@ SUBROUTINE DCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,`
`716`	`716`	`*`
`717`	`717`	`DO J = 1, KFACT-1, 1`
`718`	`718`
`719`		`- DTEMP = (( ABS( A( (J-1)*M+J ) ) -`
`720`		`- $ ABS( A( (J)*M+J+1 ) ) ) /`
	`719`	`+ DTEMP = (( ABS( A( (J-1)*LDA+J ) ) -`
	`720`	`+ $ ABS( A( (J)*LDA+J+1 ) ) ) /`
`721`	`721`	`$ ABS( A(1) ) )`
`722`	`722`	`*`
`723`	`723`	`IF( DTEMP.LT.ZERO ) THEN`
Original file line number	Diff line number	Diff line change
`@@ -715,8 +715,8 @@ SUBROUTINE SCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,`
`715`	`715`	`*`
`716`	`716`	`DO J = 1, KFACT-1, 1`
`717`	`717`
`718`		`- DTEMP = (( ABS( A( (J-1)*M+J ) ) -`
`719`		`- $ ABS( A( (J)*M+J+1 ) ) ) /`
	`718`	`+ DTEMP = (( ABS( A( (J-1)*LDA+J ) ) -`
	`719`	`+ $ ABS( A( (J)*LDA+J+1 ) ) ) /`
`720`	`720`	`$ ABS( A(1) ) )`
`721`	`721`	`*`
`722`	`722`	`IF( DTEMP.LT.ZERO ) THEN`
Original file line number	Diff line number	Diff line change
`@@ -720,8 +720,8 @@ SUBROUTINE ZCHKQP3RK( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL,`
`720`	`720`	`*`
`721`	`721`	`DO J = 1, KFACT-1, 1`
`722`	`722`	`*`
`723`		`- DTEMP = (( ABS( A( (J-1)*M+J ) ) -`
`724`		`- $ ABS( A( (J)*M+J+1 ) ) ) /`
	`723`	`+ DTEMP = (( ABS( A( (J-1)*LDA+J ) ) -`
	`724`	`+ $ ABS( A( (J)*LDA+J+1 ) ) ) /`
`725`	`725`	`$ ABS( A(1) ) )`
`726`	`726`	`*`
`727`	`727`	`IF( DTEMP.LT.ZERO ) THEN`