np.linalg.inv

np.linalg.inv(A)[source]

Computes the inverse of the matrix.

References

Examples

In [1]: GF = galois.GF(31)

# Ensure A is full rank and invertible
In [2]: while True:
   ...:     A = GF.Random((3,3))
   ...:     if np.linalg.matrix_rank(A) == 3:
   ...:         break
   ...: 

In [3]: A
Out[3]: 
GF([[13, 21,  5],
    [ 1, 25, 22],
    [10,  5, 21]], order=31)

In [4]: A_inv = np.linalg.inv(A); A_inv
Out[4]: 
GF([[26,  8, 12],
    [23, 13,  6],
    [22, 30, 18]], order=31)

In [5]: A_inv @ A
Out[5]: 
GF([[1, 0, 0],
    [0, 1, 0],
    [0, 0, 1]], order=31)