diff --git a/src/runner_doiact.h b/src/runner_doiact.h
index 3ea935b950a3989fbb79499a81a54f1bf8aca432..ba788865da4924ebcca6f46599d612afeb7e7954 100644
--- a/src/runner_doiact.h
+++ b/src/runner_doiact.h
@@ -2383,8 +2383,11 @@ void DOSUB_SUBSET(struct runner *r, struct cell *ci, struct part *parts,
 
     /* Otherwise, compute self-interaction. */
     else
+#ifdef WITH_VECTORIZATION
+      runner_doself_subset_density_vec(r, ci, parts, ind, count);
+#else
       DOSELF_SUBSET(r, ci, parts, ind, count);
-
+#endif
   } /* self-interaction. */
 
   /* Otherwise, it's a pair interaction. */