@property<EOL><INDENT>def edges(self):<DEDENT>
<EOL>return [self.delta[d, d] * np.arange(self.shape[d] + <NUM_LIT:1>) +<EOL>self.origin[d] - <NUM_LIT:0.5> * self.delta[d, d]<EOL>for d in range(self.rank)]<EOL>
Edges of the grid cells, origin at centre of 0,0,..,0 grid cell.

        Only works for regular, orthonormal grids.