Call vectorised version of doself_subset_density when vectorisation enabled.

src/runner.c

@@ -776,7 +776,11 @@ void runner_do_ghost(struct runner *r, struct cell *c, int timer) {
 
             /* Self-interaction? */
             if (l->t->type == task_type_self)
+#ifdef WITH_VECTORIZATION
+              runner_doself_subset_density_vec(r, finger, parts, pid, count);
+#else
               runner_doself_subset_density(r, finger, parts, pid, count);
+#endif
 
             /* Otherwise, pair interaction? */
             else if (l->t->type == task_type_pair) {