We compute the `r_ij` anyhow, we could store it for the edges, being smaller than the full `r_ij`? - make sure, changing positions, erases it! - does this make sence? Isn't the performance fast enough to compute it otf, if we know the connectivity.