Skip to content
Snippets Groups Projects
Commit 622df1ee authored by James Willis's avatar James Willis
Browse files

Moved kernel independent code outside of #ifdef in kernel_deval_1_vec.

parent ee70f42f
No related branches found
No related tags found
1 merge request!320Dopair1 vectorisation merge
...@@ -413,11 +413,6 @@ __attribute__((always_inline)) INLINE static void kernel_deval_1_vec( ...@@ -413,11 +413,6 @@ __attribute__((always_inline)) INLINE static void kernel_deval_1_vec(
dw_dx->v = vec_fma(dw_dx->v, x.v, w->v); dw_dx->v = vec_fma(dw_dx->v, x.v, w->v);
w->v = vec_fma(x.v, w->v, wendland_const_c5.v); w->v = vec_fma(x.v, w->v, wendland_const_c5.v);
/* Return everything */
w->v =
vec_mul(w->v, vec_mul(kernel_constant_vec.v, kernel_gamma_inv_dim_vec.v));
dw_dx->v = vec_mul(dw_dx->v, vec_mul(kernel_constant_vec.v,
kernel_gamma_inv_dim_plus_one_vec.v));
#else #else
/* Load x and get the interval id. */ /* Load x and get the interval id. */
...@@ -440,12 +435,15 @@ __attribute__((always_inline)) INLINE static void kernel_deval_1_vec( ...@@ -440,12 +435,15 @@ __attribute__((always_inline)) INLINE static void kernel_deval_1_vec(
dw_dx->v = (dw_dx->v * x.v) + w->v; dw_dx->v = (dw_dx->v * x.v) + w->v;
w->v = (x.v * w->v) + c[k].v; w->v = (x.v * w->v) + c[k].v;
} }
#endif
/* Return everything */ /* Return everything */
w->v = w->v * kernel_constant_vec.v * kernel_gamma_inv_dim_vec.v; w->v =
dw_dx->v = vec_mul(w->v, vec_mul(kernel_constant_vec.v, kernel_gamma_inv_dim_vec.v));
dw_dx->v * kernel_constant_vec.v * kernel_gamma_inv_dim_plus_one_vec.v; dw_dx->v = vec_mul(dw_dx->v, vec_mul(kernel_constant_vec.v,
kernel_gamma_inv_dim_plus_one_vec.v));
#endif
} }
/** /**
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment