In this tutorial, we will look at how to check if a numpy matrix (a 2d numpy array) is invertible or not with the help of some examples.
When is a matrix invertible?
A square matrix, for example,
M is said to be invertible, if there exists a matrix
N such that,
MN = NM = I where
I is an identity matrix. The following image shows an invertible matrix.
Alternatively, we can say that
N is equal to the inverse of the matrix
M. That is, a matrix is invertible if you can take its inverse.
A square matrix that is not invertible is called a singular matrix.
How to check if a matrix is invertible in Numpy?
To check if a matrix is invertible or not in Numpy, check if it has a non-zero determinant.
If a matrix has a non-zero determinant (the determinant is not zero), we can say that the matrix is invertible. Use the
numpy.linalg.det() function to compute the determinant of a matrix. The following is the syntax –
import numpy as np # check if matrix ar is invertible not np.isclose(np.linalg.det(ar), 0)
Let’s now look at the above method with the help of some examples. First, we will create two matrices that we will use throughout this tutorial, one invertible and the other non-invertible.
import numpy as np # create matrices ar1 = np.array([ [1, 2, 3], [0, 4, 5], [0, 0, 6] ]) ar2 = np.array([ [2, 4, 6], [2, 0, 2], [6, 8, 14] ])
Here, the matrix
ar1 is invertible and the matrix
ar2 is non-invertible.
Example – Check if the determinant is non-zero
In this method, we calculate the determinant of the matrix using the
numpy.linalg.det() function and check whether it is non-zero or not. If the determinant is non-zero, we say the matrix is invertible.
Let’s apply this method to the matrices created above.
# check if ar1 is invertible print(np.linalg.det(ar1) != 0) # check if ar2 is invertible print(np.linalg.det(ar2) != 0)
True for both matrices, which is not the correct answer. The matrix
ar1 is invertible, so the first
True is okay, but the matrix
ar2 is singular and this method should give
False but we get
Why is this happening? Let’s print out the determinant for
We get a very small value which is not exactly zero.
To work around this, you can use the
numpy.isclose() function which lets you define a tolerance to determine how close two values need to be to be considered equal. By default, the tolerance is set to 1e-05, which means that two values are considered equal if they are within 0.00001 of each other. You can adjust this tolerance by passing a different value to the
rtol (relative tolerance) or
atol (absolute tolerance) parameters of
numpy.isclose() function to compare the results from
numpy.linalg.det() to zero helps resolve the above issue.
# check if ar1 is invertible print(not np.isclose(np.linalg.det(ar1), 0)) # check if ar2 is invertible print(not np.isclose(np.linalg.det(ar2), 0))
We now get the correct answer. In the above example, we first check whether the determinant is zero (very close to zero) or not, if it isn’t, we say that the matric is invertible.
You might also be interested in –
- How to check if a matrix is symmetric in Numpy?
- How to check if a matrix is a square matrix in Numpy?
- How to check if a matrix is a diagonal matrix in Numpy?
Subscribe to our newsletter for more informative guides and tutorials.
We do not spam and you can opt out any time.