diff --git a/lapack-netlib/SRC/DEPRECATED/cgelqs.f b/lapack-netlib/SRC/DEPRECATED/cgelqs.f index 47e17a5830..aba3632a74 100644 --- a/lapack-netlib/SRC/DEPRECATED/cgelqs.f +++ b/lapack-netlib/SRC/DEPRECATED/cgelqs.f @@ -16,7 +16,7 @@ * .. * .. Array Arguments .. * COMPLEX A( LDA, * ), B( LDB, * ), TAU( * ), -* $ WORK( LWORK ) +* $ WORK( * ) * .. * * @@ -128,7 +128,7 @@ SUBROUTINE CGELQS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, * .. * .. Array Arguments .. COMPLEX A( LDA, * ), B( LDB, * ), TAU( * ), - $ WORK( LWORK ) + $ WORK( * ) * .. * * ===================================================================== diff --git a/lapack-netlib/SRC/DEPRECATED/cgeqrs.f b/lapack-netlib/SRC/DEPRECATED/cgeqrs.f index 13ac7f74fd..9d0527283d 100644 --- a/lapack-netlib/SRC/DEPRECATED/cgeqrs.f +++ b/lapack-netlib/SRC/DEPRECATED/cgeqrs.f @@ -16,7 +16,7 @@ * .. * .. Array Arguments .. * COMPLEX A( LDA, * ), B( LDB, * ), TAU( * ), -* $ WORK( LWORK ) +* $ WORK( * ) * .. * * @@ -128,7 +128,7 @@ SUBROUTINE CGEQRS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, * .. * .. Array Arguments .. COMPLEX A( LDA, * ), B( LDB, * ), TAU( * ), - $ WORK( LWORK ) + $ WORK( * ) * .. * * ===================================================================== diff --git a/lapack-netlib/SRC/DEPRECATED/dgelqs.f b/lapack-netlib/SRC/DEPRECATED/dgelqs.f index ecbb5893c3..1bab678901 100644 --- a/lapack-netlib/SRC/DEPRECATED/dgelqs.f +++ b/lapack-netlib/SRC/DEPRECATED/dgelqs.f @@ -16,7 +16,7 @@ * .. * .. Array Arguments .. * DOUBLE PRECISION A( LDA, * ), B( LDB, * ), TAU( * ), -* $ WORK( LWORK ) +* $ WORK( * ) * .. * * @@ -128,7 +128,7 @@ SUBROUTINE DGELQS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, * .. * .. Array Arguments .. DOUBLE PRECISION A( LDA, * ), B( LDB, * ), TAU( * ), - $ WORK( LWORK ) + $ WORK( * ) * .. * * ===================================================================== @@ -174,18 +174,19 @@ SUBROUTINE DGELQS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, * * Solve L*X = B(1:m,:) * - CALL DTRSM( 'Left', 'Lower', 'No transpose', 'Non-unit', M, NRHS, - $ ONE, A, LDA, B, LDB ) + CALL DTRSM( 'Left', 'Lower', 'No transpose', 'Non-unit', M, + $ NRHS, ONE, A, LDA, B, LDB ) * * Set B(m+1:n,:) to zero * IF( M.LT.N ) - $ CALL DLASET( 'Full', N-M, NRHS, ZERO, ZERO, B( M+1, 1 ), LDB ) + $ CALL DLASET( 'Full', N-M, NRHS, ZERO, ZERO, B( M+1, 1 ), + $ LDB ) * * B := Q' * B * - CALL DORMLQ( 'Left', 'Transpose', N, NRHS, M, A, LDA, TAU, B, LDB, - $ WORK, LWORK, INFO ) + CALL DORMLQ( 'Left', 'Transpose', N, NRHS, M, A, LDA, TAU, B, + $ LDB, WORK, LWORK, INFO ) * RETURN * diff --git a/lapack-netlib/SRC/DEPRECATED/dgeqrs.f b/lapack-netlib/SRC/DEPRECATED/dgeqrs.f index bfb7bd8bb8..e3e6c4048e 100644 --- a/lapack-netlib/SRC/DEPRECATED/dgeqrs.f +++ b/lapack-netlib/SRC/DEPRECATED/dgeqrs.f @@ -16,7 +16,7 @@ * .. * .. Array Arguments .. * DOUBLE PRECISION A( LDA, * ), B( LDB, * ), TAU( * ), -* $ WORK( LWORK ) +* $ WORK( * ) * .. * * @@ -128,7 +128,7 @@ SUBROUTINE DGEQRS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, * .. * .. Array Arguments .. DOUBLE PRECISION A( LDA, * ), B( LDB, * ), TAU( * ), - $ WORK( LWORK ) + $ WORK( * ) * .. * * ===================================================================== diff --git a/lapack-netlib/SRC/DEPRECATED/sgelqs.f b/lapack-netlib/SRC/DEPRECATED/sgelqs.f index 83afb4690b..2b1dd44b71 100644 --- a/lapack-netlib/SRC/DEPRECATED/sgelqs.f +++ b/lapack-netlib/SRC/DEPRECATED/sgelqs.f @@ -16,7 +16,7 @@ * .. * .. Array Arguments .. * REAL A( LDA, * ), B( LDB, * ), TAU( * ), -* $ WORK( LWORK ) +* $ WORK( * ) * .. * * @@ -128,7 +128,7 @@ SUBROUTINE SGELQS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, * .. * .. Array Arguments .. REAL A( LDA, * ), B( LDB, * ), TAU( * ), - $ WORK( LWORK ) + $ WORK( * ) * .. * * ===================================================================== @@ -174,18 +174,19 @@ SUBROUTINE SGELQS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, * * Solve L*X = B(1:m,:) * - CALL STRSM( 'Left', 'Lower', 'No transpose', 'Non-unit', M, NRHS, - $ ONE, A, LDA, B, LDB ) + CALL STRSM( 'Left', 'Lower', 'No transpose', 'Non-unit', M, + $ NRHS, ONE, A, LDA, B, LDB ) * * Set B(m+1:n,:) to zero * IF( M.LT.N ) - $ CALL SLASET( 'Full', N-M, NRHS, ZERO, ZERO, B( M+1, 1 ), LDB ) + $ CALL SLASET( 'Full', N-M, NRHS, ZERO, ZERO, B( M+1, 1 ), + $ LDB ) * * B := Q' * B * - CALL SORMLQ( 'Left', 'Transpose', N, NRHS, M, A, LDA, TAU, B, LDB, - $ WORK, LWORK, INFO ) + CALL SORMLQ( 'Left', 'Transpose', N, NRHS, M, A, LDA, TAU, B, + $ LDB, WORK, LWORK, INFO ) * RETURN * diff --git a/lapack-netlib/SRC/DEPRECATED/sgeqrs.f b/lapack-netlib/SRC/DEPRECATED/sgeqrs.f index ed11489104..bdbad5dcb4 100644 --- a/lapack-netlib/SRC/DEPRECATED/sgeqrs.f +++ b/lapack-netlib/SRC/DEPRECATED/sgeqrs.f @@ -16,7 +16,7 @@ * .. * .. Array Arguments .. * REAL A( LDA, * ), B( LDB, * ), TAU( * ), -* $ WORK( LWORK ) +* $ WORK( * ) * .. * * @@ -128,7 +128,7 @@ SUBROUTINE SGEQRS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, * .. * .. Array Arguments .. REAL A( LDA, * ), B( LDB, * ), TAU( * ), - $ WORK( LWORK ) + $ WORK( * ) * .. * * ===================================================================== diff --git a/lapack-netlib/SRC/DEPRECATED/zgelqs.f b/lapack-netlib/SRC/DEPRECATED/zgelqs.f index 5f629f8c7e..772165dfd9 100644 --- a/lapack-netlib/SRC/DEPRECATED/zgelqs.f +++ b/lapack-netlib/SRC/DEPRECATED/zgelqs.f @@ -16,7 +16,7 @@ * .. * .. Array Arguments .. * COMPLEX*16 A( LDA, * ), B( LDB, * ), TAU( * ), -* $ WORK( LWORK ) +* $ WORK( * ) * .. * * @@ -128,7 +128,7 @@ SUBROUTINE ZGELQS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, * .. * .. Array Arguments .. COMPLEX*16 A( LDA, * ), B( LDB, * ), TAU( * ), - $ WORK( LWORK ) + $ WORK( * ) * .. * * ===================================================================== diff --git a/lapack-netlib/SRC/DEPRECATED/zgeqrs.f b/lapack-netlib/SRC/DEPRECATED/zgeqrs.f index 6583e38591..cc33a45fc1 100644 --- a/lapack-netlib/SRC/DEPRECATED/zgeqrs.f +++ b/lapack-netlib/SRC/DEPRECATED/zgeqrs.f @@ -16,7 +16,7 @@ * .. * .. Array Arguments .. * COMPLEX*16 A( LDA, * ), B( LDB, * ), TAU( * ), -* $ WORK( LWORK ) +* $ WORK( * ) * .. * * @@ -128,7 +128,7 @@ SUBROUTINE ZGEQRS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, * .. * .. Array Arguments .. COMPLEX*16 A( LDA, * ), B( LDB, * ), TAU( * ), - $ WORK( LWORK ) + $ WORK( * ) * .. * * =====================================================================