https://github.com/simonsobs/pixell/blob/73d3f2b36ceeee9987d4835838831e959074c09e/pixell/enmap.py#L724 would be great to have op be an argument so that, for example, one could do something like op = lambda : a,b: np.maximum(a,b)