Hope no one minds, but I thought I’d share my “creative solution”.
def multi_max(input, axes, keepdim=False):
    '''
    Performs `torch.max` over multiple dimensions of `input`
    '''
    axes = sorted(axes)
    maxed = input
    for axis in reversed(axes):
        maxed, _ = maxed.max(axis, keepdim)
    return maxed

This style should work for sum as well.