Стратегии отладки числовых проблем устойчивости?
Быстрое и грязное решение состоит в том, чтобы использовать шарик
@files = glob ('/path/to/dir/*');
2 ответа
Я бы рекомендовал задать этот вопрос в списке рассылки scipy-user , возможно, с примером вашего кода. Обычно люди в списке кажутся очень опытными в числовых вычислениях и действительно полезны, просто следование списку - само по себе образование.
В остальном, боюсь, у меня нет никаких идей ... Если вы думаю, что это проблема числовой точности / округления с плавающей запятой, первое, что вы можете попробовать, - это увеличить все типы dtypes до float128 и посмотреть, есть ли разница.
Интервальная арифметика может помочь, но я не уверен, будет ли производительность достаточной для полноценной отладки при интересующих вас размерах матрицы (вы должны рассчитывать на пару порядков значительного замедления, что между заменой «скалярных» операций с плавающей запятой с высокой степенью поддержки HW на «интервальные» операции с тяжелым SW и добавлением проверки того, какие интервалы становятся слишком широкими, когда, где и почему). 1112465]