Segfault while removing virtual site

[pstaerk]

import espressomd
from espressomd.virtual_sites import VirtualSitesRelative

box_l = 10

system = espressomd.System(box_l=[box_l,box_l,box_l])
system.virtual_sites = VirtualSitesRelative(have_velocity=True, 
        have_quaternion=False)

skin = .6
dt = .001
system.time_step = dt
system.cell_system.skin = skin

p0 = system.part.add(type=1, pos=[0.0316-1,3.7522-1,6.4854-1])
p0.v = [1,1,1]
p2 = system.part.add(type=2, pos=[0.0300-1,3.7528-1,6.4855-1])

p2.vs_auto_relate_to(0)
system.integrator.run(971)
p2.remove()

[fweik]

Thank you for reporting this issue. It's not actually related to virtual sites, but was a general bug in the particle removal code that you happened to trigger. This is fixed in #3288.