OSchip
/
OpenBLAS
Not watched
1
0
Fork 0
Code
Releases 66 Wiki evaluate Activity
Issues 0 Pull Requests 0 Datasets Model Cloudbrain HPC
Browse Source

Merge pull request #39 from xianyi/develop 

rebase
tags/v0.3.10^2
Martin Kroeker GitHub 6 years ago
parent
c0da205412 b6a948fbee
commit
a986d42ea6
No known key found for this signature in database GPG Key ID: 4AEE18F83AFDEB23
9 changed files with 425 additions and 343 deletions
Split View
Diff Options
Show Stats Download Patch File Download Diff File
  1. +4
    -0
      Makefile
  2. +4
    -1
      cmake/fc.cmake
  3. +3
    -3
      cmake/system_check.cmake
  4. +7
    -0
      ctest/c_dblas1.c
  5. +194
    -170
      ctest/c_dblat1.f
  6. +7
    -0
      ctest/c_sblas1.c
  7. +202
    -168
      ctest/c_sblat1.f
  8. +1
    -1
      getarch.c
  9. +3
    -0
      interface/gemm.c

+ 4
- 0
Makefile View File

@@ -84,6 +84,10 @@ else
@echo " (Multi-threading; Max num-threads is $(NUM_THREADS))"
endif

ifeq ($(DYNAMIC_ARCH), 1)
@echo " Supporting multiple $(ARCH) cpu models with minimum requirement for the common code being $(CORE)"
endif

ifeq ($(USE_OPENMP), 1)
@echo
@echo " Use OpenMP in the multithreading. Because of ignoring OPENBLAS_NUM_THREADS and GOTO_NUM_THREADS flags, "


+ 4
- 1
cmake/fc.cmake View File

@@ -76,8 +76,11 @@ if (${F_COMPILER} STREQUAL "GFORTRAN")
endif ()
endif ()

if (${F_COMPILER} STREQUAL "INTEL")
if (${F_COMPILER} STREQUAL "IFORT")
set(CCOMMON_OPT "${CCOMMON_OPT} -DF_INTERFACE_INTEL")
if (MSVC)
set(FCOMMON_OPT "${FCOMMON_OPT} -names:lowercase -assume:underscore")
endif ()
if (INTERFACE64)
set(FCOMMON_OPT "${FCOMMON_OPT} -i8")
endif ()


+ 3
- 3
cmake/system_check.cmake View File

@@ -109,10 +109,10 @@ else()
endif()

if (X86_64 OR X86)
file(WRITE ${PROJECT_BINARY_DIR}/avx512.tmp "#include <immintrin.h>\n\nint main(void){ __asm__ volatile(\"vbroadcastss -4 * 4(%rsi), %zmm2\"); }")
execute_process(COMMAND ${CMAKE_C_COMPILER} -march=skylake-avx512 -c -v -o ${PROJECT_BINARY_DIR}/avx512.o -x c ${PROJECT_BINARY_DIR}/avx512.tmp OUTPUT_QUIET ERROR_QUIET RESULT_VARIABLE NO_AVX512)
file(WRITE ${PROJECT_BINARY_DIR}/avx512.c "#include <immintrin.h>\n\nint main(void){ __asm__ volatile(\"vbroadcastss -4 * 4(%rsi), %zmm2\"); }")
execute_process(COMMAND ${CMAKE_C_COMPILER} -march=skylake-avx512 -c -v -o ${PROJECT_BINARY_DIR}/avx512.o ${PROJECT_BINARY_DIR}/avx512.c OUTPUT_QUIET ERROR_QUIET RESULT_VARIABLE NO_AVX512)
if (NO_AVX512 EQUAL 1)
set (CCOMMON_OPT "${CCOMMON_OPT} -DNO_AVX512")
endif()
file(REMOVE "avx512.tmp" "avx512.o")
file(REMOVE "avx512.c" "avx512.o")
endif()

+ 7
- 0
ctest/c_dblas1.c View File

@@ -53,6 +53,13 @@ void F77_drot( const int *N, double *X, const int *incX, double *Y,
return;
}

void F77_drotm(const int *N, double *X, const int *incX, double *Y,
const int *incY, const double *dparam)
{
cblas_drotm(*N, X, *incX, Y, *incY, dparam);
return;
}

