ネットブックを落とす

物理的に高いところから落としたんじゃありません。

３時間ぐらい？　ずっと計算回していいたら、突然プツンっと逝きまし....た。

再起動したからよかったけど、ひやひやもんだった。負荷かけすぎたか。

やっぱり大きな計算は計算ノードでやろう。

三次元拡散方程式を解く 追記

断面を表示する方法があった。可視化するにしても軽い。

あと可視化はみんなParavirewでやってます。 追記：　二次元だけど陰解法も追加　http://a-daily-life-in-the-office.blogspot.jp/2012/04/blog-post_3867.html

イギリスの訛り

自分の研究室の人じゃないのだけれど、３つの研究室で大部屋を分けているのでイギリスの方がひとり同じで部屋にいます。初めて話したとき、liverpool周辺の人だと思って聞いたらやっぱりそうでした。

イギリス英語でも方言がいくつか存在します。代表的なものを紹介

標準イギリス英語

標準イギリス英語。ロンドン周辺など、南の広い地域の人は日常これで話します。BBCとか

ウェールズ（南ウェールズ）

自分が住んでた地域なのでたまり訛りが出ると言われました。北はかなり違うらしいです。

ニューキャッスル周辺

かなりはっきり訛りがわかるよ。この女性がニューキャスルの訛り

リバプール

ビートルズの故郷、これも結構簡単に分かる。このふくよかな青年がそう。
http://youtu.be/mHxhjhBsWi4

マンチェスター

バーミンガムとマンチェスター周辺の方言の違いは、いまいち違いがわからない

スコットランド

聞けば、はっきりわかる

アイルランド

なんとなくわかる。ゲストの人の方

拡散方程式を三次元で解いてみた

ボックス領域を陽解法の差分で解いてみた。これはさっくりできたけど、VTKへのIO書く方に時間を費やした。32x32x32がメモリーの限界でした。

三次元の可視化ってめっちゃ重いんだね。atom mechineではフレーム更新に20秒近くかかる。ネットブックは止まる->kill -9コース。

09/feb/2012そういやコードあげといた。Paraview用のvtk(vtu)I/O付き、三次元のvtkを結果ファイルとして出力するよ。 追記：　二次元だけど陰解法も追加　http://a-daily-life-in-the-office.blogspot.jp/2012/04/blog-post_3867.html
http://www.geocities.jp/yosuke_3141592/index.html

AVSからVTKに追記

統計を見ているとAVSからVTKフォーマットへの変換法を探している方が多いようなので書いておきます。

AVSのフォーマットはよく知らないけど（セレンディピティ serendipityとは扱えるの？）

でもメッシュのデータはconnectivityと座標と各ノードやセルでの値を与えてやれば良いので変換プログラムならば簡単にかけると思う。

ただしVTKではAVSの他にもセルあたりのノード数（オフセット）とセルのIDを与える必要がある。

今度書いてみます。

簡単に擬似的なフォーマットを書くと

AVS

ノード数　エレメント数 　変数の数
ID コネクティビティ
ID 座標
変数の数　変数の成分
スカラー データ
ベクトルデータ

そしてこれが例厳密には応力は2nd rankのテンソルなんだが

AVS

6 2 1 0 0
1 0.0 1.0 0.0
2 1.0 1.0 0.0
3 2.0 1.0 0.0
4 0.0 0.0 0.0
5 1.0 0.0 0.0
6 2.0 0.0 0.0
1 1 quad 4 5 2 1
2 1 quad 5 6 3 2
1 1
stress, si
1 0.0
2 50.0
3 20.0
4 0.0
5 0.0
6 0.0

VTK

ノード数　エレメント数
スカラーデータ
ベクトルデータ
コネクティビティ
オフセット値
セル　ID

参考までに
<UnstructuredGrid>を<PUnstructuredGrid>

などというように変更して開ければ良い。詳しくはwww.vtk.org/VTK/img/file-formats.pdf　のp16以降にに変更すべきタグについて書いてある。

VTK　旧形式

# vtk DataFile Version 3.0
vtk output
ASCII
DATASET UNSTRUCTURED_GRID
POINTS 4 float
0 1 1 1 1 1 0 0 1
1 0 1
CELLS 1 5
4 2 3 1 0

CELL_TYPES 1
9

CELL_DATA 1
SCALARS Material%20Id int
LOOKUP_TABLE default
1
POINT_DATA 4
SCALARS pressure float
LOOKUP_TABLE default
0 19 0 0
FIELD FieldData 1
velocity 3 4 float
1 0 0 0.4767 0.5721 0.6674 0 0 0
0 0 0

VTK XMLフォーマット

<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
<UnstructuredGrid>
<Piece NumberOfPoints=" 4 " NumberOfCells=" 1 ">
<PointData Scalars="pressure">
<DataArray type="Float32" Name="pressure" format="ascii">
0.000000
19.000000
0.000000
0.000000
</DataArray>
<DataArray type="Float32" Name="velocity" NumberOfComponents="3" format="ascii">
1.000000 0.000000 0.000000
0.4767 0.5721 0.6674
0.000000 0.000000 0.000000
0.000000 0.000000 0.000000
</DataArray>
</PointData>
<CellData Scalars="Material Id">
<DataArray type="Int32" Name="Material Id" format="ascii">
1
</DataArray>
</CellData>
<Points>
<DataArray type="Float32" Name="Points" NumberOfComponents="3" format="ascii">
0.000000 1.000000 1.000000
1.000000 1.000000 1.000000
0.000000 0.000000 1.000000
1.000000 0.000000 1.000000
</DataArray>
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="ascii">
2 3 1 0
</DataArray>
<DataArray type="Int64" Name="offsets" format="ascii" >
4
</DataArray>
<DataArray type="UInt8" Name="types" format="ascii">
9
</DataArray>
</Cells>
</Piece>
</UnstructuredGrid>
</VTKFile>

命名 マイコー

ヤフオク見ていたらかなり高性能な計算機が破格の価格で出品されていました。

もちろん落札。

Sun Microsystems Ultra 40

ですよ。opteron x 2, 8Gのやつです。世代は古いですが。研究室に運ぶのが大変だった。ありがとkさん,mさん。

人生は

隣の席の人（ゲーマー）の今日の格言

人生はレベル上げ。

悟ってるよあんた。