100 lines
3.6 KiB
Diff
100 lines
3.6 KiB
Diff
|
diff -up lapack-3.5.0/BLAS/SRC/dgbmv.f.R lapack-3.5.0/BLAS/SRC/dgbmv.f
|
||
|
|
--- lapack-3.5.0/BLAS/SRC/dgbmv.f.R 2011-11-03 16:32:56.000000000 -0400
|
||
|
|
+++ lapack-3.5.0/BLAS/SRC/dgbmv.f 2014-07-07 11:50:48.257175853 -0400
|
||
|
|
@@ -312,18 +312,18 @@
|
||
|
|
JX = KX
|
||
|
|
IF (INCY.EQ.1) THEN
|
||
|
|
DO 60 J = 1,N
|
||
|
|
- IF (X(JX).NE.ZERO) THEN
|
||
|
|
+c IF (X(JX).NE.ZERO) THEN
|
||
|
|
TEMP = ALPHA*X(JX)
|
||
|
|
K = KUP1 - J
|
||
|
|
DO 50 I = MAX(1,J-KU),MIN(M,J+KL)
|
||
|
|
Y(I) = Y(I) + TEMP*A(K+I,J)
|
||
|
|
50 CONTINUE
|
||
|
|
- END IF
|
||
|
|
+c END IF
|
||
|
|
JX = JX + INCX
|
||
|
|
60 CONTINUE
|
||
|
|
ELSE
|
||
|
|
DO 80 J = 1,N
|
||
|
|
- IF (X(JX).NE.ZERO) THEN
|
||
|
|
+c IF (X(JX).NE.ZERO) THEN
|
||
|
|
TEMP = ALPHA*X(JX)
|
||
|
|
IY = KY
|
||
|
|
K = KUP1 - J
|
||
|
|
@@ -331,7 +331,7 @@
|
||
|
|
Y(IY) = Y(IY) + TEMP*A(K+I,J)
|
||
|
|
IY = IY + INCY
|
||
|
|
70 CONTINUE
|
||
|
|
- END IF
|
||
|
|
+c END IF
|
||
|
|
JX = JX + INCX
|
||
|
|
IF (J.GT.KU) KY = KY + INCY
|
||
|
|
80 CONTINUE
|
||
|
|
diff -up lapack-3.5.0/BLAS/SRC/dgemm.f.R lapack-3.5.0/BLAS/SRC/dgemm.f
|
||
|
|
--- lapack-3.5.0/BLAS/SRC/dgemm.f.R 2011-11-03 16:32:56.000000000 -0400
|
||
|
|
+++ lapack-3.5.0/BLAS/SRC/dgemm.f 2014-07-07 11:50:48.257175853 -0400
|
||
|
|
@@ -311,12 +311,12 @@
|
||
|
|
60 CONTINUE
|
||
|
|
END IF
|
||
|
|
DO 80 L = 1,K
|
||
|
|
- IF (B(L,J).NE.ZERO) THEN
|
||
|
|
+c IF (B(L,J).NE.ZERO) THEN
|
||
|
|
TEMP = ALPHA*B(L,J)
|
||
|
|
DO 70 I = 1,M
|
||
|
|
C(I,J) = C(I,J) + TEMP*A(I,L)
|
||
|
|
70 CONTINUE
|
||
|
|
- END IF
|
||
|
|
+c END IF
|
||
|
|
80 CONTINUE
|
||
|
|
90 CONTINUE
|
||
|
|
ELSE
|
||
|
|
@@ -353,12 +353,12 @@
|
||
|
|
140 CONTINUE
|
||
|
|
END IF
|
||
|
|
DO 160 L = 1,K
|
||
|
|
- IF (B(J,L).NE.ZERO) THEN
|
||
|
|
+c IF (B(J,L).NE.ZERO) THEN
|
||
|
|
TEMP = ALPHA*B(J,L)
|
||
|
|
DO 150 I = 1,M
|
||
|
|
C(I,J) = C(I,J) + TEMP*A(I,L)
|
||
|
|
150 CONTINUE
|
||
|
|
- END IF
|
||
|
|
+c END IF
|
||
|
|
160 CONTINUE
|
||
|
|
170 CONTINUE
|
||
|
|
ELSE
|
||
|
|
diff -up lapack-3.5.0/BLAS/SRC/dgemv.f.R lapack-3.5.0/BLAS/SRC/dgemv.f
|
||
|
|
--- lapack-3.5.0/BLAS/SRC/dgemv.f.R 2011-11-03 16:32:56.000000000 -0400
|
||
|
|
+++ lapack-3.5.0/BLAS/SRC/dgemv.f 2014-07-07 11:50:48.257175853 -0400
|
||
|
|
@@ -278,24 +278,24 @@
|
||
|
|
JX = KX
|
||
|
|
IF (INCY.EQ.1) THEN
|
||
|
|
DO 60 J = 1,N
|
||
|
|
- IF (X(JX).NE.ZERO) THEN
|
||
|
|
+c IF (X(JX).NE.ZERO) THEN
|
||
|
|
TEMP = ALPHA*X(JX)
|
||
|
|
DO 50 I = 1,M
|
||
|
|
Y(I) = Y(I) + TEMP*A(I,J)
|
||
|
|
50 CONTINUE
|
||
|
|
- END IF
|
||
|
|
+c END IF
|
||
|
|
JX = JX + INCX
|
||
|
|
60 CONTINUE
|
||
|
|
ELSE
|
||
|
|
DO 80 J = 1,N
|
||
|
|
- IF (X(JX).NE.ZERO) THEN
|
||
|
|
+c IF (X(JX).NE.ZERO) THEN
|
||
|
|
TEMP = ALPHA*X(JX)
|
||
|
|
IY = KY
|
||
|
|
DO 70 I = 1,M
|
||
|
|
Y(IY) = Y(IY) + TEMP*A(I,J)
|
||
|
|
IY = IY + INCY
|
||
|
|
70 CONTINUE
|
||
|
|
- END IF
|
||
|
|
+c END IF
|
||
|
|
JX = JX + INCX
|
||
|
|
80 CONTINUE
|
||
|
|
END IF
|