Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
elpa
elpa
Commits
e13b6e7f
Commit
e13b6e7f
authored
Nov 20, 2017
by
Andreas Marek
Browse files
Fix error in Power8 block2 and 4 kernel
parent
db81a315
Changes
2
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
src/elpa2/kernels/real_vsx_2hv_template.c
View file @
e13b6e7f
...
...
@@ -58,7 +58,6 @@
#define _SSE_LOAD (__vector double) vec_ld
#define _SSE_ADD vec_add
#define _SSE_MUL vec_mul
#define _SSE_NEG vec_neg
#define _SSE_STORE vec_st
#define offset 2
#endif
...
...
@@ -68,7 +67,6 @@
#define _SSE_LOAD (__vector float) vec_ld
#define _SSE_ADD vec_add
#define _SSE_MUL vec_mul
#define _SSE_NEG vec_neg
#define _SSE_STORE vec_st
#define offset 4
#endif
...
...
@@ -468,7 +466,7 @@ void double_hh_trafo_real_vsx_2hv_single(float* q, float* hh, int* pnb, int* pnq
#ifdef HAVE_SSE_INTRINSICS
h1
=
_SSE_XOR
(
tau2
,
sign
);
#endif
#ifdef HAVE_
SPARC64
_SSE
#ifdef HAVE_
VSX
_SSE
//h1 = vec_neg(tau2);
h1
=
vec_mul
(
vec_splats
(
mone
),
tau2
);
#endif
...
...
src/elpa2/kernels/real_vsx_4hv_template.c
View file @
e13b6e7f
...
...
@@ -111,7 +111,7 @@ void quad_hh_trafo_real_vsx_4hv_single(float* q, float* hh, int* pnb, int* pnq,
!f>#ifdef HAVE_VSX_SSE
!f> interface
!f> subroutine quad_hh_trafo_real_vsx_4hv_double(q, hh, pnb, pnq, pldq, pldh) &
!f> bind(C, name="quad_hh_trafo_real_
sse
_4hv_double")
!f> bind(C, name="quad_hh_trafo_real_
vsx
_4hv_double")
!f> use, intrinsic :: iso_c_binding
!f> integer(kind=c_int) :: pnb, pnq, pldq, pldh
!f> type(c_ptr), value :: q
...
...
@@ -125,7 +125,7 @@ void quad_hh_trafo_real_vsx_4hv_single(float* q, float* hh, int* pnb, int* pnq,
!f>#ifdef HAVE_VSX_SSE
!f> interface
!f> subroutine quad_hh_trafo_real_vsx_4hv_single(q, hh, pnb, pnq, pldq, pldh) &
!f> bind(C, name="quad_hh_trafo_real_
sse
_4hv_single")
!f> bind(C, name="quad_hh_trafo_real_
vsx
_4hv_single")
!f> use, intrinsic :: iso_c_binding
!f> integer(kind=c_int) :: pnb, pnq, pldq, pldh
!f> type(c_ptr), value :: q
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment