My Notebook

XP PC に UTC を表示

noreply@blogger.com (TETRA) — Sat, 13 Feb 2010 11:16:00 +0000

かずぼんのホームページに存在する「Tclock Light」はタスクバー右側の時計表示を自由なフォーマットにするソフトウェア。

インストールして、書式指定を次の文字列に換えることで、JSTとUTCを並列して表示することができるようになる。

JST yyyy/mm/dd ddd hh:nn:ss\nUTC td-09:00yyyy/mm/dd dde hh:nn:ss

フォントは適当な等幅フォントに変更するとよい。

最終的に得られるのはこんなイメージ。

とりあえず雰囲気は大切。

参考

VistaとTClock light

飛騨天文台NTP

noreply@blogger.com (TETRA) — Sat, 13 Feb 2010 10:43:00 +0000

plasma が NTP サーバーらしい。

SSW (Solar SoftWare) IDL Desktop Environment

noreply@blogger.com (TETRA) — Fri, 12 Feb 2010 14:15:00 +0000

SSWIDL でデスクトップ環境を立ち上げるためのメモ。

idlh も良いけど Desktop Environment は状況によっては便利。

飛騨とかの環境では、idlh で「環境設定→sswidl起動」のスクリプトが立ち上がるようになっている。

$ whereis idlh

で場所がわかる。飛騨の場合は /usr/local/bin/idlhにあった。このスクリプトをコピーして使いたい。ホームディレクトリ下に bin という名前でディレクトリを作成して、そこに /usr/local/bin/idlh を適当に名前を変えて（idlhde とか）コピーする。

該当するスクリプトファイルを開いて、一番下の行を書き換える。

sswidl $argv[*]

となっている。（ctshの方言で、sswidlというエイリアスで指定されたコマンドに、最初に呼び出されたidlhに指定されたコマンドをそのまま引き渡している）。これを

sswidlde

に書き換える。保存後は idlhde で起動するはず。

IDL Codestyle Guide?

noreply@blogger.com (TETRA) — Fri, 12 Feb 2010 13:57:00 +0000

ないのか？

NASAの “The IDL Astronomy User’s Library” を参考にしすつ、個人的にこうやって書いたらいいんじゃないか的なのをメモ。

一般的なこと

長すぎる行は $ で分割して書く。
インデントは適切に行なう。

関数・プロシージャ定義

FUNCTION/PRO FUNCTIONNAME END は大文字で書く。引数は小文字で記入する。なるべくキーワード付き引数を使うと後の拡張が楽なんだろう。キーワードは大文字で、同名変数を右側に小文字で書く。なるべく同名にするように心がける。ちなみに IDL は配列を常に裏側で参照渡ししているらしい。気持ちわるー

定義内では必ずインデントする。関数の中身はすぐ次の行から書き始める。一般的には中身の前にコメントを入れる。関数の役割や引数に関してのコメントを関数定義の中に書くか、その前に書くかは流儀が分かれる。中に書くことにする。

関数名は全て大文字で書く。

制御構文

制御構文は必ず大文字で書く。明らかに使わない方が簡潔で美しい場合を除き、 BEGIN END ブロックを使用する。

制御構文の条件には単変数（フラグ）を使用する場合を除き必ず括弧を使用する。

ラベルは必ず大文字で書く。

変数

変数は小文字で書く。ただし、一文字だけの変数、フラグは大文字で書く。

変数は命名規則を元に命名しておくことにする。具体的には Camel 記法、Pascal記法、アンダーバーセパレーションのどれかを使う。自分はアンダーバーセパレーションを利用することにした。

ポインタとグローバル変数に関しては（のみ）ハンガリアン記法 p_ および g_ を使うことにする。配列に関しては a_ を使わない。そもそも組み込みの文字列型があるのと、配列になりそうなものは別途変数名にそう書いてある（***_image とか）なので問題はないだろう。

終わり

コメントはたくさん書いた方がよさげ。

計算機が応答しない時

noreply@blogger.com (TETRA) — Thu, 04 Feb 2010 09:25:00 +0000

計算機室の計算機が応答しない時、具体的には電源を入れても初期画面（BIOS）さえ表示されない時は、一度PC本体をコンセントから抜いてみると復旧することが！

なんでここの PC 起動がおかしいんだろう。

IDL でマルチコアを使い切る

noreply@blogger.com (TETRA) — Fri, 29 Jan 2010 07:40:00 +0000

方針的な物がわかりました。たぶん。

IDL では FORK は使えなくて、 Thread はユーザーが見えない範囲で一部の関数の計算に使用されています。

まぁでも特定関数に対する Thread の利用はそんなに効率があがらないっぽいので（少なくとも FFT の場合）、もうちょっと荒粒度でマルチコアを使う機能がないかなぁと思っていました。最悪IDLを二つ立ち上げてもいいんですが、ライセンスの問題をその分消費するはずです。

IDL_IDLBridge が Fork に一番近いんですが、FORKというよりは IDL を複数立ち上げるのに近くて、それぞれのプロセスに関して別途初期化が必要、とのことでした。しかもカレントディレクトリすら引き継がないし。

というわけでそれらを解決する方法。

fork したい直前で SAVE, /all, FILENAME=’なんとか’
p2 = OBJ_NEW(‘IDL_IDLBRIDGE’) で新しいプロセスを作成
p2->Execute, ‘CD, “’ + GETENV(‘PWD’) + ‘”’ でカレントディレクトリに写る
p2->Execute, ‘RESTORE, FILENAME=”なんたら”’ で子プロセスのメモリの内容をまとめる
p2->Execute, ‘なんたら’, /NOWAIT で子プロセスに仕事をさせる。

うーん綺麗じゃないですね。子プロセスの生成コストとしては一般的なプログラミングにおける fork よりも相当悪い（copy on write じゃないし、というかそれ以前の問題で）ですが、新規 IDL（SSWIdl）を立ち上げるよりはよい、ぐらいです。

一応起動した後はパラメーター渡して色々みたいなこともできるので、まぁ使いようが無いわけではない、ぐらいで。

世界最大の太陽専用望遠鏡、ハワイに

noreply@blogger.com (TETRA) — Sun, 24 Jan 2010 13:33:00 +0000

http://www.sorae.jp/031009/3579.html

IDL de Multithread

noreply@blogger.com (TETRA) — Mon, 18 Jan 2010 08:30:00 +0000

公式テキストはこちら『Multi-Threading in IDL』

そして公式テキストからわかること

Multithread 環境下でも FFT はぜんぜんはやくならない
具体的には2Core使っても1.5倍程度にしかならない

えーと。どうしろと。

Multithread.pro なんてものを見つけました。こちら。

えーとなになに。

MULTITHREAD.PRO provides a means to make use of multiple CPUs on a single computer. While some of the IDL routines are already able to parallel process over several CPUs, it is not possible with all IDL commands, or even user written routines.

This procedure uses the IDL_Bridge objects to create multiple IDL session instances on the executing machine (one per CPU), copies the necessary variables to each session, executes a command that use those variables and then harvests and re-assembles the data for output.

For example, to calculate Z = X^2 + Y^2

X = FINDGEN( 10000L)
Y = X
multithread,'Z=x[i:j]^2 + y[i:j]^2', ['x','y'], ['z'], /WAIT

will return Z at the calling scope.

This routine is particularly useful for processing multiple files:

filenames = FILE_SEARCH('*.txt')
multithread,'process_files,filenames[i:j]', ['filenames'], /WAIT

where the user routine PROCESS_FILES would be executed for each of the files in the filenames variable.

Multithread.pro が使えるのかどうかについてはわかりませんが、どうやら IDL_Bridge が IDL における Multithread の鍵を握っているらしい。しかし実際にはプロセスモデルっぽいですね。スレッドは？

感動した！

noreply@blogger.com (TETRA) — Sun, 10 Jan 2010 13:50:00 +0000

※このブログは Kazuki Matsuda （仮称）が彼のまじめで綺麗で純粋な部分だけを書くBLOGです。Rebuild of TETЯA-BЯOADCДST ですよ。どーでもいいけどキリル文字って便利！ Google もこれは補正できまい。

みなさん、HAYABUSAはもう見られましたか？僕はまだ見てないです。そんなわけで感動したものリスト。ここらへんがよくまとまっててお勧めです。まとめてくださった誰かに感謝です (>_<)

http://www.youtube.com/view_play_list?p=0CFEDB32F57C55D3

京都市付近だと大阪府科学館がやってるらしいですよ。京阪中之島線でいけます。よろしく！

Wiki 欲しい

noreply@blogger.com (TETRA) — Wed, 06 Jan 2010 16:45:00 +0000

学部の同回生はみんな Wiki を作ってそこに自分の Tips を書き込んでいるという事実。

MediaWiki が数式サポート付きで動くサーバーが欲しい……

Image Shift

noreply@blogger.com (TETRA) — Wed, 23 Dec 2009 15:48:00 +0000

pro correl_optimize

を使うとよさげ。＿〆(。。)メモメモ…

IDL de Fitting

noreply@blogger.com (TETRA) — Tue, 22 Dec 2009 13:36:00 +0000

IDL で既存データを多項式フィッティング（最小二乗法）する→poly_fit

IDL 非線形フィッティング（？）する → MPFIT

自分でプログラミングしてた。まぁ最初はいい経験。

Object Oriented Programing on IDL

noreply@blogger.com (TETRA) — Sun, 20 Dec 2009 04:57:00 +0000

ムービーのみ。

http://www.ittvis.com/portals/0/webinar/oop-in-idl/lib/playback.html

American Museum of Natural History

noreply@blogger.com (TETRA) — Fri, 18 Dec 2009 17:53:00 +0000

GIGAZINE というニュースサイトでおもしろいムービーが紹介されていました。

『いろいろな星の大きさを比較してこの地球がどの程度のものなのかがよくわかるムービー「Star Size Comparison HD」』

ついつい制作環境について考えてしまうのですが、途中からどう見ても Uniview でした。そしてなんとクレジットに Director: Carter ってそうですかこれも American Museum of Natural History の作品ですか。

Youtube にはAMNH の公式チャネルがあるようです。こちら。

先日のセミナーで流れてたムービーは “Journey to the Stars” というタイトルなんですが、それに関しての Carter 氏のコメントも掲載されていました。ところどころ挿入されている Uniview のムービーがなんとも素敵です。（本物はもっと綺麗でしたが）

やはり力の入れ方とか方向性が日本とは全然違いますね。

ちなみに Uniview の映像は日本では科学技術館（東京）のシンラドームで見えるらしいです。ぜひぜひ。

浮動小数点計算が遅い遅い遅い

noreply@blogger.com (TETRA) — Fri, 18 Dec 2009 08:56:00 +0000

2048 x 2048 の要素でちょっとした浮動小数点の計算をやらせてるのですが、こんなに遅かったっけなもうどうでもいいや。

（youtubeにアップロードする）クイックルック作ろうと思って ffmpeg 立ち上げたら segmentation fault するし。

コンパイルすればなおるのかしら。

screen コマンド！

noreply@blogger.com (TETRA) — Thu, 17 Dec 2009 13:54:00 +0000

これを探していた！

べ、別に最初っから気づいてたんだからNE!

noreply@blogger.com (TETRA) — Tue, 15 Dec 2009 14:14:00 +0000

2週間前に見た California Academy of Science で見たプラネタリウムはそのまんま Uniview で作成されていたことが発覚！

いや、最初から気づいてましたって。ほんと。ほんと。ほんと orz

学術研究における映像実践の最前線

noreply@blogger.com (TETRA) — Sun, 13 Dec 2009 07:06:00 +0000

第13回京都大学国際シンポジウム「学術研究における映像実践の最前線」に参加してきました。面白そうなところをつまみ食いするかんじで。Prof Shibata が実行委員をなされているということで特に１日目の第３セッションは期待して行ったのですが。大満足で期待以上でした。

#3の中にも4つの発表題目があって、そのうち後半の2つしか見れないというていたらくでした。だって最初はひのでの映像とかだからいいと思ったんですよ。でもやっぱりちょっと後悔。

そんなわけで、1日目のうちまともに見て、さらに印象に残ったのは最後の最後の発表だけだったんですが、それが抜群にすごかったのでメモです。まぁ一般向けなんですが。

発表は Dr. Carter Emmart (Rose Center for Earth and Space at the American Museum of Natural History, New York, USA) によるもので、『デジタル宇宙3 次元アトラスから共有コンテンツの世界へ』ということでした。

日本でいうところの 4D2U プロジェクトにあたるものとして、米国では Digital Universe というものがあって、そこでもちろん Mitaka みたいなソフトウェアを作っているという話でした。さらにちょっとソフトウェアを書き換えて商用利用（プラネタリウム施設への供給）もされているらしいです。

このソフトウェアを Uniview (ja) というらしくて。ここでは特別（？）に時計台で Uniview で作られた（？）動画を見ることが出来ました。太陽圏磁気とかいろいろ。

