Attachment 916539 Details for Bug 1117491 – Fix issue with zero value handling in BLAS functions

[patch] Fix issue with zero value handling in BLAS functions

lapack-3.5.0-R-blas-fixes.patch (text/plain), 3.64 KB, created by Tom "spot" Callaway on 2014-07-08 19:32:11 UTC

(hide)

Description:

Filename:

MIME Type:

Creator: Tom "spot" Callaway

Created: 2014-07-08 19:32:11 UTC

Size: 3.64 KB

patch

obsolete

>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

Actions: View | Diff

Attachments on bug 1117491: 916539