void F77_dscal(const int *N, const double *alpha, double *X,
const int *incX)
{


+ 194
- 170
ctest/c_dblat1.f View File

@@ -19,7 +19,7 @@
DATA SFAC/9.765625D-4/
* .. Executable Statements ..
WRITE (NOUT,99999)
DO 20 IC = 1, 10
DO 20 IC = 1, 11
ICASE = IC
CALL HEADER
*
@@ -40,7 +40,7 @@
ELSE IF (ICASE.EQ.1 .OR. ICASE.EQ.2 .OR. ICASE.EQ.5 .OR.
+ ICASE.EQ.6) THEN
CALL CHECK2(SFAC)
ELSE IF (ICASE.EQ.4) THEN
ELSE IF (ICASE.EQ.4 .OR. ICASE.EQ.11) THEN
CALL CHECK3(SFAC)
END IF
* -- Print
@@ -59,7 +59,7 @@
INTEGER ICASE, INCX, INCY, MODE, N
LOGICAL PASS
* .. Local Arrays ..
CHARACTER*15 L(10)
CHARACTER*15 L(11)
* .. Common blocks ..
COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS
* .. Data statements ..
@@ -73,6 +73,7 @@
DATA L(8)/'CBLAS_DASUM '/
DATA L(9)/'CBLAS_DSCAL '/
DATA L(10)/'CBLAS_IDAMAX'/
DATA L(11)/'CBLAS_DROTM'/
* .. Executable Statements ..
WRITE (NOUT,99999) ICASE, L(ICASE)
RETURN
@@ -400,199 +401,81 @@
LOGICAL PASS
* .. Local Scalars ..
DOUBLE PRECISION SC, SS
INTEGER I, K, KI, KN, KSIZE, LENX, LENY, MX, MY
INTEGER I, KI, KN, KSIZE, LEN
* .. Local Arrays ..
DOUBLE PRECISION COPYX(5), COPYY(5), DT9X(7,4,4), DT9Y(7,4,4),
+ DX1(7), DY1(7), MWPC(11), MWPS(11), MWPSTX(5),
+ MWPSTY(5), MWPTX(11,5), MWPTY(11,5), MWPX(5),
+ MWPY(5), SSIZE2(14,2), STX(7), STY(7), SX(7),
+ SY(7)
INTEGER INCXS(4), INCYS(4), LENS(4,2), MWPINX(11),
+ MWPINY(11), MWPN(11), NS(4)
DOUBLE PRECISION DX(10), DY(10), SSIZE2(10,2), STX(10),
+ STY(10), SX(10), SY(10),
+ PARAM(5, 4), DPARAM(5)
INTEGER INCXS(7), INCYS(7), NS(5)
* .. External Subroutines ..
EXTERNAL STEST,DROTTEST
EXTERNAL STEST, DROTTEST, DROT
* .. Intrinsic Functions ..
INTRINSIC ABS, MIN
INTRINSIC MIN
* .. Common blocks ..
COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS
* .. Data statements ..
DATA INCXS/1, 2, -2, -1/
DATA INCYS/1, -2, 1, -2/
DATA LENS/1, 1, 2, 4, 1, 1, 3, 7/
DATA NS/0, 1, 2, 4/
DATA DX1/0.6D0, 0.1D0, -0.5D0, 0.8D0, 0.9D0, -0.3D0,
+ -0.4D0/
DATA DY1/0.5D0, -0.9D0, 0.3D0, 0.7D0, -0.6D0, 0.2D0,
+ 0.8D0/
DATA INCXS/1, 1, 2, 2, -2, -1, -2/
DATA INCYS/1, 2, 2, -2, 1, -2, -2/
DATA NS/0, 1, 2, 4, 5/
DATA DX/0.6D0, 0.1D0, -0.5D0, 0.8D0, 0.9D0, -0.3D0,
+ -0.4D0, 0.7D0, 0.5D0, 0.2D0/
DATA DY/0.5D0, -0.9D0, 0.3D0, 0.7D0, -0.6D0, 0.2D0,
+ 0.8D0, -0.5D0, 0.1D0, -0.3D0/
DATA SC, SS/0.8D0, 0.6D0/
DATA DT9X/0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.78D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.78D0, -0.46D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.78D0, -0.46D0, -0.22D0,
+ 1.06D0, 0.0D0, 0.0D0, 0.0D0, 0.6D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.78D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.66D0, 0.1D0, -0.1D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.96D0, 0.1D0, -0.76D0, 0.8D0, 0.90D0,
+ -0.3D0, -0.02D0, 0.6D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.78D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.0D0, -0.06D0, 0.1D0,
+ -0.1D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.90D0,
+ 0.1D0, -0.22D0, 0.8D0, 0.18D0, -0.3D0, -0.02D0,
+ 0.6D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.78D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.78D0, 0.26D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.78D0, 0.26D0, -0.76D0, 1.12D0,
+ 0.0D0, 0.0D0, 0.0D0/
DATA DT9Y/0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.04D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.04D0, -0.78D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.04D0, -0.78D0, 0.54D0,
+ 0.08D0, 0.0D0, 0.0D0, 0.0D0, 0.5D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.04D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.7D0,
+ -0.9D0, -0.12D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.64D0, -0.9D0, -0.30D0, 0.7D0, -0.18D0, 0.2D0,
+ 0.28D0, 0.5D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.04D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.7D0, -1.08D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.64D0, -1.26D0,
+ 0.54D0, 0.20D0, 0.0D0, 0.0D0, 0.0D0, 0.5D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.04D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.04D0, -0.9D0, 0.18D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.04D0, -0.9D0, 0.18D0, 0.7D0,
+ -0.18D0, 0.2D0, 0.16D0/
DATA LEN/10/
DATA PARAM/-2.0D0, 1.0D0, 0.0D0, 0.0D0, 1.0D0,
+ -1.0D0, 0.2D0, 0.3D0, 0.4D0, 0.5D0,
+ 0.0D0, 1.0D0, 0.3D0, 0.4D0, 1.0D0,
+ 1.0D0, 0.2D0, -1.0D0, 1.0D0, 0.5D0/
DATA SSIZE2/0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 0.0D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0,
+ 0.0D0, 0.0D0, 0.0D0, 0.0D0, 1.17D0, 1.17D0,
+ 1.17D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0, 1.17D0,
+ 1.17D0, 1.17D0, 1.17D0/
+ 1.17D0, 1.17D0/
* .. Executable Statements ..
*
DO 60 KI = 1, 4
DO 60 KI = 1, 7
INCX = INCXS(KI)
INCY = INCYS(KI)
MX = ABS(INCX)
MY = ABS(INCY)
*
DO 40 KN = 1, 4
DO 40 KN = 1, 5
N = NS(KN)
KSIZE = MIN(2,KN)
LENX = LENS(KN,MX)
LENY = LENS(KN,MY)
*
IF (ICASE.EQ.4) THEN
* .. DROTTEST ..
DO 20 I = 1, 7
SX(I) = DX1(I)
SY(I) = DY1(I)
STX(I) = DT9X(I,KN,KI)
STY(I) = DT9Y(I,KN,KI)
DO 20 I = 1, 10
SX(I) = DX(I)
SY(I) = DY(I)
STX(I) = DX(I)
STY(I) = DY(I)
20 CONTINUE
CALL DROTTEST(N,SX,INCX,SY,INCY,SC,SS)
CALL STEST(LENX,SX,STX,SSIZE2(1,KSIZE),SFAC)
CALL STEST(LENY,SY,STY,SSIZE2(1,KSIZE),SFAC)
ELSE
CALL DROT(N,STX,INCX,STY,INCY,SC,SS)
CALL STEST(LEN,SX,STX,SSIZE2(1,KSIZE),SFAC)
CALL STEST(LEN,SY,STY,SSIZE2(1,KSIZE),SFAC)
ELSE IF (ICASE.EQ.11) THEN
* .. DROTMTEST ..
DO 90 I = 1, 10
SX(I) = DX(I)
SY(I) = DY(I)
STX(I) = DX(I)
STY(I) = DY(I)
90 CONTINUE
DO 70 I = 1, 4
DO 80 K = 1, 5
DPARAM(K) = PARAM(K,I)
80 CONTINUE
CALL DROTMTEST(N,SX,INCX,SY,INCY,DPARAM)
CALL DROTM(N,STX,INCX,STY,INCY,DPARAM)
CALL STEST(LEN,SX,STX,SSIZE2(1,KSIZE),SFAC)
CALL STEST(LEN,SY,STY,SSIZE2(1,KSIZE),SFAC)
70 CONTINUE
ELSE
WRITE (NOUT,*) ' Shouldn''t be here in CHECK3'
STOP
END IF
40 CONTINUE
60 CONTINUE
*
MWPC(1) = 1
DO 80 I = 2, 11
MWPC(I) = 0
80 CONTINUE
MWPS(1) = 0.0
DO 100 I = 2, 6
MWPS(I) = 1.0
100 CONTINUE
DO 120 I = 7, 11
MWPS(I) = -1.0
120 CONTINUE
MWPINX(1) = 1
MWPINX(2) = 1
MWPINX(3) = 1
MWPINX(4) = -1
MWPINX(5) = 1
MWPINX(6) = -1
MWPINX(7) = 1
MWPINX(8) = 1
MWPINX(9) = -1
MWPINX(10) = 1
MWPINX(11) = -1
MWPINY(1) = 1
MWPINY(2) = 1
MWPINY(3) = -1
MWPINY(4) = -1
MWPINY(5) = 2
MWPINY(6) = 1
MWPINY(7) = 1
MWPINY(8) = -1
MWPINY(9) = -1
MWPINY(10) = 2
MWPINY(11) = 1
DO 140 I = 1, 11
MWPN(I) = 5
140 CONTINUE
MWPN(5) = 3
MWPN(10) = 3
DO 160 I = 1, 5
MWPX(I) = I
MWPY(I) = I
MWPTX(1,I) = I
MWPTY(1,I) = I
MWPTX(2,I) = I
MWPTY(2,I) = -I
MWPTX(3,I) = 6 - I
MWPTY(3,I) = I - 6
MWPTX(4,I) = I
MWPTY(4,I) = -I
MWPTX(6,I) = 6 - I
MWPTY(6,I) = I - 6
MWPTX(7,I) = -I
MWPTY(7,I) = I
MWPTX(8,I) = I - 6
MWPTY(8,I) = 6 - I
MWPTX(9,I) = -I
MWPTY(9,I) = I
MWPTX(11,I) = I - 6
MWPTY(11,I) = 6 - I
160 CONTINUE
MWPTX(5,1) = 1
MWPTX(5,2) = 3
MWPTX(5,3) = 5
MWPTX(5,4) = 4
MWPTX(5,5) = 5
MWPTY(5,1) = -1
MWPTY(5,2) = 2
MWPTY(5,3) = -2
MWPTY(5,4) = 4
MWPTY(5,5) = -3
MWPTX(10,1) = -1
MWPTX(10,2) = -3
MWPTX(10,3) = -5
MWPTX(10,4) = 4
MWPTX(10,5) = 5
MWPTY(10,1) = 1
MWPTY(10,2) = 2
MWPTY(10,3) = 2
MWPTY(10,4) = 4
MWPTY(10,5) = 3
DO 200 I = 1, 11
INCX = MWPINX(I)
INCY = MWPINY(I)
DO 180 K = 1, 5
COPYX(K) = MWPX(K)
COPYY(K) = MWPY(K)
MWPSTX(K) = MWPTX(I,K)
MWPSTY(K) = MWPTY(I,K)
180 CONTINUE
CALL DROTTEST(MWPN(I),COPYX,INCX,COPYY,INCY,MWPC(I),MWPS(I))
CALL STEST(5,COPYX,MWPSTX,MWPSTX,SFAC)
CALL STEST(5,COPYY,MWPSTY,MWPSTY,SFAC)
200 CONTINUE
RETURN
END
SUBROUTINE STEST(LEN,SCOMP,STRUE,SSIZE,SFAC)
@@ -726,3 +609,144 @@
+ /1X)
99997 FORMAT (1X,I4,I3,3I5,2I36,I12)
END
SUBROUTINE DROT(N,DX,INCX,DY,INCY,C,S)
* .. Scalar Arguments ..
DOUBLE PRECISION C,S
INTEGER INCX,INCY,N
* ..
* .. Array Arguments ..
DOUBLE PRECISION DX(*),DY(*)
* ..
* applies a plane rotation.
* jack dongarra, linpack, 3/11/78.
* modified 12/3/93, array(1) declarations changed to array(*)
*
* .. Local Scalars ..
DOUBLE PRECISION DTEMP
INTEGER I,IX,IY
* ..
IF (N.LE.0) RETURN
IF (INCX.EQ.1 .AND. INCY.EQ.1) GO TO 20
IX = 1
IY = 1
IF (INCX.LT.0) IX = (-N+1)*INCX + 1
IF (INCY.LT.0) IY = (-N+1)*INCY + 1
DO 10 I = 1,N
DTEMP = C*DX(IX) + S*DY(IY)
DY(IY) = C*DY(IY) - S*DX(IX)
DX(IX) = DTEMP
IX = IX + INCX
IY = IY + INCY
10 CONTINUE
RETURN
20 DO 30 I = 1,N
DTEMP = C*DX(I) + S*DY(I)
DY(I) = C*DY(I) - S*DX(I)
DX(I) = DTEMP
30 CONTINUE
RETURN
END
SUBROUTINE drotm(N,DX,INCX,DY,INCY,DPARAM)
*
* -- Reference BLAS level1 routine (version 3.8.0) --
* -- Reference BLAS is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* November 2017
*
* .. Scalar Arguments ..
INTEGER INCX,INCY,N
* ..
* .. Array Arguments ..
DOUBLE PRECISION DPARAM(5),DX(*),DY(*)
* ..
*
* =====================================================================
*
* .. Local Scalars ..
DOUBLE PRECISION DFLAG,DH11,DH12,DH21,DH22,TWO,W,Z,ZERO
INTEGER I,KX,KY,NSTEPS
* ..
* .. Data statements ..
DATA zero,two/0.d0,2.d0/
* ..
*
dflag = dparam(1)
IF (n.LE.0 .OR. (dflag+two.EQ.zero)) RETURN
IF (incx.EQ.incy.AND.incx.GT.0) THEN
*
nsteps = n*incx
IF (dflag.LT.zero) THEN
dh11 = dparam(2)
dh12 = dparam(4)
dh21 = dparam(3)
dh22 = dparam(5)
DO i = 1,nsteps,incx
w = dx(i)
z = dy(i)
dx(i) = w*dh11 + z*dh12
dy(i) = w*dh21 + z*dh22
END DO
ELSE IF (dflag.EQ.zero) THEN
dh12 = dparam(4)
dh21 = dparam(3)
DO i = 1,nsteps,incx
w = dx(i)
z = dy(i)
dx(i) = w + z*dh12
dy(i) = w*dh21 + z
END DO
ELSE
dh11 = dparam(2)
dh22 = dparam(5)
DO i = 1,nsteps,incx
w = dx(i)
z = dy(i)
dx(i) = w*dh11 + z
dy(i) = -w + dh22*z
END DO
END IF
ELSE
kx = 1
ky = 1
IF (incx.LT.0) kx = 1 + (1-n)*incx
IF (incy.LT.0) ky = 1 + (1-n)*incy
*
IF (dflag.LT.zero) THEN
dh11 = dparam(2)
dh12 = dparam(4)
dh21 = dparam(3)
dh22 = dparam(5)
DO i = 1,n
w = dx(kx)
z = dy(ky)
dx(kx) = w*dh11 + z*dh12
dy(ky) = w*dh21 + z*dh22
kx = kx + incx
ky = ky + incy
END DO
ELSE IF (dflag.EQ.zero) THEN
dh12 = dparam(4)
dh21 = dparam(3)
DO i = 1,n
w = dx(kx)
z = dy(ky)
dx(kx) = w + z*dh12
dy(ky) = w*dh21 + z
kx = kx + incx
ky = ky + incy
END DO
ELSE
dh11 = dparam(2)
dh22 = dparam(5)
DO i = 1,n
w = dx(kx)
z = dy(ky)
dx(kx) = w*dh11 + z
dy(ky) = -w + dh22*z
kx = kx + incx
ky = ky + incy
END DO
END IF
END IF
RETURN
END

