Come provare - tranne una operazione di matrice illegale a causa della singolarità in NumPy

Question

In NumPy, sto cercando di usare linalg per calcolare gli inversi di matrice in ogni fase di uno schema di Newton-Raphson (la dimensione del problema è piccola intenzionalmente in modo da può invertire matrici hessiane calcolate analiticamente). Tuttavia, dopo essere andato molto avanti verso la convergenza, l'Iuta si avvicina al singolare.

C'è qualche metodo in NumPy che mi consente di verificare se una matrice è considerata singolare (il determinante del calcolo non è abbastanza robusto)? Idealmente, sarebbe bello se ci fosse un modo per usare un blocco tryexcept per catturare l'errore di un singolo array di NumPy.

Come faccio? L'errore NumPy data al terminale è:

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

Answer 1

La sintassi sarebbe come questo:

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 risposta non lavora più per le versioni attuali di NumPy (sto usando 1.13 presso il tempo di scrittura). Invece:

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