この記事を Notebook で開く計算速度よりもプログラミングの正確さや改変のしやすさにより重点がある数値計算の場合、Mathematica を用いたプログラミングに大きなアドバンテージがあります。ここでは、Mathematica で可能なより効率的かつ間違いの少ないプログラミング方法として、連想を用いた方法を紹介します。
今回のプログラムの内容は四次関数の極値を求める簡単なものですが、より複雑なプログラムにも応用できます。
プログラムで一貫して使われ、変更されない定数や数式などはグローバルに定義しておきます。また、入力する変数によらず計算できるものはここで計算しておくことで実行時間を短縮できます。
注意すべきことは、変数など実行時に変更したいものはグローバルに定義しないことです。変数がグローバルに漏れると予期しないバグの原因になり得ます。
（このプログラムでは使用しないので）余談ですが、値が十分定まっている物理定数などは手入力しなくても直接欲しい単位で取得できます。数値の打ち間違えや単位変換の間違いを防ぐことができます。誤差が大きい定数や、定義に幾つかの流儀があるものなどは注意が必要です。
例えば自然単位系への変換定数は以下のようになります。
作成したいプログラムを小さな作業のまとまりに分け、それぞれをモジュールとしてコーディングしていきます。
モジュールの入出力は連想を用いて統一し、入力と出力の両方が名前付きで格納されます。これにより、数値の取り出し間違いが減り、デバッグもしやすくなります。
ここでは、ポテンシャル V の極値を求める最初のモジュールで
を入力とし、
を出力するものを作成します。
異なるアルゴリズムを比較したい場合、出力が同じになるようにモジュールを作成することで、いつでもアルゴリズムを変更してプログラムを実行できます。
次に、二つ目のモジュールを作成しましょう。
ここでは、ポテンシャルの2階微分の値を極値の一つで求めるモジュールを作成します。
前のモジュールの出力が次のモジュールの入力となり、すべての途中結果が保存されます。これによって、何かおかしな結果が出たときに遡ってデバッグが可能になります。
また、モジュール作成時には前のモジュールまでの計算を一時的に保存し、テスト入力として使用することでデバッグ時間を短縮できます。
作成したモジュールは /@ や @ などを用いて合成し、全体のプログラムを構築します。
もし結果がおかしかった場合、途中結果を使ってデバッグすることも可能です。
異なるアルゴリズムで計算し、結果を比較することも簡単にできます。
ここでは簡単な例を用いて連想を用いたプログラミングについて紹介しました。特に結果の正確さや計算方法の比較などに重点がある数値計算を行う場合、今回紹介した方法が参考になればと思います。また、速度が重要な数値計算においても、最初に Mathematica で簡単なプロトタイプを作っておくことにより、Rust や C 言語などで書いたプログラムのデバッグとして役に立ちます。
※編集者注：導入事例をお伺いした際の内容に補足したものを寄稿いただきました。ご協力いただきありがとうございました。
