def _remove_closest_particle(self, p):
<EOL>dp = self.pos - p<EOL>dist2 = (dp*dp).sum(axis=<NUM_LIT:1>)<EOL>ind = dist2.argmin()<EOL>rp = self.pos[ind].copy()<EOL>self.pos = np.delete(self.pos, ind, axis=<NUM_LIT:0>)<EOL>return rp<EOL>
removes the closest particle in self.pos to ``p``