def center(system):
    """
    Center particles in the simulation cell.

    It won't check if that is done multiple times.
    """
    for p in system.particle:
        p.position -= system.cell.side / 2.0
    return system
