計算屋さんの日常

2013/09/30

太った

ヤバい、ヤバい、ヤバい、ヤバい、ヤバい、ヤバい、ヤバい、ヤバい、ヤバい、
ふとった、ふとった、ふとった、ふとった、ふとった、ふとった、ふとった

百円ショップのニッケル水素電池

百均で充電電池買ってきた. Eneloopとかではなくなんかいかにもパクリっぽい名前のやつ. 純正と比べてみると軽い.

必要なものが入ってないんじゃ？

2013/09/27

だんだんわかってきた。このリファレンスがかなり使える

http://sourceforge.net/projects/openfvm/files/Reference%20Manual/

OpenFVMはcollocate grid実装なのか、しかも壁面のfluxの計算が線形になってる。高次スキームにするにはこの辺から弄るしかないか。

collocate実装ならfull implicitの実装も実現できるな。これならば数値的な安定性もいいはずだ。

2013/09/26

金沢カレー

なんか最近金沢カレーの店に行くことが多い

ゴーゴーカレー　今一だった、学食より微妙
チャンピョンカレー　普通

寝過ごしてしまった

起きたら11:00だったぜ。社会人になったら完全にアウトだな。

2013/09/25

これで問題解決できる

この手があったか（棒）

http://northwood.blog60.fc2.com/blog-entry-6914.html

プロのコード

旧帝大のせんせーの書いた差分のノードとか読んでるけど、みんな係数を時間更新のループ前に配列に入れて処理してる。このほうがいいのか？　いつも計算するときに読んできて係数計算しているのだけれど。

大規模問題とか記憶容量がかかる仕様になってしまうのでは？

VOFコードのsurfer

VOFコードのsurferのウェブサイトが復活してたので、作ってみた。一部F77か...

2001から作って何故F77を選んだ？　

コード使いまわしたのか ?

コンパイルはさっくり完了, serialしか作ってない .

GetParCをsurferのdirに置けばとりあえずできる。ただsub dirにあるMakefileの中のインストール先を手動で全部なおせとREADMEに書いてある。そのままつくったら$HOMEにできた。HAVE_FORTRANフラグはoffにした。gcc-3.4.5以降はoffにしたほうがいいと書いてあるし。

sub dir 'test' 内を見ると正常に動いているようだ。

フランスの機関で作っているコードはどうも一般的なbuildの手順と異なることが多いような気がする。scotchとか

2013/09/24

指導教官のプログラミングスキルの腕

かなり疑問に思うことが最近立て続けに起きているのだが、しかもf90とc書けないし。今時,情報系でf77だけってひどくないか？

パンチカードでも打ってろ

2013/09/23

同僚

キャンプ用品を買ってオフィスで使っている同僚がいる。

災害時は迷わず彼についていこう

2013/09/22

lenovo G 530

amazonのtime saleにlenovo G530が連日の陳列されている。つい数カ月ほど前ならば、30分もすれば売り切れだったのにこの頃は在庫がある状態いが続いている。

さすがに需要を食い尽くしたか、amazonのtime saleの割引きは大して安くないのに消費者が気づいたのか。

はじめてみた

日本で初めてBUTを使っている人見た。イギリスでも相当古い熱力の本でも見ない限り出会えないのに。

BTU = British Thermal Unit

約　1 BUT = 1 J

2013/09/21

Surfer

VOFコードのSurferの開発もとのサーバーが落ちてるな。確かフランスの研究機関だったな。名前は忘れたけど。

Appleの新製品の発売といえば...

S. jobsの死後もブレないな

http://blog.esuteru.com/archives/7323842.html

2013/09/20

コードハック中

FEMのコードをhack中やっぱり人の書いたのをどんどん読まないとダメだ。

論文通った

論文が通りました。

これで卒業できる。

’俺、論文通ったら新しいマシン買うんだ’

とフラグを立てましたが破りました。

OpenFVM

OpenFVMいいな

OpenFOAMはソースコードが巨大すぎて読みづらい。あと自分はC++はあんまり得意じゃないので。

OpenFVMはcで書かれているし、コンパクトにまとまっている。ソースコードを見るとx-flowという会社のコードらしい元商用なのかな？　アルゴリズムはsimple系ですね。速度3成分+pを陰的に解いてる。

非構造格子への有限堆積法の実装を勉強中。

sunのマシン

個人でsunのwork station( ultra 40, ultra 20, w100 )を3台保有している.

同僚によると私はヘンタイらしい.

だってsun好きだもん.

sparcのマシンを持っていないのでヤフオクで狙ってる。t1000あたりを安価で落とせないかな？

2013/09/18

kindleは今のところ文系向きだと思う

理系の専門書のkindleのサンプルをいくつかダウンロードして読んでみたけど、式の形がぐちゃぐちゃになってしまって読めたものではない。文庫本とかはストレスなく読める。

PDFで取り込めば型は崩れないけど字が小さくて読みづらい。もっと大画面のkindleが早いうちに出て欲しい。

環境

能力に恵まれた人が回りにいるのは幸いなことだ。自分との差ばかり見えて自己嫌悪に墜ちるということもあるが。

2013/09/17

有限体積法の非構造格子における実装

どうやるのかやっと見えてきた。

今考えている方法

gmshメッシュの作成したら、有限要素メッシュへ変換すればいいはずなんだけど、セルの境界面に対して垂直なベクトルと圧力と速度の格納位置で悩む。あと隣り合うセル同士のconnectivityの計算が必要。これは簡単に実装できるはずだ。

http://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/20070030207_2007030377.pdf

あと境界の設定とかゴーストセルいるよね？　
ゴーストセルどうやて作ったらいいんだろう？　メッシャー必要じゃないのか？

あんまり参考書が出ていないんだよね。論文から起こすのは結構きつい。

　