powered by
a version of onBPBoundary to use with neighbours
onBPBoundary(se, so, j)
The sum of even bits, e.g. sum( f$even )
The sum of odd bits, e.g. sum( f$odd )
The resolution parameter nside = 2^j