698a2e
From ad278554860b0da7d5848262a7bf35e058266cb1 Mon Sep 17 00:00:00 2001
698a2e
From: Andreas Arnez <arnez@linux.ibm.com>
698a2e
Date: Wed, 12 Dec 2018 20:06:27 +0100
698a2e
Subject: [PATCH 5/8] Optimizations for IBM z13
698a2e
698a2e
Perform some optimizations for IBM z13:
698a2e
- Compile with -O2 instead of -O.
698a2e
- Streamline vector loads/stores.
698a2e
- Define the vvrsum2 macro.
698a2e
698a2e
Also, use the compile option -march=z13 instead of -march=native.
698a2e
---
698a2e
 CONFIG/src/atlcomp.txt |  8 +++-----
698a2e
 include/atlas_simd.h   | 11 +++++------
698a2e
 2 files changed, 8 insertions(+), 11 deletions(-)
698a2e
698a2e
diff --git a/CONFIG/src/atlcomp.txt b/CONFIG/src/atlcomp.txt
698a2e
index aa31604..2ac71cf 100644
698a2e
--- a/CONFIG/src/atlcomp.txt
698a2e
+++ b/CONFIG/src/atlcomp.txt
698a2e
@@ -246,12 +246,10 @@ MACH=IBMz9,IBMz10,IBMz196 OS=ALL LVL=500 COMPS=f77
698a2e
    'gfortran' '-O3 -funroll-loops'
698a2e
 MACH=IBMz9,IBMz10,IBMz196,IBMz12 OS=ALL LVL=500 COMPS=smc,dmc,skc,dkc,icc,xcc,gcc
698a2e
    'gcc' '-O3 -funroll-loops'
698a2e
-MACH=IBMz13 OS=ALL LVL=1000 COMPS=dmc,skc,dkc,icc,xcc,gcc
698a2e
-   'gcc' '-march=native -O -mvx -mzvector'
698a2e
-MACH=IBMz13 OS=ALL LVL=1000 COMPS=smc
698a2e
-   'gcc' '-march=native -O -mvx -mzvector -fno-peephole -fno-peephole2'
698a2e
+MACH=IBMz13 OS=ALL LVL=1000 COMPS=smc,dmc,skc,dkc,icc,xcc,gcc
698a2e
+   'gcc' '-march=z13 -mtune=z13 -O2'
698a2e
 MACH=IBMz13 OS=ALL LVL=1000 COMPS=f77
698a2e
-   'gfortran' '-march=native -O -mvx -mzvector'
698a2e
+   'gfortran' '-march=z13 -mtune=z13 -O2'
698a2e
 #
698a2e
 # Windows defaults ; need to make SSE/SSE2 arch dep.
698a2e
 #
698a2e
diff --git a/include/atlas_simd.h b/include/atlas_simd.h
698a2e
index 68daf79..f171933 100644
698a2e
--- a/include/atlas_simd.h
698a2e
+++ b/include/atlas_simd.h
698a2e
@@ -384,8 +384,8 @@
698a2e
    #endif
698a2e
    #define ATL_VTYPE vector double
698a2e
    #if (defined(DREAL) || defined(DCPLX))
698a2e
-      #define ATL_vld(v_, p_) {v_[0] = *(p_); v_[1] = (p_)[1]; }
698a2e
-      #define ATL_vst(p_, v_)  {*(p_) = v_[0]; (p_)[1] = v_[1];}
698a2e
+      #define ATL_vld(v_, p_) v_ = *(ATL_VTYPE *)(p_)
698a2e
+      #define ATL_vst(p_, v_) *(ATL_VTYPE *)(p_) = v_
698a2e
    #else
698a2e
       #define ATL_vld(v_, p_) v_ = vec_ld2f(p_);
698a2e
       #define ATL_vst(p_, v_) vec_st2f(v_, p_);
698a2e
@@ -400,10 +400,9 @@
698a2e
    #define ATL_vmul(d_, s1_, s2_) d_ =  s1_ * s2_
698a2e
    #define ATL_vmac(d_, s1_, s2_) d_ = __builtin_s390_vec_madd(s1_, s2_, d_)
698a2e
    #define ATL_vvrsum1(s0_) \
698a2e
-   {  ATL_VTYPE t_;\
698a2e
-      t_ = vec_splat(s0_, 1); \
698a2e
-      s0_ += t_; \
698a2e
-   }
698a2e
+   {  s0_ = vec_mergeh(s0_, s0_) + vec_mergel(s0_, s0_); }
698a2e
+   #define ATL_vvrsum2(s0_, s1_) \
698a2e
+   {  s0_ = vec_mergeh(s0_, s1_) + vec_mergel(s0_, s1_); }
698a2e
    #define ATL_vsplat0(d_, s_) d_ = vec_splat(s_, 0)
698a2e
    #define ATL_vsplat1(d_, s_) d_ = vec_splat(s_, 1)
698a2e
 #elif defined(ATL_NEON) && (defined(SREAL) || defined(SCPLX))
698a2e
-- 
698a2e
2.23.0
698a2e