Skip to content

Commit 0e4c25b

Browse files
hjjvandamhjjvandam
committed
Adding extra space for extended API
1 parent 37e4061 commit 0e4c25b

25 files changed

+229
-171
lines changed

‎SRC/DEPRECATED/cggsvp.f

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -401,8 +401,8 @@ SUBROUTINE CGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
401401
*
402402
* Update A := A*Z**H
403403
*
404-
CALL CUNMR2( 'Right', 'Conjugate transpose', M, N, L, B, LDB,
405-
$ TAU, A, LDA, WORK, INFO )
404+
CALL CUNMR2( 'Right', 'Conjugate transpose', M, N, L, B,
405+
$ LDB, TAU, A, LDA, WORK, INFO )
406406
IF( WANTQ ) THEN
407407
*
408408
* Update Q := Q*Z**H
@@ -489,8 +489,8 @@ SUBROUTINE CGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
489489
*
490490
* Update Q( 1:N,1:N-L ) = Q( 1:N,1:N-L )*Z1**H
491491
*
492-
CALL CUNMR2( 'Right', 'Conjugate transpose', N, N-L, K, A,
493-
$ LDA, TAU, Q, LDQ, WORK, INFO )
492+
CALL CUNMR2( 'Right', 'Conjugate transpose', N, N-L, K,
493+
$ A, LDA, TAU, Q, LDQ, WORK, INFO )
494494
END IF
495495
*
496496
* Clean up A

‎SRC/DEPRECATED/dgegs.f