+ 7
- 0
ctest/c_sblas1.c View File

@@ -55,6 +55,13 @@ void F77_srotg( float *a, float *b, float *c, float *s)
return;
}

void F77_srotm(blasint *N, float *X, blasint *incX, float *Y, blasint *incY,
float *param)
{
cblas_srotm(*N, X, *incX, Y, *incY, param);
return;
}

void F77_srot( blasint *N, float *X, blasint *incX, float *Y,
blasint *incY, const float *c, const float *s)
{


+ 202
- 168
ctest/c_sblat1.f View File

@@ -19,7 +19,7 @@
DATA SFAC/9.765625E-4/
* .. Executable Statements ..
WRITE (NOUT,99999)
DO 20 IC = 1, 10
DO 20 IC = 1, 11
ICASE = IC
CALL HEADER
*
@@ -40,7 +40,7 @@
ELSE IF (ICASE.EQ.1 .OR. ICASE.EQ.2 .OR. ICASE.EQ.5 .OR.
+ ICASE.EQ.6) THEN
CALL CHECK2(SFAC)
ELSE IF (ICASE.EQ.4) THEN
ELSE IF (ICASE.EQ.4 .OR. ICASE.EQ.11) THEN
CALL CHECK3(SFAC)
END IF
* -- Print
@@ -59,7 +59,7 @@
INTEGER ICASE, INCX, INCY, MODE, N
LOGICAL PASS
* .. Local Arrays ..
CHARACTER*15 L(10)
CHARACTER*15 L(11)
* .. Common blocks ..
COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS
* .. Data statements ..
@@ -73,6 +73,7 @@
DATA L(8)/'CBLAS_SASUM '/
DATA L(9)/'CBLAS_SSCAL '/
DATA L(10)/'CBLAS_ISAMAX'/
DATA L(11)/'CBLAS_SROTM'/
* .. Executable Statements ..
WRITE (NOUT,99999) ICASE, L(ICASE)
RETURN
@@ -396,203 +397,92 @@
* .. Scalar Arguments ..
REAL SFAC
* .. Scalars in Common ..
INTEGER ICASE, INCX, INCY, MODE, N
INTEGER ICASE, INCX, INCY, N
LOGICAL PASS
* .. Local Scalars ..
REAL SC, SS
INTEGER I, K, KI, KN, KSIZE, LENX, LENY, MX, MY
INTEGER I, K, KI, KN, KSIZE, LEN
* .. Local Arrays ..
REAL COPYX(5), COPYY(5), DT9X(7,4,4), DT9Y(7,4,4),
+ DX1(7), DY1(7), MWPC(11), MWPS(11), MWPSTX(5),
+ MWPSTY(5), MWPTX(11,5), MWPTY(11,5), MWPX(5),
+ MWPY(5), SSIZE2(14,2), STX(7), STY(7), SX(7),
+ SY(7)
INTEGER INCXS(4), INCYS(4), LENS(4,2), MWPINX(11),
+ MWPINY(11), MWPN(11), NS(4)
REAL DX(19), DY(19),
+ SSIZE2(19,2), STX(19), STY(19), SX(19), SY(19),
+ PARAM(5, 4), SPARAM(5)
INTEGER INCXS(7), INCYS(7), NS(7)
* .. External Subroutines ..
EXTERNAL SROTTEST, STEST
EXTERNAL SROTMTEST, SROTM
* .. Intrinsic Functions ..
INTRINSIC ABS, MIN
INTRINSIC MIN
* .. Common blocks ..
COMMON /COMBLA/ICASE, N, INCX, INCY, MODE, PASS
* .. Data statements ..
DATA INCXS/1, 2, -2, -1/
DATA INCYS/1, -2, 1, -2/
DATA LENS/1, 1, 2, 4, 1, 1, 3, 7/
DATA NS/0, 1, 2, 4/
DATA DX1/0.6E0, 0.1E0, -0.5E0, 0.8E0, 0.9E0, -0.3E0,
+ -0.4E0/
DATA DY1/0.5E0, -0.9E0, 0.3E0, 0.7E0, -0.6E0, 0.2E0,
DATA INCXS/1, 1, 2, 2, -2, -1, -2/
DATA INCYS/1, 2, 2, -2, 1, -2, -2/
DATA NS/0, 1, 2, 4, 5, 8, 9/
DATA DX/0.6E0, 0.1E0, -0.5E0, 0.8E0, 0.9E0, -0.3E0,
+ -0.4E0, 0.5E0, -0.9E0, 0.3E0, 0.7E0, -0.6E0,
+ 0.2E0, 0.8E0, -0.46E0, 0.78E0, -0.46E0, -0.22E0,
+ 1.06E0/
DATA DY/0.5E0, -0.9E0, 0.3E0, 0.7E0, -0.6E0, 0.2E0,
+ 0.6E0, 0.1E0, -0.5E0, 0.8E0, 0.9E0, -0.3E0,
+ 0.96E0, 0.1E0, -0.76E0, 0.8E0, 0.90E0, 0.66E0,
+ 0.8E0/
DATA SC, SS/0.8E0, 0.6E0/
DATA DT9X/0.6E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.78E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.78E0, -0.46E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.78E0, -0.46E0, -0.22E0,
+ 1.06E0, 0.0E0, 0.0E0, 0.0E0, 0.6E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.78E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.66E0, 0.1E0, -0.1E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.96E0, 0.1E0, -0.76E0, 0.8E0, 0.90E0,
+ -0.3E0, -0.02E0, 0.6E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.78E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.0E0, -0.06E0, 0.1E0,
+ -0.1E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.90E0,
+ 0.1E0, -0.22E0, 0.8E0, 0.18E0, -0.3E0, -0.02E0,
+ 0.6E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.78E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.78E0, 0.26E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.78E0, 0.26E0, -0.76E0, 1.12E0,
+ 0.0E0, 0.0E0, 0.0E0/
DATA DT9Y/0.5E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.04E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.04E0, -0.78E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.04E0, -0.78E0, 0.54E0,
+ 0.08E0, 0.0E0, 0.0E0, 0.0E0, 0.5E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.04E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.7E0,
+ -0.9E0, -0.12E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.64E0, -0.9E0, -0.30E0, 0.7E0, -0.18E0, 0.2E0,
+ 0.28E0, 0.5E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.04E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.7E0, -1.08E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.64E0, -1.26E0,
+ 0.54E0, 0.20E0, 0.0E0, 0.0E0, 0.0E0, 0.5E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.04E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.04E0, -0.9E0, 0.18E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.04E0, -0.9E0, 0.18E0, 0.7E0,
+ -0.18E0, 0.2E0, 0.16E0/
DATA PARAM/-2.0E0, 1.0E0, 0.0E0, 0.0E0, 1.0E0,
+ -1.0E0, 0.2E0, 0.3E0, 0.4E0, 0.5E0,
+ 0.0E0, 1.0E0, 0.3E0, 0.4E0, 1.0E0,
+ 1.0E0, 0.2E0, -1.0E0, 1.0E0, 0.5E0/
DATA LEN/19/
DATA SSIZE2/0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 0.0E0, 1.17E0, 1.17E0, 1.17E0, 1.17E0, 1.17E0,
+ 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0, 0.0E0,
+ 1.17E0, 1.17E0, 1.17E0, 1.17E0, 1.17E0, 1.17E0,
+ 1.17E0, 1.17E0, 1.17E0/
+ 1.17E0, 1.17E0, 1.17E0, 1.17E0, 1.17E0, 1.17E0,
+ 1.17E0, 1.17E0, 1.17E0, 1.17E0, 1.17E0, 1.17E0,
+ 1.17E0/
* .. Executable Statements ..
*
DO 60 KI = 1, 4
DO 60 KI = 1, 7
INCX = INCXS(KI)
INCY = INCYS(KI)
MX = ABS(INCX)
MY = ABS(INCY)
*
DO 40 KN = 1, 4
DO 40 KN = 1, 7
N = NS(KN)
KSIZE = MIN(2,KN)
LENX = LENS(KN,MX)
LENY = LENS(KN,MY)
*
IF (ICASE.EQ.4) THEN
* .. SROTTEST ..
DO 20 I = 1, 7
SX(I) = DX1(I)
SY(I) = DY1(I)
STX(I) = DT9X(I,KN,KI)
STY(I) = DT9Y(I,KN,KI)
DO 20 I = 1, 19
SX(I) = DX(I)
SY(I) = DY(I)
STX(I) = DX(I)
STY(I) = DY(I)
20 CONTINUE
CALL SROTTEST(N,SX,INCX,SY,INCY,SC,SS)
CALL STEST(LENX,SX,STX,SSIZE2(1,KSIZE),SFAC)
CALL STEST(LENY,SY,STY,SSIZE2(1,KSIZE),SFAC)
CALL SROT(N,STX,INCX,STY,INCY,SC,SS)
CALL STEST(LEN,SX,STX,SSIZE2(1,KSIZE),SFAC)
CALL STEST(LEN,SY,STY,SSIZE2(1,KSIZE),SFAC)
ELSE IF (ICASE.EQ.11) THEN
* .. SROTMTEST ..
DO 90 I = 1, 19
SX(I) = DX(I)
SY(I) = DY(I)
STX(I) = DX(I)
STY(I) = DY(I)
90 CONTINUE
DO 70 I = 1, 4
DO 80 K = 1, 5
SPARAM(K) = PARAM(K,I)
80 CONTINUE
CALL SROTMTEST(N,SX,INCX,SY,INCY,SPARAM)
CALL SROTM(N,STX,INCX,STY,INCY,SPARAM)
CALL STEST(LEN,SX,STX,SSIZE2(1,KSIZE),SFAC)
CALL STEST(LEN,SY,STY,SSIZE2(1,KSIZE),SFAC)
70 CONTINUE
ELSE
WRITE (NOUT,*) ' Shouldn''t be here in CHECK3'
STOP
END IF
40 CONTINUE
60 CONTINUE
*
MWPC(1) = 1
DO 80 I = 2, 11
MWPC(I) = 0
80 CONTINUE
MWPS(1) = 0
DO 100 I = 2, 6
MWPS(I) = 1
100 CONTINUE
DO 120 I = 7, 11
MWPS(I) = -1
120 CONTINUE
MWPINX(1) = 1
MWPINX(2) = 1
MWPINX(3) = 1
MWPINX(4) = -1
MWPINX(5) = 1
MWPINX(6) = -1
MWPINX(7) = 1
MWPINX(8) = 1
MWPINX(9) = -1
MWPINX(10) = 1
MWPINX(11) = -1
MWPINY(1) = 1
MWPINY(2) = 1
MWPINY(3) = -1
MWPINY(4) = -1
MWPINY(5) = 2
MWPINY(6) = 1
MWPINY(7) = 1
MWPINY(8) = -1
MWPINY(9) = -1
MWPINY(10) = 2
MWPINY(11) = 1
DO 140 I = 1, 11
MWPN(I) = 5
140 CONTINUE
MWPN(5) = 3
MWPN(10) = 3
DO 160 I = 1, 5
MWPX(I) = I
MWPY(I) = I
MWPTX(1,I) = I
MWPTY(1,I) = I
MWPTX(2,I) = I
MWPTY(2,I) = -I
MWPTX(3,I) = 6 - I
MWPTY(3,I) = I - 6
MWPTX(4,I) = I
MWPTY(4,I) = -I
MWPTX(6,I) = 6 - I
MWPTY(6,I) = I - 6
MWPTX(7,I) = -I
MWPTY(7,I) = I
MWPTX(8,I) = I - 6
MWPTY(8,I) = 6 - I
MWPTX(9,I) = -I
MWPTY(9,I) = I
MWPTX(11,I) = I - 6
MWPTY(11,I) = 6 - I
160 CONTINUE
MWPTX(5,1) = 1
MWPTX(5,2) = 3
MWPTX(5,3) = 5
MWPTX(5,4) = 4
MWPTX(5,5) = 5
MWPTY(5,1) = -1
MWPTY(5,2) = 2
MWPTY(5,3) = -2
MWPTY(5,4) = 4
MWPTY(5,5) = -3
MWPTX(10,1) = -1
MWPTX(10,2) = -3
MWPTX(10,3) = -5
MWPTX(10,4) = 4
MWPTX(10,5) = 5
MWPTY(10,1) = 1
MWPTY(10,2) = 2
MWPTY(10,3) = 2
MWPTY(10,4) = 4
MWPTY(10,5) = 3
DO 200 I = 1, 11
INCX = MWPINX(I)
INCY = MWPINY(I)
DO 180 K = 1, 5
COPYX(K) = MWPX(K)
COPYY(K) = MWPY(K)
MWPSTX(K) = MWPTX(I,K)
MWPSTY(K) = MWPTY(I,K)
180 CONTINUE
CALL SROTTEST(MWPN(I),COPYX,INCX,COPYY,INCY,MWPC(I),MWPS(I))
CALL STEST(5,COPYX,MWPSTX,MWPSTX,SFAC)
CALL STEST(5,COPYY,MWPSTY,MWPSTY,SFAC)
200 CONTINUE
RETURN
END
SUBROUTINE STEST(LEN,SCOMP,STRUE,SSIZE,SFAC)
@@ -726,3 +616,147 @@
+ /1X)
99997 FORMAT (1X,I4,I3,3I5,2I36,I12)
END
SUBROUTINE SROT(N,SX,INCX,SY,INCY,C,S)
*
* --Reference BLAS level1 routine (version 3.8.0) --
* --Reference BLAS is a software package provided by Univ. of Tennessee, --
* --Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* November 2017
*
* .. Scalar Arguments ..
REAL C,S
INTEGER INCX,INCY,N
* ..
* .. Array Arguments ..
REAL SX(*),SY(*)
* ..
* .. Local Scalars ..
REAL STEMP
INTEGER I,IX,IY
* ..
IF (n.LE.0) RETURN
IF (incx.EQ.1 .AND. incy.EQ.1) THEN
DO i = 1,n
stemp = c*sx(i) + s*sy(i)
sy(i) = c*sy(i) - s*sx(i)
sx(i) = stemp
END DO
ELSE
ix = 1
iy = 1
IF (incx.LT.0) ix = (-n+1)*incx + 1
IF (incy.LT.0) iy = (-n+1)*incy + 1
DO i = 1,n
stemp = c*sx(ix) + s*sy(iy)
sy(iy) = c*sy(iy) - s*sx(ix)
sx(ix) = stemp
ix = ix + incx
iy = iy + incy
END DO
END IF
RETURN
END
SUBROUTINE srotm(N,SX,INCX,SY,INCY,SPARAM)
*
* --Reference BLAS level1 routine (version 3.8.0) --
* --Reference BLAS is a software package provided by Univ. of Tennessee, --
* --Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* November 2017
*
* .. Scalar Arguments ..
INTEGER INCX,INCY,N
* ..
* .. Array Arguments ..
REAL SPARAM(5),SX(*),SY(*)
* ..
*
* ====================================================================
*
* .. Local Scalars ..
REAL SFLAG,SH11,SH12,SH21,SH22,TWO,W,Z,ZERO
INTEGER I,KX,KY,NSTEPS
* ..
* .. Data statements ..
DATA zero,two/0.e0,2.e0/
* ..
*
sflag = sparam(1)
IF (n.LE.0 .OR. (sflag+two.EQ.zero)) RETURN
IF (incx.EQ.incy.AND.incx.GT.0) THEN
*
nsteps = n*incx
IF (sflag.LT.zero) THEN
sh11 = sparam(2)
sh12 = sparam(4)
sh21 = sparam(3)
sh22 = sparam(5)
DO i = 1,nsteps,incx
w = sx(i)
z = sy(i)
sx(i) = w*sh11 + z*sh12
sy(i) = w*sh21 + z*sh22
END DO
ELSE IF (sflag.EQ.zero) THEN
sh12 = sparam(4)
sh21 = sparam(3)
DO i = 1,nsteps,incx
w = sx(i)
z = sy(i)
sx(i) = w + z*sh12
sy(i) = w*sh21 + z
END DO
ELSE
sh11 = sparam(2)
sh22 = sparam(5)
DO i = 1,nsteps,incx
w = sx(i)
z = sy(i)
sx(i) = w*sh11 + z
sy(i) = -w + sh22*z
END DO
END IF
ELSE
kx = 1
ky = 1
IF (incx.LT.0) kx = 1 + (1-n)*incx
IF (incy.LT.0) ky = 1 + (1-n)*incy
*
IF (sflag.LT.zero) THEN
sh11 = sparam(2)
sh12 = sparam(4)
sh21 = sparam(3)
sh22 = sparam(5)
DO i = 1,n
w = sx(kx)
z = sy(ky)
sx(kx) = w*sh11 + z*sh12
sy(ky) = w*sh21 + z*sh22
kx = kx + incx
ky = ky + incy
END DO
ELSE IF (sflag.EQ.zero) THEN
sh12 = sparam(4)
sh21 = sparam(3)
DO i = 1,n
w = sx(kx)
z = sy(ky)
sx(kx) = w + z*sh12
sy(ky) = w*sh21 + z
kx = kx + incx
ky = ky + incy
END DO
ELSE
sh11 = sparam(2)
sh22 = sparam(5)
DO i = 1,n
w = sx(kx)
z = sy(ky)
sx(kx) = w*sh11 + z
sy(ky) = -w + sh22*z
kx = kx + incx
ky = ky + incy
END DO
END IF
END IF
RETURN
END

