Есть ли способ эффективно инвертировать массив матриц с помощью numpy?

Обычно я бы инвертировал массив матриц 3x3 в цикле for, как в примере ниже. К сожалению, циклы forработают медленно. Есть ли более быстрый и эффективный способ сделать это?

import numpy as np
A = np.random.rand(3,3,100)
Ainv = np.zeros_like(A)
for i in range(100):
    Ainv[:,:,i] = np.linalg.inv(A[:,:,i])