手を入れてみました
オリジナル版: 25.02s user 4.89s system 101% cpu 29.414 total
修正版: 5.24s user 0.47s system 109% cpu 5.220 total
†SIMD†: 1.50s user 0.41s system 134% cpu 1.412 total
i7-8700K @ 3.70GHz
https://qiita.com/omi_UT/items/a436ac1c77c20efe4e1d
以下オリジナル版
二次元音響FDTD法（Two-dimensional finite-difference time-domain method）の各種プログラミング言語の速度比較デモです。用意しているのは C言語（.c），Fortran 90（.f90），Julia（.jl），MATLAB（.m），Processing 3（.pde），Perl（.pl），Python 3（.py），Scilab（.sci） の各言語のプログラムです。
モデルは均一な気体で，モデル端面の吸収境界条件は実装していません。FDTD 法については https://ultrasonics.jp/nagatani/fdtd/ なども参照して下さい。
計算結果の音場は Visualization_with_Scilab.sci で可視化することができます（Scilab 5 以降対応）。
Speed_Comparison.txt は手元の環境で測ってみた結果です。ファイル出力および画面出力をオフにした状態で10000ステップ分を演算するのに掛かった時間を人力で記録しました。
Contents are licensed under CC BY-SA 3.0 (Creative Commons Attribution-ShareAlike 3.0 Unported License).
これらのコンテンツはクリエイティブコモンズ CC BY-SA 3.0 ライセンスの元で公開しています。
長谷 芳樹 （神戸市立工業高等専門学校）
Yoshiki NAGATANI, Kobe City College of Technology, Japan
https://ultrasonics.jp/nagatani/
https://twitter.com/nagataniyoshiki