numpy.linalg.solve(a, b)

numpy.linalg.matrix_power(M, n)

numpy.linalg.matrix_rank(M, tol=None)

numpy.linalg.svd(a, full_matrices=1, compute_uv=1)

numpy.inner(a, b) Inner product of two arrays. Ordinary inner product of vectors for 1-D arrays (without complex conjugation)

numpy.einsum(subscripts, *operands, out=None, dtype=None, order='K', casting='safe') Evaluates the Einstein summation convention

numpy.kron(a, b)

numpy.linalg.inv(a)

numpy.linalg.qr(a, mode='reduced')