быстрый способ инвертировать или расставить точки по матрице kxnxn

Есть ли быстрый способ вычислить обратную матрицу kxnxn с помощью numpy (обратная величина вычисляется на каждом k-слое)? Другими словами, есть ли способ векторизовать следующий код:

>>>from numpy.linalg import inv
>>>a-random(4*2*2).reshape(4,2,2)
>>>b=a.copy()
>>>for k in range(len(a)):
>>>    b[k,:,:] = inv(a[k,:,:])