Numpy: Как попасть внутрь "переполнения, обнаруженного в exp"?

Question

Первое сообщение, пожалуйста, будьте осторожны!

Я пытаюсь root найти экспоненциальную функцию 1/(np.exp((j-w)/x)-1), и я обнаружил, что получаю RuntimeWarning: overflow encountered in exp.

Я понимаю, что это связано с ограничениями на float64, но я не могу вставить код, чтобы понять, откуда это происходит. Я попытался включить несколько try...except OverflowError блоков, чтобы я мог заглянуть внутрь, где это происходит, но, кажется, никогда не поймал это.

У меня есть идеи, как справиться с этим, пока я могу получить мой код, чтобы сделать что-то другое, когда он сталкивается с переполнением.

Заранее спасибо!

Answer 1

Я предполагаю, что проблема в том, что вы не можете понять, когда происходит предупреждение. Для этого используйте try / except, но в слегка измененном виде:

import warnings
warnings.filterwarnings("error")
...
try:
    1/(np.exp((j-w)/x)-1)
except RuntimeWarning:
    print(j)
    print(w)
    print(x)
    print((j-w)/x)

Это должно дать подсказку, в этом проблема