Fix argument in SLASET call to zero S

martin-frbg · web-flow · commit a6a52a73f7b0 · 2018-11-10T17:16:53.000+01:00
fixes #1859 in accordance with https://github.com/LAPACK-Reference/issue/296
diff --git a/lapack-netlib/SRC/sgelss.f b/lapack-netlib/SRC/sgelss.f
@@ -407,7 +407,7 @@ SUBROUTINE SGELSS( M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK,
 *        Matrix all zero. Return zero solution.
 *
          CALL SLASET( 'F', MAX( M, N ), NRHS, ZERO, ZERO, B, LDB )
-         CALL SLASET( 'F', MINMN, 1, ZERO, ZERO, S, 1 )
+         CALL SLASET( 'F', MINMN, 1, ZERO, ZERO, S, MINMN )
          RANK = 0
          GO TO 70
       END IF

Original file line number	Diff line number	Diff line change
`@@ -407,7 +407,7 @@ SUBROUTINE SGELSS( M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK,`
`407`	`407`	`* Matrix all zero. Return zero solution.`
`408`	`408`	`*`
`409`	`409`	`CALL SLASET( 'F', MAX( M, N ), NRHS, ZERO, ZERO, B, LDB )`
`410`		`- CALL SLASET( 'F', MINMN, 1, ZERO, ZERO, S, 1 )`
	`410`	`+ CALL SLASET( 'F', MINMN, 1, ZERO, ZERO, S, MINMN )`
`411`	`411`	`RANK = 0`
`412`	`412`	`GO TO 70`
`413`	`413`	`END IF`