portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-04 19:25:56 +01:00
Code Releases Activity

Fix NEON SIMD

Browse Source
This commit is contained in:
Antonin Portelli 2023-04-06 11:30:48 +01:00
parent ffa7fe0cc2
commit 18ce23aa75
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Grid/simd/Grid_neon.h
View File

@ -320,7 +320,7 @@ struct Conj{
struct TimesMinusI{ struct TimesMinusI{
//Complex single //Complex single
inline float32x4_t operator()(float32x4_t in, float32x4_t ret){ inline float32x4_t operator()(float32x4_t in){
// ar ai br bi -> ai -ar ai -br // ar ai br bi -> ai -ar ai -br
float32x4_t r0, r1; float32x4_t r0, r1;
r0 = vnegq_f32(in); // -ar -ai -br -bi r0 = vnegq_f32(in); // -ar -ai -br -bi
@ -328,7 +328,7 @@ struct TimesMinusI{
return vtrn1q_f32(r1, r0); // ar -ai br -bi return vtrn1q_f32(r1, r0); // ar -ai br -bi
} }
//Complex double //Complex double
inline float64x2_t operator()(float64x2_t in, float64x2_t ret){ inline float64x2_t operator()(float64x2_t in){
// a ib -> b -ia // a ib -> b -ia
float64x2_t tmp; float64x2_t tmp;
tmp = vnegq_f64(in); tmp = vnegq_f64(in);
@ -338,7 +338,7 @@ struct TimesMinusI{
struct TimesI{ struct TimesI{
//Complex single //Complex single
inline float32x4_t operator()(float32x4_t in, float32x4_t ret){ inline float32x4_t operator()(float32x4_t in){
// ar ai br bi -> -ai ar -bi br // ar ai br bi -> -ai ar -bi br
float32x4_t r0, r1; float32x4_t r0, r1;
r0 = vnegq_f32(in); // -ar -ai -br -bi r0 = vnegq_f32(in); // -ar -ai -br -bi
@ -346,7 +346,7 @@ struct TimesI{
return vtrn1q_f32(r1, in); // -ai ar -bi br return vtrn1q_f32(r1, in); // -ai ar -bi br
} }
//Complex double //Complex double
inline float64x2_t operator()(float64x2_t in, float64x2_t ret){ inline float64x2_t operator()(float64x2_t in){
// a ib -> -b ia // a ib -> -b ia
float64x2_t tmp; float64x2_t tmp;
tmp = vnegq_f64(in); tmp = vnegq_f64(in);