code: def max(self, **kwargs):
          
          out = kwargs.get('out', None)
          if out is not None:
              out[:] = self.max_pt
              return out
          else:
              return self.max_pt
docs: Return `max_pt`.
  
          Parameters
          ----------
          kwargs
              For duck-typing with `numpy.amax`
  
          See Also
          --------
          min
          odl.set.domain.IntervalProd.max
  
          Examples
          --------
          >>> g = RectGrid([1, 2, 5], [-2, 1.5, 2])
          >>> g.max()
          array([ 5.,  2.])
  
          Also works with Numpy:
  
          >>> np.max(g)
          array([ 5.,  2.])