Làm thế nào để thử, ngoại trừ một hoạt động ma trận bất hợp pháp do số ít trong NumPy

Question

Trong NumPy, tôi đang cố gắng sử dụng linalg để tính toán nghịch đảo ma trận ở mỗi bước của một chương trình Newton-Raphson (kích thước vấn đề là nhỏ cố ý để chúng tôi có thể đảo ngược ma trận Hessian được tính toán phân tích). Tuy nhiên, sau khi tôi tiến xa về phía sự hội tụ, Hessian trở nên gần kỳ dị.

Có phương pháp nào trong NumPy cho phép tôi kiểm tra xem liệu ma trận có được coi là số ít hay không (yếu tố tính toán không đủ mạnh)? Lý tưởng nhất, nó sẽ là tốt đẹp nếu có một cách để sử dụng một khối tryexcept để bắt lỗi mảng đơn số của NumPy.

Tôi làm như thế nào? Các lỗi NumPy trao tại nhà ga là:

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

Answer 1

Cú pháp sẽ là như thế này:

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 câu trả lời không còn hoạt động cho các phiên bản hiện tại của NumPy (Tôi đang sử dụng 1,13 tại thời gian viết). Thay vào đó, hãy thực hiện:

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