Недостатък в Python скрипт е могъл да доведе до неверни резултати в повече от 100 публикации по химия

Аспирант на хавайския университет откри проблема в Python-скрипта, използван за изчисления на химическия сдвиг, определящ химическата структура на изследваното вещество при спектралния анализ на сигналите с метода на ядрено-магнитен резонанс. При проверка на резултатите от изследванията на един от своите професори, аспирантът забелязал, че при изпълнение на скрипта в различни операционни системи над един и същ набор от данни е получавал различен резултат.

Например, при стартиране в macOS 10.14 и Ubuntu 16.04 за тествания набор от данни, скриптът е давал неправилна стойност 172.4 вместо 173.2. Скриптът съдържа около 1000 реда код и се използва от химици от 2014 година. Проучването на кода показа, че изходът на неправилната стойност е обусловен от различията при сортиране на файловете в различните операционни системи. Авторите на скрипта са смятали, че функцията “glob()“ винаги връща файлове, които са сортирани по име, докато в документацията на glob е посочено, че редът на изхода не е гарантиран. Поправката се състои в добавянето на list_of_files.sort() след извикването на glob().

Откритият проблем поставя под съмнение коректността на над 100 публикации по химия, чийто изводи са направени на база на изчисления химически сдвиг, получен от скрипта. Точният брой на изследванията, в които е използван скрипта, не е известен, но публикацията с неговия код е цитирана в 158 труда. На авторите на тези трудове е препоръчано да оценят коректността на работата на скрипта в използваните за изчисление операционни системи и да направят повторни изчисления, за да се уверят в верността на получените стойности. Инцидентът е отличен пример за това, че не само качеството на провежданите експерименти, но и коректността на обработката на получените данни в програмите, които до
това широко се използваха, може да повлияе на крайния резултат.

Източник: opennet.ru