Skip to content

Commit a1e8cb4

Browse files
mkrainiukmkrainiuk
committed
LAPACK SRC long lines updated by script
1 parent b891145 commit a1e8cb4

File tree

1,414 files changed

+14437
-7213
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

1,414 files changed

+14437
-7213
lines changed

‎SRC/cbbcsd.f

Lines changed: 24 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -326,7 +326,8 @@
326326
*> \ingroup bbcsd
327327
*
328328
* =====================================================================
329-
SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
329+
SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P,
330+
$ Q,
330331
$ THETA, PHI, U1, LDU1, U2, LDU2, V1T, LDV1T,
331332
$ V2T, LDV2T, B11D, B11E, B12D, B12E, B21D, B21E,
332333
$ B22D, B22E, RWORK, LRWORK, INFO )
@@ -373,7 +374,8 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
373374
$ UNFL, X1, X2, Y1, Y2
374375
*
375376
* .. External Subroutines ..
376-
EXTERNAL CLASR, CSCAL, CSWAP, SLARTGP, SLARTGS, SLAS2,
377+
EXTERNAL CLASR, CSCAL, CSWAP, SLARTGP, SLARTGS,
378+
$ SLAS2,
377379
$ XERBLA
378380
* ..
379381
* .. External Functions ..
@@ -560,9 +562,11 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
560562
*
561563
* Compute shifts for B11 and B21 and use the lesser
562564
*
563-
CALL SLAS2( B11D(IMAX-1), B11E(IMAX-1), B11D(IMAX), SIGMA11,
565+
CALL SLAS2( B11D(IMAX-1), B11E(IMAX-1), B11D(IMAX),
566+
$ SIGMA11,
564567
$ DUMMY )
565-
CALL SLAS2( B21D(IMAX-1), B21E(IMAX-1), B21D(IMAX), SIGMA21,
568+
CALL SLAS2( B21D(IMAX-1), B21E(IMAX-1), B21D(IMAX),
569+
$ SIGMA21,
566570
$ DUMMY )
567571
*
568572
IF( SIGMA11 .LE. SIGMA21 ) THEN
@@ -719,10 +723,12 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
719723
CALL SLARTGP( Y2, Y1, RWORK(IV2TSN+I-1-1),
720724
$ RWORK(IV2TCS+I-1-1), R )
721725
ELSE IF( .NOT. RESTART12 .AND. RESTART22 ) THEN
722-
CALL SLARTGP( B12BULGE, B12D(I-1), RWORK(IV2TSN+I-1-1),
726+
CALL SLARTGP( B12BULGE, B12D(I-1),
727+
$ RWORK(IV2TSN+I-1-1),
723728
$ RWORK(IV2TCS+I-1-1), R )
724729
ELSE IF( RESTART12 .AND. .NOT. RESTART22 ) THEN
725-
CALL SLARTGP( B22BULGE, B22D(I-1), RWORK(IV2TSN+I-1-1),
730+
CALL SLARTGP( B22BULGE, B22D(I-1),
731+
$ RWORK(IV2TSN+I-1-1),
726732
$ RWORK(IV2TCS+I-1-1), R )
727733
ELSE IF( NU .LT. MU ) THEN
728734
CALL SLARTGS( B12E(I-1), B12D(I), NU,
@@ -781,7 +787,8 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
781787
* chasing by applying the original shift again.
782788
*
783789
IF( .NOT. RESTART11 .AND. .NOT. RESTART12 ) THEN
784-
CALL SLARTGP( X2, X1, RWORK(IU1SN+I-1), RWORK(IU1CS+I-1),
790+
CALL SLARTGP( X2, X1, RWORK(IU1SN+I-1),
791+
$ RWORK(IU1CS+I-1),
785792
$ R )
786793
ELSE IF( .NOT. RESTART11 .AND. RESTART12 ) THEN
787794
CALL SLARTGP( B11BULGE, B11D(I), RWORK(IU1SN+I-1),
@@ -790,14 +797,16 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
790797
CALL SLARTGP( B12BULGE, B12E(I-1), RWORK(IU1SN+I-1),
791798
$ RWORK(IU1CS+I-1), R )
792799
ELSE IF( MU .LE. NU ) THEN
793-
CALL SLARTGS( B11E(I), B11D(I+1), MU, RWORK(IU1CS+I-1),
800+
CALL SLARTGS( B11E(I), B11D(I+1), MU,
801+
$ RWORK(IU1CS+I-1),
794802
$ RWORK(IU1SN+I-1) )
795803
ELSE
796804
CALL SLARTGS( B12D(I), B12E(I), NU, RWORK(IU1CS+I-1),
797805
$ RWORK(IU1SN+I-1) )
798806
END IF
799807
IF( .NOT. RESTART21 .AND. .NOT. RESTART22 ) THEN
800-
CALL SLARTGP( Y2, Y1, RWORK(IU2SN+I-1), RWORK(IU2CS+I-1),
808+
CALL SLARTGP( Y2, Y1, RWORK(IU2SN+I-1),
809+
$ RWORK(IU2CS+I-1),
801810
$ R )
802811
ELSE IF( .NOT. RESTART21 .AND. RESTART22 ) THEN
803812
CALL SLARTGP( B21BULGE, B21D(I), RWORK(IU2SN+I-1),
@@ -806,7 +815,8 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
806815
CALL SLARTGP( B22BULGE, B22E(I-1), RWORK(IU2SN+I-1),
807816
$ RWORK(IU2CS+I-1), R )
808817
ELSE IF( NU .LT. MU ) THEN
809-
CALL SLARTGS( B21E(I), B21E(I+1), NU, RWORK(IU2CS+I-1),
818+
CALL SLARTGS( B21E(I), B21E(I+1), NU,
819+
$ RWORK(IU2CS+I-1),
810820
$ RWORK(IU2SN+I-1) )
811821
ELSE
812822
CALL SLARTGS( B22D(I), B22E(I), MU, RWORK(IU2CS+I-1),
@@ -992,7 +1002,8 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
9921002
IF( B12D(IMAX)+B22D(IMAX) .LT. 0 ) THEN
9931003
IF( WANTV2T ) THEN
9941004
IF( COLMAJOR ) THEN
995-
CALL CSCAL( M-Q, NEGONECOMPLEX, V2T(IMAX,1), LDV2T )
1005+
CALL CSCAL( M-Q, NEGONECOMPLEX, V2T(IMAX,1),
1006+
$ LDV2T )
9961007
ELSE
9971008
CALL CSCAL( M-Q, NEGONECOMPLEX, V2T(1,IMAX), 1 )
9981009
END IF
@@ -1059,7 +1070,8 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P, Q,
10591070
IF( WANTU2 )
10601071
$ CALL CSWAP( M-P, U2(1,I), 1, U2(1,MINI), 1 )
10611072
IF( WANTV1T )
1062-
$ CALL CSWAP( Q, V1T(I,1), LDV1T, V1T(MINI,1), LDV1T )
1073+
$ CALL CSWAP( Q, V1T(I,1), LDV1T, V1T(MINI,1),
1074+
$ LDV1T )
10631075
IF( WANTV2T )
10641076
$ CALL CSWAP( M-Q, V2T(I,1), LDV2T, V2T(MINI,1),
10651077
$ LDV2T )

‎SRC/cbdsqr.f

Lines changed: 26 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -280,7 +280,8 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
280280
EXTERNAL LSAME, SLAMCH
281281
* ..
282282
* .. External Subroutines ..
283-
EXTERNAL CLASR, CSROT, CSSCAL, CSWAP, SLARTG, SLAS2,
283+
EXTERNAL CLASR, CSROT, CSSCAL, CSWAP, SLARTG,
284+
$ SLAS2,
284285
$ SLASQ1, SLASV2, XERBLA
285286
* ..
286287
* .. Intrinsic Functions ..
@@ -361,10 +362,12 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
361362
* Update singular vectors if desired
362363
*
363364
IF( NRU.GT.0 )
364-
$ CALL CLASR( 'R', 'V', 'F', NRU, N, RWORK( 1 ), RWORK( N ),
365+
$ CALL CLASR( 'R', 'V', 'F', NRU, N, RWORK( 1 ),
366+
$ RWORK( N ),
365367
$ U, LDU )
366368
IF( NCC.GT.0 )
367-
$ CALL CLASR( 'L', 'V', 'F', N, NCC, RWORK( 1 ), RWORK( N ),
369+
$ CALL CLASR( 'L', 'V', 'F', N, NCC, RWORK( 1 ),
370+
$ RWORK( N ),
368371
$ C, LDC )
369372
END IF
370373
*
@@ -488,7 +491,8 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
488491
$ CALL CSROT( NCVT, VT( M-1, 1 ), LDVT, VT( M, 1 ), LDVT,
489492
$ COSR, SINR )
490493
IF( NRU.GT.0 )
491-
$ CALL CSROT( NRU, U( 1, M-1 ), 1, U( 1, M ), 1, COSL, SINL )
494+
$ CALL CSROT( NRU, U( 1, M-1 ), 1, U( 1, M ), 1, COSL,
495+
$ SINL )
492496
IF( NCC.GT.0 )
493497
$ CALL CSROT( NCC, C( M-1, 1 ), LDC, C( M, 1 ), LDC, COSL,
494498
$ SINL )
@@ -621,7 +625,8 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
621625
CALL SLARTG( D( I )*CS, E( I ), CS, SN, R )
622626
IF( I.GT.LL )
623627
$ E( I-1 ) = OLDSN*R
624-
CALL SLARTG( OLDCS*R, D( I+1 )*SN, OLDCS, OLDSN, D( I ) )
628+
CALL SLARTG( OLDCS*R, D( I+1 )*SN, OLDCS, OLDSN,
629+
$ D( I ) )
625630
RWORK( I-LL+1 ) = CS
626631
RWORK( I-LL+1+NM1 ) = SN
627632
RWORK( I-LL+1+NM12 ) = OLDCS
@@ -637,10 +642,12 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
637642
$ CALL CLASR( 'L', 'V', 'F', M-LL+1, NCVT, RWORK( 1 ),
638643
$ RWORK( N ), VT( LL, 1 ), LDVT )
639644
IF( NRU.GT.0 )
640-
$ CALL CLASR( 'R', 'V', 'F', NRU, M-LL+1, RWORK( NM12+1 ),
645+
$ CALL CLASR( 'R', 'V', 'F', NRU, M-LL+1,
646+
$ RWORK( NM12+1 ),
641647
$ RWORK( NM13+1 ), U( 1, LL ), LDU )
642648
IF( NCC.GT.0 )
643-
$ CALL CLASR( 'L', 'V', 'F', M-LL+1, NCC, RWORK( NM12+1 ),
649+
$ CALL CLASR( 'L', 'V', 'F', M-LL+1, NCC,
650+
$ RWORK( NM12+1 ),
644651
$ RWORK( NM13+1 ), C( LL, 1 ), LDC )
645652
*
646653
* Test convergence
@@ -659,7 +666,8 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
659666
CALL SLARTG( D( I )*CS, E( I-1 ), CS, SN, R )
660667
IF( I.LT.M )
661668
$ E( I ) = OLDSN*R
662-
CALL SLARTG( OLDCS*R, D( I-1 )*SN, OLDCS, OLDSN, D( I ) )
669+
CALL SLARTG( OLDCS*R, D( I-1 )*SN, OLDCS, OLDSN,
670+
$ D( I ) )
663671
RWORK( I-LL ) = CS
664672
RWORK( I-LL+NM1 ) = -SN
665673
RWORK( I-LL+NM12 ) = OLDCS
@@ -672,7 +680,8 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
672680
* Update singular vectors
673681
*
674682
IF( NCVT.GT.0 )
675-
$ CALL CLASR( 'L', 'V', 'B', M-LL+1, NCVT, RWORK( NM12+1 ),
683+
$ CALL CLASR( 'L', 'V', 'B', M-LL+1, NCVT,
684+
$ RWORK( NM12+1 ),
676685
$ RWORK( NM13+1 ), VT( LL, 1 ), LDVT )
677686
IF( NRU.GT.0 )
678687
$ CALL CLASR( 'R', 'V', 'B', NRU, M-LL+1, RWORK( 1 ),
@@ -727,10 +736,12 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
727736
$ CALL CLASR( 'L', 'V', 'F', M-LL+1, NCVT, RWORK( 1 ),
728737
$ RWORK( N ), VT( LL, 1 ), LDVT )
729738
IF( NRU.GT.0 )
730-
$ CALL CLASR( 'R', 'V', 'F', NRU, M-LL+1, RWORK( NM12+1 ),
739+
$ CALL CLASR( 'R', 'V', 'F', NRU, M-LL+1,
740+
$ RWORK( NM12+1 ),
731741
$ RWORK( NM13+1 ), U( 1, LL ), LDU )
732742
IF( NCC.GT.0 )
733-
$ CALL CLASR( 'L', 'V', 'F', M-LL+1, NCC, RWORK( NM12+1 ),
743+
$ CALL CLASR( 'L', 'V', 'F', M-LL+1, NCC,
744+
$ RWORK( NM12+1 ),
734745
$ RWORK( NM13+1 ), C( LL, 1 ), LDC )
735746
*
736747
* Test convergence
@@ -777,7 +788,8 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
777788
* Update singular vectors if desired
778789
*
779790
IF( NCVT.GT.0 )
780-
$ CALL CLASR( 'L', 'V', 'B', M-LL+1, NCVT, RWORK( NM12+1 ),
791+
$ CALL CLASR( 'L', 'V', 'B', M-LL+1, NCVT,
792+
$ RWORK( NM12+1 ),
781793
$ RWORK( NM13+1 ), VT( LL, 1 ), LDVT )
782794
IF( NRU.GT.0 )
783795
$ CALL CLASR( 'R', 'V', 'B', NRU, M-LL+1, RWORK( 1 ),
@@ -833,7 +845,8 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
833845
IF( NRU.GT.0 )
834846
$ CALL CSWAP( NRU, U( 1, ISUB ), 1, U( 1, N+1-I ), 1 )
835847
IF( NCC.GT.0 )
836-
$ CALL CSWAP( NCC, C( ISUB, 1 ), LDC, C( N+1-I, 1 ), LDC )
848+
$ CALL CSWAP( NCC, C( ISUB, 1 ), LDC, C( N+1-I, 1 ),
849+
$ LDC )
837850
END IF
838851
190 CONTINUE
839852
GO TO 220

‎SRC/cgbbrd.f

Lines changed: 9 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -223,7 +223,8 @@ SUBROUTINE CGBBRD( VECT, M, N, NCC, KL, KU, AB, LDAB, D, E, Q,
223223
COMPLEX RA, RB, RS, T
224224
* ..
225225
* .. External Subroutines ..
226-
EXTERNAL CLARGV, CLARTG, CLARTV, CLASET, CROT, CSCAL,
226+
EXTERNAL CLARGV, CLARTG, CLARTV, CLASET, CROT,
227+
$ CSCAL,
227228
$ XERBLA
228229
* ..
229230
* .. Intrinsic Functions ..
@@ -243,7 +244,9 @@ SUBROUTINE CGBBRD( VECT, M, N, NCC, KL, KU, AB, LDAB, D, E, Q,
243244
WANTC = NCC.GT.0
244245
KLU1 = KL + KU + 1
245246
INFO = 0
246-
IF( .NOT.WANTQ .AND. .NOT.WANTPT .AND. .NOT.LSAME( VECT, 'N' ) )
247+
IF( .NOT.WANTQ .AND.
248+
$ .NOT.WANTPT .AND.
249+
$ .NOT.LSAME( VECT, 'N' ) )
247250
$ THEN
248251
INFO = -1
249252
ELSE IF( M.LT.0 ) THEN
@@ -339,7 +342,8 @@ SUBROUTINE CGBBRD( VECT, M, N, NCC, KL, KU, AB, LDAB, D, E, Q,
339342
NRT = NR
340343
END IF
341344
IF( NRT.GT.0 )
342-
$ CALL CLARTV( NRT, AB( KLU1-L, J1-KLM+L-1 ), INCA,
345+
$ CALL CLARTV( NRT, AB( KLU1-L, J1-KLM+L-1 ),
346+
$ INCA,
343347
$ AB( KLU1-L+1, J1-KLM+L-1 ), INCA,
344348
$ RWORK( J1 ), WORK( J1 ), KB1 )
345349
10 CONTINUE
@@ -378,7 +382,8 @@ SUBROUTINE CGBBRD( VECT, M, N, NCC, KL, KU, AB, LDAB, D, E, Q,
378382
* apply plane rotations to C
379383
*
380384
DO 30 J = J1, J2, KB1
381-
CALL CROT( NCC, C( J-1, 1 ), LDC, C( J, 1 ), LDC,
385+
CALL CROT( NCC, C( J-1, 1 ), LDC, C( J, 1 ),
386+
$ LDC,
382387
$ RWORK( J ), WORK( J ) )
383388
30 CONTINUE
384389
END IF

‎SRC/cgbcon.f

Lines changed: 10 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -143,7 +143,8 @@
143143
*> \ingroup gbcon
144144
*
145145
* =====================================================================
146-
SUBROUTINE CGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND,
146+
SUBROUTINE CGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM,
147+
$ RCOND,
147148
$ WORK, RWORK, INFO )
148149
*
149150
* -- LAPACK computational routine --
@@ -185,7 +186,8 @@ SUBROUTINE CGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND,
185186
EXTERNAL LSAME, ICAMAX, SLAMCH, CDOTC
186187
* ..
187188
* .. External Subroutines ..
188-
EXTERNAL CAXPY, CLACN2, CLATBS, CSRSCL, XERBLA
189+
EXTERNAL CAXPY, CLACN2, CLATBS, CSRSCL,
190+
$ XERBLA
189191
* ..
190192
* .. Intrinsic Functions ..
191193
INTRINSIC ABS, AIMAG, MIN, REAL
@@ -260,13 +262,15 @@ SUBROUTINE CGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND,
260262
WORK( JP ) = WORK( J )
261263
WORK( J ) = T
262264
END IF
263-
CALL CAXPY( LM, -T, AB( KD+1, J ), 1, WORK( J+1 ), 1 )
265+
CALL CAXPY( LM, -T, AB( KD+1, J ), 1, WORK( J+1 ),
266+
$ 1 )
264267
20 CONTINUE
265268
END IF
266269
*
267270
* Multiply by inv(U).
268271
*
269-
CALL CLATBS( 'Upper', 'No transpose', 'Non-unit', NORMIN, N,
272+
CALL CLATBS( 'Upper', 'No transpose', 'Non-unit', NORMIN,
273+
$ N,
270274
$ KL+KU, AB, LDAB, WORK, SCALE, RWORK, INFO )
271275
ELSE
272276
*
@@ -281,7 +285,8 @@ SUBROUTINE CGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND,
281285
IF( LNOTI ) THEN
282286
DO 30 J = N - 1, 1, -1
283287
LM = MIN( KL, N-J )
284-
WORK( J ) = WORK( J ) - CDOTC( LM, AB( KD+1, J ), 1,
288+
WORK( J ) = WORK( J ) - CDOTC( LM, AB( KD+1, J ),
289+
$ 1,
285290
$ WORK( J+1 ), 1 )
286291
JP = IPIV( J )
287292
IF( JP.NE.J ) THEN

‎SRC/cgbequ.f

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -150,7 +150,8 @@
150150
*> \ingroup gbequ
151151
*
152152
* =====================================================================
153-
SUBROUTINE CGBEQU( M, N, KL, KU, AB, LDAB, R, C, ROWCND, COLCND,
153+
SUBROUTINE CGBEQU( M, N, KL, KU, AB, LDAB, R, C, ROWCND,
154+
$ COLCND,
154155
$ AMAX, INFO )
155156
*
156157
* -- LAPACK computational routine --

‎SRC/cgbequb.f

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -157,7 +157,8 @@
157157
*> \ingroup gbequb
158158
*
159159
* =====================================================================
160-
SUBROUTINE CGBEQUB( M, N, KL, KU, AB, LDAB, R, C, ROWCND, COLCND,
160+
SUBROUTINE CGBEQUB( M, N, KL, KU, AB, LDAB, R, C, ROWCND,
161+
$ COLCND,
161162
$ AMAX, INFO )
162163
*
163164
* -- LAPACK computational routine --

‎SRC/cgbrfs.f

Lines changed: 8 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -201,7 +201,8 @@
201201
*> \ingroup gbrfs
202202
*
203203
* =====================================================================
204-
SUBROUTINE CGBRFS( TRANS, N, KL, KU, NRHS, AB, LDAB, AFB, LDAFB,
204+
SUBROUTINE CGBRFS( TRANS, N, KL, KU, NRHS, AB, LDAB, AFB,
205+
$ LDAFB,
205206
$ IPIV, B, LDB, X, LDX, FERR, BERR, WORK, RWORK,
206207
$ INFO )
207208
*
@@ -245,7 +246,8 @@ SUBROUTINE CGBRFS( TRANS, N, KL, KU, NRHS, AB, LDAB, AFB, LDAFB,
245246
INTEGER ISAVE( 3 )
246247
* ..
247248
* .. External Subroutines ..
248-
EXTERNAL CAXPY, CCOPY, CGBMV, CGBTRS, CLACN2, XERBLA
249+
EXTERNAL CAXPY, CCOPY, CGBMV, CGBTRS, CLACN2,
250+
$ XERBLA
249251
* ..
250252
* .. Intrinsic Functions ..
251253
INTRINSIC ABS, AIMAG, MAX, MIN, REAL
@@ -332,7 +334,8 @@ SUBROUTINE CGBRFS( TRANS, N, KL, KU, NRHS, AB, LDAB, AFB, LDAFB,
332334
* where op(A) = A, A**T, or A**H, depending on TRANS.
333335
*
334336
CALL CCOPY( N, B( 1, J ), 1, WORK, 1 )
335-
CALL CGBMV( TRANS, N, N, KL, KU, -CONE, AB, LDAB, X( 1, J ), 1,
337+
CALL CGBMV( TRANS, N, N, KL, KU, -CONE, AB, LDAB, X( 1, J ),
338+
$ 1,
336339
$ CONE, WORK, 1 )
337340
*
338341
* Compute componentwise relative backward error from formula
@@ -390,7 +393,8 @@ SUBROUTINE CGBRFS( TRANS, N, KL, KU, NRHS, AB, LDAB, AFB, LDAFB,
390393
*
391394
* Update solution and try again.
392395
*
393-
CALL CGBTRS( TRANS, N, KL, KU, 1, AFB, LDAFB, IPIV, WORK, N,
396+
CALL CGBTRS( TRANS, N, KL, KU, 1, AFB, LDAFB, IPIV, WORK,
397+
$ N,
394398
$ INFO )
395399
CALL CAXPY( N, CONE, WORK, 1, X( 1, J ), 1 )
396400
LSTRES = BERR( J )

0 commit comments

Comments
 (0)