The R Project SVN R

Rev

Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
83470 ripley 1 
--- orig/dlatrs.f	2022-11-11 17:49:54
2 
+++ dlatrs.f	2022-12-16 14:50:10
3 
@@ -259,7 +259,7 @@
4 
       LOGICAL            NOTRAN, NOUNIT, UPPER
5 
       INTEGER            I, IMAX, J, JFIRST, JINC, JLAST
6 
       DOUBLE PRECISION   BIGNUM, GROW, REC, SMLNUM, SUMJ, TJJ, TJJS,
7 
-     $                   TMAX, TSCAL, USCAL, XBND, XJ, XMAX
8 
+     $                   TMAX, TSCAL, USCAL, XBND, XJ, XMAX, WK(1)
9 
 *     ..
10 
 *     .. External Functions ..
11 
       LOGICAL            LSAME
12 
@@ -361,17 +361,19 @@
13 
 *
14 
 *              A is upper triangular.
15 
 *
16 
+*     R change: SUMJ -> WK to avoid rank mismatch.
17 
                DO J = 2, N
18 
-                  TMAX = MAX( DLANGE( 'M', J-1, 1, A( 1, J ), 1, SUMJ ),
19 
+                  TMAX = MAX( DLANGE( 'M', J-1, 1, A( 1, J ), 1, WK ),
20 
      $                        TMAX )
21 
                END DO
22 
             ELSE
23 
 *
24 
 *              A is lower triangular.
25 
 *
26 
+*     R change: SUMJ -> WK to avoid rank mismatch.
27 
                DO J = 1, N - 1
28 
                   TMAX = MAX( DLANGE( 'M', N-J, 1, A( J+1, J ), 1,
29 
-     $                        SUMJ ), TMAX )
30 
+     $                        WK ), TMAX )
31 
                END DO
32 
             END IF
33 
 *