Solved by Pietro_Cicalese in post #2 
              
              
                This appears to have been solved here. This is however not in-place as far as I can tell. 
The solution: 
def batched_index_select(input, dim, index):
    for ii in range(1, len(input.shape)):
        if ii != dim:
            index = index.unsqueeze(ii)
    expanse = list(input.shape)
    expanse[0…