Removed last change as Sedov Blast 3D fails to run.

src/cache.h

@@ -172,11 +172,6 @@ __attribute__((always_inline)) INLINE void cache_read_particles(
     ci_cache->vz[i] = ci->parts[i].v[2];
   }
 
-  /* Pad cache with fake particles that exist outside the cell so will not interact.*/
-  float fake_pix = 2.0f * ci->width[0] * ci->parts[ci->count - 1].x[0];
-  for (int i = ci->count; i < ci->count + (2 * VEC_SIZE); i++)
-    ci_cache->x[i] = fake_pix;
-
 #endif
 }
 

src/runner_doiact_vec.c

@@ -445,6 +445,14 @@ __attribute__((always_inline)) INLINE void runner_doself1_density_vec(
       int pad = (num_vec_proc * VEC_SIZE) - rem;
 
       count_align += pad;
+
+      /* Set positions to the same as particle pi so when the r2 > 0 mask is
+       * applied these extra contributions are masked out.*/
+      for (int i = count; i < count_align; i++) {
+        cell_cache->x[i] = pix.f[0];
+        cell_cache->y[i] = piy.f[0];
+        cell_cache->z[i] = piz.f[0];
+      }
     }
 
     vector pjx, pjy, pjz;