La inversa ordinaria

Para una matriz cuadrada, A⁻¹ cumple A⁻¹A = I. Es una operación exacta de deshacer.

Pero solo existe cuando las columnas son linealmente independientes:

S = np.array([[2., 1.], [1., 3.]])
np.linalg.inv(S) @ S            # ≈ identidad, bien

Singular = np.array([[1., 2.], [2., 4.]])   # columna 2 = 2 x columna 1
np.linalg.inv(Singular)          # lanza LinAlgError