[dennou-ruby:000138] Re: reading binary files

[ken-suke@xxxxxx]

中島です。

堀之内> これって、ある (ix,iy) という一点を考えてるような気がしますがそ
堀之内> うですか。この点でのラプラシアンをとるなら、おっしゃるような行列
堀之内> を掛けて足算すればいいんですが、勘違いしてるでしょうか。各点での
堀之内> 値なら
堀之内> 
堀之内>  ( f[2:-1,1:-2] + f[0:-3,1:-2] + f[1:-2,2:-1] + f[1:-2,0:-3]
堀之内>   - 4*f[1:-2,1:-2] ) / dx^2
堀之内> 
堀之内> でいいですし、メモリー上の並びは気にしないでいいはずなんですが。
堀之内> (といっても f[2:-1,1:-2] みたいなオペレーションは f を再配分した
堀之内> 配列を作り直すので、効率重視なら嫌われるかもしれませが...)
堀之内> 
堀之内> ＃ なんかひどく誤解してるかな？

第一には、まさに効率の問題です。
あとは、単純な座標でない場合の問題ですね。

ともあれ、「後々考える」ということで OK だと思います。

-- 
 中島健介＠九州大学 大学院理学研究科 地球惑星科学専攻 
           流体圏科学大講座  地球流体力学
 〒812-8581 福岡市 東区 箱崎 6-10-1
 ken-suke@xxxxxx
 Tel   : 092-642-2682    fax : 092-642-2685