Cómo probar, excepto una operación de matriz ilegal debido a la singularidad en NumPy

Question

En NumPy, estoy tratando de usar linalg para calcular inversos de matriz en cada paso de un esquema de Newton-Raphson (el tamaño del problema es pequeño intencionalmente para que podamos puede invertir matrices de Hessian calculadas analíticamente). Sin embargo, después de llegar lejos hacia la convergencia, el Hessian se acerca al singular.

¿Hay algún método dentro de NumPy que me permita probar si una matriz se considera singular (el determinante informático no es lo suficientemente robusto)? Idealmente, sería bueno si hay una manera de utilizar un bloque tryexcept para capturar el error de matriz singular de NumPy.

¿Cómo voy a hacer esto? El error NumPy dada en la terminal es:

raise LinAlgError, 'Singular matrix' 
numpy.linalg.linalg.LinAlgError: Singular matrix 

Answer 1

La sintaxis sería así:

import numpy 

try: 
    # your code that will (maybe) throw 
except numpy.linalg.linalg.LinAlgError as err: 
    if 'Singular matrix' in err.message: 
     # your error handling block 
    else: 
     raise 

Answer 2

wim 's respuesta ya no funciona para las versiones actuales de NumPy (estoy usando 1,13 en el tiempo de escritura). En su lugar, haga:

import numpy as np 

try: 
    # your code that will (maybe) throw 
except np.linalg.LinAlgError as e: 
    if 'Singular matrix' in str(e): 
     # your error handling block 
    else: 
     raise