def get_mean_probe_position(mean_masks, site_positions):
    return (np.sum(site_positions - mean_masks[:, np.newaxis], axis=0) /
            max(1, np.sum(mean_masks)))
dataset/ETHPy150Open kwikteam/phy/phy/stats/clusters.py/get_mean_probe_position/operator `*` -> `-`