プログラミング

物理量の空間分布を VTK ファイルに出力する [Julia]

数値計算で得られる結果は数値の列でしかないので,何かしらの方法で可視化しないといけないことが多い.2次元の簡単なものなら Excel でもいいが,3次元のデータとなると結構厄介である.幸いなことに ParaView という非常に高機能な可視化ソフトウェアがフ…

プロセス並列化+リダクション演算による配列の和の計算 [Julia]

配列要素の総和を求めるような演算は,数値計算においてしばしば表れる.たとえば Julia で書かれた次のようなコードがあるとする. # 配列 a, b の確保 sum_a = zero(eltype(a)); sum_b = zero(eltype(b)); for j in 1:N for i in 1:M sum_a += a(i,j); sum…

シルベスターの方法で生成したアダマール行列からウォルシュ型アダマール行列を得る方法 [C++]

信号処理や統計学で使われるアダマール行列は,以下のようにシルベスターの方法で再帰的に生成するのが(おそらく)一般的. ここで はクロネッカー積を表す.例として を示す. さて,この方法で生成したアダマール行列について,各行の符号反転の回数(周…

The Story of Mel

(以下の記事は2006年3月5日に書いたものです)ハッカー界に伝わる英雄叙事詩「The Story of Mel」の翻訳です。有名な割に Web 上で読める日本語訳がなかったので訳してみました。誤訳や間違いはコメント欄で指摘してくれれば幸いです。 - プログラミングの…