def make_edges(rcen): edges = np.empty(len(rcen)+1) spacing = 0.5*np.diff(rcen) edges[0], edges[-1] = rcen[0]-spacing[0], rcen[-1]+spacing[-1] edges[1:-1] = rcen[1:]-spacing return edges
def make_kedges(data): k = data[:, 0] kedges = np.empty(len(k) + 1) dk = 0.5 * np.diff(k) kedges[0], kedges[-1] = k[0] - dk[0], k[-1] + dk[-1] kedges[1:-1] = k[1:] - dk return kedges