+ 1
- 1
getarch.c View File

@@ -1038,7 +1038,7 @@ USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
"-DL1_DATA_SIZE=65536 -DL1_DATA_LINESIZE=64 -DL1_DATA_ASSOCIATIVE=4 " \
"-DL2_SIZE=1048576 -DL2_LINESIZE=64 -DL2_ASSOCIATIVE=16 " \
"-DDTB_DEFAULT_ENTRIES=64 -DDTB_SIZE=4096 " \
"-DHAVE_VFPV4 -DHAVE_VFPV3 -DHAVE_VFP -DHAVE_NEON -DARMV8" \
"-DHAVE_VFPV4 -DHAVE_VFPV3 -DHAVE_VFP -DHAVE_NEON -DARMV8 " \
"-march=armv8.2-a -mtune=cortex-a72"
#define LIBNAME "neoversen1"
#define CORENAME "NEOVERSEN1"


+ 3
- 0
interface/gemm.c View File

@@ -272,6 +272,9 @@ void CNAME(enum CBLAS_ORDER order, enum CBLAS_TRANSPOSE TransA, enum CBLAS_TRANS
PRINT_DEBUG_CNAME;

#if !defined(COMPLEX) && !defined(DOUBLE) && defined(USE_SGEMM_KERNEL_DIRECT)
#ifdef DYNAMIC_ARCH
if (gotoblas == &gotoblas_SKYLAKEX)
#endif
if (beta == 0 && alpha == 1.0 && order == CblasRowMajor && TransA == CblasNoTrans && TransB == CblasNoTrans && sgemm_kernel_direct_performant(m,n,k)) {
sgemm_kernel_direct(m, n, k, a, lda, b, ldb, c, ldc);
return;


Write Preview
Loading…
Cancel
Save