Pythonの数値計算の保守 - XICA Tech blog

はじめに はじめまして、サイカのソフトウェアエンジニアの松嶋です。 サイカでは数理・統計処理を Python、Julia、R で実装しています。 今回、保守の観点から、数値計算の数値ズレについてお話ししようと思います。 皆さんは、数値計算の結果が、パッケージの更新前後や実行環境の違いで微妙に変わったことはないでしょうか? numpy のバージョンを上げるとテストが失敗するようになったり、ローカル環境ではテストが通るのにCI環境ではテストが失敗したり。 この記事では、Python を使って以下の数値変化をご紹介したいと思います。 入力変数の順序 数値解析ライブラリのバグ修正 加算アルゴリズム 尚、ここで扱うのは、情報落ちなどの浮動小数点演算で発生する計算誤差そのものの話ではないことにご注意ください。 入力変数の順序 最初に入力変数の順序による数値変化をご紹介します。 Python の dic