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