Lines changed: 10 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -358,7 +358,8 @@ SUBROUTINE DGEGS( JOBVSL, JOBVSR, N, A, LDA, B, LDB, ALPHAR,
358358
END IF
359359
*
360360
IF( ILASCL ) THEN
361-
CALL DLASCL( 'G', -1, -1, ANRM, ANRMTO, N, N, A, LDA, IINFO )
361+
CALL DLASCL( 'G', -1, -1, ANRM, ANRMTO, N, N, A, LDA,
362+
$ IINFO )
362363
IF( IINFO.NE.0 ) THEN
363364
INFO = N + 9
364365
RETURN
@@ -378,7 +379,8 @@ SUBROUTINE DGEGS( JOBVSL, JOBVSR, N, A, LDA, B, LDB, ALPHAR,
378379
END IF
379380
*
380381
IF( ILBSCL ) THEN
381-
CALL DLASCL( 'G', -1, -1, BNRM, BNRMTO, N, N, B, LDB, IINFO )
382+
CALL DLASCL( 'G', -1, -1, BNRM, BNRMTO, N, N, B, LDB,
383+
$ IINFO )
382384
IF( IINFO.NE.0 ) THEN
383385
INFO = N + 9
384386
RETURN
@@ -496,7 +498,8 @@ SUBROUTINE DGEGS( JOBVSL, JOBVSR, N, A, LDA, B, LDB, ALPHAR,
496498
* Undo scaling
497499
*
498500
IF( ILASCL ) THEN
499-
CALL DLASCL( 'H', -1, -1, ANRMTO, ANRM, N, N, A, LDA, IINFO )
501+
CALL DLASCL( 'H', -1, -1, ANRMTO, ANRM, N, N, A, LDA,
502+
$ IINFO )
500503
IF( IINFO.NE.0 ) THEN
501504
INFO = N + 9
502505
RETURN
@@ -516,12 +519,14 @@ SUBROUTINE DGEGS( JOBVSL, JOBVSR, N, A, LDA, B, LDB, ALPHAR,
516519
END IF
517520
*
518521
IF( ILBSCL ) THEN
519-
CALL DLASCL( 'U', -1, -1, BNRMTO, BNRM, N, N, B, LDB, IINFO )
522+
CALL DLASCL( 'U', -1, -1, BNRMTO, BNRM, N, N, B, LDB,
523+
$ IINFO )
520524
IF( IINFO.NE.0 ) THEN
521525
INFO = N + 9
522526
RETURN
523527
END IF
524-
CALL DLASCL( 'G', -1, -1, BNRMTO, BNRM, N, 1, BETA, N, IINFO )
528+
CALL DLASCL( 'G', -1, -1, BNRMTO, BNRM, N, 1, BETA, N,
529+
$ IINFO )
525530
IF( IINFO.NE.0 ) THEN
526531
INFO = N + 9
527532
RETURN

‎SRC/DEPRECATED/dgegv.f

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -301,8 +301,9 @@
301301
*> \endverbatim
302302
*>
303303
* =====================================================================
304-
SUBROUTINE DGEGV( JOBVL, JOBVR, N, A, LDA, B, LDB, ALPHAR, ALPHAI,
305-
$ BETA, VL, LDVL, VR, LDVR, WORK, LWORK, INFO )
304+
SUBROUTINE DGEGV( JOBVL, JOBVR, N, A, LDA, B, LDB, ALPHAR,
305+
$ ALPHAI, BETA, VL, LDVL, VR, LDVR, WORK, LWORK,
306+
$ INFO )
306307
*
307308
* -- LAPACK driver routine --
308309
* -- LAPACK is a software package provided by Univ. of Tennessee, --

‎SRC/DEPRECATED/dgelqs.f

Lines changed: 6 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -174,18 +174,19 @@ SUBROUTINE DGELQS( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK,
174174
*
175175
* Solve L*X = B(1:m,:)
176176
*
177-
CALL DTRSM( 'Left', 'Lower', 'No transpose', 'Non-unit', M, NRHS,
178-
$ ONE, A, LDA, B, LDB )
177+
CALL DTRSM( 'Left', 'Lower', 'No transpose', 'Non-unit', M,
178+
$ NRHS, ONE, A, LDA, B, LDB )
179179
*
180180
* Set B(m+1:n,:) to zero
181181
*
182182
IF( M.LT.N )
183-
$ CALL DLASET( 'Full', N-M, NRHS, ZERO, ZERO, B( M+1, 1 ), LDB )
183+
$ CALL DLASET( 'Full', N-M, NRHS, ZERO, ZERO, B( M+1, 1 ),
184+
$ LDB )
184185
*
185186
* B := Q' * B
186187
*
187-
CALL DORMLQ( 'Left', 'Transpose', N, NRHS, M, A, LDA, TAU, B, LDB,
188-
$ WORK, LWORK, INFO )
188+
CALL DORMLQ( 'Left', 'Transpose', N, NRHS, M, A, LDA, TAU, B,
189+
$ LDB, WORK, LWORK, INFO )
189190
*
190191
RETURN
191192
*

‎SRC/DEPRECATED/dgelsx.f

Lines changed: 18 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -173,8 +173,8 @@
173173
*> \ingroup doubleGEsolve
174174
*
175175
* =====================================================================
176-
SUBROUTINE DGELSX( M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK,
177-
$ WORK, INFO )
176+
SUBROUTINE DGELSX( M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND,
177+
$ RANK, WORK, INFO )
178178
*
179179
* -- LAPACK driver routine --
180180
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -283,20 +283,23 @@ SUBROUTINE DGELSX( M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK,
283283
*
284284
* Scale matrix norm up to SMLNUM
285285
*
286-
CALL DLASCL( 'G', 0, 0, BNRM, SMLNUM, M, NRHS, B, LDB, INFO )
286+
CALL DLASCL( 'G', 0, 0, BNRM, SMLNUM, M, NRHS, B, LDB,
287+
$ INFO )
287288
IBSCL = 1
288289
ELSE IF( BNRM.GT.BIGNUM ) THEN
289290
*
290291
* Scale matrix norm down to BIGNUM
291292
*
292-
CALL DLASCL( 'G', 0, 0, BNRM, BIGNUM, M, NRHS, B, LDB, INFO )
293+
CALL DLASCL( 'G', 0, 0, BNRM, BIGNUM, M, NRHS, B, LDB,
294+
$ INFO )
293295
IBSCL = 2
294296
END IF
295297
*
296298
* Compute QR factorization with column pivoting of A:
297299
* A * P = Q * R
298300
*
299-
CALL DGEQPF( M, N, A, LDA, JPVT, WORK( 1 ), WORK( MN+1 ), INFO )
301+
CALL DGEQPF( M, N, A, LDA, JPVT, WORK( 1 ), WORK( MN+1 ),
302+
$ INFO )
300303
*
301304
* workspace 3*N. Details of Householder rotations stored
302305
* in WORK(1:MN).
@@ -350,8 +353,8 @@ SUBROUTINE DGELSX( M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK,
350353
*
351354
* B(1:M,1:NRHS) := Q**T * B(1:M,1:NRHS)
352355
*
353-
CALL DORM2R( 'Left', 'Transpose', M, NRHS, MN, A, LDA, WORK( 1 ),
354-
$ B, LDB, WORK( 2*MN+1 ), INFO )
356+
CALL DORM2R( 'Left', 'Transpose', M, NRHS, MN, A, LDA,
357+
$ WORK( 1 ), B, LDB, WORK( 2*MN+1 ), INFO )
355358
*
356359
* workspace NRHS
357360
*
@@ -408,18 +411,22 @@ SUBROUTINE DGELSX( M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK,
408411
* Undo scaling
409412
*
410413
IF( IASCL.EQ.1 ) THEN
411-
CALL DLASCL( 'G', 0, 0, ANRM, SMLNUM, N, NRHS, B, LDB, INFO )
414+
CALL DLASCL( 'G', 0, 0, ANRM, SMLNUM, N, NRHS, B, LDB,
415+
$ INFO )
412416
CALL DLASCL( 'U', 0, 0, SMLNUM, ANRM, RANK, RANK, A, LDA,
413417
$ INFO )
414418
ELSE IF( IASCL.EQ.2 ) THEN
415-
CALL DLASCL( 'G', 0, 0, ANRM, BIGNUM, N, NRHS, B, LDB, INFO )
419+
CALL DLASCL( 'G', 0, 0, ANRM, BIGNUM, N, NRHS, B, LDB,
420+
$ INFO )
416421
CALL DLASCL( 'U', 0, 0, BIGNUM, ANRM, RANK, RANK, A, LDA,
417422
$ INFO )
418423
END IF
419424
IF( IBSCL.EQ.1 ) THEN
420-
CALL DLASCL( 'G', 0, 0, SMLNUM, BNRM, N, NRHS, B, LDB, INFO )
425+
CALL DLASCL( 'G', 0, 0, SMLNUM, BNRM, N, NRHS, B, LDB,
426+
$ INFO )
421427
ELSE IF( IBSCL.EQ.2 ) THEN
422-
CALL DLASCL( 'G', 0, 0, BIGNUM, BNRM, N, NRHS, B, LDB, INFO )
428+
CALL DLASCL( 'G', 0, 0, BIGNUM, BNRM, N, NRHS, B, LDB,
429+
$ INFO )
423430
END IF
424431
*
425432
100 CONTINUE

‎SRC/DEPRECATED/dgeqpf.f

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -253,7 +253,8 @@ SUBROUTINE DGEQPF( M, N, A, LDA, JPVT, TAU, WORK, INFO )
253253
* Generate elementary reflector H(i)
254254
*
255255
IF( I.LT.M ) THEN
256-
CALL DLARFG( M-I+1, A( I, I ), A( I+1, I ), 1, TAU( I ) )
256+
CALL DLARFG( M-I+1, A( I, I ), A( I+1, I ), 1,
257+
$ TAU( I ) )
257258
ELSE
258259
CALL DLARFG( 1, A( M, M ), A( M, M ), 1, TAU( M ) )
259260
END IF
@@ -264,8 +265,8 @@ SUBROUTINE DGEQPF( M, N, A, LDA, JPVT, TAU, WORK, INFO )
264265
*
265266
AII = A( I, I )
266267
A( I, I ) = ONE
267-
CALL DLARF( 'LEFT', M-I+1, N-I, A( I, I ), 1, TAU( I ),
268-
$ A( I, I+1 ), LDA, WORK( 2*N+1 ) )
268+
CALL DLARF( 'LEFT', M-I+1, N-I, A( I, I ), 1,
269+
$ TAU( I ), A( I, I+1 ), LDA, WORK( 2*N+1 ) )
269270
A( I, I ) = AII
270271
END IF
271272
*

‎SRC/DEPRECATED/dggsvp.f

Lines changed: 13 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -392,8 +392,8 @@ SUBROUTINE DGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
392392
*
393393
* Update Q := Q*Z**T
394394
*
395-
CALL DORMR2( 'Right', 'Transpose', N, N, L, B, LDB, TAU, Q,
396-
$ LDQ, WORK, INFO )
395+
CALL DORMR2( 'Right', 'Transpose', N, N, L, B, LDB, TAU,
396+
$ Q, LDQ, WORK, INFO )
397397
END IF
398398
*
399399
* Clean up B
@@ -439,9 +439,10 @@ SUBROUTINE DGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
439439
*
440440
CALL DLASET( 'Full', M, M, ZERO, ZERO, U, LDU )
441441
IF( M.GT.1 )
442-
$ CALL DLACPY( 'Lower', M-1, N-L, A( 2, 1 ), LDA, U( 2, 1 ),
443-
$ LDU )
444-
CALL DORG2R( M, M, MIN( M, N-L ), U, LDU, TAU, WORK, INFO )
442+
$ CALL DLACPY( 'Lower', M-1, N-L, A( 2, 1 ), LDA,
443+
$ U( 2, 1 ), LDU )
444+
CALL DORG2R( M, M, MIN( M, N-L ), U, LDU, TAU, WORK,
445+
$ INFO )
445446
END IF
446447
*
447448
IF( WANTQ ) THEN
@@ -460,7 +461,8 @@ SUBROUTINE DGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
460461
90 CONTINUE
461462
100 CONTINUE
462463
IF( M.GT.K )
463-
$ CALL DLASET( 'Full', M-K, N-L, ZERO, ZERO, A( K+1, 1 ), LDA )
464+
$ CALL DLASET( 'Full', M-K, N-L, ZERO, ZERO, A( K+1, 1 ),
465+
$ LDA )
464466
*
465467
IF( N-L.GT.K ) THEN
466468
*
@@ -472,8 +474,8 @@ SUBROUTINE DGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
472474
*
473475
* Update Q( 1:N,1:N-L ) = Q( 1:N,1:N-L )*Z1**T
474476
*
475-
CALL DORMR2( 'Right', 'Transpose', N, N-L, K, A, LDA, TAU,
476-
$ Q, LDQ, WORK, INFO )
477+
CALL DORMR2( 'Right', 'Transpose', N, N-L, K, A, LDA,
478+
$ TAU, Q, LDQ, WORK, INFO )
477479
END IF
478480
*
479481
* Clean up A
@@ -497,9 +499,9 @@ SUBROUTINE DGGSVP( JOBU, JOBV, JOBQ, M, P, N, A, LDA, B, LDB,
497499
*
498500
* Update U(:,K+1:M) := U(:,K+1:M)*U1
499501
*
500-
CALL DORM2R( 'Right', 'No transpose', M, M-K, MIN( M-K, L ),
501-
$ A( K+1, N-L+1 ), LDA, TAU, U( 1, K+1 ), LDU,
502-
$ WORK, INFO )
502+
CALL DORM2R( 'Right', 'No transpose', M, M-K,
503+
$ MIN( M-K, L ), A( K+1, N-L+1 ), LDA, TAU,
504+
$ U( 1, K+1 ), LDU, WORK, INFO )
503505
END IF
504506
*
505507
* Clean up

‎SRC/DEPRECATED/dlahrd.f

Lines changed: 11 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -221,8 +221,8 @@ SUBROUTINE DLAHRD( N, K, NB, A, LDA, TAU, T, LDT, Y, LDY )
221221
* w := V1**T * b1
222222
*
223223
CALL DCOPY( I-1, A( K+1, I ), 1, T( 1, NB ), 1 )
224-
CALL DTRMV( 'Lower', 'Transpose', 'Unit', I-1, A( K+1, 1 ),
225-
$ LDA, T( 1, NB ), 1 )
224+
CALL DTRMV( 'Lower', 'Transpose', 'Unit', I-1,
225+
$ A( K+1, 1 ), LDA, T( 1, NB ), 1 )
226226
*
227227
* w := w + V2**T *b2
228228
*
@@ -236,8 +236,9 @@ SUBROUTINE DLAHRD( N, K, NB, A, LDA, TAU, T, LDT, Y, LDY )
236236
*
237237
* b2 := b2 - V2*w
238238
*
239-
CALL DGEMV( 'No transpose', N-K-I+1, I-1, -ONE, A( K+I, 1 ),
240-
$ LDA, T( 1, NB ), 1, ONE, A( K+I, I ), 1 )
239+
CALL DGEMV( 'No transpose', N-K-I+1, I-1, -ONE,
240+
$ A( K+I, 1 ), LDA, T( 1, NB ), 1, ONE,
241+
$ A( K+I, I ), 1 )
241242
*
242243
* b1 := b1 - V1*w
243244
*
@@ -258,12 +259,12 @@ SUBROUTINE DLAHRD( N, K, NB, A, LDA, TAU, T, LDT, Y, LDY )
258259
*
259260
* Compute Y(1:n,i)
260261
*
261-
CALL DGEMV( 'No transpose', N, N-K-I+1, ONE, A( 1, I+1 ), LDA,
262-
$ A( K+I, I ), 1, ZERO, Y( 1, I ), 1 )
263-
CALL DGEMV( 'Transpose', N-K-I+1, I-1, ONE, A( K+I, 1 ), LDA,
264-
$ A( K+I, I ), 1, ZERO, T( 1, I ), 1 )
265-
CALL DGEMV( 'No transpose', N, I-1, -ONE, Y, LDY, T( 1, I ), 1,
266-
$ ONE, Y( 1, I ), 1 )
262+
CALL DGEMV( 'No transpose', N, N-K-I+1, ONE, A( 1, I+1 ),
263+
$ LDA, A( K+I, I ), 1, ZERO, Y( 1, I ), 1 )
264+
CALL DGEMV( 'Transpose', N-K-I+1, I-1, ONE, A( K+I, 1 ),
265+
$ LDA, A( K+I, I ), 1, ZERO, T( 1, I ), 1 )
266+
CALL DGEMV( 'No transpose', N, I-1, -ONE, Y, LDY, T( 1, I ),
267+
$ 1, ONE, Y( 1, I ), 1 )
267268
CALL DSCAL( N, TAU( I ), Y( 1, I ), 1 )
268269
*
269270
* Compute T(1:i,i)

‎SRC/DEPRECATED/dlatzm.f

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -147,7 +147,8 @@
147147
*> \ingroup doubleOTHERcomputational
148148
*
149149
* =====================================================================
150-
SUBROUTINE DLATZM( SIDE, M, N, V, INCV, TAU, C1, C2, LDC, WORK )
150+
SUBROUTINE DLATZM( SIDE, M, N, V, INCV, TAU, C1, C2, LDC,
151+
$ WORK )
151152
*
152153
* -- LAPACK computational routine --
153154
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -202,8 +203,8 @@ SUBROUTINE DLATZM( SIDE, M, N, V, INCV, TAU, C1, C2, LDC, WORK )
202203
* w := C1 + C2 * v
203204
*
204205
CALL DCOPY( M, C1, 1, WORK, 1 )
205-
CALL DGEMV( 'No transpose', M, N-1, ONE, C2, LDC, V, INCV, ONE,
206-
$ WORK, 1 )
206+
CALL DGEMV( 'No transpose', M, N-1, ONE, C2, LDC, V, INCV,
207+
$ ONE, WORK, 1 )
207208
*
208209
* [ C1, C2 ] := [ C1, C2 ] - tau* w * [ 1 , v**T]
209210
*

‎SRC/DEPRECATED/dtzrqf.f

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -194,7 +194,8 @@ SUBROUTINE DTZRQF( M, N, A, LDA, TAU, INFO )
194194
* Use a Householder reflection to zero the kth row of A.
195195
* First set up the reflection.
196196
*
197-
CALL DLARFG( N-M+1, A( K, K ), A( K, M1 ), LDA, TAU( K ) )
197+
CALL DLARFG( N-M+1, A( K, K ), A( K, M1 ), LDA,
198+
$ TAU( K ) )
198199
*
199200
IF( ( TAU( K ).NE.ZERO ) .AND. ( K.GT.1 ) ) THEN
200201
*
@@ -216,8 +217,8 @@ SUBROUTINE DTZRQF( M, N, A, LDA, TAU, INFO )
216217
* and B := B - tau*w*z( k )**T.
217218
*
218219
CALL DAXPY( K-1, -TAU( K ), TAU, 1, A( 1, K ), 1 )
219-
CALL DGER( K-1, N-M, -TAU( K ), TAU, 1, A( K, M1 ), LDA,
220-
$ A( 1, M1 ), LDA )
220+
CALL DGER( K-1, N-M, -TAU( K ), TAU, 1, A( K, M1 ),
221+
$ LDA, A( 1, M1 ), LDA )
221222
END IF
222223
20 CONTINUE
223224
END IF

0 commit comments

Comments
 (0)