全体として行ってしまえば、Mitaka にシミュレーションで得られた磁力線データとかを重ねた感じなんですが、実際にサーベイで得られた画像とシミュレーションの3D（磁力線とかバウショックとか）が非常にマッチしてました。「色とか描き方の面で違和感のあるヴィジュアライジングは失敗作だよねー」みたいな話がこの後に合ったんですが、その点でいうとこれは非常に完成度が高い物だったなーと。

いやはや、まとめるといい物見れました。

あとは

場の空気を可視化するシステム！

のデモとかありました。超未来！

Prosillica Camera のデータ保存

noreply@blogger.com (TETRA) — Sun, 22 Nov 2009 14:47:00 +0000

現在の卒業研究には光学系はともかく、太陽像の取得には Prosillica という研機会社（？）の High Speed Shutter Camera を利用しています。ちなみにこれ、Smart にも搭載してそこでも Phase Diversity を行なう予定らしいです。

Prosillica のカメラのうち、こちらで現在利用しているのは GE1650 です。スペックとしては最小シャッタースピードが 50 μs （カメラ的表現で行くと 1/20,000）で、それ以降 60s = 60,000,000μs まで 1μs 単位で調整できます。シーイングさえよければ DST と適当な倍率で接続したときに 50μsで適当な intensity のイメージが撮影できます。まぁここらへんはどうでもいいや。

で、インターフェースとしては 1G の Ethernet がついてます。最近はなんでも PC 化していますね。この前、オシロスコープに USB 端子が付いてて唖然としました。

GE1650 は各ピクセルが 12bit の深さをもってます（Canon とか Nikon のセンサーといっしょ）。画素数は1600x1200 なんで 1,920,000 pix（≒ 2Mpix）です。なので画像1枚あたりの情報量は 2.88e6 bits です。で、毎秒30枚の連射ができます。

というか実はデータシートを直接見たことがないので何もわかっていないのですが、1G ether を全力で使って、12bit が 2 byte になってPC のメモリに保存されるころには 100Mbyte/s ぐらいになります。一秒で100Mbyteです。

Phase Diversity なので最低でも2台のカメラを同時に利用します。

毎秒毎秒これを確実にデータとして保存し続けるというのはあまり簡単なことではないのです。なぜかっていうとHDDの典型的な書き込み速度って最大で80MByte/sぐらいなんですね。（外側の一番おいしいところだと、 Max 100MByte 出る HDD もありますが、内側と平均したら80MByte/sぐらいがせいぜい、だと思う……）

で、 PCI-e 接続の SSD ストレージを知っていたので、G-Monster というそのうちの一つ（民生品）の話とかをしていたんですが、その後色々調べてみました。やっぱりやっぱり大手もこういうのを出していて、例えば HP のやつなんかが見つかりました。

まぁそれはいいんですけど、Phase Diversity を実用でやろうとしたら大量のデータをどう保存するのかが重要ですよね。まぁ適当に考えてみよー

独自シェル？

noreply@blogger.com (TETRA) — Thu, 12 Nov 2009 13:40:00 +0000

ssh でログアウトされた後もプログラムを走らせ続ける一番簡単な方法は、 Linux 風にnohup を使うことです。まぁシミュレーションだったら batch とか使ってもいいんじゃないでしょうか。それぞれ使い方は。

echo "コマンドを文章で" | batch midnight とやるとコマンドが勝手に夜中に実行されます。

nohup command はその場ですぐコマンドを実行しますが、ログアウト後（HUPシグナルが送られる）に自動終了しません。

プログラムが自動終了するものでなかった場合、後で自分で kill しないと root に殺されるのでご注意。（プロセスも自分も）

あ、ちなみに標準入力が閉じると終了してしまうプログラムは nohup していても落ちます。 Perlが止まったのはこのせいですね。標準入力と標準出力をファイルに結びつけておくとおもしろいことができます。

適当な名前で標準入力・標準出力用のファイルを作る ex. touch hup.in hup.out
コマンドの標準入力と標準出力をそれぞれファイルに結びつけて実行します。 nohup tail -n 0 -f hup.in | idlh >> hup.out
あとは tail -f hup.out として出力 hup.out を見ながら、片手間で hup.in にコマンドを流し込んでやります。 perl -pe'open(STDOUT,">>hup.in")
とかやって。これをあわせると疑似シェルとか作れますね。というかこんな稚拙なやり方ではない何かがすでにある気がする。

まとめ

例えば IDL の場合、
nohup tail –n 0 –f hup.in | idlh >> hup.out
とするとログアウト後もずっと存在するプロセスが生成できます。もちろん自分で書いたプログラムでも一緒。
そしてそれは以下のコマンドを入力するといつでも操ることができます。
tail -f hup.out & perl -pe'open(STDOUT,">>hup.in");'

