1 ：DNS厨：2012/04/23(月) 01:28:05.51
CFD全般、圧縮・非圧縮、単相流・多相流何でも議論しましょう
ただし、OpenFOAMは専用スレかGoogleフォーラムで議論しましょう
役に立つリンクは>>2から
11 ：名無しさん＠５周年：2012/04/27(金) 22:20:30.00
> 10
試しに書き込んでみたがダメだった。1スレ512kまでらしい。
久しぶりに来た有限要素屋だが
続・流れの有限要素法のシミュレーションのソース間違ってるのを見つけた。
GPBiCG solverの係数betaの計算、そのままだと誤差が出る。
除算は分母、分子に分けて計算した後、逆数を掛けなきゃ桁落ちするよ。
一時変数 neum, denm, inv_denm 作ってこう書き直せ！
neum = alph * r0_r
denm = zeta * r0_r0
inv_denm = 1.0d0 / denm
beta = neum * inv_denm
他の係数計算も全部直したけれど、他は大丈夫みたい。結構危険な書き方してるな。結果には、小数第四位ぐらいからしか影響しない。10反復ぐらい余計にかかるけど。
12 ：名無しさん＠５周年：2012/04/28(土) 00:03:05.73
いつのまに次スレがｗ
あのテンポでAA貼ってると500もいかないのか
しかしDNS厨なら許す
13 ：名無しさん＠５周年：2012/04/29(日) 07:14:43.44
＜世界版ＳＰＥＥＤＩ（ＷＳＰＥＥＤＩ）＞
国内だけでなく世界の原発事故などによって放出される放射性物質の拡散状況を
気象データなどを基に計算して予測するシステム。旧ソ連チェルノブイリ原発
事故を受け、日本原子力研究所（現日本原子力研究開発機構）が１９８７年に
開発に着手し、２００９年に完成した。１００キロ～地球の半分程度まで広域に
試算できる。ＳＰＥＥＤＩの試算範囲は最大１００キロ。
14 ：DNS厨：2012/05/04(金) 18:37:40.51
､___ﾉヽ.＿,､_ﾉ､__, ＿ﾊ＿
） ゝ, ヽ /
ヽ / ,,....:'"￣:＼ ﾉ も ヽ._
/ | -|-iヽ 〉.:.:.:.:.:.:.:.:.:.:.:.:.＼ ） .う /
〉 ゝ‐' /_ﾉ （:;.-､:.:.:.:.:.:.:.:.:.:.:.:.:＼ヽ っ ヽ
､ ーｧ / 〉´ﾉ.:.:.:.:.;:ｨ'"｀ヽ、.:.:ｿ ! ! （
ﾉ （__ /^Lﾉ ､:.:.;.ｨ''､ /| ﾊ .＼:.) /
ヽ ー､ /＼・・ （', /l Xヽl. /,__Ｘli ヽ}-､,.--v'⌒
〈 ‐'ﾞ / ＼ （./ / >l/ ＜｀ lﾉ.lヽ.:l
〉 | | | / 〉〇‐'＿＿__ .ヽOﾉｲ/
〈 ・.・.・ /^Lﾉ {ヽ_〈 ｀ヾ! ,,.ｨ':／
. ) rﾞ'.:.:.｀''''===ｨ'.:.:／
},..-､,.---､,..--､fﾞ'.:.:.:.:.:fi.:.:.:.<､了
/.:.:.:.:.:.:.:.:.:.:.:fjﾆ).:.:.:.:ヽ'ｿ
./.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.:.､
Ubuntu12.04 人柱結果です
1.Unity が訳分からないんだけど～
入れたアプリをどう呼び出せばいいか分からないw
2.デフォルトコンパイラ gcc4.6(gfortran4.6) が10.04の時の Makefile でエラーになる
ライブラリの順番入れ替えると通ったけど
3.BSDライセンスになったGotoblas2がエラーでコンパイルできない
ttp://www.tacc.utexas.edu/tacc-projects/gotoblas2
=>これはOpenblasで解決したけど
4.OpenMPが1コア70%程度しか使わない
ほかにもあるけど、どゆこと？
15 ：DNS厨：2012/05/07(月) 22:48:09.93
.,, -. . . . . . . . . . . . . . ,ヽ
／. . . . . . . . λ. . . . . . . . . . ..ヽ λ
/. . . . . . . . . ../. .i. . . . . . . . . . . . ,V...）
./. . . . . . . /|. ./ .|. .ﾄ. ..λ. . . . . . . . ,k
/. . . . n. . / .| { . |::/ り い. . . . . . . . ゝ
. {. . . . /弋T￣ﾚ, ,ﾚ￣ ＼ .ゝ . . . . . . .}
. .. i. . . ./ ／￣ヽ r二 ヽ 〉 . . . . ..i
. ﾚ|. ./i { r''し . { . し } i. . . . . . .,|
. V | { .i } .i .} ./ ..i. .ｒヽ. . ...|
i | 乂ツ 乂_ツ .|...ﾄ 〉. . .|
. 〉.i xx xx |....} ./. . ..|
. i〈 __,--― 、 .|..ｲｲ. . . ..|
⊂ . i へ 「/ ） ﾚ. . . . . ...|
⊂ . |. . . . .ヽ ヽ―--- イ. . . . . . . . ...|
⊂ .|. . . . . . . .イコ二二_/ / ＼. . . . . . . .|
. |. . . . .／/ .i /＼ / /へ. . . . . . |
┌─────────────―――――――――┐
│今回の竜巻は被害がおっきかったみたいで .｜
│心からのお見舞いを申し上げます .｜
│ .｜
│米国なんかだとOptimal perturbationをハリケーンや .｜
│竜巻のシミュレーションに取り入れてるみたい .｜
│ttp://ams.confex.com/ams/pdfpapers/83940.pdf .｜
│日本もはかせの記憶が正しければ地球流体電脳 .｜
│倶楽部が1980年代にはOptimalの勉強会始めてる .｜
│みたいだけど .｜
│ .｜
│災害に対する攻性予防のためにもこの手のシミュ .｜
│レーションは予算がいっぱい付いた方がいいと思う .｜
└─────────────────――――─┘
16 ：DNS厨：2012/05/07(月) 23:16:10.15
. - ￣ 、
／ ヽ
/ /| / ヽ ｀ 、
/ /-―-| / ヽ― ,
｜ / _, -- 、' _,--、 l
｜ _| ノ ｒ'::::::ヽ ヽ r'::::ヽヽ り
｜ / .| | l::::::::::ｌ |:::::::::ｌ | /
｜ ｌ | | ゝ::::::ﾉ ヽ::::ﾉ {.ノ
｜ ＼| |" " 〉|
｜ | |. ν ------- ､ / |
｜ ゝ （ l ／ |
｜ /｀ ｰ ー―― |
｜ /ゝ＿_ 人＿ヽ |
｜ / ヽ |
┌─────────────――――――――――――――――┐
│あれ？挙げる物間違えた？Optimal perturbationの例で挙げるなら .｜
│ttp://www.fas.harvard.edu/~epsas/dynamics/vortex/generalized.pdf .｜
│の方が良いみたいだけど .｜
│参考文献的な意味でも .｜
└────────────────────────――――─┘
17 ：DNS厨：2012/05/10(木) 02:03:49.97
.． ｚ==ﾆﾍ￣ ' ､
.: ´二＞=､ /V//,､
/ ／ 孑z .ﾊ ヽY ;
/ / 斗γ⌒ヽ V
､ i,..イ ,:'入_l{ ﾘﾊ
!ー | .〃{ﾟ::!ﾜ __.z:ｒ! !
`ー=ｲiﾄ. ::i` ￣ (ﾟ'ﾊ/i! i !
ﾘﾍ i !:､_ _，r＝1ﾘ 人 ）
.::ﾊト.r!>― <ーイzイゝ-′
/イ .:ﾘ >.::∀:.L==ｲ.ﾘ
/.ノ .:/! 込/_Y_| !1
..ノ/ｲｒ-ミ_＿_..｡!!__L__ /.:.､
≧=z' ( ( !:::::::..｀.＜::::::..,.´.:::::::イ_ヽ､_
￣￣≧=-=:::::..＞ﾟ.:. ＜⌒¨´
了三三ﾆ）
┌─────────────――――――――――┐
│ふふん>>14のOpenMPについて遊びという名の調査中 .｜
│-fopenmp を付ければ -lgomp -lpthreadを付けなくても .｜
│OpenMPとして動いてくれるが -lgomp -lpthread を .｜
│付けた方が己のコードだと一割ほど速いじゃないか .｜
│ .｜
│1コア70～80%の表示は気持ち悪いけど計算時間は .｜
│ちょっと速い位なのでこのままでも良いのかな？ .｜
└────────────────────―――┘
18 ：名無しさん＠５周年：2012/05/11(金) 01:18:56.20
> -fopenmp を付ければ -lgomp -lpthreadを付けなくても
> OpenMPとして動いてくれるが -lgomp -lpthread を
> 付けた方が己のコードだと一割ほど速いじゃないか
なにそれ謎々。コア辺り平均負荷7~8割って結構同期しているっぽいが
19 ：DNS厨：2012/05/12(土) 19:31:14.25
／≧≦ゝ"7ｰ〈:⌒: . ＜ ノ ｜ｌ｜ｌ i
价'⌒／: : : :! : : : ､: : : ー: ＜ ヽ ｜ｌ ｌ ｌ
/ Ⅵ／: : : : : :}:i ;'へ: :ミ: : : : : : ＼} f￣ ﾌ -〉 ヽ
{ ／: : : : ;': : 从:{ ＼`＜: : :＼: :ヽ| Т
`ー ７ :/: : :i: :./ ､ ヾ／ヽ: : ＼:.＼ ヽ ﾄ
/ V: :/: : : :‐/- ,, ＼､ ／}! j: : : : i: :ヽ } ＼､ ＿ -‐-
,' i: : : : : : ﾚ ≧x ゝ:zｨﾃ弐ト}: !: : :| ﾘ ﾟ - へ ノi ヾ
， {: : :/: : : { rﾃ弐ﾒヾ 弋:ｿ j:从: :i i＼ ＿j }! 〉
{i : : ,'{: : : : ､ 弋:り ､ ∧: :.∧ ＿,j ＞ j ﾉ
i! V:{ ､: : : ､:_ゝ v ｱ ｲ: :レ(⌒ ﾉ＜ ﾉ ／ ―-＝ー ＞
i: : ＼: : : :ゝ_ ゝ - ／ {: :i: :､ { ﾆ-‐￣
､:ゝ ＼:ゝ: : :≧ ｰ- ,/∧: : : :＼ ﾆ-‐￣
ヾﾐ=-＞}: :i: i_イ ∧=《》 / ゝ: ＼: : :ﾆ=‐ ／: : : : ﾆ=‐-
┌─────────────―――――――――――――――――┐
│>>18 gcc v4.2からgompの成果が取り込まれているらしいので . ｜
│フラグfopenmpだけでもOpenMP計算が出来ると考えているが .｜
｜ttp://gcc.gnu.org/projects/gomp/ .|
｜-lgomp -lpthreadでどれくらい速くなったかちゃんと調べたら .|
｜Ubuntu10.04では4.5%、Ubuntu12.04で6.3%だったよ .|
｜ .|
｜>コア辺り平均負荷7~8割って結構同期しているっぽいが .|
｜Ubuntu10.04の時はコアあたり9割超えていたので違和感があるんだよな .|
｜でもUbuntu12.04の方が計算が速いんだけどw .|
｜カーネル、GCC、モニタのどれが原因なんだろう？ .|
└───────────――――――――――――――――───┘
20 ：名無しさん＠５周年：2012/05/12(土) 22:38:27.42
エスパーするしかないが、コード自体にOpenMPディレクティブが無く
リンク要件として -lgomp が必要な場合には変わるかもしらん
実処理（CPUがアクティブ）の効率が上がって
全実行時間に対する同期（CPUが暇）の割合が増えたか？
まぁ並列は１コアvs複コアで処理時間が1/コア数に近づいたかってのが
重要なんで、速くなって結果が壊れてなければいんじゃね
スレッドを全部読む
