Как поднять массив numpy в степень? (соответствует повторяющимся умножениям матриц, а не поэлементно)

Я хочу возвести двумерный массив numpy , назовем его A , в степень некоторого числа n , но мне пока не удалось найти функцию или оператор для этого.

Я знаю, что могу привести его к типу matrix и использовать тот факт, что тогда (аналогично тому, как будет действовать в Matlab), A ** n делает именно то, что я хочу (для array то же выражение означает поэлементное возведение в степень). Однако преобразование в матрицу и обратно кажется довольно уродливым решением.

Конечно, должен быть хороший способ выполнить это вычисление, сохранив формат array ?