なんか perl のワンライナーとか自分で見返してもカオスですけどまぁいいことに。シェルスクリプトだったらもっと短く書けそう。でも Perl でできればそれでいいや。

計算機室プリンタ使用法

noreply@blogger.com (TETRA) — Thu, 12 Nov 2009 08:53:00 +0000

学生控え室からプリンタを使用する方法メモ。

2009年11月12日時点で Phaser0 の IP アドレスは hoge.ho.ge.203 でした。http://10.226.87.203/ からプリンタの状態まで見ることが出来ます。

機種は Phaser 850 なのでドライバとかは富士ゼロックス Phaser 850 にあります。

ここまで調べて Windows 7 用のドライバが無いことに気づきましたとさ。めでたしめでたし…… orz

kips 性能一覧

noreply@blogger.com (TETRA) — Thu, 05 Nov 2009 06:03:00 +0000

KIPSカードは、クレジットカード会社が近鉄グループKIPS参加企業と提携し発行するクレジットカード (日本)である。カード名は、「Kintetsu International Passport System」の頭文字と、鉄道の切符に由来している。（from Wikipedia）

何かと話題になる？kipsの性能一覧です。先輩が性能を取得するコマンド（聞いたこと無かった）を教えてくれたのですが忘れてしまいました。シンプルに /proc/cpuinfo とか /proc/meminfo を見てちまちまやってます。

セキュリティ上の問題はないはず。ホスト名の頭の文字は各自で補完して下さい。

Host Name	CPU	Core Clock	Cores	Total Clock	Memory
ub	Xeon X3220	2.4 G	4	9.6 G	8 G
uf	Pentium 4 無印 2.8G	2.8 G	(2)	2.8 G	2 G
uk	Athlon 64 3200+	1.0 G	1	1.0 G	3 G
um	Opteron 246	1.0 G	2	2.0 G	2 G
un	Core 2 6700	2.66 G	2	5.3 G	3 G
uo	Core 2 Quad Q9450	2.66 G	4	10.6 G	8 G
up	Core 2 Quad Q9650	3.00 G	4	12.0 G	8 G

Pentium 4 がデュアルコアに見えるのは見せかけなのでだまされてはいけませんって感じですね。ub がなかなか高速なんですが、ub はウェブサーバーやメールサーバーの役割もあるっぽい（？）ので、ub 上でジョブを走らせるのはやめたほうがよさげです。

間違いなどあれば matsuda@kwasan まで。

Streampix...

noreply@blogger.com (TETRA) — Tue, 03 Nov 2009 06:39:00 +0000

Streampix から取得した CCD のデータがおかしかった。

Prosilica の 1600x1200 のカメラから、2x2 binding （800x600 になる）で画像を取得してファイルに保存したつもりが、どうやらハードウェアのほうでは binding がきいていなかったみたい。

おそらく 1600x1200 の画像がPCまで転送されて先頭の 800x600 分（領域では上側 1600x300に相当）が 800x600 の画像であるかのように保存されてる…

原因がわかるまでに飛騨での1日と今日の半日を費やしました orz

以下画像の例

タイムススタンプも間延びしているので気づいてよかった。

10/30の観測

noreply@blogger.com (TETRA) — Sun, 01 Nov 2009 02:37:00 +0000

ほとんど忘れかけてるけどメモメモ

10/30 に晴れていたので観測を行なった。たまたま目視で黒点が一つ観測できたのでそれを主なターゲットにすることに。後で調べところでは黒点番号1029で10月24日から発生していたらしい。 > 宇宙天気ニュース

500μm で大量に取る取る。デフォーカス量は 15mm と 9.3 mm。でも15mmのセッティングをしたときのほうが使えてる気がする。

ダークを取った後、太陽中心を使ってフラットフィールドをとろうとしたらはまった、500μs ではサチレーションしてしまうので 250μs にしてフラットフィールドを取得することに。250μs でダークを取り直し。

たぶん前半のデータはほとんど使えない。11時周辺のが一番シーイングがいいはず。

後は解析ー

SSWidl

noreply@blogger.com (TETRA) — Sun, 01 Nov 2009 02:23:00 +0000

へーそんなものがあったのか orz

後でぜったいインストールしてやる。

IDL の Ver 7.1 がなぜか共同研究室に眠っていた。7.0から7.1になるとどう変わるんだろう？