From f151cbca0d1e873a248d59fe672e277ba119d493 Mon Sep 17 00:00:00 2001
From: James Willis <james.s.willis@durham.ac.uk>
Date: Wed, 26 Apr 2017 14:10:48 +0100
Subject: [PATCH] Values not set for pOrhoj2 vector.

---
 src/runner_doiact_vec.c | 7 ++-----
 1 file changed, 2 insertions(+), 5 deletions(-)

diff --git a/src/runner_doiact_vec.c b/src/runner_doiact_vec.c
index c37d2f52c6..a531802959 100644
--- a/src/runner_doiact_vec.c
+++ b/src/runner_doiact_vec.c
@@ -400,6 +400,7 @@ __attribute__((always_inline)) INLINE static void calcRemForceInteractions(
     
     rhoj.v = vec_load(&int_cache->rhoq[*icount_align]);
     grad_hj.v = vec_load(&int_cache->grad_hq[*icount_align]);
+    pOrhoj2.v = vec_load(&int_cache->pOrho2q[*icount_align]);
     balsara_j.v = vec_load(&int_cache->balsaraq[*icount_align]);
     cj.v = vec_load(&int_cache->soundspeedq[*icount_align]);
     hj_inv.v = vec_load(&int_cache->h_invq[*icount_align]);
@@ -568,6 +569,7 @@ __attribute__((always_inline)) INLINE static void storeForceInteractions(
       
       rhoj.v = vec_load(&int_cache->rhoq[pjd]);
       grad_hj.v = vec_load(&int_cache->grad_hq[pjd]);
+      pOrhoj2.v = vec_load(&int_cache->pOrho2q[pjd]);
       balsara_j.v = vec_load(&int_cache->balsaraq[pjd]);
       cj.v = vec_load(&int_cache->soundspeedq[pjd]);
       hj_inv.v = vec_load(&int_cache->h_invq[pjd]);
@@ -1774,7 +1776,6 @@ __attribute__((always_inline)) INLINE void runner_doself2_force_vec_3(
     struct runner *r, struct cell *restrict c) {
 
 #ifdef WITH_VECTORIZATION
-  static int intCount = 0;
   const struct engine *e = r->e;
   int doi_mask;
   struct part *restrict pi;
@@ -2050,14 +2051,10 @@ __attribute__((always_inline)) INLINE void runner_doself2_force_vec_3(
       pi->force.v_sig = max(pi->force.v_sig, v_sigSum.f[k]);
     VEC_HADD(entropy_dtSum, pi->entropy_dt);
 
-    intCount += icount;
-
-    message("No. of interactions for particle %lld: %d", pi->id, icount);
     /* Reset interaction count. */
     icount = 0;
   } /* loop over all particles. */
 
-  message("No. of force interactions: %d", intCount);
   //TIMER_TOC(timer_doself_force);
 #endif /* WITH_VECTORIZATION */
 }
-- 
GitLab