Monologue at midnight

2010年開始

2010-01-01T20:52:00.003+09:00

あけましておめでとうございます。
もう2010年か。早いもんです。

去年は本当に忙しくて結構限界チックだったので、今年はもうちょっとゆっくりしたいです。周りから聞いてはいたけど、博士の３年ってまじでしんどいですわ。正直病みました。

体調壊したり病んだりした結果、個人の生活を大事にするって重要だなと思いました。ついついもっと結果を、もっと早くとか焦り、常に研究について考えてしまいがちでした。休みの日も家で大学の端末に繋いで仕事したり、家帰ってからも夜中に作業するのがずっと続いていたのでそりゃ病むってもんですな。研究したり知らんことを勉強するのがもともと好きでやってたはずが、いつのまにか「せねばならないこと」になってたり、常に「無駄な時間を過ごしてはいかん」という気分が強くて、ちょっとした用事を頼まれても時間がつぶされたように感じてイライラしたり、楽しくなく日々を過ごすようになったりしてました。博士論文が忙しい以外にも、指導教官が忙しい時期によりによって超しんどい仕事を回してきたり、公聴会前のピリピリしてるときにトンデモなメールを送付してきたり、就職に関して（実家の事情がいろいろあるので）意見をいうと、全否定してきたり、研究室の生活において精神的にこたえることが多かったのもキツかったです。いろんな意味で一生忘れないと思います。

実家がすぐ近所なので晩ごはんを大抵親と一緒に食べてたこと、あとツレが延々と話を聞いてくれたのでなんとか乗り切れたようなもんです。ツレの勧めで専門医に相談したのもよかったかな。それらがなかったら、どっかで心が折れててこのお正月迎えてないでしょう。親とツレには本気で感謝です。

去年一年の経験から、どんだけ仕事を効率よくこなし、結果を残そうとも、体を壊したり、日々の生活が楽しくなければ何の意味もないなと思いました。頑張るのはいいことだけど、苦痛になるレベルでやってはいかんなという教訓を得た気がします。あと、他人の言うことをあまり気にしないようにしようとも思います。参考や反面教師にするのは良いけど、振り回されて無理して倒れるのはあまりにもばかばかしいです。指導教官とか上司にあたる存在は、それらしいこと言うだけ言うけれど、その意見にそって失敗した場合、結果やこっちの人生に責任持ってくれるわけじゃないですし。なので、今年はもうちょっと気楽に、自分や家族の生活を楽しくすることを主体に生きようと思います。さぼらないけど無理しない感じで。

というわけで、今年も皆様よろしくお願いします。
初詣のおみくじは２か所ともに大吉だったので幸先は良さそうです。
幸多い一年となりますように。

パブリックコメントfor文科省

2009-11-19T20:52:00.003+09:00

文科省には1200通のパブリックコメントが届いたとか。
きのう出したパブリックコメントを一応アップ。身バレはさすがにアレなので一部伏字。文章もアレだし議論の詰め方も甘い…。書いてて自分は科学技術政策とかほとんど知らないなぁというのがよく分かった。

なんとなく今まで、自分が頑張っていることを外部にアピールすることや、政治的なニオイのする行動にかかわることは、若干気恥かしく、あまりするもんじゃないという意識があった。たぶん同業者のなかにはそういう人も多いんじゃないかと思う。昔、知り合いに「今日はこれとこれと頑張った」というようなことを話すと「頑張ったことはわざわざ人に言うもんじゃない、だまって評価されるのを待つべき」というようなことをいわれたこともあったし。なので口を閉ざし、黙って論文を書いていればそれなりに評価されるだろうと言う風に進んできたらこういうことが起きた。今回も「そんなことは研究者の仕事じゃない」と黙って耐え忍ぶのもありかもしれないけど、あとから「あぁーあのときになにかしておけば…」てのは嫌だから。

文部科学省
中村正春　副大臣殿
後藤　斎　政務官殿

事業番号13　競争的資金（若手研究育成）について意見を送らせていただきます。

私は○○大学○○研究科にて博士課程に在学し、○○○○の研究を行っております○○と申します。民主党政権が目指す科学技術政策「科学技術のフロントランナーを目指して」を拝見し、「モノから人へ」という理念とそれに基づく研究者への資金直接配分や若手への支援等の方針に非常に期待しております。しかし、今回の行政刷新会議の議論を拝聴し、若手研究者育成に対する予算の減額要求に関して強い危機感を覚えました。

特に行政刷新会議にて予算減額が要請された、
１）日本学術振興会による、特別研究員制度、科学研究費
２）テニュアトラック制度
について意見があります。

私自身、特別研究員に本年度より採用され、その研究費によって数編の論文を出版し、国内外の学会に参加し発表を行ってきました。研究費を税金から頂いている以上、科学技術の発展に寄与せねばならないと日々真摯に研究を行っております。その研究費や月々の奨励費に関する「生活保護のようである」とのコメントや、博士取得者について「社会的に需要がないのに供給過多になっている」というワークグループ委員の発言は大変残念に思います。

日本の研究機関において、博士課程学生・ポストドクターといった私たち20代～30代の若手研究者の寄与は非常に大きいです。名称上、博士課程学生は単なる学生、ポストドクターは卒業後進路が決まらず研究室に残っている存在のように思われがちですが、実際は自ら研究テーマを設定し遂行する自立した研究者として活動しています。ですが、研究費とポストに恵まれておらず、非常に不安定な立場に苦しみつつ研究活動を行っています。

研究費に関して申し上げれば、特別研究員制度は博士課程学生が申請できる唯一の公募研究費であり、科学研究費のなかの若手研究はポストドクターが採択される可能性がある数少ない種別の科学研究費です。この2つの科学研究費を減らすことは、若手研究者の海外流出や自由な発想に基づくイノベーション創出の機会を奪うことにつながります。

またポストに関しては、常勤教員採用数が博士取得者数に比べて明らかに少ない現状では、テニュアトラック制度が最も現実的な制度だと考えます。大学における３７歳以下の若手教員比率が年々減少している（平成１９年度では21.3パーセント）現状を見ましても、テニュアトラック制度のように若手を積極的に採用する仕組みを整備・維持することが必要です。博士取得者がキャリアパスとしてアカデミックのみを追求することが慢性的なポスト不足を招いているかのような委員の指摘もございましたが、周囲の博士学生は企業への就職にも意欲的です。ですが、博士卒学生に関して「狭い専門領域ばかりを勉強している」、「コミュニケーション能力を欠くのでは」等の偏見を持っている人事担当者が多いため、企業は博士学生をあまり採用しない傾向があります。企業の40パーセント以上が博士課程卒の学生を、60パーセント以上がポストドクターを、そもそも採用する予定がないのが現状です。また日本企業の新卒最重視の姿勢が中途採用として扱われることの多い博士取得者とりわけポストドクターの民間企業就職を困難にしています。博士課程卒業後の民間企業へのキャリアパス形成、それに向けた博士課程在学中に企業が求める能力を備えることができるような大学カリキュラムの改革などは今後必要でしょうし、現場の学生としても望むところであります。しかし、この不況下では、テニュアトラック制度を縮小したからといって、民間企業への就職者数が増えるとは考えづらく、単なる若手研究者の切り捨てとなるでしょう。

科学研究費制度・大学制度に数々の旧弊があることは事実です。しかし、代替制度を運用する前に若手研究者育成関連の予算を削減することは、研究者の人材喪失につながります。知識・技術の引き継ぎが重要となる科学分野において1世代の人材を欠くことは、今後長きにわたり悪影響を及ぼすことでしょう。

また、民主党政権の掲げるCO2 25%削減等の実現におきましても、科学技術力の更なる増進とそれを担う若手人材育成は必須であると考えます。予算に限りがあること、削減が必要であることは理解しておりますが、とくに若手研究者人材育成に関連する予算は保守していただくよう切にお願いする次第であります。

事業仕分けの結果（科学技術関連）に対して思うこと

2009-11-18T23:25:00.006+09:00

ここ２～３日、研究室では科学技術関連の事業仕分けの話題でもちきり。新聞では主にスパコンであるとか、ＧＸロケットの話題がメインですけど研究室での話題は若手育成関連の予算削減。だって切実なんですもの。

私たち博士課程の学生が唯一申請できる公的な研究費は日本学術振興会というところが取り仕切っている「特別研究員」。この制度は、審査を経て採用された博士課程の学生に月々20万（税込）の奨励費と、年間上限を150万とした科学研究費（最近、理論系だと大体50~80万？）が支給されるというものです。博士課程の学生は身分でこそ学生ですが、自分で学会発表を行い、論文を書き投稿するなど自立した研究者としての働きを求められます。支給される科学研究費は主に学会への出張旅費や論文投稿料、必要な器材を買うことに費やされます。実験に比べてかかるお金が少ない理論系ですら国際会議に参加すると一回で30万、計算機として使うためのパソコンが1台40万以上（メモリ12Gまで増設とか数値計算用にいろいろカスタマイズするから）、論文一本投稿すると2万とか普通にかかるので、埋蔵金とか作りようがありません。ちなみに審査は結構キビシいです。審査委員会が開かれて、大勢の委員の採点の上で採否が決まるので（噂はいろいろあれど）公平な方だと思います。決して事業仕分け会議で委員の方がおっしゃったような「ドクターの生活保護」のような制度ではありません。

さてこの特別研究員の予算も削減される方向で議論されているようです。さらには任期付正職員として若手研究者を採用する制度（テニュアトラック制度）の予算も縮減するそうです。若手研究者育成関連の予算を切って、今後の科学技術政策をどういった方向に持って行くつもりなんでしょう？

応用的、すぐに実用につながるような研究内容であれば、仕分け会議の委員が言うように、企業から直接必要十分なだけ資金をいただくことも可能でしょう。しかし、そうではない基礎的な研究についてはある程度、政府が支援する必要があるかと思います。基礎研究に携わる若手研究者を支援することが、国民の生活向上にどう反映されるかは非常に見えづらいとは思います。しかし、基礎研究というのは種のようなもので、数十年後に幅広く応用され花開く可能性を秘めています。新薬の開発や高性能デバイスも、数十年前の基礎研究が基となっています。また、最近では企業がコストダウンのため、基礎研究を大学に委託することも多く、得られた成果や特許が企業の製品を通じ国民へ還元されています。そうした研究活動に大きく寄与し、下支えしている博士課程学生、ポストドクターへの予算を削ることは、日本の科学技術力全体を下げることになりかねません

Web上では危機意識を持った方々が動き始めています
Twitterハッシュタグ#f_o_s
http://twitter.com/#search?q=%23f_o_s

科学研究費補助金の一部の執行停止に対する反対署名
http://www.shomei.tv/project-1343.html

大隅典子先生のブログ
http://nosumi.exblog.jp/10452121/#10452121_1

仕分けまとめwiki
http://mercury.dbcls.jp/w/index.php?%BB%B2%B9%CD%BB%F1%CE%C1%A5%EA%A5%F3%A5%AF#bccc7185

じゃあ私たち学生でもできることは何でしょう？
sivadさんのブログならびに科学政策ニュースクリップにて研究者ができるロビイング活動についてまとめられております。
http://d.hatena.ne.jp/sivad/20090903/p1
http://d.hatena.ne.jp/scicom/20091115/p1

また文部科学省はパブリックコメントを募集しています
http://www.mext.go.jp/a_menu/kaikei/sassin/1286925.htm

パブリックコメントを投稿された方がブログにアップされておられます
http://d.hatena.ne.jp/roadman2005/

これを参考にさせてもらいつつ、私も本当に稚拙な文章ながらパブリックコメントを作成し、文科省に送信しました。民主党本部の意見コーナーにも送りました。

「メール送ったところでなにか変わるんかいな」とは思わないでもありません。でも黙っていては変わる可能性は０です。短くても、文章構成が微妙でもメールを送り意見を言った方が言わないよりは断然いいと思います。

鹿と学会と宝物

2009-10-31T22:16:00.006+09:00

気がつけば10月も終わり。今月はハイパー忙しかった。なぜならば科研費申請書シーズンですから。なんで学生なのに助教とかのポスト付きじゃねーと申請できない科研費で忙しくなるかは察してください。（下請けってやつですね）博士論文かかなあかんＤ３に計16ページの空白を埋めないといけない科研費のしごとふるとか某研究室はブラックにもほどがありますね

しかも国際会議があったりしました。場所は奈良です。おもいっきり台風来てる日にかぶってました。会場に行くと、能楽堂で講演していてかなりびっくりしました。松が描かれた金屏風を背景に光電子分光の結果を解説とかかなりシュールです。そしてコーヒーブレイクのおやつが毎回ゴージャスという。

で、会場でひろったのか、風邪こじらせてまたしばらく寝込みました。9月に引き続き毎月寝込んでますな。

最近、忙しい中でもわりと遊びにく時間はキープするように心がけてます。気分転換しないとマジで病みます。今月は水都大阪のフィナーレと、正倉院展に行ってきました。

ラッキードラゴンが火を吹いているのが見れて満足です。
正倉院展は相変わらずスゴかったです。今年はやっぱり琵琶かなぁ。伎楽？の面とかユニークなものもありました。帰りに興福寺の阿修羅像も見てきました。

というわけで今月は奈良行きまくりでした。最近奈良のガイドもよく出てるし、店も増えてるっぽいし、奈良いい感じです。近所なんだからもっとちょくちょく行きたいなぁ

熊本行き

2009-09-29T23:21:00.003+09:00

世の中は先週シルバーウィークとやらで賑やかだったようですが、私は高熱を出した上におなかを壊し、連休をベットで死んだように寝て過ごしました。しかも連休明けから熊本で学会でした。ここのところポスター発表で2時間たっているのがしんどいなぁという理由で口頭発表を選ぶことが多く、今回も口頭にしてたので、風邪ひいていよーがなんとかしていかなくては、ってことで急遽、一日遅れていくことにしました。飛行機を格安でとっていたので変更がきかないし、直前にとると高いので電車でGoです。向こうで「電車で大阪から来た」周りに言うと唖然とされましたが4時間半でつきました。しかもJR西日本のひかりレイルスターは超快適なのでゆうほど大変じゃないなぁと。まぁかなり暇ですが。通学で電車に乗りなれている身なので電車にのってぼんやりするのは結構すきです。特急リレーつばめにも乗れて満足なり。

病み上がりでふらふらしながら発表もお仕事も終わり、2日目にはがっつり名産馬刺しを食べてました。高級馬専門店「菅乃家」でランチ食べました。まじウマーー。馬刺しって高級店で食べるとこんなにもうまいんだなと。石焼セットにして焼いた肉も食べたんですが、焼いたのと刺身で味が違うのにも驚き。とりあえず熊本いったら菅乃家いっとくべきっす。

学会いったらがっつり観光する主義なのですが、今回は病み上がりなのでセーブして熊本城と水前寺公園だけ行きました。熊本城は結構見ごたえありました。なぞのキャラ加藤清正君と忍者とかもいて面白いです。

あひる@水都大阪

2009-09-12T20:09:00.003+09:00

今、大阪中之島では水都大阪なるイベントを行っているらしい。
そのイベントがあることは、例の府知事と市長の顔が水面から出てくるポスターで知っていたが、中身がどんなものかは知らなかった。
が、一昨日、某巨大掲示板で水都大阪のイベントの一環で巨大なアヒルのオブジェが川に浮かんでいるのを知って、猛烈に見たくなった。巨大なアヒルが堂島川に浮かんでるとか銀河ヒッチハイクガイドなシュールさだ。シュールなもの好きとしては見なければ。

水都大阪の会場にはアヒル以外にもいろいろ前衛的オブジェがあった。結構会場が広いので一気に見るのは無理そうなので、一番気になったラッキードラゴンを見に行く。

私が見た時はおとなしくしてたが、これ巡航したり火を噴いたりするらしい。
http://www.yanobe.com/LD_info
すごい見たいなー。いけそうなのは10月12日かぁ、覚えとこ。
で、お目当てのアヒルちゃん。

でかいwww
やっぱでかいわ。なんか満足しました。
このあひるちゃんは大阪以外の都市でこれまでにも浮かべられてきて、国境の無意味さを訴えているそうです。以外に意味深なメッセージが含まれているのだなぁ。
http://www.namura.cc/ahiru/concept.html
なんにせよ、この光景は気に入ったので、近場だしあと何回か見に行こうかなぁと思ってたり。大阪に期間中来るなら話のネタに見てみることをお勧めします。

夏休み（貴船）

2009-09-12T19:59:00.003+09:00

ちょい前の話になりますが、8月末に避暑も兼ねて貴船に母と行ってきた。
京阪が鞍馬・貴船１dayチケットというのも出しているので、それを利用。
母は体力がある方なので、鞍馬から山を越えて貴船に向かうルートにてガチで鞍馬・貴船周遊。
学部生のころにも貴船に流しそうめんを食べるべく、同じルートで行ったことがある。その時は観光地だし～とか舐めてて、サンダルで行ったら本気の登山道で涙目になった。鞍馬の山ルートはかなりきついっす。

この山を越えると、貴船。涼しい清流の上の川床はマジ天国。

アユの塩焼きウマー。
市内で食べるのと違ってはらわたのえぐみがなく、ほろにがいって感じで美味しいです。
いやー貴船いいですな。紅葉もいいらしいので、秋にも行ってみたいです

なんとも胡散臭い

2009-08-18T22:38:00.003+09:00

月末は選挙らしく、ポスターがあちこちに貼られている。

今回の選挙ってどうなるんだろう。TVはひたすら「政権交代」とわめいているけれど、政権が交代したからってよくなるわけで無いだろうに。マスコミ関係者の、聞こえが良い4文字熟語を作って自己満足してるようすがとても気持ち悪い。

現政権が不満だ→ちがうのに変えてみようというのは、とても幼稚だと思う。ひとのやっていることの上げ足をとって批判するのはとても簡単だ。批判することと、良くすることは全然違う。大体が、民主党の目玉にしている「子ども手当」と「高速道路無料化」にどれぐらいの効果があるんだろうか。子供手当は少子化対策らしいが、はたして月2万いくらかもらえるから子供をつくるかって気分になるだろうか。個人的な意見を言わせてもらうと、全くそんな気にならないな。何で、若年層が子供もたないかていうと、雇用が安定しない＆子育てに魅力を感じないからだと思うぞ。ベースとなる賃金がしっかりないと、手当もらったからって助けにならない。逆に仕事している女性の場合、子供できちゃうとほぼイコール現役引退になってしまうので、やっぱりそこでも躊躇するとおもう。研究職にしたって、いろんな公募で若手研究者＝35歳までとかいう枠があるので、それ考えると…ねぇ。

だから、お金を配るというよりも、30～40歳代の人間が、家庭があっても無理なく働けるような体制であるとか（「社畜」のように働かなくていい、給料は抑えめになるけど6時間程度勤務というのも可とか）、再雇用機会を増やすとか、働き方をどうするかって方が大事だとおもうんだけどな。

高速道路はもっと効果のほどが謎である。CO2削減って言っておいて高速道無料は本当に訳が分からない。しかも阪神とか首都高は有料。通行料の少ない田舎道を無料にして、そこの補修であるとか狸しか通らん高速の建設を税金から持って行って、挙句の果てにCO2が増加した分、産業界や電力業界に押し付けてその分こっちの支出に回ってくるのが目に見えてる。生活道路としては国道があるんだから高速道路は受益者負担でいいだろう…。

このなんとも怪しい公約はしかも、どっかから突っ込みを受けるたびに細部がコロコロ変わる。マニフェストとやらもver.1.12ぐらいになってるんじゃなかろうか。マイクロソフトも目じゃないね。一体どんなセキュリティホールがあることやら末恐ろしい。

そして国旗切り貼りして党の旗つくちゃったってか。国旗・国歌を敬えとかは押し付けないが、常識的に考えてそりゃいかんだろ。知ったら大半の日本人が不快に思うんじゃなかろうか。
選挙のある月に自国の国旗を刻んじゃうような、常識も危機管理能力０の政党。こんなのが支援されるのが本当に不思議でしょうがない。これがこのまま行くと与党とか思うと、夏バテも進みやる気も失せるってもんです。

インストールVASP（ifort11.0+mkl10.2）

2009-08-07T22:56:00.005+09:00

ちとわけあって、研究室で4年ほど稼働していたPCクラスターをばらしたものにOSとコンパイラつんでvasp動かそうというプロジェクトに今週は追われていた。学生が運用できる計算機（のちのちはMPIつかえるクラスターへ）を整備しようという計画＋ちと早急に数台計算機がいるイベントがあったので。研究室として保有しているものは、現在主力のマルチコアで早いのあるけど、いろいろインストールしてみたり、ちょっといじったりができないから学生が運用できるのがあればいいなぁと思った次第。私はVASP使わないんだけどな。
OSはCentOS5を使用。CPUはいまどきPen4。VASPつかうのでFortranコンパイラがいるのだが、もちろんifort（非商用）をもらってくる。(11.1)
あと、blas.lapackを最適化されたので使いたいのでintel math kernel library (mkl)も非商用フリーをもらってくる
アドレスはこちら。メール送るだけでもらえます
ありがたい。
http://software.intel.com/en-us/articles/non-commercial-software-download/

pen4なのでifortは32ビット版をもらってきました。(ia32)
ifortとmklのインストールはinstall.shを実行するだけなんだが、あとめんどくさいのは、pathの設定。
ユーザーのhome/aaa/.bashrcに

. /opt/intel/Compiler/11.1/046/bin/ia32/ifortvars_ia32.sh
(これは. /opt/intel/Compiler/11.1/046/bin/ifortvars.sh　ia32 とCPUの種類を引数で与えてもいいっぽい)
. /opt/intel/mkl/10.2.1.017/tools/environment/mklvars32.sh

LD_LIBRARY_PATH=/opt/intel/mkl/10.2.1.017/lib/32:$LD_LIBRARY_PATH
export LD_LIBRARY_PATH

を書き足しました。

vaspのインストールはMakefileのチューニングで盛大にこけました。ifortのバージョンが新しく変わったせいで、オプションを書き変えないといけなかったり。pen4仕様（MPI無）にするには、

FC=ifort

OFLAG=-O3 -xW -mtune=pentium4

BLAS= -L/opt/intel/mkl/10.2.1.017/lib/32 -lmkl_intel -lmkl_sequential -lmkl_core -lpthread -liomp5
LAPACK= -L/opt/intel/mkl/10.2.1.017/lib/32 -lmkl_intel -lmkl_sequential -lmkl_core -lpthread -liomp5

fft3dlib.o : fft3dlib.F
$(CPP)
$(FC) -FR -lowercase -O1 -xSSE4.1 -unroll0 -w95 -vec_report3 -c $*$(SUFFIX)

あたりを書き換えればいいみたいです。
（しかし、自分のデスクにあるCore2機はx86_64なせいかpen4とは同じようにはいかなかった）
さて、次は並列化してみるか。いやそれよりD論進めろよと。

人はさして合理的ではないらしい

2009-08-01T20:50:00.003+09:00

最近読んだ本：
「世界は分けてもわからない」
「人は原子、世界は物理法則で動く-社会物理学で読み解く人間行動」

上のは生物、下のは物理よりの内容なんだが、なんだかわりに共通の話題を扱ってる部分もあり、並行して読んでたらおもしろかった。一言で言うとMany is different　(P. W. Anderson)。集合体は部分の和以上のモノになる。しかしなんだ、下のはタイトルがトンデモ本ちっくではあるな。生協で買ってたら後輩に「なんかそれヤバい本ぽいタイトル」みたいなコメントをもらった。内容はまともだし面白いんだけどなー。めずらしく、物性物理よりの話題にからませて書いてあるから近しく感じた。物理よりの一般書といえば宇宙（ひも理論）関連な内容多いからなぁ。本文の最後の文がまたイイ引用句だし。よけいなお世話だけど原題のThe Social Atomをそのまま「ソーシャル・アトム」にしたほうが手に取りやすいような…

内容は、社会に現れる現象（人種差別や少子化や格差）なんかは、相転移とか物理で出てくる現象に近いということ。人の行動が総じて合理的ではないこと、それゆえ人の合理性を第一原理として作られた経済学ではうまくいかない例がいろいろ挙げてあって面白い。規制緩和による自由競争でつねに価格が下がるわけではないとか（カリフォルニアの電気料金とか）。規制緩和や民営化が万能薬のように言われているが、そう単純ではないみたいだな。

一番わかりやすいのは最初にでてくる「なぜ、アメリカでは白人が住んでいる区域と黒人が住んでいるエリアが分かれるか」という話題。普通に考えればそりゃ何らかの差別があるからでは？ってなるが、差別とかなんもなくても、ただ「自分が極端なマイノリティにはなりたくない」という条件を与えるだけでクラスターが出来上がる。磁性体とかと全く同じだ。シュミレーション結果の図もイジング模型のモンテカルロの結果と似てる気が。

「極端なマイノリティになりたくない→近しい人間のまねをしとこう」ってのは人類が古くから持っている性質みたいで、理性を軽く圧倒してしまうそうだ。”サクラ”が凄く有効に働くのもそんな理由だろう。さて、そんな性質があるならばいわゆる”民意”ってのは軽く操作されてしまうんではなかろうか。マスコミはいわば大規模なサクラみたいなもので、マスコミが大声で意見を垂れ流せば、世間の数割はそれに影響されるだろう。その割合がある程度大きければ”近しい人間のまねをしておく”相互作用によってその影響は伝搬され、なんとなく社会全体に広まり、マスコミ発の意見＝”民意”になっちまう可能性がある。”民意”とかそれによって選ばれた政治家は一体何者ってことになる。ならば２大政党制はあぶないよなぁと思う。upとdownしかないスピン系で相転移がおこり強磁性体になるように、２つの選択肢しかなければ、どちらかの圧勝が起こりやすくなるんじゃなかろうか。そのパラメーターがマスコミの影響で変わるとかだったら実におそろしい。

液漏れ

2009-07-22T23:17:00.002+09:00

エアコンのリモコンの電池をかえたら液漏れしていた。
謎の液体Xが手につくわ、まさか液漏れしてると思ってなかったのでベットの上で作業してたもんでシーツにもつくわでげんなり。実家で電池液漏れ勃発のときに「電池の液はヤバいから」ってことを言ってたので調べてみたら、水酸化カリウムって強アルカリ性の物体らしいので焦る。化学実験のときに強アルカリのやけどはこわいでー肉まで溶けるでーと脅されまくった学部時代の記憶がよみがえる。（そんなふうに脅されまくったから化学系に行かなかった。）
掃除用においてあったクエン酸の水溶液作ってそれで仕上げにゆすいどいた。中和されてますように。リモコンは接点が緑青吹いてたので捨てた。ダイキンはリモコンだけ買うのOKだったので。

ググったら電池の液漏れは頻発しているみたいで、それに伴うやけど事故も頻発してるみたいだ。KOHみたいな強アルカリな液体が、ちょっとリモコンとか玩具のなかに長時間入れてたからって即効漏れ出たり（自宅エアコンのは２－３年しか使ってないはず）、しかも「電池って液漏れするから皮膚に付いたら洗ってね」てな対応でいいのかと思う。説明書にかいてあるじゃんっていっても皮膚に付いたら洗え程度じゃ洗剤の説明とかわらんし。アルカリのやけどは塩酸とか硫酸かかるより怖いっていうし。液漏れする前提なら、もうちょっと安全な物質使えよ。身近なくせに危険な物質が漏れ出てくる可能性大てのはかなりどうかと思う。代替材料を研究しないのかな。大容量な電池より安全なのが欲しい。

アルカリ乾電池以外にも、携帯のバッテリーに使われてるリチウム電池も実はあまり安定性がよろしくない物質つかってるので、強打したり火にほりこむと爆裂する可能性があるらしい。電池怖い。

それを予知してかどうかしらんが、今日交換するために買ってきてた電池はマクセルの液漏れ０設計とうたってるやつだった。はやく液漏れ撲滅してくれ。マンガン電池は液漏れしてもわりかし安全な物質みたいだ。もうマンガン電池かマクセルのしか買わない。

誕生日

2009-07-16T21:09:00.002+09:00

今日は誕生日で、しかし平日なので大学に行った。
大学にいったら、研究室の秘書さんたちに保湿クリームをプレゼントしてもらい、さらに学内にあるフレンチなレストランでお昼を奢ってもらった。ありがたやー。そしてサプライズで寄せ書きをもらった。
誕生日に寄せ書きをもらうというのはレア体験で驚いた。寄せ書きをもらうこと自体、ひたすら大学に残留してたので、送るばっかりであまりなかったからなぁ。後輩がいろいろ書いてて、ネタな文章書いてあったりで面白い。
というわけで、年をとるのも悪くないなと思った。

来週は日食があるらしい。日食グラスを買ってみた。家族が、自分もほしいというので調べたら売り切れててた。宇宙関係は人気ありすぎだ。

経済物理学て実のところどうよ

2009-07-12T21:33:00.007+09:00

経済のことは本当によく分からない。
一番謎なのが、おそらく経済学が、市場経済は無限に成長すると思っているところだ。
ニュースでもなんでも、GDPなり、とある指標の伸び率が前の年とくらべて何パーセント上下したかで騒いでいる。しかし毎年数パーセントの成長が永遠に続くはずなかろうと思う。熱力学の第一法則的に考えて。
同様にして、日本の人口が永遠に増加するというモデルで年金なり社会インフラのあれこれを設定しているのが意味不明すぎる。モデルがまちがってるんだからそこにどんだけ高度な統計解析したところでだめだろ。

一番ダメなのがまともなモデルもなく「とりあえず最低賃金は1000円、消費税は上げません、年金ダイジョーブ（財源とかこまけぇことはいいんだよ）」とか言っちゃうどっかの政治家だが。

ぎりぎりわかるのが物理なので熱力学とのアナロジーを考えてみる。外部からなんの働きかけもない閉じた系では、いずれ動きのない熱的に死んだ状態になってしまう。これはどう考えても芳しくない。
が、外部とエネルギーなり粒子なりのやりとりがある開放平衡系での定常状態というのは対流に代表されるように内部での動きがある。外部とのやりとりを貿易、内部での動きを税金とかに置き換えて、非平衡系の力学を応用したロバストで成長しなくても定常な経済状態を目指せるモデルなんてできないんかねぇという妄想。

まぁこんどこそ物理学会で経済物理のセッションを覗いてこよう

R で描画しよう

2009-06-28T21:56:00.004+09:00

MathematicaがUpdateしてから調子悪い。
描画するとラベルの字がつぶれるつぶれる。25万もするくせになにこれ、うぇーん。
EPSで出力するとフォントつぶれは回避できるんですが、MathematicaのEPSはクソ重くて、それだけで10Mとか論文に貼り付けるにはあまりに激しい容量です…。
そこで、３Dプロットを美しく書けるフリーソフトを探していたわけです。
要求することはただ一つ。離散的な3次元のデータをデータ間の補間を行いつつ３Dプロットすること。
まぁGnuplotでもそれなりのものはかけて、
set pm3d
でpm3dを読み込んでやれば、splotで
http://t16web.lanl.gov/Kawano/gnuplot/plotpm3d.html
ここにあるような３Dのカラーグラフはかけます。しかし、いかんせんGnuplotは３Dデータの補間がないようなのですよ。csplineは基本2次元plot用だし…。
そこで手をだしてみたのが"R"です。
"R"の存在は前々から知っていてその異様にシンプルな名前にびびって手を出せていなかった代物です。統計解析ツールというかそれ専用の言語＋統合開発環境です。言語なのでガチでプログラミングぽいです(ToT)
さて、Rの強みの一つはグラフィックが美しいことです。そして統計解析専用に作られているんで、補間用のパッケージとかもごろごろ転がっています。といいことづくめなようなのですが、めっちゃ躓きました。
グラフィックを出力したいのに、デフォルトのグラフィック用の関数が、等間隔のグリッドにしか対応してねぇ…。
私の場合、計算のデータポイントは基本、変動が大きそーなとこは細かく、どーでもよさそうなところは粗く、対数的にとっているのでこれはかなり困りました。
4時間ほど苦闘した結果、
akima（Akima先生？とにかくありがてぇ）っていうライブラリ中にあるinterpという関数を使って補間をつかって等間隔グリッドを作成、データをリストに格納→描画
っていう流れでなんとかなりそうな予感です。
備忘録として、
CSVからデータを取り込み→それをリストに格納→interp→描画
のテストコードをあげておきます。


book1<-read.csv("book1.csv") 
library(akima)
rlist<-as.list(NULL)
rlist$x <- c(book1[,1])
rlist$y <- c(book1[,2])
rlist$z <- c(book1[,3])
rlist.li<- interp(rlist$x,rlist$y,rlist$z, 
xo=seq(min(rlist$x),max(rlist$x),length=50),
yo=seq(min(rlist$y),max(rlist$y),length=50))
image(rlist.li)
persp(rlist.li)

read.csv()：CSVから読み込み
c()：リストへの要素の追加や結合
persp()：俯瞰図の描画
interp：補間。lengthでメッシュの数を決めることができます

CSVファイルとして

こんなのを用意しとくと

こんなかんじのができます。
これはx^2+y^2を描画してみたものです。

旧友とランチ＠牡丹園

2009-06-27T21:49:00.002+09:00

久しぶりに高校時代の友人に会った。なつかしい。メールのやり取りはちらほらあったのだが、会うのはずいぶん久しぶりである。
関東からほぼ日帰り、帰りは夜行バスらしい。母は強しだな～。

昼ごはんを一緒に食べようということで、EST内にある神戸別館牡丹園に行った。本店は名前の通り神戸にあって、昔は千里中央の阪急にも支店が入っていた。ここの焼きそばやひき肉を甘辛く炒めた物をレタスで包む料理が好物で、何度か頼んで祖父母に連れて行ってもらった記憶がある。牡蠣のお好み焼きはとくにハイパーおいしかった記憶がある。十数年ほど前？に千里中央店がなくなってしまい、それ以降神戸はさすがに遠くてなかなかいけなかったのだが、最近ESTに入ったときき、行く機会をうかがっていたのだ。梅田飲食店激戦区なだけあって、ランチメニューやらオンラインクーポンなどあってお得である。

いろいろ話を聞いていると、女の集団のなかに入るってのは大変そうだなぁと思った。女性がある程度多いとなぜか「女王蜂」が誕生してそれ以外が（男女を問わず）「働き蜂」的ポジションでへいこらして持ち上げるてな構図になりがちっぽいなぁ。それが集団を平和に維持するための最適なモデルなのかもしれんがな、個人的にはちょっとなぁ。ネットワーク理論でこういうのも解けるのかしら。

あと、母親に関する思い出で共通する部分があってものすごくウケた。あぁみんなやっぱりそんな体験があるんだなと。

Tex用エディタ

2009-06-15T22:41:00.004+09:00

texのエディタでかなり高機能なものを見つけた。
ここ
http://www.juen.ac.jp/math/nakagawa/texguide.html
にある、EasyTex。数学の先生が作られただけあって、数式関係の機能が凄い。
数式の行数を聞いてくれて、左寄せするテンプレートまで作ってくれるとは。腐るほど数式かかなあかん理論系にはありがたや。

Winshellも高機能なんだけど、Vistaだとどうもうまくいかないみたいで、後輩が困ってた。こればVista対応とあるし、勧めてみよう。サクラエディタにマクロで十分便利なんだが、ギリシア記号の入力支援とかラベル参照のときの支援があると断然楽だと思う。

もうひとつ、EclipseのTex用プラグインがなかなか秀逸そうなので、そっちも試そうと思う（タブで開けないと気持ち悪い病にすっかり侵されている）

f77とf90の混在

2009-06-07T20:13:00.006+09:00

さて、最近良くも悪くも忙しい。論文なぞ書いていて、ふと、「この数値計算ほんとに大丈夫なんだろうか」という不安が頭をよぎりだしたりして、D論を終わらせる前に計算の妥当性をちゃんと調べておこうという気になった。まーすでに何本かジャーナルに出てるんで妥当でないと困るんだが。

なんでそんなことが気になるかというと、NRGにおけるちょとしたテクがいまの系ではうまく働かないからだ。NRG（数値くりこみ群）は伝導電子バンドを対数的に離散化することで、問題を有限サイズ系の対角化にはめこんでいる。離散化の影響を最低限に抑えるために、いろんなメッシュとろうぜってのを試すとちょっとアレ？ってな感じだったりする。いまの自分のコードでは一般的な場合は影響がほとんどないんだが、伝導バンドとの混成にエネルギー依存をつけると、あやしい。一般的な場合に影響がでてないので、NRGの本体ではなくそこにほりこむパラメータ類の計算が怪しい。

パラメータの類はMathematicaで計算していた。Mathematicaは高性能かもしれんが、ブラックボックスなので不安である。なんでMathematicaかというと、計算に特殊関数のBessel関数とその数値積分がいるからだ。ニューメリカルレシピは単精度だったのでつかえねぇ…。というわけでMathematicaつかってたんだが、つい最近、高精度のBessel関数のフリーなライブラリを見つけた。
http://www.kurims.kyoto-u.ac.jp/~ooura/index-j.html
オオウラ先生ありがとぅぅぅ。Bessel関数だけじゃなく、高精度数値積分のパッケージまで作っておられる。まさにネ申。というかここのパッケージ使えば大半の物理計算ができてしまうと思われ。

恥ずかしながら、DE（二重指数関数型）公式という数値積分法は知らなかった。もう一個のIMT公式もだけど、日本人が発見した方法らしい。日本ってやっぱ数学つよいんだな。変数変換をもちいて区分求積することで猛烈に精度がよくなるらしい。なんでうちの専攻の授業ではガウス積分とかばっかでもっとつかえそうなこっち教えないんだ？

さて、ありがたくもらってきたんだが、ソースはf77で書かれている。私の普段使ってるのはf90。でも大丈夫。基本的にはおなじコンパイラ(g95,ifort)でコンパイルするので、-cオプションでオブジェクトファイルにしておいてリンクすれば.f90のファイルから呼び出すことができる。

ただし、問題が一つ。f90は結構型宣言について厳しい。私は、先頭にimplicit noneつけないと気持ち悪いんだが、これをするとfunctionを使うときにやっかい。たとえば、外部のファイルにfunction hogeがあるとして


program main

implicit none

integer :: i
double precision :: x,res
external :: hoge

res=hoge(x)

end program main

は、hogeが型宣言されてないとしてはじかれる。関数を　double precision :: hoge みたいに型宣言してもあいてがf77なので仕様が違う→結構ヘンなことになりそうで怖い。


program main

implicit none

integer :: i
double precision :: x,res

interface
    double precision function hoge(x)
    double precision , intent (in) :: x
    end function
end interface


res=hoge(x)

end program main

とinterface文でhogeの属性を指定してやるのがおそらくベスト。うぅ。いちいち使うたびにinterfaceかくのめんどくさいので、



subroutine ldbesj0(x,res)
implicit none

double precision :: x,res
interface
    double precision function dbesj0(x)
    double precision , intent (in) :: x
    end function
end interface

res=dbesj0(x)

end subroutine ldbesj0

というのを各ベッセル関数について作った。こうすると、ライブラリにしたときに、ldbesj0とかを呼べばinterfaceについて考えなくてもよいしなぁ。というわけで、明日はライブラリ作ろうと思います。

多世界解釈

2009-05-26T22:22:00.003+09:00

「宇宙を復号する」読み終わった。おもしろかったんだが、著者が多世界解釈持ち上げまくっているのが気になる。エヴァレットの多世界解釈はなんだか考えているとむずむずする概念なので。

量子力学の世界では、とある粒子がいくつかの状態をとりえる場合、それらが重ね合わさった状態をとっている。確率に応じて観測したときにどう見えるかが決まるだけであって観測されるまで粒子自体も自分がどんな状態か分かってない。（粒子がわかるてな言い方は変だけど）そんなわけで、例の半分生きてて半分死んでいるシュレディンガーの猫が出てくるわけなんだが、世の中そんな不気味なことは起こらない。私の場合、寝起きは半分死んでるけどな。なぜかというと、周囲の環境によって常にものは観測されているので（光があたるとか、熱とか）、自分たちが感知できるようなタイムスケールでは普通重ね合わせがおこらないからだ。さてその重ね合わせ→普段感知している古典世界へ、どうやって移行してるのかってのが結構な問題だったりする。なぜAでもありBでもある波動関数が、観測された瞬間AもしくはBのどちらかに決まるのか。ふつう授業とかでは、「波動関数が収縮するのです」と習う。収縮ってぱっと聞き都合がいいけどよく分からない。この「収縮するのです」理論はコペンハーゲン解釈というらしい。

これに真っ向対立するのがエヴァレットというひとが院生のときにぶちたてた多世界解釈。なんと観測された瞬間に、われわれの住んでる宇宙が状態Aをとる宇宙と状態Bをとる宇宙に分裂するらしい。観測って、真空のゆらぎやら光があたることまで考えると、毎秒何億ってレベルじゃない規模で宇宙が分裂していくわけだろ。宇宙すべての原子についてとりうる状態の個数っていくらなんだか。で自分が存在する宇宙はその無限の中の一個で…とか考えると自分がいま生存してなんかやってるってことがすごい危ういというか足元がぐらつく気分になるな。

なので個人的には多世界解釈は嫌いです。無限個の宇宙とかそういうのは物理的じゃないと思うなー。こういう観測問題とかわかりやすく書いた解説書ないかなぁ。

「宇宙を復号する」ってすごいタイトルだな

2009-05-17T21:57:00.003+09:00

最近、量子情報に興味深々。もともとはD論やら論文のイントロに書く「今やってる研究が～～の役に立つよん」の～～のところに量子コンピュータを入れれるのかどうか調べるためにあれこれ巡回してたんだが、それ自体がおもしろいなぁと。量子論ベースのことやってながらいまさらすぎる。

なんでこれまで量子コンピュータやら量子暗号やら量子テレポーテーションの話を避けていた感があるかというと、一方的に「あれって光学のほうがメインっしょ」ってイメージがあったからだ。学会とかでも量子情報な話は大抵レーザーの偏光を使っていて偏光状態のポアンカレ球うんぬんいわれても固体物理な私にゃ？？だったのだ。なんで光なのかっつうと特殊な結晶にレーザー光を当てるとうまいこと対になった光子ができて、それを量子ビット、キュビットとして使えるから便利らしい。
この対になるものが重なりあってる状態（エンタングルメント）を作れれば、それを量子情報に使える可能性があるわけで、光以外には、分子振動とか分子回転で計算しようって試みもあるらしい。光以外で有力とされてるのがスピンを使うってやつで、それは知ってたんだけど、「あーNMRとかイオントラップわかんね」とかこれもまたスルーしてた。なんか最近ダイヤモンド使うと室温でエンタングル状態が作れたとか固体のほうでも熱いっぽい。エンタングルメントの基が不純物核スピン＋電子スピンの相互作用、うは、なんか理論でも面白そう。

情報理論は興味あるけどさっぱり知らないし、しょっぱなに専門書読むのもしんどいよなぁとおもって「宇宙を復号する」ってタイトルの本を買ってみた。なーんか量子力学も情報理論を通じて考えれば、観測問題も情報を移転する操作として見るとわかりやすいんだと。大きな物体ではなんで古典的なふるまいなのかってことについていろいろ書かれててちょっとわかった気になるな。フラーレンですら、波動性をもって２重スリット実験で干渉するらしいてのはびっくりした。ウイルスで実験しようとしている人もいるらしいけど波動性があったらトンネル効果で感染するのか？

ウイルスつながりで。新型インフルで大学は休講っぽい。まぁ研究室の活動には関係ないんですがね。

インフル対策

2009-05-04T21:12:00.002+09:00

新型インフルは一体どうなるんだろう。連休明けのパンデミックとかがいやなのでとりあえずマスクと食料を備蓄。Amazonでクレべリンゲルも買った。二酸化塩素ってすごいな。インフルエンザウイルスまでノックアウトとは。しかもウイルス表面のタンパク質を修飾するって仕組みがなんかすげぇ。タンパク質を修飾しちゃうんなら長期使用するとなんか肌荒れしそうでこわいが。

Amazonで防護服まで売ってるのにはびっくりした。ゴーグルは化学実験のときに買わされたのがあるのでいざとなったらそれを転用だ。今回のは一度は収束すると思うけど、秋が怖い。スペイン風邪も春の第一波、秋の第二波とあったらしく、しかも後半のほうがやばかったらしい。南半球でまぜまぜされてパワーアップとかやめてくれよ…。夏中にマスクと殺菌用品を備蓄せんとなぁ。

秋は国際会議のシーズンなのでそのあたりも悩ましい。ことしはイタリア北部であるらしいんだが、10数時間飛行機乗って時差ボケでふらふらで、各国から大勢人が来るところに行くってパンデミック対策としては最悪だと思う。海外で病院いかなあかんとか絶対嫌だ…。なので国内のをがんばって探そうと思います。まぁそうしたほうが旅費をi7搭載したPCとでかいディスプレイ買うのにつかえそうだから個人的にはうれしいかも。ただ、今年を逃すと、来年の国際会議は北京ばっかり（なんとなく避けたい）だし、2010年以降数年間、猛烈な太陽風→電磁嵐が吹きつけて飛行機やばいかもねーシーズンの到来らしく、しばらく海外行く機会なさそうだ。

まぁそんなわけで、連休だけど家で仕事な日々。

Mathematicaメモ

2009-04-24T22:32:00.002+09:00

Mathematica覚書

論文用の図をキレーに作るため、Mathematicaと格闘して一日終わる…
Mathematicaは７にアップグレードしてからなんか挙動が怪しい。グラフを書かせると、字化けするようになった。EPSなら大丈夫なんだがMathematicaのEPSとイラレが相性超悪なのでEPSで出されても…。なんでだ…。たぶんデフォルトのフォントがうちの機会と相性悪いんだろう。
そんなわけで、グラフのラベル等のオプションを覚書

たとえば、

Plot[Sin[x], {x, 0, 10},
LabelStyle ->
Directive[Orange, FontFamily -> "Courier", FontSize -> 20]]

こんなふうにすると、ラベルの文字がオレンジ色で、FontはCourierでサイズは２０ptの図ができる。
基本的に、Labelの文字の設定はLabelStyleってのを使う。そのなかでFontのオプションはFontFamilyやらFontSizeやらいろいろあるんだが、組み合わせて使うには、Directiveってのを使う必要があるそうな。それにしてもMathematicaの解説サイトのわかりにくさはありえねー。

へたれが群論を勉強してみるテスト

2009-04-11T21:13:00.002+09:00

群論を勉強せねばならなくなった。前々から嫌な予感はしてたんだがSU(2)だのそんな用語を駆使しないと結論部分がしまりのないかんじになってしまふ。まぁそんなわけでとりあえずSU(N)だの直積表現だのがわかるようになるのを目標としてぼちぼちと勉強することにする。

せっかくなので勉強したことはチラシの裏に書きとめておく。間違ってそーな予感プンプンするな。

とりあえず、群ってなんやねんてとこから。群というのはヒジョーに抽象的な概念であるらしく、とりあえずどんなものの集合でもいいからその任意の要素AとかBとかCに対して

１．積A・Bも集合の要素
２．(AB)C=A(BC)
３．単位元がある（作用しても相手を変化させない要素、１とか単位行列みたいな）
４．逆元が要素のなかにある（元に戻すような要素が定義できる）

という4つの条件を満たせばそれらは「群」ですよと言えるっぽい。オプショナルな性質が加わると、リー群とか巡回群とか可換群とかちょっとややこしい名前になってくる。要素は具体的なモノや数字じゃなくてもよくてたとえば、右方向に90度回転させるであるとか、ｘ軸方向に30pt進めるとかのなんらかの操作でもいい。考えてみればカードをカットするというのも群、しかも巡回群を作るんじゃねーの。

群がそんな概念ですよと言われても、はぁさいでっかって感じなのだが、物理でありがたいのは対称性を群論で扱えるというところ。系の対称性はその系がどんな状態をとるかというのにダイレクトに効いてくる。量子力学だったら固有状態の縮重度が対称性で決まってくる。

対称性は、どの系がどんな対称操作（回転とか反転とか）に対して不変であるかで決まってくる。対称操作といっても、いちいち「z軸を中心とした90度回転」とかいうわけにはいかないので、基本的に行列で表わすことを考える。操作をうけるブツとして関数（ある条件でどんな値をもつかを具体的に記述したモノ、量子力学への応用なら波動関数と思っておけばいいぽい）を考えておく。こうした関数を基底という。ターゲットとする、群をなす対称操作の集合に対して閉じた（つまり、ある操作の結果として出てくる関数も集合の要素に含まれる）関数系を用いれば、対称操作の集合を行列をつかって表記できる。量子力学で出てくる完全系によるハミルトニアンの行列表記と似たよーな感じ？抽象的な対称操作を、「こんな関数に作用させるとこんな結果が出てくるよ～」として行列で見える化してるわけか。こうして得られた行列の集合を、群の表現という。関数系の選び方には任意性があるので、表現は何パターンもある可能性がある。そのなかでも適切な基底関数系を選んで、うまいことなってる状態（対角化され度合が一番高いってのが妥当なイメージか？）を既約表現という。

とりあえず、こんなところで。続くかも…？

プロキシ立てた

2009-04-10T22:57:00.004+09:00

SquidでCentOS機にプロキシサーバーを立てた。
ここの設定にあるようにsquid.confをいじる。
http://centossrv.com/squid.shtml

chkconfig squid onが「コマンドが見当たんねーよ」てなエラーがでて使えないのは、/sbin/chkconfig とフルパス書くことで回避できた。そしてこれはrootになってないと使えないコマンドらしい。

さてプロキシ立てたは良いけど、使うとき、使わないときでいちいちツール→オプション→…はめんどくさい。Firefoxを使っているので、プロキシ経由用のプロファイルを作ることにした。
プロファイルを作るには
http://support.mozilla.com/ja/kb/%E3%83%97%E3%83%AD%E3%83%95%E3%82%A1%E3%82%A4%E3%83%AB%E3%81%AE%E7%AE%A1%E7%90%86
にあるように、
ファイル名を指定して実行から
firefox.exe -ProfileManager
としてできるぽい。とりあえずproxyという名前でプロファイル作ってみる。
次にfirefox起動するとどのプロファイル使うか聞かれるので、proxyを選んで起動、ツール→オプション→詳細でプロキシの設定をする。
プロファイルが違うものを同時起動することもできて、その時は　-no-remote -p proxyのようにオプションを付ける。バッチをつくるなら

"C:\Program Files\Mozilla Firefox\firefox.exe" -no-remote -p proxy

のようにしておくとよいぽい。

Xming+putty＋ポートフォワーディング

2009-04-09T21:38:00.005+09:00

Xmingの存在は前々から知ってたけど、なんかこないだ使おうとしたときにはうまくいかなかったのだな。Xmingってのはなにかってゆーと、CygwinのXサーバー部分をとってきたようなもんで、Windows上でLinuxのXwindowつかえるというなんだか便利な代物。CygwinもVMwareも中々重たい代物なのでXwindow使うためだけには立ち上げたくないよね～てな向きには最適ぽい。研究室の自分のデスクにLinux機＋windows機の２台が置いてあるってな状態になったのであきらめずに試してみる。
通常の場合のXming＋puttyの使い方はココに書いてある通り。
http://zodiac30.cse.kyutech.ac.jp/~fujiwara/setup/xming-putty/
XmingつかったXDMCPもできます。その場合はLinux側の設定をちょっといじる必要あり。
CentOS5の場合、/etc/gdm/custom.confをちょろっといじって

[xdmcp]
Enable=True

として、rebootすればOK。あとはWindows側でXLaunch起動、One Windowを選択、IP入れるとログイン画面が出るはず。

さて問題は、自宅から研究室のLinux機を使うとき。セキュリティ上、Linux機にはグローバルIPふってないし、LAN外からはアクセスできない。なわけでポートフォワーディングかませて、大学のSSH接続用の鯖→グローバルIP付きの研究室の鯖→自分とこのマシン、ていうややこしいトンネルを掘らなければいけない。Xmingでそれがつかえるのか？？というのが今回の問題。
結論から言うと、Puttyのコマンドラインオプションで解決できた。
たとえば、目的のマシンの２２番ポートに自分とこの5100番をポートフォワードしてるとき、
Xlaunchで設定するときに
start a programを選択 →Using PuTTYを選択
Connect to computer ：localhost
で
次のページで

Additional parameters for PuTTY or SSHのところに
-P 5100

と入れておけばよい。
-Pはポート番号を指定するputtyのオプション。
puttyはそれ以外にもコマンドラインでオプションを与えて実行できるので、バッチも作れる。いつも使うのは-loadでセッションをロードするバッチを作っておいてもいいかも。
http://the.earth.li/~sgtatham/putty/0.57/htmldoc/Chapter3.html#S3.7

てなわけで、MultiWindowのは起動できた。mathematica使ってみたがおっそい…。まぁいろいろ経由してきてるんだからしゃぁないか。LAN外からのXDMCPはうまくいかず。これだけ遅いならVNCを試したほうがいいかもなぁ。

Mercurial始めました

2009-04-05T21:56:00.006+09:00

書けども書けどもD論と投稿論文が埋まらないので、ぢつとディスプレイを見て、我泣きぬれてWinshellとたわむれる土日。せっかく書いた文章を研究室の端末とやり取りしててうっかり上書きしてしまったりすると数日ガチへこみしそうなのでウマいバックアップ管理の仕方を模索してた。そんなわけで以前いいよとの噂を聞いてたMercurial。MercurialはSubversionとかCVSとかに代表されるバージョン管理システムの一種らしいんだが、Subversionとかとちがって、どっかのサーバーに変更履歴（リポジトリというらしい？）を溜めるのではなく、各端末に履歴がたまっていく。サーバー立てるとかプロキシ云々のめんどくさい設定なしで、ローカルで動くバックアップツールみたいに気軽に使うこともできるぽい。Mercurial使うときは作業してるディレクトリの中に、.hgというフォルダができてそこに変更履歴がためられる。まぁなによりMercurialという名前がええよな。Mercuryから派生してるからコマンドもHgとか無駄に凝ってる。
参考文献
http://www.02.246.ne.jp/~torutk/mercurial/intro.html
http://www.lares.dti.ne.jp/~foozy/fujiguruma/scm/mercurial.html
http://www.lares.dti.ne.jp/~foozy/fujiguruma/scm/mercurial-tortoise.html#encoding

MercurialはWindowsでもLinuxでも使える。基本はコマンドラインで動かすんだけど、Windows用にはエクスプローラーのアドオンとして使えるTortoiseHGというのがある。カメさんかわゆい。一応両方入れておく。ついでに前々からエクスプローラーからコマンドプロンプトを呼べればいいのにと思っていてそれが実現できるソフトFolderPromptも見つけたので入れておく。いやーこれ便利だわ。なんでいままで探そうとしなかったんだろう。
で、あとMercurialで差分をとってそれを表示するときに、コマンドプロンプトだとあまりにわかりにくいし、TortoiseHGのデフォで入っているKdiffは日本語字化けしたので、Winmergeをインストール。これをMercurialで使うには、Mercurial.iniをちょっといじる必要があるぽい。
TortoiseHG用


[extdiff]
;cmd.vdiff = C:\Program Files\TortoiseHg\kdiff3
;cmd.vdiff = C:\Progra~1\TortoiseSVN\bin\TortoiseMerge.exe
;cmd.vimdiff = gvim.exe
;opts.vimdiff = -f '+next' '+execute "DirDiff ".argv(0)." ".argv(1)'
cmd.vdiff = C:\Program Files\WinMerge\WinMergeU.exe
opts.vdiff = /r /e /x /ub

Mercurial本体の場合


[extensions]
; The win32text extension is available and installed by default.  It
; provides built-in Python hooks to perform line ending conversions.
; This is normally much faster than running an external program.
hgext.win32text =
hgext.extdiff =
[extdiff]
cmd.wmdiff = C:\Program Files\WinMerge\WinMergeU.exe
opts.wmdiff = /r /e /x /ub

はーMercurial関連いじってPCが全般的に使いやすくなった。

とりあえず、自分が作業してるフォルダでリポジトリ作るには


hg init

とする。（TortoiseHGだと右クリック→Create Repository Here）
そうすっと　.hgというフォルダができる。
変更履歴をとりたいファイル、たとえばAAA.txtを


hg add AAA.txt

として追加する。
とりあえずは、最初の状態を保存したいのでこの時点で


hg commit -m "initial state（なんかコメントを書く）"

とするとこれが最初の”チェンジセット”（変更履歴のまとまり）として保存される。
そのあとファイルになんか変更すると、


hg status

とやると


M AAA.txt

みたいに、変更されたことがMで表示される。


hg diff (もしくはエクステンションで定義したwmdiff等のコマンド)

で、差分を表示。


hg log

で変更履歴を表示。ある程度変更したらまたcommitしてチェンジセットを作る。
まぁしばらくこんな感じで。のちのち端末間のをうまいことまとめたりする方法をしらべていこうかなぁ。

追記１
WinMergeのxdocdiffてゆープラグインを入れると、なんとOffice系統のファイルの差分も見れるようになる。いやぁ便利便利。注意すべきはWinMergeの設定でプラグインの検出を自動にしておくこと。

あとMercurialをRedmineと連動させてみたら予想以上に便利そう。RedMine1.6.3使ってるけど、プロジェクトの設定のとこから、リポジトリのタブでSCMをMercurialで選択,ルートディレクトリを.hgがあるディレクトリへのパス(/home/aaa/workspace/test/とか)って入れるだけでおけ。差分見れたりできます。他にも色々使い方ありそうなんで、解説本に興味が…

追記２
CentOS5へのインストールでは、yum install mercurial てな呪文が使えない…。かわりに、easy_install mercurial でpythonのお力を借りるのが妥当っぽい

気疲れ

2009-04-03T22:13:00.003+09:00

新年度だ。研究室にも新人が来た。
ハイパー人見知りなもので新人の相手は猛烈に気疲れする、いやまぁほっとんど絡んでないのだけどそれでも疲れるってどんだけー。相手のキャラも不明だし、傾向と対策は謎いので話すときに必要な演算がその分多いんだよ。この先数年顔つき合わさなあかんかもしれん相手やし地雷は踏めん。むしろ全然違う研究室で初対面とかのほうがまだ話しやすいな、その場だけ猫かぶって乗り切れるもんな。とか思っておそらく露骨なほど避けまくっている。向こうから見たら話しにくい先輩ってな感じなんだろな、やれやれ。

珍しく同性の後輩が研究室に入ってきた。男の人からみたら、女性同士のほうが話しやすいんでないのとか思うかも知れんが、女同士のほうがいろいろ難しいんだよな、と女子高出身者としては思うのだ。つーか長年工学部にいるもんで、同年代の女の人と話すネタをわすれました、何話していいかわからん。絶対「痛いニュース」とか/.しらんだろうしな…。最近自分の話すネタがオタっぽいものしかない気がする。

女の人のほうが細かいことに気づきやすい分、些細な発言を気にしたりする気がする。下手な発言かませて凹まれたり恨まれたりしても嫌だしなぁ。女の後輩とかどう扱ったらいいかわからんや。いかん、思考が後ろ向きだ。

そんなわけで黙々とD論書き進めている、しばらく研究室では引き籠ってひたすらひっそり作業してるとおもう。はー、人見知りなんとかせんとなぁ。

池袋出張

2009-03-30T21:45:00.003+09:00

物理学会に行ってきました。
場所は立教大学。池袋です。すんげーレトロで美しい建物の大学でよかったです。そして立地良すぎ。池袋、個人的には気に入ったかも。でかいデパート、ジュンク堂、ハンズとかでなんでもそろうし全く不便ない。Echikaの甘物などウマーなものもあるしなぁ。ビール半額につられて行ったのに、夢の庭なる居酒屋が料理おいしかった。から揚げが予想外にうまくてびっくりした。その前に東京在住の友人が連れて行ってくれた「音音」なるところもご飯おいしかったので、池袋の印象はイイ。

今回は、はりこんでいつもより高い目のホテルに泊まった。なんか元々不眠気味で、物音が気になると眠れない。ちょっと安い目だと廊下の音声が猛烈に聞こえたり、隣の鼾が聞こえたりで安眠できない。これまでは出張で足でるとかあほらしいてやってられへんとギリッギリのところに泊ってた、そしてそのたびに疲労困憊。健康は金に代え難いと最近思うようになってきたので快適さ重視で選んだ。ストリックス東京というとこにした。コンフォートシングルにしたら広くて快適だった。ベット寝心地いいし、クーラーから変な音するのとか言ったらちゃんとチェックしてくれたし加湿器貸してくれるし、朝ごはんバイキングじゃなしに各人に出してくれるし、ネットの調子がちと悪かった以外文句なし。快適に過ごしたんでかなりまじめに学会出た。快適さはモチベーションに影響するな。

学会も今回は自分の研究にかなりリンクするのが多くて業界の流行とか探れてよかった。レビューセッションで近藤効果の話もあったので聞きに行く。教科書とかだと書いてない、歴史というか流れみたいなものをちょっと窺い知れてよかったかも。大抵の教科書はS-d模型とAnderson模型の話を全く別物のようにぶったぎって書いていて、初見だとものすごく混乱する。S-d模型のほうは、スピンありきの領域だけが有効であるとしている。その結果、Anderson模型の部分空間での有効ハミルトニアンとなっている。そのせいで、摂動展開では局在スピン領域と近藤芳田一重項領域がつながらない。Anderson模型の数値繰り込み群の論文に書いてある謎の図の見方が初めてわかったわ。模型をシンプルにしたせいでシンプルな手法では解けなくなったとかその辺が教訓的であると思う。

sakuraエディタでTeX

2009-03-20T22:00:00.005+09:00

こないだ、後輩にTeXの使い方を教える機会があって、フリーでいかにTeX環境を構築するかちゅうのをいろいろ考えた。

TeX+GS+GSViewとか本体
http://www.ms.u-tokyo.ac.jp/~abenori/mycreate/index.html
からインストーラーを頂いて、楽にインストールさせてもらう。

エディタ
sakuraでキーワード補完、色分け
マクロでplatex実行→dviout表示

がフリーでは一番楽なんじゃないかと思った。Winshellはなんか個人的に使いにくい…
そんなわけでplatex実行用のマクロを探す。
http://doter.seesaa.net/article/38297286.htmlからもらってきたものの、自分とこの環境ではMyDocumentとかにおいてあるTexファイルをコンパイルしようとするとエラーが出る。原因は、Document and Settingの空白っぽい。これだからWindowsでTexは…。というわけで、C直下に標準でインストールされるtexフォルダに一時的にコピーして、コンパイルするように変更する。
こんな感じ。あらかじめC:\tex\tempを作っとかないとエラーでる。終わった後にコピーしたtemp内のは消す仕様。JavaScript


var objShell = new ActiveXObject("WScript.Shell");
var myObject = new ActiveXObject("Scripting.FileSystemObject");
var i;
var pathlength;
var FileName;    //原型ファイル名格納用
var TeXFileName;    //TeXファイル名格納用
var DviFileName;    //DVIファイル名格納用
var PathName ; //path名格納
var FileName2; //コピー先の名前を格納
Editor.FileSave();        //ファイル上書き保存
FileName=Editor.GetFileName();        //ファイルパスの取得

TeXFileName="";        //初期化
DviFileName="";        //初期化
PathName="";
FileName2="";


for(i=0;i<FileName.length;i++){
        if(FileName.charAt(i)=="\\"){    //ファイルパスの「\」を「/」に変換
                TeXFileName=TeXFileName+"/";
        }else{
                TeXFileName=TeXFileName+FileName.charAt(i);
        }
        if(i<FileName.length-3){    //拡張子をdviに変換
                DviFileName=DviFileName+TeXFileName.charAt(i);
        }else{
                if(i==FileName.length-3)DviFileName=DviFileName+"d";
                if(i==FileName.length-2)DviFileName=DviFileName+"v";
                if(i==FileName.length-1)DviFileName=DviFileName+"i";
        }
}

for(i=0;i<FileName.length;i++){
        if (FileName.charAt(i)=="\\"){
        pathlength=i+1;
        }
}

for (i=0; i<pathlength;i++){
        PathName=PathName+TeXFileName.charAt(i);
}

for (i=0;i<FileName.length-pathlength;i++){
        FileName2=FileName2+TeXFileName.charAt(i+pathlength);
        
}
//objShell.Popup('FileNameは"' + FileName2 + '"です', 0, "ファイル名", 0); //debug用
myObject.CopyFile(TeXFileName, "C:/tex/temp/"+FileName2);        //tempフォルダをあらかじめ作っておくこと
objShell.Run("C:/tex/bin/platex "+"C:/tex/temp/"+FileName2,1,1);        //TeXコンパイラ実行(各自でパス変更)
objShell.Run("C:/tex/dviout/dviout.exe -1 "+DviFileName,1,0);    //DVIビューアの実行(各自でパス変更)
myObject.DeleteFile("C:/tex/temp/"+FileName2);

CUDA+Fortran

2009-03-13T22:07:00.009+09:00

昨日から今日にかけてホワイトデーのお返しで菓子をあちこちからいっぱいもらったのでなんか嬉しい。もともと甘いもんスキーなので。先生方からもお返しもらってしまった。若干恐縮するなぁ。
研究室の工事があって、夕方から猛烈に接着剤のにおいがしだして、頭痛＋アレルギー症状が劇化したので早々に帰ってきた。昔から溶剤なにおいとか防虫剤なにおいとか芳香剤なにおいが嫌いというか嗅ぐと気分悪くなったり頭痛くなる。なわけで今日はあんまり仕事してない。まぁ実は数値計算の新バージョンの計算コードの動作チェック（発狂しそうにめんどくさかった）も完了してあとはひたすら計算機のパワーに任せるだけの段階になったので気分的にのんびりしてたりする。

そんなわけで昨日に引き続きCUDAで遊ぶ。
計算物理のぎょーかいでは速度命だったり昔々の偉人がつくったコードを使いまわしたり、線形代数の神ライブラリLapackを使うことが多いのでいまだにFortranが幅を利かせている。他の分野ではFortranなにそれオイシイの？だろうけどなぁ…。自分もメインはFortran90でありこれまで作ってきたのも大半がそれベース。なのでその一部をCUDAで高速化しようと思うと、FortranからCUDAコードを呼び出す方法をお勉強しとかなあかん。
FortranとCはわりかし互換性がよくて、もともと相互にFortran中でCの関数を呼ぶとかC中でFortranのサブルーチンを呼び出すとかができる。CUDAはCがベースなので、Fortranから呼び出すことも実は比較的容易っぽい。

CUDAに限らず、FortranとCを混在させるときの注意点
・FortranからCの関数を呼ぶときには、呼び出される側のCの関数にexternを付けておく、最後に_を付けておく。
・Fortranはそもそもサブルーチン等の引数が参照渡しなのでC側では引数をポインタとして扱わないといけない。
・Fortranは配列が１から始まる、Cでは0から始まる。Fortranでは列優先なので(a,b)の次は(a+1,b)、Cでは行優先なので(a,b)の次は(a,b+1)。

とりあえず昨日の配列の和のコードをベースに作る
Fortranで配列を準備→CUDA側に引数で渡す、GPUで計算→Fortranサイドで出力チェックてなことをやってみようと。環境はLinux(CentOS5)+ifort

CUDA側(testf.cu)


#include <stdio.h>
#include <stdlib.h>
#include <cuda.h>
#include <cuda_runtime.h>


__global__ void wa ( int *array1, int *array2, int len)
{
    int i;
    for (i=0; i<len ;i++)
     array1[i]= array1[i]+array2[i];
    return;
}

extern "C" void cudafunc_ (int* arrayH1, int* arrayH2, int* Np)
{
    int i;
    int N;
    int* arrayD1;
    int* arrayD2;
     
    size_t array_size;

    N=*Np;
    
    printf("inside the C code \n");
    printf("input array \n");    
    
    printf("input of H1 \n");
    for (i=0;i<N;i++)    
    printf("%d\n",arrayH1[i]);
    
    printf("input of H2\n");
    for (i=0;i<N;i++)
    printf("%d\n",arrayH2[i]);


    array_size = sizeof(int) * N;
    
    cudaMalloc( (void **) &arrayD1, array_size);
    cudaMalloc((void **) &arrayD2, array_size);
    
    cudaMemcpy(arrayD1,arrayH1,array_size,cudaMemcpyHostToDevice);
    cudaMemcpy(arrayD2,arrayH2,array_size,cudaMemcpyHostToDevice);
    
    wa <<<dim3(1,1),dim3(1,1,1)>>>(arrayD1,arrayD2,N);
    
    cudaMemcpy(arrayH1,arrayD1,array_size,cudaMemcpyDeviceToHost);
    
    printf("output\n");
    for (i=0;i<N;i++)
    printf("%d\n",arrayH1[i]);
    
    
    return ;
    
}

Fortran側(testcuda.f90)


program main
!test for CUDA+Fortran

!the length of array
integer,parameter :: N=8
integer :: arrayH1(1:N)
integer :: arrayH2(1:N)
integer :: i


!setting the input array
    do i=1,N
     arrayH1(i)=i
     arrayH2(i)=i*2
    end do

    call cudafunc(arrayH1,arrayH2,N)

print*,"output in fortran code"

do i=1,N
  print*,arrayH1(i)
end do

end program main

Makefileの例


all: CUDAfortran

# Define Fortran compiler
FC= ifort

CUDAfortran: testcuda.f90 testf.o
 $(FC) -o CUDAfortran testcuda.f90 testf.o -L/usr/local/cuda/lib  -lcudart

testf.o: testf.cu
 nvcc -c testf.cu

clean: 
 rm testf.o CUDAfortran

うむ、なんとか動いた。これであれやこれやをGPUで計算できるわけだがなんと10秒ルールとやらがあるらしい。GPU上のジョブの時間制限10秒ってマジっすか？10秒で計算終わってもデータ転送に時間かかってるっぽいしなんだか用途の制限多いんじゃ…？まぁ面白そうではあるので引き続きいろいろ調べると思う。

CUDAその３

2009-03-12T22:31:00.004+09:00

研究室にNIVIDA Tesla　C1060が入ったやつが導入されたのでぼちぼちCUDAをいじってみる。

TeslaはGPUのくせに、ビデオ出力の機能がない。ひたすら数値演算に特化している。とりあえず、C1060だとマルチプロセッサが240あってしかも各プロセッサに8つのコアが含まれているらしい。超並列処理とかいっちゃってるわけだ。
GPGPUのウリはこの並列化らしい。各プロセッサの8つのコアはSIMT(single instruction multi thread)という形式で同じ内容の手続きを複数のスレッドに分けて実行する（SIMDと似ているがこっちのほうがややこしい処理を並列化できるんか）。スレッドは8*4=32が基本単位になっていて、32スレッドを1ワープと数えるらしい。変なの。

母体となるコードはCPU上で実行される。コードのなかでも並列化が肝になってくるところをGPU上で行うのが普通。GPU側で処理される手続きをカーネルと呼ぶ。
ふつーカーネルで行われる手続きは並列化される。並列化の基本単位はスレッドなんだが、CUDAではスレッドの塊とさらにその塊の集まったモノを定義する。スレッドの塊はスレッドブロック、スレッドブロックの集まりをグリッドという。ひとつのカーネルに対して、ひとつのグリッドが割り当てられる。スレッドブロックのインデックスは2次元まで、スレッドそのもののインデックスは3次元までとることができる。こんなに入れ子な仕組みにするのはたぶんGPUのハード側の特性（マルチプロセッサが並列＋そのなかのコアが並列）になってるからだろう。それがメモリの階層性にも反映されている。メモリはグローバルメモリとシェアードメモリの2つに大別される。グローバルメモリはグリッド内すべてで共有であり、シェアードメモリはブロック内のみで共有される。シェアードメモリのほうが圧倒的にアクセスが早い。

CUDAのプログラムの大まかな流れは

CPUサイドでデータ処理
↓
GPU上でメモリ確保（cudaMalloc）
↓
GPUへデータ転送　(cudaMemorycpy)
↓
GPUで計算（aaa<<>>(x,y,z)）
↓
GPUからデータ転送(cudaMemorycpy)

こんな感じ。めっちゃ簡単な一次元の配列を足し算するだけのをつくってみた。バージョンはCUDA v2.0。2.0になってちょこちょこ変更があったらしく、Web上の旧バージョンのサンプルとはヘッダファイルの名前が変わってたりというトラップあり。要注意。


#include <stdio.h>
#include <stdlib.h>
#include <cuda.h>
#include <cuda_runtime.h>

#define N (32)

__global__ void wa ( int *array1, int *array2, int len)
{
    int i;
    for (i=0; i<len ;i++)
     array1[i]= array1[i]+array2[i];
    return;
}

int main (int argc,char *argv[])
{
    int i;
    int arrayH1[N];
    int arrayH2[N];
    int arrayoutput[N];
    
    int* arrayD1;
    int* arrayD2;
  
    size_t array_size;
        printf("input array \n"); 
    for (i=0; i<N;i++) arrayH1[i]=i;
    for (i=0; i<N;i++) arrayH2[i]=i*2;
    
    printf("input of H1 \n");
    for (i=0;i<N;i++)
    
    printf("%d\n",arrayH1[i]);
    
    printf("input of H2\n");
    for (i=0;i<N;i++)
    
    printf("%d\n",arrayH2[i]);


    array_size = sizeof(int) * N;
    
    cudaMalloc( (void **) &arrayD1, array_size);
    cudaMalloc((void **) &arrayD2, array_size);
    
    cudaMemcpy(arrayD1,arrayH1,array_size,cudaMemcpyHostToDevice);
        cudaMemcpy(arrayD2,arrayH2,array_size,cudaMemcpyHostToDevice);  
    wa <<<dim3(1,1),dim3(1,1,1)>>>(arrayD1,arrayD2,N);
    
    cudaMemcpy(arrayoutput,arrayD1,array_size,cudaMemcpyDeviceToHost);
    
    printf("output\n");
    for (i=0;i<N;i++)
    
    printf("%d\n",arrayoutput[i]);
    
    
    return 0;
    
}

なんか(void**)が？？だったけど、CUDAではこういう型もあるっぽい。void*はvoid型ポインタだからその拡張か？C言語自体あんまりしらんのでなかなか目新しいのが満載だ。今回つかったのは単純のためまったく並列化もなんもせんと1ブロック1スレッド。さて次は並列化を学習しよう。

面談

2009-03-12T00:29:00.002+09:00

気がつけば3月も中旬になりつつある。早い。
ここんとこなんだか忙しかった。毎回忙しいって言ってるな。ワークショップの手伝いがあったり、面談うけたり、論文の再投稿したり、OPローズダスト読み出したりしてたからだ。OPローズダストは亡国のイージスと同じ著者。文庫本で上中下なのでなかなか時間かかっている。ここんとこ夜更かししてしまう原因。

今週は月曜に、他の学科の先生に研究内容を説明して評価をうけるっていうなかなか重い仕事があったので、一週間が長く感じる。なんでそんなお仕事があったかというとGCOE関連。うちの研究室がかんでいるGCOEでは学生を研究者に養成するのに力点が置かれているらしく、プロジェクトは学生主体でチームリーダーも学生。一応数少ない日本人Dの一人なので必然的にチームリーダーとやらになってた。まぁとりたてなにもしていないんだけど（予算はおいしく使わせてもらいました）。養成事業の一環？として、外部評価を受けるというのがあるんだなー。ワークショップに駆り出されたついでに、そこに参加している研究室の教授の方にお願いして、面談してもらった。前々から知っている先生とはいえやっぱその道のプロに説明して評価されるっちゅうのは緊張する。どこがいい点なのか、今度どういう方向にすればいいかとかいろいろコメントもらえたのでよかった。なかなか学生にしてみればしんどいタスクだけどこの試みっていいんじゃないのかなぁと。自分の研究てのは意外なほど自分ではウリが分からない。外部の先生に見てもらって、どこが興味深いと思ってもらえるのかとか、業界のトレンドにどう絡めていけばいいのかとかをコメントしてもらえるのはすごくいいことだと思う。他の研究室の先生に「研究内容を聞いてください！」なんて、こんな大義名分がなきゃ言いづらいしな。研究はそりゃ自分がやってて楽しいのが何よりだけど、ひとから興味持ってもらってナンボなとこある。他人が興味持ってくれる切り口がどこかを知るのは強みになるだろう。来年度は学振あるのでリーダーはできない、補佐にまわることになる。さて誰を後釜にするんだ、勧誘せんと。

まぁそんなわけで今後の方向なんぞも模索しつつある今日この頃。GPGPUも導入されたことだしいろいろやること多いわ。

道具

2009-03-01T20:44:00.002+09:00

気がつけば三月。2月はかなり忙しかったので、なんか本当にあっというまだった。

私のどーでもいい持論の一つは、「人のやる気は道具の性能に大きく左右される」ってやつだ。PCとかソフトウェアの起動が遅いとその時点でなんかやる気失せるし、入力補完ないエディタとか打つのめんどくさくなってくる。とまぁやる気のないのを環境のせいにしてみてるわけだが。
研究はともかく、日常生活でどうもやる気が起きないことの一つが掃除だったりする。さすがに最低限掃除機かけたりはしているので、汚宅とかゴミ屋敷レベルには達してはいないが、なんとはなしに埃が多い。掃除機をエルゴラピードにしてコードレスにしてから前より頻繁に掃除機つかってる割には埃レヴェルが低下しない。絨毯が敷いてあるのが主な原因と思われる。これが掃除機の限界というものかとか勝手に納得して、埃じゃ死なんと達観してた。
が、こないだ、電球が切れたのでコーナンに行ったら、粘着テープをコロコロころがすタイプのブツが目に入ったので、買って試してみた。埃めっちゃとれるやんか…。もしかして、みんな掃除機とこれの合わせ技をつかってたのか？というかそれが常識だったのか…？家事全般の方法がよくわからん上に常識がないのでコロコロが進化しているとは知らなかった。フローリング、畳、タイル、カーペット全部これ一個でOKタイプ、しかも持ち手が伸びるタイプまである。ニトムズとかいうメーカーのは、切れ目やら粘着剤が溝上についているとか、ワンタッチでケースが開くとかなかなか工学的にイイ感じ満載だ。うーん持ってコロコロしてるだけでなんとかなるなら、掃除してやろうではないかという気にもなるってもんだ。やっぱりやる気を改善するにはいい道具を見つけることだな。
で、休みの日に仕事に飽きた合間に最近ちょくちょくコロコロしてたら、なんか家にいる時の花粉症の症状がすごくマシである。花粉症対策の一歩はまず掃除からってな文章を見るたびに掃除苦手なもんで、「けっ、原因はヒスタミンの分泌っしょ」とか思ってたけど、本当だったんだな。3月中はマメにコロコロ転がすことにしよう。

一安心

2009-02-24T22:37:00.004+09:00

これまでで最もしんどかった発表が終わった。
英語で口頭、持ち時間トータル50分。学内とはいえインバイテッドの外国人研究者が聴衆に数名。これってぺーぺーの学生にはなかなかきついシチュだとおもう。まぁ英語口頭プレゼン初体験としては、いきなり国際会議とかに比べるとハードルがだいぶひくい機会をもらえてありがたい。いつも大抵ぶっつけなんだけど、さすがに今回はがっちり原稿つくった。
30分ぐらい発表して10分ぐらい質疑応答?キンチョーのあまり時間の感覚があんまりなかったり。外国の人が質問タイムに手をあげて涙目。ものすごい気合いでなんとか聞き取ってなんとか返答。はぁーやれやれ、今年入って最も集中した瞬間だったわ。すんごい疲れました。バンケットで学内の比較的ウマーなレストラン（リーガロイヤルが経営？）でそこの名物と伝え聞く、えびふりゃーだのケーキだのをがっつりいけたのでちょっと癒されました。私が個人的に萌えを感じてひたすら実装に勤しんでいる計算手法NRGの説明に興味もってくれた人もいてうれしい限り。

なにがしんどいって開始が9時だったのが一番しんどかったなぁ…。うーむ一般企業につとめるのは無理ですな。

しかしつくづく私は人の話を聞くのがにがてだなぁと。会話とかなら聞く側に回るの好きなんだけど、学習に関しては、聞いて学ぶというのが苦手だ。学習のパターンとして、「見て学ぶ」、「読み書きして学ぶ」、「聞いて学ぶ」、「動いて学ぶ」の4パターンほどがあるらしく、脳のタイプによって違うとかなんとか。私は完璧に読み書き派で自分でちまちま参考文献読んでノート作るのが好きだし一番頭に入る。わからんくてもとりあえず写経のように書きまくるとそのうちわかってくるてなことも多。なので学会とか授業とか座ってひたすら人の話聞くのはあんまり頭にはいらないし苦痛なんだよな…。なんとかならんもんかな。

粉の季節

2009-02-18T22:00:00.003+09:00

年度末が迫ってきている。
大学の年度末は3月なんだが、予算を使うにあたっていろいろ手続きせにゃぁならんので2月中に使い切ろうっていう風潮があるので年度末感が強い。そのせいでこの時期、学会だのワークショップだのが連発で、なんか自分の事を腰据えてできないのがつらい。まぁ業績になるし、たまにいいこともあるんだけどな。

花粉の時期らしい。先週は喉のイガイガとか目のしょぼしょぼ感がすごかったんだが、週明けからなんだかちょっと楽になった。抗ヒスタミン剤が効いてきたんかなぁ。周りは逆にこの２－３日でひどくなったって言っているけれど。花粉症とはもう長らく戦っていて、毎年乳酸菌だのトマトの力だのなんとか茶を試しているがあんまり劇的な効果があるものはないな。花粉症でひどいときには熱が出たり、口が腫れたりしたことがあるので対策は結構本気だ。花粉症もアレルギーの一種なのでアナフィラキシ―にならんとも限らんしなぁ。ちなみに今年はつくし飴とやらを実験している。つくしが花粉症に効くらしい。つくし＝スギナ？なんか薬効はあるっぽいが…。意外にフリスクが喉のイガイガとか鼻のムズムズに効く。のどあめと違って糖分入ってないから歯に影響も少なそう。とりあえず花粉症には医者の薬とフリスクが最強なんじゃないかと思っている。

前衛

2009-02-11T20:05:00.003+09:00

そういえばこないだの日曜は久し振りに出かけた。後輩が大阪に帰ってきていたので、遊びに行くことに。遊びに行くといっても自分が超インドアなもので勝手がよくわからん。さてどうしよう、とりあえず中之島エリアに行くことにする。
中之島線に初めて乗った。地下鉄といいながら、表示がすべて「京阪」なので戸惑う。初乗り運賃も市営地下鉄と違うし。新型車両は一人掛けがあった。それにしても空いてる。大丈夫なんかな。

中之島は国際国立美術館やら科学館やらインドア派人間にやさしそうな施設が多い。なんだか、大学がキャンパスなんちゃらメンバーになったとかで国立の美術館・博物館に格安で入れるらしいし。ノリで国立国際美術館に行くとアヴァンギャルド・チャイナなる展覧会をやってた。アヴァンギャルド自体は死語だそーです。
もともと美術芸術の類を解さない上に前衛となるとさっぱり？？本を洗濯機にほりこんでくたくたにした塊で作品とかいっちゃったりするんだもんなぁ。嫌いじゃないけど。わけわからんなりに面白い。芸術家って変態と紙一重なんかしら（スンマセン）なキワドイのもあってびびる。
腐っても工学部としては蚊の精巧な模型みたいなん（腹部が膨らんだりへこんだりする）と超リアルじーちゃんon全自動車椅子がインパクト大だったな。ランダムネスの応用とかオートマトンみたいな手続きで図形描くみたいな結構理系ライクなのもあった。それならバンド図も芸術とかいえるんじゃないかと思う。

そういや中之島バーガーがこの近所なことを発見した。松坂牛100パーセントのハンバーガーか。気になる。

迷惑メール

2009-02-07T22:11:00.002+09:00

PCのメールにはよく迷惑メールと思しきものが来るが、これまで携帯メールにはあんまり来なかった。ここんとこしょっちゅうくる。大半が「パチンコ・スロットの打ち子バイト勧誘」いやーそれしょっちゅう詐欺特集でとりあげられてる「打ち子詐欺」ってやつやないっすか。

速効消して拒否リストに登録してるけどそれでも来る。うぜ～。メールアドレスが流出したんかと疑ったが、ToでなくBccに入っている模様＆Toのアドレスが自分とよく似てちょっと違うってパターンなんでランダム生成して自動で送りつけているんだろう。はぁ～。サルがキーボードをウン万年叩いたらシェークスピアの全作品が出力されるっちゅう話があるけどランダム生成でも結構存在するメールアドレスにヒットするもんなんだな、ランダムすげー。まぁ携帯のは単純なアドレスだしなぁ、そろそろ変え時か。

研究室のほかの人のところにも、「サイトに登録され、その後退会手続きを取られていないため＊＊万円（ｒｙ」が来たらしい。不況やと迷惑メールが活気づくのかねぇ。

そういえば、研究室のOBの方で迷惑メール対策のビジネスを立ち上げたという人がいるらしい。その人が月曜日にくるらしい。何がすごいって役職名が総帥らしい。え、代表取締役とかじゃなくて総帥なんですか、まじで？と思ったけど公的書類にも書いてきはったしマジで総帥らしい。なんかわくわくするな～

おわらん仕事

2009-02-04T23:02:00.003+09:00

いちおう講義するってのは終わった。90分はさすがに話がもたんし、相手もうんざりだろうから一時間ちょいで切り上げた。すげぇ疲れる。これ毎週何個とかこなす教官ってすごいな。教壇に立って見ると意外なほど生徒がなにしてるかよくわかる。うーん自分が聞いてる側のころは始終爆睡しまくってたので、悪いことしたなぁと。（まぁ熱心に聞かれても自分の場合はなんか怖いというか、あー間違ってるかも知れんし、ごまかしたいとこあるからあんまり見んといてって気分になるけど）なんかもうちょっと学生さんが興味持てるような手持ちのネタがあったらよかったのになぁと反省。もうちょい知ってること他人に説明できることの幅広げんとなぁ

あぁこれでしばらく人前で話すのとはおさらばじゃと思ってたら、口頭発表が回ってきた。なんでも例のG●OEで国際シンポジウムを開く必要ありとかで、それの内部の招待講演枠がのこってたとかで教授が推薦してくれはったとかなんとか経緯がよくわからんうちに決まってた、あまりにも謎い。というかいくら身内のシンポジウムとはいえその枠って学生が使っていいのか？履歴書的にはうれしいけど…。英語なんでさすがにわりとしっかりめに準備せんとやばいなぁ。英語の口頭発表を学生のうちに一回はするってのはひそかな目標ではあったんだが、なんせ話が急なもんだしあと20日（しかもその前にポスター発表も一件ある）ってことで若干焦る…。うぁぁ…

昼ごはん

2009-02-02T21:45:00.003+09:00

たまには理系な話以外を

ちょっくら母に連れられリッツカールトンの一階のスプレンディードなるレストラン行ってきた。リッツカールトンは内装が本気でヨーロッパ調なので行くだけでなんかキャラちゃうけど優雅な気分になるな。日本で無い度が大阪市内最強だと思う。はー普段いる場所とはえらい違いだわ。暖炉で火燃えてるし。なんつうか普段があまりにガサツな日々を送っているのでたまにはこういうとこ来んとなぁ。

ランチに行ったんだが、ここの凄いところはランチにデザートビュッフェがついてることだ。ランチビュッフェとかデザートビュッフェに軽食メニューがあるのはちょくちょく見るけど、ここはランチ本体は一品づつちゃんとサーブしてくれて、デザートだけがビュッフェ。ランチがうまー＆しっかりしてる分デザートビュッフェゆーてもそんなに食べられへんのが悩みやけど。リッツカールトンのケーキ有名やし、ビュッフェやのにかなり凝ったモノが食べられて大満足。レアチーズとマロンクリーム？がウマーだった。どうも私はグラスにクリームやらが詰まった系が好きらしい。二時以降はティータイムになるぽく、デザートビュッフェ単品でも行けるっぽい。

とりあえず2日分ぐらいのエネルギーを摂取してもうた気がするのでしばらく運動と粗食生活かねー。まぁ明日の仕事でエネルギー使い果たしそうやけど。

仕事が降ってくる

2009-02-01T21:42:00.003+09:00

あまり心穏やかでない週末。

金曜日にいろいろくらったもので。
ひとつには論文査読結果が返ってきた。
内容とか計算についてはさほど突っ込まれなかったので、書き方を改めればいいっぽいんだが「英語が論文としてプッ…」ってなコメントが。うーん中学生英会話な単語しか知らないので、あまりに表現が論文として拙かったらしい。英作文とか会話と違って意味がとおりゃぁいいじゃんだといかんのだなぁ。We think は論文ではつかっちゃイカンのだな、基本かもorz。とすると、これまで書いてきたやつもかなりひどかったろうによく通してくれたな、あーはずかし。それにしても、いろいろ直し方まで書いてくれてたので査読してくれた人は親切だ。自分もたまに査読を代わりにすることあるけど、スペルミスとか文法ミスとか全然気づかないわ…。日本語も英語も硬い文章書くのは苦手。論文は特に苦手っす。業績＝論文数なところあるから苦手苦手いうてたらあかんのやけどなー。結果でるのはうれしいけど論文書くのがユーウツや。

もう一個のユーウツは、2日後に講義を教官のかわりにせにゃいかんこと。90分なに話せと。ていうか金曜日に言って次の火曜にやれってそりゃきつくないっすかセンセー？とりあえず、こないだのレビュートークで準備不足のツラさが身にしみたので、一応何話すかを紙にぐだぐだ書きつけてる。研究内容話せばいいやんって感じだろうがいかんせん自分の研究内容はそのジャンルでも、おたっきーすぎるので専門外の学生に話すのはこっちも相手も苦痛すぎる。そんなわけで近藤効果を時間を稼ぐために式書きまくって紹介とかで済まそうかと。黒板に式書くのは時間稼ぎにはもってこいだ。服がチョークで汚れそうでやだなー。

ifortでのベクトル化

2009-01-28T21:44:00.004+09:00

今後、なにかと計算規模をスケールアップさせたいんだが、何が問題って計算時間である。今ので一個あたり16時間～20時間かかっていて、扱う状態の数を2倍にすれば、行列要素は4倍になるので単純に考えて４日ぐらいかかるのか。４～５日ですめばいいけど、基本的にあちこち入れ子になったループが存在しまくっているので実際にはもっと時間かかるんじゃねーかと。さすがに処理の高速化をぼちぼち図らないといけないか。

処理を早くするには並列化が基本みたいだ。どんな作業も一人ですべてやるよか何人かで手分けしたほうが早く終わるにきまってる。でも手分けするときに作業内容とかを説明したり、終わったあとに作業内容を統合する必要が出てくる。PCでも同じことで、スケジューリングやらバリア同期をしなあかんのでよっぽど並列処理に向いている部分がないと、並列化させたところでさほど恩恵があるわけではないらしい。（アムダールの法則ってやつか）並列計算に向いてるのは検索、大規模行列の線形代数な演算とからしい。

並列化で最も手っ取り早いというかなんとかつかえそうなのが、ベクトル化とマルチスレッドやろうと思ったので調べてみる。
最近ifortでコンパイルしてるとしょっちゅうdo ループがベクトル化されました　っていうメッセージが出てくる。ベクトル化ってスパコンしか関係ない話じゃねーのと思ってたのでなんか不思議に思ってたが、IA-32なCPUでは疑似ベクトル化ができるらしい。
ベクトル化ってのがそもそもなんやねんって話だが、要はレジスタを塊で使って一気に大量の数値を演算しようっていうことみたいだ。


do i=1,10
 x(i)=a(i)+b(i)
end do

みたいに演算の順番を変えようがまったく影響ないループだと、aの1～１０、bの１～１０を一次元配列に格納してその配列の和を計算するのと同等なので、レジスタの集合を配列のように使えば10回計算しなあかんところを一回で済む。10回程度だと大してなんも変わらないけど一万回が一回に化けるとだいぶ違う。地球シミュレータとかSXシリーズとかのスパコンはそれ専用のベクトルレジスタを大量に積んでてそれゆえ行列演算が異様に早いっぽい。（とはいえベクトル型のスパコンを作ってるのはいまやNECのみらしいが）こういう並列化の仕方は、SIMD（Single　Instruction Multi Data）というらしい。

実はIA-32なアレにはこのSIMDな処理ができる部分が存在するらしい。CPUのレジスタというと32ビットの汎用レジスタeaxなどなどが使われてるイメージだが、実は128ビットのレジスタxmm0～xmm7というのが存在するそうな。128ビットなので倍精度実数なら2個格納できる。浮動小数点演算は通常、名前も浮動小数点計算ユニットであるFPUとかゆーst(0)～st(7)レジスタを持つところで行うらしいが、xmmレジスタをつかうこともできる。xmmレジスタをまとめて使ってベクトル化したような処理をおこなうための拡張命令がSSE/SSE2/SSE3（SSEはStreaming SIMD Extensionsらしい）
SSEはもともとインテルがつくったのでインテルのコンパイラはこれを使えるのが多い。ifortで出てくるベクトル化しましたってのは実はこのSSE2やらSSE3やらをつかってxmmレジスタで疑似ベクトル化してるということみたいだ。実際、－Sオプションで、アセンブラコードを吐かせてみると、movsd だの %xmm0だのそれつかってるぽい形跡がちらほら。ちなみにmovsdはmovのSSE2バージョンぽい。やはりインテルのコンパイラは頭ええんだな。ところで、DOループ中にprint,read,writeがあるとこのベクトル化ができなかったりする。ベクトル化の詳細レポートがほしい時には-vec_report3とするとよし。

長年の疑問

2009-01-26T21:23:00.003+09:00

前々から思うことなんだが、なんで日本の大学は留学生超優待なんだろう。
こないだも円高とかリーマンショックで留学生の生活が大変、救いの手をって新聞記事になってたけど、救わなあかんのは留学生だけちゃうやろうと。いっとくが国費で来てる人々はかなりいろいろ手当もらってて余裕アリだぞ。

まぁこういうこと外で言うとなにかと叩かれるし、自分が困ってるわけでもないのに嫌な奴と思われそうなのでなにもいえないけどな。

自分たちによこせとか彼らにやるなってわけじゃなくて、なんだかリソースの配分にあまりにも偏りがあると、国内の学生のやる気を削ぐよ、それって結果的に損なんじゃないのってことが言いたいのだ。

国内の学生への補助は手薄だと思う。奨学金はほとんどが返済義務アリだし。企業や国が出す、返済義務なしの奨学金はほとんどが留学生対象。そりゃDに進学する日本人が少ないわけだわ。博士課程の日本人率が少ないって問題だと思うんだが。進学する人数が少ないからって留学生補充というのは本末転倒だろ、たいてい皆本国に帰っちゃうし。卒業後に日本の企業や研究機関に定着する率とか考えても、日本人の学生への補助も増やして進学しやすくするのは、投資に見合うリターンが国にあるかと。優秀な子でも、金かかるし、早いこと稼ぎたいですって皆就職してしまう。ある程度の数が基礎研究に居ないと、何十年後かに困ると思うぞ。

私費留学生で大層苦学してる人もいるのでその人たちには助成が必要なのはわかる。だが留学生＝まじめ・苦労人、助けてあげなきゃ一辺倒ではないと思う。なにがほんとうんざりするって、大学外の少し年長な方々は「日本の若い学生はさぼってばっかりで留学生のほうがまじめなんでしょ」（ソースはうちの親族）てな発言をされる。それ実情見てから言えよ、まじでやる気そがれるわ。政府高官の方々同じような意見やからのいまの政策なんやろな。頼むから日本人の困ってる子も助けたってくれよ。

たしかに留学生のなかにゃスゲぇ優秀で成果バンバンあげてしかも性格もイイてな超人もいるが、そういう人ばっかでは無いんだよな。なんか長年、あの研究室にいると嫌な部分もいろいろ見てしまって「国際」ていう言葉にアレルギー持ちになる。日本の大学にいるはずなのに聞こえてくるのがよくわからん言語のほうが多いってのは本当に憂鬱。大学院行って外国嫌いになりました。だめじゃん。個人単位ではいい人なんはわかるが、あの環境にいるとフィルターが頭の中にできちゃって狭量になってる自分がいて嫌だ。

寒すぎ

2009-01-25T21:53:00.003+09:00

久しぶりに日記的な文章
とにかく寒すぎる。風邪をひいた。
基本的に寒いの嫌いで、寒い時期は体調がよろしくない。こーなんちゅうかヤバくはないけどうっとおしい不調に悩まされることが多い。しもやけとかなんとなく風邪とか。血行不良とか冷え症がすべての原因なんだろうなぁ。起きてる時間大半座って作業してんだからしょうがないけど。
そろそろ学会の準備を始めないとあかんかなと。データ処理に追われる。はよ終わらんかなぁー、遊びてー、データ処理は単調作業ですぐ飽きるので嫌いだ。プログラム作って計算結果出た時点で満足してしまってさて次のことしたいなぁってなる。でもデータ処理をせんと「研究成果」として扱われる発表とか論文とかのネタにならんわけで。うぅジレンマ。シェルスクリプトで大半なんとかしようと思ったけど、細かいチューニングとかExcel使う部分やグラフ・図の作成はどーしても手作業が必要。

定例の3月の以外に2月にも身内な学会とはいえ一応国際な学会がある。2月の会場が自分とこの大学でほんとよかった。出張はめんどい。
3月の場所が東京なのでホテルを確保した。東京でホテルを安く上げようとケチるとロクなことがないので、今回はわりといい目のを予約した。地方やと一泊7000円でだいぶ快適に過ごせるのに…。やっぱ学会は地方がいいです。そして、こないだ使った丸の内メトロポリタンがちょー快適で「あぁやっぱええホテルってええわぁ…」と味をしめてしまったのもある。さすがにメトロポリタンクラスのは無理やけど、3泊せにゃならんので健康と睡眠には代えがたい。睡眠と水分をしっかりとらんとダメな体質ぽい。

あぁとりあえずさっさとグラフつくらんと。Excleいじるだけってつまんねー。データ処理よかOpenMPで計算の並列化うまいことできるか試してみたいんやけどなー。せっかく自由にいじれるクアッドコアが手元にきたんだから並列化を試さん手はないと思うのだ。

Fortran備忘録：ファイル名とかいろいろ

2009-01-24T22:19:00.011+09:00

数値計算でよくある「繰り返し計算の各回での結果を、第何回目であるかの数字をつけたファイル名にして出力する」とか「あるパラメータでの計算結果をパラメータの数値をファイル名に入れて出力させる」とかのFortranでのやり方まとめ。Fortranは文字列と数値間のキャストの組み込み関数がないのでメンドイ。

1.まず、ファイル名は文字列変数で指定することができる。
2.文字列変数の結合は//で可能
なのでたとえば
character(30) ::param,filename
ていうふうに文字列変数を用意しておいて


open(8,file=param//filename)

なんてのも可能。

3.integer,real,double型の変数を文字列変数に変換やその逆をする場合には内部ファイルを使うとよい。
内部ファイルというのは、入出力先をプログラム中の変数にしたものである。
たとえば


program main
implicit none
integer :: i,j
character(10) :: a,b
i=3
a="5"
!b に　i の値を文字列として代入
write(b,*) i
!aの値を数値としてjに代入
read(a,*) j

print*,b
print*,j

end program main

とすれば、３　５てな出力が出てくる。

4.ただ出力を見てもわかるように、bでは3の前にいらん空白が入っているし、これをそのままファイル名で使うと、１０文字分の領域を確保しているせいで後ろにも空白が入ってきてなんとも間抜けたファイル名ができる。そんなときには文字列の組み込み関数のTRIMとADJUSTLを組み合わせるといい感じ？TRIMは右側の空白を取り除く、ADJUSTLは文字列を左詰めにする関数。


open(10,file="a"//TRIM(ADJUSTL(dis))//"total")

的な。

5.ファイルに埋め込みたいのが整数ならこれでいいんだが、実数型だと1.000000.txtみたいにいらん長さのファイル名になってしまう。
対策例１
・文字列の部分配列を取り出してくる
文字列は配列のように、1文字目から3文字目みたいに切り出してくることができる。
たとえば
a="abcdefg"
なら
a(1:3)
は"abc"である。
実数を文字列に変換した場合には最初に符号とかを入れる関係で頭に2文字ほど確保されているらしく、a=1.00000000から1.0だけを取り出すという風に切り出す場合にはa(1:5)という風に指定すればいいみたいだ。

対策例２
・何文字目までが有効桁かを探る→切り出し
例
valはcharacter(30)
nonzero,iはinteger


 nonzero=0
 do i=1,30
 if ( val(i:i) > "0") then
  nonzero=i
 end if
 end do
open(15,file="pot"//TRIM(ADJUSTL(val(1:nonzero))))

CUDAをコマンドラインから使えるようにしよう

2009-01-20T22:31:00.004+09:00

VisualStudioが意味フなのでWindowsでもコマンドラインから使えるようにしたいなぁと。どうせ本番はLinux機の可能性100パーですし。
ちなみにOSはWindowsXPです。
コマンドラインから使えるようにするには環境変数をいろいろいじらんとあかんみたいだ。
先人の知恵がつまった参考文献3つで言っていることが違うぜ、どうしよう
http://tech.ckme.co.jp/cuda_inst.shtml
http://exth.net/~ohshima/wordpress/2007/10/30/5/
http://chihara.naist.jp/people/STAFF/imura/computer/OpenGL/cuda1/disp_content

環境変数INCLUDEなるものが見当たらないので、まず真ん中の人がやっているのと同様に、システムのPATHをいじることにする。
とりあえずPATHに

C:\Program Files\Microsoft Visual Studio 8\VC\bin;C:\Program Files\Microsoft Visual Studio 8\Common7\IDE;C:\Program Files\NVIDA Corporation\NVIDA CUDA SDK\common\lib

を足してみる。
で、昨日ためして動いたbandWidthTest.cuをコンパイルして実行すべく

nvcc -o test.exe bandWidthTest.cu -I"C:\Program Files\Microsoft Visual Studio 8\VC\include" -I"C:\Program Files\NVIDIA Corporation\NVIDIA CUDA SDK\common\inc" -L"C:\Program Files\Microsoft Visual Studio 8\VC\lib" -L"C:\Program Files\NVIDIA Corporation\NVIDIA CUDA SDK\common\lib" -lcutil32

てなバッチファイルを作って投下。
test.exe無事生成
が、cutil32.dll見つからないというエラーメッセージ出て動かず

そこで最後の文献を参考にして
C:\Program Files\NVIDIA Corporation\NVIDIA CUDA SDK\common\lib\cutil32.dllをC:\Windows\system32へコピーすると動いた…いいのかこれで？
たぶんファイル構造をつらつら見るに、
C:\Program Files\NVIDIA Corporation\NVIDIA CUDA SDK\bin\win32\Release
にPATH通しといても同様に動くはず。このwin32\Releaseやwin32\DebugはVisualStudioのコンパイラ経由でDebugなりReleaseなりするときに必要なDLLとかを全部ほりこんであるぽいのでPATH通しとくのがいいんじゃないかなぁ。やっとこ

なんにしろ無事うごいたっぽいのでCUDAが出始めのころに試して記録を残してくれた先人たちに感謝感謝。

CUDAを試す

2009-01-19T22:18:00.004+09:00

どうやらCFSは順調に動いたようで磁場下での1粒子励起スペクトルが綺麗に出た。今までの方法と段違いにきれいな結果だ。ひじょーに満足。作るの苦労した分、上手くいくとめちゃうれしいわぁ～。これで夏の国際会議のネタが確保できる見通しもついた。

ほっとしたので、自宅機で噂のCUDAを試す環境を作ってみる。CUDAって何て読むんだろ。きゅーだ？くだ？
研究室でも最新版Teslaを試そうぜっていう話がでてたんだが、取引先の業者さんがまだ性能評価終わってないらしい。Teslaは結構熱がでるとかでファンをどうするとかGPUコンピューティングするなら帯域がどーとかそれなりに相性とか考えなあかんぽいとか。評価終わるの待ってて予算執行がぎりぎりになるんもあれだし、導入はたぶん来年度かねぇ。それまで自宅のにはせっかくCUDAが動かせるGeForce9800GTが入っているので試しにいじってみようかと。GPUつかったら画面が出力できないとかOTLなことになるんじゃないのとかいろいろ無駄におそれたがべつにそんなことはないらしい。
http://tech.ckme.co.jp/cuda_inst.shtml
に詳しくインストールの仕方が書いてあるのでとりあえず試す。
ドライバは対応してるのがもともと入っているのでいじらない（グラボのドライバってなんかいじってえらい目にあった知人の話聞いて以来怖い）
NVIDAのページhttp://www.nvidia.com/object/cuda_get.htmlからcuda toolkit cuda sdkのバージョン2.0をダウンロード。バージョン2.0を動かすためにはWindowsではVisual Studio 2005C++が必要らしいのでそれもMicrosoftのページからダウンロード。最新版は2008らしいがCUDAの2.0は2005じゃないといかんらしい。あとDirectXつかうサンプルもあるらしく、DirectX SDKも入れといたほうがいいらしいが、入れ忘れた。まぁとりあえず放置。

CUDAについてはtoolkit →　sdkの順でインストールすべしとのこと。インストールはデフォルト通りにしといた。後でみたらいろいろパスが増えとる。参考資料によると、CUIでコンパイル、ビルドするにはパスをいろいろいじる必要があるらしいが、今日のところはとりあえず動くかどうかをチェックしたいのでまたとりあえず放置。
デフォルト通りにすると
C:\Program Files\NVIDIA Corporation\NVIDIA CUDA SDK\projects\simpleTemplates
下にsimpleTempletes.slnなるものあり。
とりあえずこれが動くかどうかでCUDAが無事インストールされたかどうかがわかる模様。
これをクリックするとVisual Studioが立ち上がる
ビルド→ソリューションのビルド
なんか一応正常終了したっぽい
Visual Studio使ってないもんでよくわからんのだが、たぶん緑の矢印が「実行」てのは共通だろうとあたりを付けてとりあえずそれを押す（これデバッグらしいが）

とりあえず動いたっていう解釈でいいのか…？謎いわ。やっぱLinuxのほうがgccで直コンパイル＋ビルドできて見通しつけやすそう。sampleにMonteCarloあんのが気になる
結果がWebにあがってたbandWidthTestも試す

出てる項目大体いっしょぽいなぁ。ビルド時に警告でたけど動いてるってことか、よーわからん。CUDAなかなか手ごわそうじゃ。

CFSバージョン起動

2009-01-18T22:03:00.003+09:00

最近このブログ内容が偏ってると周囲の知人からは不評を頂いてたりする。しゃーねーやん、日々そーゆーんしか扱ってないんだもん。ここんとこ毎日家と研究室を点と線やし。google analyticsをみると、Fortran関連の検索で来ている人が多い模様。Fortran関係者の役に少しでもたってるとよいのぅ。Fortranのちょっと知っとくと便利情報とかエラーメッセージの解読方法とかあんまりWebで載ってないもんなぁ。わかったらまたぼちぼち書いていこう。

さてそんな無機質な日々の甲斐あって、なんだか数値計算の新バージョンがなんだかそれっぽい答えを返し始めた。ぬかよろこびでないことを祈る。数値計算ってテストケースがうまくいったーと思ったら本番投入したらわけわからんエラー出ることよくあるし。

新バージョンてのは12月の頭ぐらいから「学振結果待ちうぇぇ辛い…」な気分をごまかすべく作り始めたComplete Fock Space略してCFSをつかったNRGの計算コードでなんつうかまたマニアックな代物である。まぁそうでもなきゃこれに手だしてないからある意味学振に感謝かな。NRG自体が非常にアレゲというか人口少ない（一不純物ソルバーとしては史上最強のはずなんだが）計算手法でその中でもまたマニアックな手法らしく、似ているDM-NRGについては結構取り上げられてるけどCFSはあんまりというかほとんど見ない。DM-NRGは状態を二重にカウントしてて誤差でるってな噂あるんだけどなぁ。NRGでは計算機で計算可能にするためにステップごとにあまり大事じゃなさそうな状態を捨てていく、つまり繰りこんでいく。これが頭イイなやり方の反面いろいろ問題を起こす。CFSはこの捨てちゃった分も拾ってきて計算できるので精度がものすげー向上する＆磁場が入っても大丈夫らしい。磁場はいったときにSTMのスペクトルがどうなるかを計算しようとおもったらこいつが一番だろうというわけでちまちま作ってたわけだ。大量のデータ扱わなあかんしめっちゃ実装するのめんどくさかった。モデル計算のくせにメモリ数G使うとか。とにかくなんとか動くもんができてやったーてな気分。なんかもうちょっと高速化できそうな方法も思いついたんで次はそれだな。

RedmineをFedoraで動かす

2009-01-17T20:33:00.007+09:00

サーバー立てても文章置いとくだけじゃあれなので、なんかサービスを入れようと。とりあえず研究室で役立ちそうなものということで、Redmineをいれたいなぁと。実機でやってミスってはまずいので家のFedora10 on VMware で試す。LinuxにはInstantRailsがないのでメンドイ。大学行くまでに忘れそうなんでメモ

まずruby on railsが動く環境にしないといけないので、それの設定。
ruby周辺のインストール


yum install ruby
yum install ruby-devel
yum install icb
yum install rdoc

RubyGemsをhttp://rubyforge.org/frs/?group_id=126からもらってくる。適当な場所で展開して


ruby setup.rb
gem install rails

Ruby on rails は最近のはデータベースにsqlite3を使うらしいので


gem install sqlite3-ruby

Redmineをhttp://rubyforge.org/frs/?group_id=1850
でもらってきて、展開。展開したフォルダ中の/config/database.yml.examle をdatabase.ymlとしてコピー。
インストールの手引きにあるように、database.ymlのproductionの部分を


production:
  adapter: sqlite3
  dbfile: db/redmine.db
  timeout: 5000

として


$ rake db:migrate RAILS_ENV=production
$ rake load_default_data RAILS_ENV=production
$ script/server -e production

でWEBrickが起動、localhost:3000で待ち受け状態になる。デフォルトランゲージの選択でjaにしていても、最初のadminでログインするとすべて英語…。管理のところの言語設定だけでなく、マイアカウントのところの言語設定も日本語にしてやる必要あり。

せっかくApache起動しているのにもう一個Webrick起動してポートも別に開放するのがなんか微妙？なのでApacheと連携させるってのを一度試してみんとなぁ。

あたらしいものは面白い

2009-01-16T22:48:00.004+09:00

今週はなんだか忙しかった。
というのも新しい計算機（Quad、i7は高くて買えず。）がきてその設定やらなんやらの作業が入ったから。クラスターではなく大きな筐体なので計算機室におけずに自分の机の横にでかいのが陣取ることになった。ディスプレーも買ったのでディスプレーに前後から挟まれる格好になっている。電磁波とか大丈夫なんやろか…。やっぱワイド液晶万歳。国産のちと高い目にしただけあってわりかしきれいで満足。いままで使ってたD●LLの付属品は比較すると涙目すぎる。

OSはCentOS5にしてもらった。業者の人が一通り設定してくれた。Linuxよくわからんからありがてぇ。念願の管理者権限げと。CentOSは扱ったことがないのでビビりながら使いだす。マウスの右ボタンで端末をディレクトリの中から呼べるのが地味に便利だ。しかし、FTPクライアントもっとええのないのか？デフォルトのgFTPはWindowsで使ってるWinSCPに比べて超使いにくい。WinSCPのいいところは転送時に新しいものだけを上書きとかミラーリングとかのオプションがある＋いつも使ってるエディタでリモートのファイルを編集できるところ。なんでLinuxにこれに匹敵するもんがないのか謎い。FileZillaがましっぽいけど、編集を直にはできないっぽいなぁ…。ファイルブラウザも付属のがビミョーに使いにくいしなんかええのないかなぁ。とまぁ不満はちょっとあれど、コンパイラとか満載やし、X-Windowやし、いろんなもんのインストールがrpmとかで楽勝なLinuxは研究用としてはWindowsより使いやすいかも。大学のサイバーメディアセンターがなぜLinuxからWindowsに切り替えるのかまったく理解不能。とりあえずWineもつかってみたいな。

Mathematicaの設定もすんだ。公費で買っている分、使い方を研究室の人間に知らせなければいけない。何度も聞かれるのが面倒なので、サーバーたててWebページで公開しとくことにした。Webページつくるん結構面白いけどいろいろよくわからん…。週末はこれと、いま作ってる計算コードいじって終わりそやなぁ。GPUコンピューティング（CUDA）にも手を出そうかって話も出つつあるしどうなるんやら。おもろいけどなんか妙に忙しくなりつつある。

BLAS+LAPACK

2009-01-12T21:35:00.005+09:00

VMwareのFedora10にBlas+Lapackを導入してみるテスト。
BlasとLapackってのが何かというと、線形代数の計算用のライブラリで、行列の対角化とかをサクサクやってくれる賢い子。数値計算には欠かせない。Fortranで書かれているのでFortran使いにはなんとも使いやすい。（Cからも比較的簡単に使える）ソースはF77だけどF90からもなんも変更なしに使える。

ここを参考に：http://www.obihiro.ac.jp/~suzukim/masuda/ifc_lapack31.html

ソースからビルドしなあかんような噂を聞いていたのできつそうやなぁとか思いつつまずはNetlibからソースを落としてくる。展開しててけとーなところにフォルダを置いて、まずはmake用のmake.incをいじる。フォルダ内にmake.inc.exampleがあるので
cp make.inc.example make.inc
としてmake.incを作ってやる。intel Fortran Compire用はこんな感じに


####################################################################
#  LAPACK make include file.                                       #
#  LAPACK, Version 3.1.1                                           #
#  February 2007                                                   #
####################################################################
#
# See the INSTALL/ directory for more examples.
#
SHELL = /bin/sh
#
#  The machine (platform) identifier to append to the library names
#
PLAT = _LINUX
#  
#  Modify the FORTRAN and OPTS definitions to refer to the
#  compiler and desired compiler options for your machine.  NOOPT
#  refers to the compiler options desired when NO OPTIMIZATION is
#  selected.  Define LOADER and LOADOPTS to refer to the loader and 
#  desired load options for your machine.
#
FORTRAN  = ifort 
OPTS     = -O3
DRVOPTS  = $(OPTS)
NOOPT    = -OO
LOADER   = ifort
LOADOPTS = -L /opt/intel/Compiler/11.0/074/lib/ia32
#
# Timer for the SECOND and DSECND routines
#
# Default : SECOND and DSECND will use a call to the EXTERNAL FUNCTION ETIME
TIMER    = EXT_ETIME
# For RS6K : SECOND and DSECND will use a call to the EXTERNAL FUNCTION ETIME_
# TIMER    = EXT_ETIME_
# For gfortran compiler: SECOND and DSECND will use a call to the INTERNAL FUNCTION ETIME
# TIMER    = INT_ETIME 
# If your Fortran compiler does not provide etime (like Nag Fortran Compiler, etc...)
# SECOND and DSECND will use a call to the INTERNAL FUNCTION CPU_TIME 
# TIMER    = INT_CPU_TIME
# If neither of this works...you can use the NONE value... In that case, SECOND and DSECND will always return 0
# TIMER     = NONE
#
#  The archiver and the flag(s) to use when building archive (library)
#  If you system has no ranlib, set RANLIB = echo.
#
ARCH     = ar
ARCHFLAGS= cr
RANLIB   = ranlib
#
#  The location of the libraries to which you will link.  (The 
#  machine-specific, optimized BLAS library should be used whenever
#  possible.)
#
BLASLIB      = ../../blas$(PLAT).a
LAPACKLIB    = lapack$(PLAT).a
TMGLIB       = tmglib$(PLAT).a
EIGSRCLIB    = eigsrc$(PLAT).a
LINSRCLIB    = linsrc$(PLAT).a

っていじるのは

FORTRAN = ifort
OPTS = -O3
DRVOPTS = $(OPTS)
NOOPT = -OO
LOADER = ifort
LOADOPTS = -L /opt/intel/Compiler/11.0/074/lib/ia32

だけですが。OPTSやNOOPTは最適化のオプションなのでOがゼロではなくオーなあたりがちょっと罠な気がする。
で、
make blaslib
make lapacklib
make tmglib
と3回いろんなブツをコンパイルする。結構大量にコンパイルするのでCore2Duo上位機種の実力とやらをみせてもらおうか、って感じだ。まぁなんとかこれで必要なblas_LINUX.a,lapack_LINUX.a,tmg_LINUX.aができたので、あとはこれをlibってなフォルダにコピーして名前変えてシンボリックリンクとやらを貼るかーってところで、Fedoraならyumでblas,lapackをインストールできることを知る。いままでの苦労は…orz
というわけで、CentOSやFedoraなら

yum install blas lapack

で一発です。まぁソースからビルドする方法は管理者権限がない場合には必要なので知っといて損はないか…。
yumでインストールすると、/usr/lib下にlibblas.so.* とか　liblapack.so.*ができる。
コンパイル時にリンクさせるのにフルパスはめんどくさいので

ln -s ~/usr/lib/libblas.so.3 ~/usr/lib/libblas.so
ln -s ~/usr/lib/liblapack.so.3 ~/usr/lib/liblapack.so

とすると、GNUのライブラリの規則ってもんが働いて、コンパイル時に
ifort -llapack -lblas test.f90
と2つのオプション -llapack と　-lblasをつけるだけでいいぽい。ちなみに/usr/lib/以外のところにライブラリを作った場合には
-L ~/lib
のように-Lオプションでライブラリの場所をしめさないといけないらしい。

さてこれでlapack使ってるNRGのコードが家のPCで動くぜ～ってmakeしてみるとコンパイル＋リンクはうまくいって実行ファイルができるのに動かした瞬間「強制終了」ってコンソールに出る。gdbのお告げによるとSIGKILLキター→強制終了らしい。SIGKILLはプロセス強制終了のときに出すシグナル…。どうやらメモリを大量に喰うプログラムなためOOM Killerに叩かれた模様。まぁ倍精度で4400＊4400の配列を10個とか用意したら１G超えてきついわな、これOKな大学の計算機ってすげーんだなぁ。配列サイズを小さくしたら無事実行された。これで家の環境だけでプログラム作る→テスト→バグ取りが完結するようになったので作業効率は上がるだろう。とりあえずDDDっていうgdbと連携するGUI付きデバッガもyumで入れてみた。今まで見たく気合いと勘だけでバグとるのもどうかと思ったので次はデバッガの使い方を…。

Fedora10+ifort

2009-01-11T20:27:00.003+09:00

VMware上Fedora10にインテルのフォートランコンパイラをインストールするテスト。
参考資料：http://www.obihiro.ac.jp/~suzukim/masuda/ifc101_install.html

インテルはFortranコンパイラも非商用の目的（研究とか）に使うのはタダでダウンロードさせてくれる。条件に"accept"で、メールアドレスを登録すると、ライセンスキーが送られてくる。
VineLinuxの場合は、ほんとお手軽にsuでルートでログインして、./install.shってするだけでできたんだが、Fedoraはこうはいかない。Fedoraの場合はセキュリティ対策の一環としてSELinuxというものを導入しているらしく、root権限だけではコンパイラのインストールとかができないらしい。インストールするためには、SELinuxをPermissiveモードとやらにする必要があるらしい。モードの切り替えは/etc/sysconfig/selinuxファイルのSELINUX=enforcingをSELINUX=permissiveにするか、setenforceコマンドで一時的に切り替えるぽい。

あとなぜか、libstdc++.so.5 ライブラリがないよという問題が発生。こーゆーC++のライブラリってFedoraは標準じゃないのか？？必要なライブラリを調べるために
yum provides libstdc++.so.5
とするとどうやら　compat-libstdc++-33-3.2.3-64.i386 ってのがいるっぽかったので
yum install compat-libstdc++-33-3.2.3-64.i386　としてインストールすると解決。

インテルのコンパイラはFedora9までが標準っぽく、それよりあたらしいFedora10ではアーキテクチャーが判別できないとかglibcのバージョンがわかんねぇとかのメッセージがでたが、スルーしてインストールしても無事っぽい。
で、インストールした現行バージョン(ファイル名：l_cprof_p_11.0.074_ia32）の場合
.bashrcの最後に
. /opt/intel/Compiler/11.0/074/bin/ia32/ifortvars_ia32.sh
とすることで、ログインするたびにパスが通ってコマンドifortでコンパイラを使えるようになる。

とりあえずFedora10にもGeanyとか入れてFortran用の環境を作ってみる。なぜかコマンドラインからGeany起動した場合はifortでのコンパイルがうまくいくんだがショートカットから起動すると無理。謎、たぶんパス関係？まぁいいや、ifortが手に入ったのであとはlapack+blasをなんとか積みたいところ。

Geany(FortranのIDEにもなるらしい)

2009-01-08T22:08:00.003+09:00

なんだか初期の選択をミスったせいで、ずーっとFortranをつかっているわけだ。古代な言語とはいえ物理の数値計算には使い勝手がかなりよいしとにかく速いのでC++とかに乗り換える理由もさしてなくそのままずーっと使っている。
唯一不満なのが、IDEがあんまりないこと。情報の実習のTAやってて「Eclipseすげー」だったのでFortran用のCDT探したが、リリースされてるPhotranはうまく動かなかった…。Silverforestが出してるPlato3も一応IDEなんだが見た目やコンパイラ自体の仕様が？なのでどうも…。
今日偶然に「おぉこりゃ満足できるレベルかも」なフリーのIDEを見つけた。Geanyという。フリーなIDEはそりゃたくさんあるんだが、マイナー言語Fortranをサポートしてくれてるのはなかなかない。Geanyはコードの色分け、サブルーチンとかファンクションとかグローバル変数とかのシンボルの表示、変数名の補間、コードの折りたたみが一通りそろっている。シンボルブラウザと折りたたみできるのがうれしい。（IDEならとーぜん持ってる機能なんだろが、Fortranでこれがあるのなかなかなかったりする…）動作画面はこんな感じ

見た目もなかなかよく、アイコンも「ジーニー」って名前だけあって魔法のランプなのでなかなか可愛かったりする。ちゃんとプロジェクトとして管理する機能もあるし。コンパイラとしてはFortranはgfortranがデフォルトのようだが、ビルド→インクルードと引数の設定　をいじればg95でも実行できた。でもgfortranつかったほうが、エラー出た時にその位置を知らせてくれるのでいいかも。ブックマークの機能がないのがちと残念だけど、使いやすい。たっかいVisual Fortran買わなくてもフリーでなんとかなりそうやん。

いまさらながら学振面接体験記

2009-01-06T21:56:00.004+09:00

ついに来た。
日本学術振興会特別研究員、採用内定。やたーっ！！
3度目の正直ってやつだ。学振はこの3年何度も鬱状態に叩き込んでくれたからなぁ。
面接までもつれ込んでなかなかにしんどかった。面接もコメントが「とにかく凄い自信だ…www」だったから正直どーなんって感じで不安だったし。（直後はすごいぐだぐだ考えてたくせに、25日に結果公表されてることを全く知らなくて、正月に教授からメール来て知った。年末はそれ以外に頭悩ますことがいっぱいあったからなぁ）

とにかく、発表練習見てくれたり、しんどいときに相談のってもらったりなにかと付き合ってくれたりした皆さんに感謝です。

DC2で面接てのはボーダーラインで、面接受けてとるのは少数派らしいので今後受ける人用に一応体験記を書いておこうかと。（書類については、面接ってこたぁ微妙だったのかってことで自信ないので書きません。）

時期：今年は11月末だった。午前と午後で合格率に差がでるとか、午後なら一番じゃないとヤヴァイとかいう噂があるけどデマだと思う。とりあえず自分は午後のしかも一番じゃなかった。順番が前後の人にリサーチしてみた結果、午後組は地方から来ている人多数っぽい。たぶん領域ごとに日が決まっていて、その中で地方からで来るのに時間かかりそうな人を午後にしてるんだと思う。同じ研究室でも領域違うと日が違ったし。

ポスター：とりあえずボードのサイズが要綱に書いてあるだろうけど、横長が絶対おすすめ。マグネットは透明で長細いのを用意してくれてるので余白はそこまで気にしなくてもいいんじゃないかなぁ。字は思ってるよりでかいほうがいいと思う。イラストレーターで作って、小塚ゴシックで、現場で貼ったときに、50ｐでもちょっと小さめかなぁとか思った。内容は理論の式とかぐだぐだ書いてもよっぽど専門近くないと分かってもらえなさそうなので、この半年でどんだけ業績作ったかを強調して見るテスト。

面接官：大人数。10人ぐらい居てはった気がする。たぶんうち数名は学振の事務の人。実は、審査官がだれかという情報はウェブに上がっている。（私はそれに気がついたのあとからで、あちゃーって感じだった。ちゃんと確認しなきゃ…）領域の中でも、何個かに分類されてそれぞれに審査官がつくので、領域あたり6人とか結構な人数になる模様。で、その細目が違うと、専門も全然違うので、自分にしてみりゃ聞きなれてて日常語な単語も向こうにしてみれば意味フということも十分考えられるので気をつけたほうがいいと思う。（知らなかったのでその辺考慮せんとはしょりまくってしまった…。）

面接：なんせ発表時間4分で原稿見れないので、セリフはかなりしっかり覚えないとやばい。時間内に説明終えると、ブザーは鳴らないみたいです。質疑応答は一番専門が近い先生2人ぐらいが質問してきはります。雰囲気とかはその先生のキャラ次第…。まぁorzなこと言われても通ってることもあります。たぶん、答えにつまったり、おどおどせんと、虚勢でもいいから「自分はこう思う」てな答えを自信ありげに返すのがいいとおもわれる？

その他：結構みんなはやくから集まっているみたいだ。結構直前に行ったら、すわる場所なかった。

終わってみたらそんなもんなんだけど、渦中にいるときは、ものすごく緊張するわ鬱はいるわしんどいよねぇ。だいぶ先だけど今年11月、このしんどいのに立ち向かう皆さんに幸ありますよう。とにかくなんらかの経験値はつくと思います。

Lyx再戦

2009-01-04T21:44:00.002+09:00

論文がいやだぁぁぁぁ。でも明日の仕事始めである程度進めてたよんてのを示さんと上からの御怒りがきそうなんでいまから焦って校正するか。
Texは仕上がりがコマンドからは想像もつかないのが問題。なのでLyxを懲りずに試す。
Lyxこないだ試したのがうまく動いてなかったのは、babelってのの設定とインストーラーの選択をまちがったからぽい。ここ→http://www.int.otaru-uc.ac.jp/lyx-howto/install/win
に書いてあるようにWindowsへの導入をおこなうと、まぁサンプルは普通に動いた。LyX-160-4-15-AltInstaller-Small.exeてのをインストールしないといけなかったのだな。

なぜか、サンプルがある場所（/example/ja）と違うところにファイル置いてdvioutでみようとするとうまくいかないのが謎すぎる。パスとかの問題？
Lyxは式とかが出力される状態で出てくるのはありがたいんだが、引用文献の番号の書き方とか細かい部分をプリアンプルでちょっといじるってのが難しそうな気がする。ソースに直で書けないし。まぁいろいろやり方はあるんだろうけど。なのでしばらくは下書きをLyx、ソース吐かせて仕上げをWinshellな感じでいいかなぁと。

どうやら明日大学行くと、とある重要書類が届いているみたいだ。うぁぁなんか胃が痛い

メモリの種類とか（いまさら）

2009-01-03T20:12:00.006+09:00

昨日は久し振りに研究室の同期と会った。iPhone持ってたので見せてもらう。iPhoneスゲー。最近ブログ更新してねーよねって聞いたらむしろTwitter派らしい。Twitterなぁ、導入したらすごい気が散りそうでなんか手を出せないわ。気晴らしになんかいいのないか聞いたら、Travianがおすすめらしい。興味もったんで調べてみたがなかなかに複雑かつ難しそうじゃ…。
ネット上のゲームなりSNS的サービスなりってどうもあんまり知らん相手とコミュニケーションとることを求められることが多い気がして苦手だなぁ。人脈広げたりさらっと知り合い作るにはいいツールなんだろうけど、携帯メールと一緒で相手の表情読めない分、文章から心理状態を把握せにゃならん気がして無駄な深読みしまくって凄い疲れそう。

さてぼちぼち本業に復帰するか
たぶんこれから実装する計算方法はモーレツにメモリを喰うのでメモリの使われ方ってのを今更ながらに調べる。
メモリは大まかに分けて、コードセグメント、データセグメント、スタックセグメントの3種類のエリアに分けて使われるらしい。アドレス指定とかに使うレジスタもそれぞれに応じてCS,DS,SSが用意されてるとな？
コードセグメントはプログラムコードが機械語に翻訳された本体を収める部分。

データセグメントはデータ、.bssとヒープの3つに大別される。データ、.bssは静的な変数とグローバルな変数が納められる。データには初期化されてるのが、.bssには初期化なしのが振り分けられる。ヒープは動的に割りつけられる変数なり配列なりインスタンスなりに使われる。

スタックセグメントは手続き中の局所的な変数や配列、関数の引数の引き渡しなんかで使われる。

だいたいのところは上に書いたようになってるらしいが、具体的にどんな変数がどの領域に確保されるかというのは言語やコンパイラやその他もろもろによって変わってくる。
Fortranの場合だと
module中の大域的な変数、save属性付きの変数→.bss
allocatable属性な配列→ヒープ
普通にサイズが宣言されている自動配列→スタック
手続き中の変数（mainでも同様）→スタック
てのがデフォみたいだ。

フラットメモリモデルとかいう一般的な方法だとメモリの0x00000000番から順にコードセグメント、データセグメント、そしてメモリ0xffffffff番（32bitなOSでは1プロセスで使えるメモリの最大値らしい）から番号が小さいほうにむけてスタックセグメントが配置される。動的割り付けが行われるとヒープがアドレスの番号の大きいほうに伸びていき、関数の呼び出しなんかをしまくるとスタックが番号の小さいほうに伸びていく。で、動的割り付けでばかでかいのをやっちゃったり再帰手続きでばかでかい配列をつかったりするとヒープとスタックが衝突することがあるらしい。ふつー、スタックのサイズには制限がかかっていて、2～8Mになっている。Linuxだとlimitっていうコマンドで調べることができるらしい。スタックのサイズ制限はunlimit stacksizeで外せるのでそうするとこういう事態が起こってしまう可能性があるらしい。OSの側でそういうの監視してないんかよって気もするが、glibcはスタックがヒープ領域を破壊するのに無頓着らしい…。

そしてFortranの処理系ではヒープ領域へのアクセスがなんか怪しい。Fortranでの配列の動的割り付けはallocateで行われるんだが、たとえば　allocate(A(2,2))と2×2の配列を割り付けたあとで、A(1,3)とか「あーそこ配列の定義外やろ…」ってとこに普通にアクセスできるし代入とかもできてしまう。コンパイルもされるし実行も普通にされるし（当然結果は変になるが）、gdbもスルーすることがある。動的な配列割り付けを使うとこのやりがちなバグを見つけるのが難しくなる。

アクセス速度自体もヒープは他に比べて遅いらしいので、数値計算なら配列は大きめに見積もってmoduleで静的に宣言しとくのが正解ってことか？

元旦

2009-01-01T19:21:00.003+09:00

あけましておめでとうございます。

うちの家では仕事の関係上、元旦は5時起きだったりする。たぶん一年で最も早く起きる日が元旦…。5時って眠すぎる＆寒すぎるだろ常識的に考えて。
なんか軽い時差ボケになりそうだ。

朝、仕事片付けて、初詣に行く。おみくじを2か所で引いて二個とも吉だったのでよかった。元旦に行く近所の神社は毎年結構まとを得たコメントが書いてある気がしておみくじひく時なんだか緊張する。なんとなく神社仏閣にいくとおみくじを引きたくなる。星占いとかとちがって自分が選択している感じがなんか面白いんだろう。今年の占い関連は何見ても自分が該当する部分は「幸運だが注意が必要」みたいなことが書いてある気がする。今年は注意深く生活しよう。

正月の仕事が終わると一気に暇になって結局学会用のデータ処理とか本業関連の作業している。元旦からそれ以外することないんか自分、なんか趣味見つけんとなぁ。

さて明日はバーゲンだ。ええもん見つかるといいなぁ

大晦日

2008-12-31T20:14:00.003+09:00

あと数時間で2008年が終わる。今年はものすごく一年過ぎるのが早かった。ここ数年一年が過ぎるのが早い。たぶん出張が多い＆研究のスパンが数か月とか半年単位だからなんだろうなぁ。

今年の漢字は「変」らしい。自分的今年の漢字は「転」かなぁ。いろいろ超展開が多かった。公募の書類の経験値上がるしまぁ無理とは思うが出して見るかーと応募した賞が予想外にもらえてしまったりまぁそのほかにも年末ギリギリまで超展開がいろいろある一年で面白かった。そういや去年ぐらいからわりと日々面白く過ごしてる気がする。ありがたや。

今年も各地を転々としたなぁ。遠出では東京4回、アムステルダム、岩手、松江。合間に京都行ったり奈良行ったり。自分ではひきこもってる傾向ある気がするけどなんだかんだいって結構出歩いているかも知れん。来年はどこ行けるんだろう？とりあえず春学会は東京らしい。また東京か。

研究も結構展開が進んだ。ちょうど一年前はNRGの計算コード作り出してて一不純物問題のテストケースがなかなかうまくいかずに四苦八苦してたんだなぁと見返して思い出した。去年も今年も大晦日似たようなことしてるなぁ。今から思えばプログラミング経験ほぼ０で一からNRG自作とかよくやる気になったもんだ。無謀というか今から思えば半分自棄だったんだな。かなり苦しかったが、選択してよかったと思う。他の人に役立つイノベーションを創出したわけではないが、自力で何か知らんことを把握して使えるようになったという経験はすこしは自信の元になるだろう。

さて学年的にも、流れ的にも今年は起承転結の転だったわけで、来年は今年いろいろいじったり拡張したものをまとめる「結」の一年となりそうだ。D論も書かないといけないし。気は抜かないけどまぁ楽しく面白く一年やっていけますように。来年もよろしく。

年賀状とか

2008-12-28T22:27:00.002+09:00

今年は年越しの準備もぼちぼちすすめ、割合順当にすべてが進んでいるはずだったんだがここにきて親戚が入院とかアクシデントかい、ということで今年もなんだか順当ではない年末年始になりそうな気配だったりする。年末年始は通常業務＋恒例行事な仕事＋たいていなんかハプニングですごいドタバタすることが多いので正直「年末年始だからって特別に構えず普通の日々として過ごせばええがな…」と思う。酢の物ぎらいとしては、おせちの具は正直苦手なもののオンパレードだしな。酢蓮根とか酢牛蒡まじ勘弁。

でもなんか大掃除はしなくてはならない気がして玄関を水で流すとかしててとうとう足にしもやけが発生する。痒くて発狂もの。大掃除は日本的しきたりなのかして研究室の大掃除のときは、留学生に「みんなやたら掃除してるけど査察がはいるのか？」って聞かれた。査察っていう発想もすごいと思うけど。むしろ向こうの国には大学にも査察があるってことか？

年賀状を作る。新PCはメーカー製じゃないので、年賀状作成ソフトとかが入っていない。年一回のことに専用ソフトを買うのもなんか嫌なので、WordとExcelでなんとか凌ごうとする。
参考サイト
http://kokoro.kir.jp/know/nenga.html

Wordの機能に「差し込み文章」というのがあって、その「作成」っていうボタンから、はがき文面印刷ウィザードなるものが立ち上がる。宛名の作成のときに、Excelで作った住所録を読み込むことができる。なのでWordとExcelで年賀状も何とかなるわけだ。まぁいろいろ使いにくいところも多々あって（Word側でウィザード使ってる間Excelの元ファイルを開けないとか、なぜかウィザード使ってる間は互換バージョンなのでしょっちゅうバグるとか）専用ソフトのほうがそりゃ使いやすいわ。まぁ試行錯誤のすえなんとか刷ることはできて問題は添える一言だったりする。こういうん本気で苦手やわぁ。気の利いた短文なんてなかなか思いつかんよなぁ。

仕事納め

2008-12-26T22:34:00.002+09:00

仕事納めしてきた。まぁ研究そのものは計算機も動いてるし年末年始かんけーなくたぶん家でぼちぼちやってるんだろうが。事務が今日で年内終りなので大学行く必要は明日以降来年5日までだいぶ薄れる。企業は29日締めなのかなぁ。大学エリアは小雪舞うほど寒くて電車から出た瞬間「やべ…さむっ」って感じで家帰りたくなる。寒いのとなんか眠いのやその他もろもろで締めなのに大半ぼけーっとしてなんもせんと終わってしまった。

いい加減26日なので年末を実感してもいい気がするんだが、あんまり年の瀬な気分しない。研究室の人間も結構ぎりぎりまで実家帰らない人多い。大学はむしろ3月が年末な感じだからなぁ。そして本業が一部現在進行中なせいもある。年内の目標は完了したけど、じゃー次の段階に進めるかって結局進行中。あと5日でキリいいところまですすむんかー？？掃除もせにゃならんしなんか忙しいのぅ。

そうだ休みだった

2008-12-23T21:39:00.002+09:00

今日が祭日で休みであることに昨日の夜に気がついた。やばい、普通に大学行くところだった。月曜に普通の生活して、また休みってなんか勘が狂うのでやだな。まぁ実際はあまり休みじゃなくて、掃除と年末恒例家内手工業的家業で暮れていったので普通に大学に行ったほうが楽だったかも知れん。大学の机は昨日大掃除で片付けたのですごいきれいで快適な状態だし。

自宅の大掃除を一日でするのは死亡確定だなぁと思ったので、分割統治方式でいくことにした。プロジェクト管理ツールを無駄に使ってみるテスト。だって本当に掃除が苦手なんだもん。管理しないと、手付けてない部分と無駄にすげー頑張った部分とができて全体としてあまりきれいじゃない状態で終わりそうだ。とりあえず今日はベランダとかサッシとかの水洗いと机の周辺＝腐海をなんとかしようとする。机の上は一応表面が出るようにはしてるんだが、その周辺のファイル入れとかに紙類をろくに考えずにほりこんでるのでわけがわからない。いちおー研究してるもんで紙類は異様にたまる。文科系だと溜まるのは本らしいが理科系の場合オンラインでPDFのファイル落としてくるのが大半なのでめんどくさがってファイリングとかしないと使えない紙類がたまりにたまる。PDFなんだから画面でみりゃいいようなもんだが書き込みとか目の疲れを考えるとやっぱ印刷してから見たくなるな。なんとなくプログラムも印刷してからチェックしたくなる。数年前にプリントアウトして未読のいつか読むかもな論文コレクションとか結構思い切っていろいろ捨てたワリにはすっきりしないなぁ。寒い中水仕事したのと腐海の毒で風邪ひいたわ…。
あしたはクリスマスイブらしいが、B4がゼミするらしいので大学に行かねば。まぁ異教徒の祭りだからいいんだけどさ…。今年の4年生まじめやなぁー。自分が4年のころってクリスマスはもちろん、ものすごい長い冬休み取ってた気がする。あのころ勉強してりゃいまごろ楽できたやろな…。工学部自体は空気読めるのか22日で授業は終わりらしい。大学すげー人口密度低そうだぁ。しゃーない粛々と研究すっか。

SYNC

2008-12-19T21:39:00.002+09:00

朝の4時に胃がモーレツに痛くて目が覚める。揚げ物の威力恐るべし。それとも風邪なのか。なんでよりによって早朝に腹痛が来るかね。そういえばこないだもこんなことがあった。

まぁそれで思い出したのが体内時計の話だったりする。内臓機能や睡眠周期は体温の上下と同期しているそうな。で、通常の生活をしていると体温の最下点が朝3時とか5時にくるらしくこのときに一番人間の認知能力が下がるらしい。（私の場合、腹痛が朝来るのも体温下がるからか？）

てな話が最近読んだ「SYNC」つー本に出てた。SYNCは別に睡眠周期についての本というわけでなく、いろんなところに出てくる同期現象についてのレビューみたいな内容になっている。こないだ読んだ蔵本先生の本で非線形科学なんかおもしろそーやんとおもってその流れで読んでみた。ホタルの点滅も睡眠周期も同じような数理モデルで扱えるそうな。物性系でもモデル扱うけど、数理工学系のモデルはいろんな要素をばさばさ切りまくっててなんか感触が違うなぁと思う。でもこの同期現象についてウィーナーって人が模式的に提案したピークが妙に近藤ピークの傾向に似てて気になる。超伝導は同期現象らしいが近藤効果も同期現象の一種とみなせるんか？

ホタルの話等もおもしろかったんだが、なんか一番興味引かれたのは「スモールワールドネットワーク」。ランダムなコネクションが入ることで、近隣に影響されるクラスターな性質を持ちつつすごい離れた部分と相互作用するようになり、系が一気に同期しやすくなるそうな。実社会の人間関係とかとよく対応してそう。近場とも遠くとも相互作用することで、「スケールフリー」になって、ってな話は相転移ともなんか共通するものを感じる。スケールフリーとかべき乗則はなんかいろんなとこに出てくるキーワードだな。そういうものをいかに扱うかっていう方法論がいろんなところで必要になるんだろう。覚えとこ。

ボトルネック

2008-12-17T22:39:00.003+09:00

使う対称性を変えて計算するバージョンも出来上がった。保存量として全スピンを採用するほうが基底数が減るし保持する基底が少なくても精度がでて断然よいだろうから計算効率としてQS対称性＞＞QSz対称性だと勝手に思っていたが、QSzのバージョンは比較にならんほど速度がでてびっくりする。5時間かかってた計算が20分弱に短縮される。なんかすごいミスしてるんじゃないかと逆に焦るな…してるんかもしれんが…。そら短時間ですめば済むほどいいんだが、何が結局前のバージョンのボトルネックだったのか探らんと気持ち悪い。

fortranで実行時間の計測は
system_clock(integer型)
という組み込み手続きを使える。たとえばcall system_clock(t1)のようにすると、t1に呼び出された時点での「処理系の時間＝1970 年 1 月 1 日の 00:00 からの経過秒数？」なる壮大な数字が返ってくる。ほんとは引数は3つあってsystem_clock(count,count_rate,count_max)らしいけど。二個目は一秒に何回時計が刻むかとか3つ目は最大値がいくらかとかなのでまぁ省略してもよかろうと。何回か呼び出してその差をとればどのサブルーチンがどんだけ時間くってるかわかるわけで、さてこれでボトルネックをちまちま探そう…

今日は忘年会で、串揚げ自分で揚げ放題で食べ放題のところに行った。揚げ物普段あんまし食べないのでまずまずおいしかったけどこたえまくり。あー忘年会終わったとはいえまだ忘年できないなぁ。とりあえず予算関連を年内に片付けないと…

フォック空間（チラシの裏レベル）

2008-12-15T21:33:00.003+09:00

参考論文の式は導出さっさとできたぜちょろいぜとナメてかかってたらそんなに甘い話じゃなかった。第二量子化された形式の式変形は数学とゆーより記号の規則に沿ってがちゃがちゃ並べ替えるだけなので慣れてしまえば追うのはそんなに難しくない。やばいのはその意味の把握。そもそも完備なフォック空間云々とかいう以前にFock空間がなにかはっきり分かってないのがよろしくない。
量子論というと、あーシュレディンガー方程式でしょ？ってイメージだが、実際使うときにはあの微分方程式の形式をそのままごりごり解くわけでは無かったりする。だって大抵の微分方程式ってそのままではうまく解けないし。演習でやるようにポテンシャル障壁が無限で閉じ込められてますよみたく極端な境界条件とかじゃないと解析的な解は出ない。ひとつの方法として微分方程式解く代わりに、例の波動関数をベクトルとみなして、固有値問題とか変分の問題に読み替えたりする。実際使う段になると量子論＝ほとんど線形代数。（微分方程式の形状で攻めていくにしろ数値的に扱うときにはほとんど差分方程式とかに落として線形代数の問題にしてしまうし）このベクトルが張る空間がヒルベルト空間に相当するらしい。

一粒子の系ならベクトルとみなすだけでええかんじになるんだけど、粒子が何個もあると、波動関数の形を決めるのに反対称性だのなんだのの規約が入ってきて波動関数ひとつ書くのもとにかくめんどくさい。粒子の入れ替えに対して波動関数の符号がどうなるか＋どの粒子が何個あるかだけで状態を書けばええやんってのが第二量子化の発想（だと思われる）。なので第二量子化された場合の基底ベクトルは各粒子の個数の情報で決まるベクトル（アップスピンのが一個とか二個とか）の直積（テンソル積？）の形状になっている。|n1>|n2>…|nN>みたいなイメージかね。そんなベクトルで張られる空間をフォック空間というらしい。なーんか裏にいろいろ数学的構造があって対称性・反対称性の性質をうまく使って縮約してるんだよとか聞くけどよくわからん。

とにかくフォック空間は生成消滅演算子を作用させたベクトルの直積で張られるっぽいので今使ってる基底（全スピンの固有状態になるようにいろいろひねっててシンプルな直積になってない）でフォック空間を定義していくのはまずい気がしてきて、結構最初からコード組み直しかよって感じだ。基底が直積で書けるようにQ,Szの対称性でするもんなんかなぁ。教えて偉い人。まぁそれだと磁場の影響とか異方性の影響も見れるからいいんだけど、一からまたかんがえるのめんどくせー。

TA関連２

2008-12-14T21:05:00.004+09:00

大学が停電してなんも作業できなくなってチョー暇暇になる。まず研究がらみのこと以外基本することないんかい…てのがとても問題なんだけど。暇なんで掃除したりうだうだしたりついでにTAで面倒見ている課題の完成形を作ってみる。
前で一応「点滅」まではなんとかなったわけで、クリスマスツリーらしきものにするためには、背景画像を読み込んで、その上に描画する＋描画する色を選択可能、点滅は色の系統ごとにばらばらに点滅する＋点滅の開始、停止、クリアができるようにするという仕様にする。gif読み込んで表示するほうが速そうなので、電球とかも全部画像で保存しておく。

javaにおいてはインターフェースってのが多用されてるみたいだ。イベント処理とかにはそれに応じたインターフェースをimplement（なにそれ）してつかうっぽい。インターフェースってそもそもなんやねんって感じだが、こういう部品をつかいますよってのだけが定義されてて実際の動作は使う側で記述しなあかんものらしい。んーメソッドを直に書いたらだめなんかいな。なにが嬉しいんかちと分かりにくいぞ。（f90にもinterface文はあるけど形状引継ぎ配列使うときとかに要る、サブルーチンとかの引用時の仕様定義ってだけでどっちかというとめんどくさーって感じかなぁ。）

背景画像をアプレットの表示サイズが変えられてもリサイズして表示するようにするためには、インターフェースComponentListenerをimplementして必要なメソッドを記述＆特にcomponentResizedにリサイズ処理を書けばよいらしい。
ボタン関連のイベント処理はActionListenerていうインターフェース＋ボタンにActionListenerを追加＋actionPerformedにボタンごとに処理をちまちま書いていけばいいぽい。
画像の読み込み待ちにMediaTrackerなるものをつかうといいらしい。
ボタンの配置をいじりたいときはPanel上にボタンを配置するようにする模様。
↓とりあえずできた物。たぶん効率悪い処理しまくり


import java.awt.*;

import java.applet.Applet;

import java.awt.Graphics;

import java.awt.event.*;

import java.net.URL;


public class test2 extends Applet
implements Runnable,MouseListener,ActionListener,ComponentListener{
     private static final long serialVersionUID = 1L;

    String status,msg;
    int change=0;
    int posX;
    int posY;
    int radius=20;
    public int col=-1;
    int radiuspink,radiusgreen,radiusblue;
    int BUTTON_NUM=6;
    Image off=null,img;
    Image pink,blue,green;
    int w,h;
    boolean flag=false;
    Graphics g;
    public final int MAX=100;    //かける最大数
    private Button[] buttons ;
    public int nump=0, numb=0,numg=0; //個数
    public Point[] pinkoval=new Point[MAX];// 座標
    public Point[] blueoval=new Point[MAX];
    public Point[] greenoval=new Point[MAX];
    Thread th;
    boolean start=false;//点滅の制御
    Color c0,c1,c2; //色の設定用
   
    //イベント処理用のメソッド
    //マウス関連（MouseListener）
    public void mouseEntered(MouseEvent e) {}
    public void mouseExited(MouseEvent e) {}

    public void mouseClicked(MouseEvent e) {
      if ( col != -1){
       if (col==0){//pink
         if(nump<MAX){        // for Max number
             pinkoval[nump]=new Point(e.getX(),e.getY());    // get coodinate¸
             nump++;            //
             repaint();        //
           }
       }
       if (col==1){//blue
             if(numb<MAX){        // for Max number
                 blueoval[numb]=new Point(e.getX(),e.getY());    // get coodinate¸
                 numb++;            //
                 repaint();        //
               }
         }
       if (col==2){//green
             if(numg<MAX){        // for Max number
                 greenoval[numg]=new Point(e.getX(),e.getY());    // get coodinate¸
                 numg++;            //
                 repaint();        //
               }
           }
      }
    }

    public void mousePressed(MouseEvent e) {}
    public void mouseReleased(MouseEvent e) {}
   
   
    //ボタンのイベント処理(ActionListener)
      public void actionPerformed (ActionEvent e){
             if (e.getSource()==buttons[0]){
                 col=0;
             }
             if (e.getSource()==buttons[1]){
                 col=1;
             }
             if (e.getSource()==buttons[2]){
                 col=2;
             }
             if (e.getSource()==buttons[3]){
                 start=true;
             }
             if (e.getSource()==buttons[4]){
                 start=false;
             }//startボタン
             if (e.getSource()==buttons[5]){
                 nump=0;numb=0;numg=0;
             }//clearボタン
         }
     
      //リサイズ処理(ComponentListener)
       public void componentResized(ComponentEvent e){
            Dimension dim = getSize();
            if( img != null)
                /* サイズを変更した画像を作成する */
                img = img.getScaledInstance(dim.width, dim.height, Image.SCALE_SMOOTH);
        }

        public void componentMoved(ComponentEvent e){}
        public void componentShown(ComponentEvent e){}
        public void componentHidden(ComponentEvent e){}
     

    public void init(){
        buttons = new Button[BUTTON_NUM] ;
        setBackground(Color.white);
        setForeground(Color.gray);
        setFont(new Font("Serif", Font.BOLD, 20));
        setLayout(new BorderLayout());
        posX=0;
        posY=0;

        status="out";
        Panel Panel1 = new Panel();
        Panel1.setLayout(new GridLayout(2,3));
        Panel1.setBackground(Color.WHITE);
        this.addMouseListener(this);
        this.addComponentListener(this);
        //画像のリサイズ
        this.componentResized(null);
        //ボタン設定
        buttons[0] = new Button("pink");
        buttons[1] = new Button("blue");
        buttons[2] = new Button("green");
        buttons[3] = new Button("start");
        buttons[4] = new Button("stop");
        buttons[5] = new Button("clear");
        for ( int i = 0 ; i < BUTTON_NUM ; ++i ){
            //ボタンに処理を追加
            buttons[i].addActionListener(this);
           //パネルにボタンを追加
            Panel1.add(buttons[i]);

          }
        //色の設定
        c1=new Color(180,210,255);
        c2=new Color(180,255,180);
        buttons[0].setBackground(Color.PINK);
        buttons[1].setBackground(c1);
        buttons[2].setBackground(c2);
        //下に配置
        add ("South",Panel1);

    }


    public void start(){
        //画像読み込み

        pink=getImage(getCodeBase(),"pink.png");
        green=getImage(getCodeBase(),"green.png");
        blue=getImage(getCodeBase(),"blue.png");

        String param0 = "tree.png";
        if(param0 != null)
        {
             img = getImage(getCodeBase(), param0);
              URL urlCode = getCodeBase();
                System.out.println("Code : " + urlCode.getPath());

            if(img == null)
            {
                String msg = "fail to read file";
                System.out.print("fail to read file");

                repaint();
                return;
            }
        }
        th = new Thread(this);
        th.start();
    }

    private void loadBackground(Image back)
    {
        msg = "loading";
        repaint();
        // 画像が読み込み終わるのを待つ (Media Tracker)
        MediaTracker mt = new MediaTracker(this);
        mt.addImage(img, 0);
        mt.addImage(pink, 1);
        mt.addImage(green, 2);
        mt.addImage(blue, 3);
        while((img.getWidth(this)) == -1 || pink.getWidth(this)==-1
                || green.getWidth(this)==-1 || blue.getWidth(this)==-1)
        {
            try
            {
                mt.waitForID(0);
                mt.waitForID(1);
                mt.waitForID(2);
                mt.waitForID(3);
                if(mt.isErrorID(0) && mt.isErrorID(1)&& mt.isErrorID(2) && mt.isErrorID(3))
                    throw new NullPointerException();
            }catch(Exception e)
            {
                msg="finish";
                repaint();
                return;
            }
        }

                System.out.println("img.getWidth="+img.getWidth(this));
                Graphics gg = back.getGraphics();
                gg.drawImage(img, 0,0,this);
                gg.dispose();
                radiuspink=pink.getWidth(this);
                radiusblue=blue.getWidth(this);
                radiusgreen=green.getWidth(this);
                // check for loading
                System.out.println("pink"+radiuspink+"blue"+radiusblue);


        }

    public void run(){

        Image back = null;
        Dimension dim = getSize();
        w = dim.width;
        h = dim.height;
        if(w <= 0 || h <= 0)
        {
            w = 500;
            h = 500;
        }
        off = createImage(w, h);
        back = createImage(w, h);
        g = off.getGraphics();
        loadBackground(back);

        if (start=false){
            repaint();
        }
    else{
      for(;;){
          repaint();
            if (change< 300){
                change=change+1;
            }
            else {change=0;}

            try{
                   Thread.sleep(500);
                }
                catch(InterruptedException e){
                }

        }
        }
    }

    public void update(Graphics g)
    {
        paint(g);
    }

    public void paint (Graphics g){

        g.drawImage(img, 0, 0, this);

       if (start==false){
           for (int i=0;i<nump ;i++){
                 g.drawImage(pink,pinkoval[i].x-radiuspink/2,pinkoval[i].y-radiuspink/2,this);
                }
           for (int i=0;i<numb ;i++){
                   g.drawImage(blue,blueoval[i].x-radiusblue/2,blueoval[i].y-radiusblue/2,this);
                   }
           for (int i=0;i<numg ;i++){
                   g.drawImage(green,greenoval[i].x-radiusgreen/2,greenoval[i].y-radiusgreen/2,this);
                   }

       }

       else{
       if ((change%3)==0 ||(change%3)==1){
        for (int i=0;i<nump ;i++){
         g.drawImage(pink,pinkoval[i].x-radiuspink/2,pinkoval[i].y-radiuspink/2,this);
        }}
       if ((change%3)==1 || (change%3)==2){
        for (int i=0;i<numb ;i++){
           g.drawImage(blue,blueoval[i].x-radiusblue/2,blueoval[i].y-radiusblue/2,this);
           }}
       if ((change%3)==2||(change%3)==0 ){
        for (int i=0;i<numg ;i++){
           g.drawImage(green,greenoval[i].x-radiusgreen/2,greenoval[i].y-radiusgreen/2,this);
           }}
       }

     }

    public void stop(){}

    public void destroy(){}

   }

さて学生さんになんて説明したらええんやら。300行ぐらいのサンプルで「え、長い（読むのやだなぁ）」って言われるしなぁ。輪講とかでも思うけどダイレクトに答え言わずほのめかして、考えさせて、正解に持っていかせるのって超むずかしい…

新規作成

2008-12-11T22:48:00.002+09:00

新しい数値計算用のコードを作り始めた。Complete Fock Spaceを用いる方法で、これを把握すれば磁場下もOK,時間発展もOKっていうなかなか強力な手法らしい。時間発展とかいわれるとなんか手を出してみてたくなるよなぁ。ダイナミクスだもん。

NRGは「くりこみ群」なので本来、あるエネルギースケールでの固有値を順に求めていく＋各段階で次のスケールにかかわってこないであろう高エネルギーの状態を捨てていくことで計算機のメモリがかなり少なくても計算できるスキームになっている、スバラシイ、個人的にこの辺も気に入っている。通常の不純物問題ならば、このエネルギースケールでの階層分け＆どう考えても低エネルギー状態に寄与してこない状態をばしばし切るのは賢いやりかたなんだけど、磁場がある場合みたいにそんな階層をぶちぎったべつのエネルギースケールがある場合には切っている部分の効果が無視できるのかという問題が出てくる。

そこでreduced density matrixですよ。このreduced density matrixをつかえば、系が環境とシステム（注目する部分系）からなっているときに、システムの部分だけにかかわってくる物理量の計算に使う密度行列に環境の影響を含ませることができ、部分系を扱うことによって生じる「端の効果」みたいなもんを抑えることができる模様。DMRGでもつかわれているっぽいな。計算の流れとしては、トータル2回NRGの計算を行って、2回目のときに捨てている部分の寄与を合成してしまうようなイメージになっているようだ。なんだかアルゴリズム考えるの難しそーだが面白そう。最初の計算のときにでてくるものすごい大量のデータを保持する必要がある。まぁ最近のメモリ容量なら余裕なんだろうが。

で、データを格納していく部分をつくってたら、久しぶりに意味フなエラーをくらった。
*** glibc detected *** double free or corruption (out): 0x000000002f5a0270 ***
forrtl: error (76): IOT trap signal
Aborted

データ入出力関係でエラーが出て、強制終了したらしい。不気味なことに、結果はちゃんとそれらしいことが返ってきているという。

*** glibc detected *** double free or corruption　はRedHat系ではデータの破損を検出したときに出るっぽい。
forrtl: error (76): IOT trap signal
入出力関係で異常終了。(通常coreダンプ出るらしいが、今回なんもでなかったんだけど。）

ちなみに原因は4列で宣言した配列に5列あるデータを代入しようとしていたから。IOT trap signalくらったら配列を疑ってかかるのがいいのかもしれない。

Texメモ

2008-12-08T22:03:00.005+09:00

学会で岩手行った時に帰りに中尊寺によったら、そこの神社から新年のご祈祷しませんかっていうDMが来た。最近は宗教関係も営業をしっかりしてるんだなぁと妙に感心する。顧客というか定期的にきはる人々にはご案内だしたりするけど、一見さんのデータまで把握して発送するとは気合い入っているなぁ。たしかにわざわざ封筒まで付けて送ってきてもらったのをスルーするのもなんだしやってみるかねぇとか思ってしまうな。普通の会社よりDMの効果は大だと思う。

論文の苦行続く…こんどは短めのletterの形式で書こうということになったんだが、letterの場合、大抵スタイルシートが配布されてないのでReferenceの形式とかを自分でいじらないといけない。アブストでもなんでもTexのスタイルシート落としてきて、ってのに慣れきっているのでTexのあれこれを忘れまくっている。ので調べたことをメモ

日本語対応のplatexとかだと、「参考文献」とか「図」とかが出てきてしまって困る
→プリアンプル部分に
\newcommand\refname{References}
\newcommand\figurename{Fig.}
\newcommand\tablename{Table}
\newcommand\appendixname{Appendix}

デフォルトの参考文献の行間でかすぎ対策
JPSJのスタイルファイルを凝視した結果プリアンプル部分にこういう感じのマクロ書いとくといいぽいのがわかった（これは　1) A.AAA et.al. みたいに番号）で出力されるタイプになってる模様）

\def\thebibliography#1{%
%改ページする仕様
\clearpage
%参考文献はボールドで出力される
\par\noindent\textbf{\refname}\par\vskip -4pt %\vskipで行間を設定
\small
\list{\arabic{enumi})}{\settowidth\labelwidth{[#1]}
%↑で出力の形式を番号）にしているぽい
\leftmargin\labelwidth
\advance\leftmargin\labelsep
\usecounter{enumi}}
\sloppy
\sfcode`\.=500\relax}
\let\endthebibliography=\endlist

Tex奥が深すぎ。やばい

筋違え

2008-12-07T21:46:00.002+09:00

どんな寝方したんだか、起きたら肩甲骨近傍の筋肉が攣ってかなり痛い。あんまり動けないのでひたすら寝る。親にそれは４０肩っていわれた。いやまだありえんから。

学部時代に所属してた部活の発表会に行った。受付とかめずらしく働いてみる。研究室じゃかなり年長の部類なんで基本働かんもので。営業スマイル＆お辞儀とか普段せん活動してみるのもいいかと。めっちゃ久し振りに会った同期に「OBがなにやってるん？？」て驚かれる。なれんことやったから肩がつったのかね。手品見に行ってメーカー研究職の後輩となぜか専門分野トーク。やはりミクロスコピックとマクロスコピックの中間をつなぐ理論がねーのが問題だなぁという結論。あぁーやっぱりそこだよなぁ。

発表会はリハーサル等を全く見ていなかったので純粋にお客さんとして楽しんだ。見て「おぉ、なんかおもしれえ」と思うのはどこがツボなんだろうかとかをちょっと考える。研究発表とかでもそうなんだろうけど、まずわかりやすさじゃないかと。それに加えてとても独創的なことをするか、もしくはある程度定型っぽいんだけど美しく仕上げているものに持っていくのがウケる王道なんだろうかなぁと。いいなと思ったのが3つほどあって、ひとつは前者、2つは後者だと思う。そして独創性の追求は両刃だと思う。独創性を生み出そうとしてわけわからんくなってしまうと見ている間中ひたすら「これはなにをしたいんだ？？？」になる。って分析してんと普通に楽しめよって感じですな。職業病ですわ。

さて今日は一日寝てマンガよんで終わってしまったので明日からまた働かんとなぁ。イムリおもしれ～。なんとなく砂の惑星ちっくなSFっぽさ漂ってていい。SFっぽいの結構好きなんだが、たまに厨ぽいのひいちゃってあちゃーになる。けどこれは面白そうな予感。続き一気買いしたいけど新幹線代のせいで限度額越えでカードがしばらくAmazonで使えない罠。CDも買えないし不便でしょうがない。

TA関連

2008-12-05T21:56:00.005+09:00

プログラミングの授業をTAでしていて、今週から最終のグループ課題が始まった。ワタシらは言語がJavaに切り替わった最初の代だったので、方針とか固まってなかったせいか、個人でする課題だけで、アプレットに適当に図形出して終わりとかだったんだが、今のグループ課題ってのは班ごとに描画ソフトウェアつくってこいとか結構いかつい課題が出るみたいだ。最近の子はたいへんやねぇ～とか思ってたら、各班に一人TAがついて面倒見るってことなので、他人事じゃなくなった。

課題は選択なので、なんとかして楽勝なのを選んでもらおうと誘導を試みる。非線形微分方程式の数値計算のが一番楽だしGUIも大していらんし超お勧めだったんだが、学部生は数値計算を避けたがるぽい。全然凝ってない差分方程式とかアルゴリズムも課題文に書いてあるようなもんだし一番楽やのに。カオスとかにも絡んできて結構面白いと思うんやけどなぁ。結局クリスマスツリー＋イルミネーションを表示するていう課題になったらしい。

Fortranの人なのでJavaで描画とか全然わかんねぇよと思いながら、片手間にちょろっと調べてみる。聞かれてわかんねーとかDの面目がね・・・（専門外だけど）。まぁそれにDMRGやらサイトパーコレーションのJavaアプレット作成例とかもあるし、物理屋も可視化の方法を知っておくのは悪くなかろうと。

まず必須の機能である「・マウスで選択した部分に・電球を示すブツを配置して・点滅させる」をどう実現するかを考える。Swingはなんかしんどそーなのでとりあえずアプレットを作る。
いちおうそれっぽいのが実現したのでメモ

import java.awt.*;
import java.applet.Applet;
import java.awt.Graphics;
import java.awt.event.*;

public class Myapplet extends Applet implements Runnable{
    String status;
    int radius=20;
    boolean flag=true;
    private final int MAX=100;        // 描ける丸の最大数
    private int num=0;            // 描かれている丸の数
    private Point[] oval=new Point[MAX];    // 座標を保持するための配列
    public static final Color[] pal = new Color[256];//色の情報を保持
    Thread th;

    static
    {
        for(int c=0;c<256;c++)
        {
            pal[c]=new Color(c,c,0);
        }
    }
    public void init(){
        //初期化

        setBackground(Color.black);
        setForeground(Color.black);
        setFont(new Font("Serif", Font.BOLD, 20));
        //リスナークラスのインスタンス化
        EventHandler eh = new EventHandler(this);
        //フィールドの初期化
        status="out";
        //リスナーの追加
        this.addMouseListener(eh);

    }

    public void start(){
        //開始
        th = new Thread(this);
        th.start();
    }

    public void run(){

      for(;;){
          repaint();
            flag=! flag;
            try{
                   Thread.sleep(300);//スレッドの停止（点滅させる）
                }
                catch(InterruptedException e){
                }

        }
    }

    public void paint (Graphics g){
        if (flag==true){
             for (int i=0;i<num ;i++){
                 g.clearRect(oval[i].x-radius/2,oval[i].y-radius/2,radius,radius);
                 //描画した円がある部分をクリアする
                 }
        }

        else{
     for (int i=0;i<num ;i++){
        for(int j=radius;j>0;j-=2)
            {
                int c=255-255*j/25;
                g.setColor(pal[c]);//重ねて描画すると電球ぽく
                g.fillOval(oval[i].x-j/2,oval[i].y-j/2,j,j);
            }
        }
        }

    }

    public void stop(){}
    public void destroy(){}

    class EventHandler implements MouseListener {
        Myapplet apl;
        EventHandler(Myapplet obj){
            apl = obj;
          }
        public void mouseEntered(MouseEvent e) {}
        public void mouseExited(MouseEvent e) {}

        public void mouseClicked(MouseEvent e) {
             if(num<MAX){        // 最大数を超えていなければ、
                 oval[num]=new Point(e.getX(),e.getY());    // 座標を保存
                 num++;            // 円の数を１つ増やす
                 repaint();        // 再描画する
               }

        }
        public void mousePressed(MouseEvent e) {}
        public void mouseReleased(MouseEvent e) {}
    }
}

まぁ最低限、点滅する電球的なものを作れることはわかった。あとは背景に画像を入れたり、描画した部分のクリアのしかたをもうちょっと何とかしたり、ボタン作って色選択したりとかかねー。それにしてもJavaって長い名前多いな。ややこしい。イベントハンドラをインスタンス化とか言われても宇宙語？って感じっす。あー学生になんて説明したらいいんだ？

バイナリ

2008-12-03T22:47:00.003+09:00

もしかしたら近々Fortran90のバイナリの入出力がいるかもしれないのでメモ

program main
integer :: i
real(8) :: a,b

open(12,file="test_unformat.dat",access="direct",&
form="unformatted",recl=64)
do i=1,10
a=dble(i)
write(12,rec=i) a
end do
close (12)

open(13,file="test_unformat.dat",access="direct",&
form="unformatted",recl=64)
do i=1,10
read(13,rec=i) b
print*,b
end do
close(13)

end

バイナリを扱うときには
openでaccess="direct"が便利（直接探査、recl=で決めたデータ長のデータを扱い、データ番号で指定して自由に読み書きできる）ていうかdirectにした時点でバイナリで書かれてしまうぽいけど。
recl指定しとかないとコンパイルエラー来た。
form="unformatted"もしくはform="binary"
データに互換性を持たせたいときはform="binary"がよいらしい
write,readでどのデータ番号に入出力するかはrec=で決定

手帳

2008-12-01T21:59:00.003+09:00

手帳買うか迷ってる。
なんちゅーか大半を家と研究室で過ごし、その大半PCの前にいるために、スケジュール管理はgoogleカレンダーに記入で間に合ってしまっている。手帳かったらそこにまたその予定を書き写すとか二度手間やんけとか思うわけで。Moleskinかっけーとか思ったけど、サイズの割に高い＋ペンさすところない＋あんまり結局使わんよなぁと思って買ってない。電車の中で思いついたことなり計算なりをすべくノートの類は鞄に入ってるけど、手帳ってそういう用途には使いにくいしなぁ。そんなわけでgoogleカレンダーのプリントアウトをいつも持ち歩いてるノートに貼っとくでいいかってのがファイナルアンサーになりそうだ。あーでもMoleskinええよなぁ。

Wordの数式エディタは使いにくさ最強だったんだがOffice2007になってめちゃめちゃ改造されたみたいだ。Tex記法が効くようになった。とういうわけで使ってみる。\sumとか打ってスペースキー打つと∑でてくるやん、すげー。\cdotもちゃんとドットに変換されるし、_で下付きに認識される。\sum_(i)でちゃんとiで総和をとるの形に変換されるとは、MSは腕をあげたな。基本的にマウスで選んで押して・・・って作業が嫌いであぁもうTexでいいよとかTex派だったんだが、これならWordでTexより早く論文かけるかも知れん。あとは参考文献引用関連をTeｘ並みにしてくれればWordに文句がなくなってしまう。
PowerPointの数式エディタはまだ旧バージョンみたいだが、旧でも一応キーボードショートカットはいくつかあるんだな。知らんかった。うーん、もうOffice主流でいいかも。さてD論どっちで書くかね。

目新しいこと

2008-11-29T20:28:00.003+09:00

結果待ちってそわそわするな。で、まあなんちゅうかやる気が減ってる。体調もわるいし。ってぐだぐだしててもだめなんだけどな。

研究については、ざっくりしたデータは得られたし、結果それっぽいし、論文も投稿したし、物理的背景もある程度把握したしで、あとは微調整とかコードをグレードアップしたりデータの補間とかできれいなグラフをプロットする方法をしらべるとか正直あんまり好みじゃない作業が待っている。微調整とか仕上げとかおおざっぱな人間にはむかないこと甚だしい。なんか試作品がうまく走ったら終わった気分になってやる気なくなんのよね。うーんそろそろ新しいことに手を出したいぞ。

しばらく廃人生活しようかと思ったんだがそれってあんまり身にも心にもよろしくないよなぁと思ったので、本業とはちょっと（だいぶ？）ずれてるけど興味あることをしばらくちょこちょこいじることにした。Google教じゃないけど、なんか本業以外にも手を出したほうが脳にはいい気がする。ひとつはBlender。新PCならまぁ動くだろうと。
Blenderはフリーの３DCG制作ツールらしい。なんか流体のシミュレーションやら物理エンジンも入っているとかいうからかなーり興味が。最新版のはPython2.5いれとかないと動かないっぽい。微妙に2.6じゃなくて2.5なのか。2.6だったら「コンパイルできないよ」みたいなエラーでた。
とりあえず、起動。ボタンが大量にあってとにかく難しそう。ありがたいことに日本語のチュートリアルがあるので、それをみながらぼちぼち立方体作ったり、それを分割したりしてみる。まずはチュートリアルの例にある人形みたいな３Dモデルつくって動かすのをやってみようとしてとりあえず人形つくるまではなんとか。こういう３Dで描画するツール使ったことないけど、なんかこれ機能豊富ですごいなぁとひたすら感心する。あー物理演算まではすごく遠い気もするけど休日にちまちまいじって遊んでみよう。

もう疲れました

2008-11-27T20:07:00.002+09:00

面接おわった。結果は不明。
2日前からひどい風邪で、昨日は水のんでものどの痛みで呻くような状態だったので、次の日面接とかさすがにもうだめかとおもった。夜中にすごい寒気とかもしてたんでたぶん熱もあったかと。今日おさまったのはほんまありがたい。あちこち詣でたご利益か。
ネットでしらべると、のどの痛みにはイプブロフェンとのことなので、家にあった鎮痛剤でイプブロフェン入っているのを飲む。なんという対症療法。でも実際きくわぁ・・・。鎮痛剤おそるべし。あとひたすらヴィックスなめて、アズレンのどスプレーしまくってたら、面接のときには普通とかわらんのど状態になってた。

たぶんプレゼン自体はとどこおりなく進んだんだが、申請した領域ちょっとずれたのにしたのはまずかったかなぁという反応だったのでなんかすごく不安。「研究の有用性はよくわからんがとにかくすごい自信だ」みたいな反応とコメントが・・・orz。普通に話しただけなのによっぽど尊大な態度と思われたのか。これが吉なのか凶なんかいつ結果くるんだろう。あぁぁーウツだ。

連休(外出編)

2008-11-24T19:49:00.008+09:00

連休中特に予定なくぼーっとしてたら、大阪に偶然来ていた後輩から連絡が入って出かけようってことになった。家にこもってたら正直気分的にやばかったのでとてもありがたい。
とりあえずどこ行くよって話になって、秋だしなぁとよくわからん文脈で、行ったらなんかありそうだしと結局京都に行く。一週間のうちに2回も京都行くとは予想外な展開だ。まるで紅葉厨。
秋の連休の京都をナメてたらやばかった。本気で人が多すぎる。計画性少なめに京都駅にとりあえず行くと、伊勢丹のごはん食べるエリアは何時間待ちって感じだった。
後輩いわく、「京都人に聞いたら高台寺がお勧めらしい」が、それは火曜日に行ったので今回は八坂神社～円山公園～知恩院に行った。地下鉄四条から祇園界隈の道が東京レヴェルに人が多かった。ありえん。都路里どんだけ人ならんでんねん。あんだけ並ぶってことはおいしいんかなぁ～。
八坂神社に詣でておみくじも引く。受験「冷静にのぞめばよし」でちょっとほっとする。勝負事は「油断禁物」らしい。公募の面接ってどっちだろ？まぁなんにせよ、あとは自分が頑張るしかない。

知恩院は行ったら通常の拝観時間が終わってたのでライトアップに行く。券売所に列ができてたのでならんどくと、最初のほうにさくっと入れるので結構空いている状態で拝観できた。方丈庭園のなかに激しく赤いもみじの木があってインパクトがあった。黒い空を背景に独特に真っ赤な葉がライトアップというのはなかなかに凄まじいものがある。紅葉はたぶん日本人のDNAレベルでなにか訴えかけるものがあるんだで。知恩院は本殿も夜間拝観していて、ろうそくが何百本も灯された本殿は怖いというか迫力があった。

券売所に並んでいる間とか近況なんぞを雑談する。頭に残ったものをメモっておくと
・「理系の人々」が面白い
・最近のギャルゲーは携帯にメールを送ってくるらしい。なんかすごいな。
・加速度センサーの仕組みはどうなってるのか（あとで調べると基本的には、おもりの移動をつかった原理なんだが、MEMS技術ですごいことになってるらしい。）

理系の人々、確かにおもろいな。「すぐ検索する」、「めんどくさいからマクロ作る」、「エクセルの基本色が気に入らない」、「無駄に裏を読む」、「知っていることは興味ない」は当てはまるかも。まーそれ以外はあんましあてはまらんからそこまで理系に精神汚染されてねぇってことですね、よかった。
私が勝手に思う理系の人々は「食品買うとき、成分を細かくチェック」、「ストレージの大容量に無駄にときめく」、「購入前にまずgoogle」かな。うーんちょくちょく観察してみよう。

胃痛

2008-11-22T21:57:00.002+09:00

土曜日って基本やる気がもっとも下がる。
月曜から金曜まで（まぁたまにさぼるけど）大抵登校してそれなりに時間と用事に追われてる分、大学行かなくてもいいとなると果てしなくぐだぐだしてしまいがち。
ぐだぐだできる暇があるってのは普段はいいんだけど、なんか次にストレスフルな用事が控えてる局面ではそれについてひたすら気に病んでしまってよろしくないな。今は面接受けに行かなあかんてのが結構なストレスで胃が痛い。これは今まで2回落ちてるやつなんで今回もあかんかったらどないしよーとか考えるとマジしんどい。それ以外の悩みもうだうだ考えちゃって、かなりぐだってる。なんてへたれな。すべてにおいて現状よりひどい状況になるのがいやで勝負をかけずに準安定状態でうだうだしているのが原因な気がしてならんけどドラスティックに動くのは怖すぎる。

まぁ悩んでてもしょうがないので発表資料を仕上げる。イラストレーターの便利な使い方をちょっと勉強。テキストにもグラデーションをつけることができるんだな。これはいい。なんか仕事してるほうが気が紛れるのであとはひたすら数値計算いじったりして一日終わった。あぁなんかもういろいろだめぽ

保養

2008-11-18T21:28:00.006+09:00

普段そこまで季節とか気にしないけど、桜と紅葉は見に行きたい。
そりゃ大学までの道も紅葉してるしそれなりにはきれいで一応名所っぽいんだけど、通学路ゆえまったりと見蕩れるという風にはなりがたい。そうだ京都行こう。
大阪の民としては若干くやしいのうなんだが、紅葉とか桜は京都のほうが断然きれいだと思う。観光客も多いし、つい敬遠するが。観光都市だけあって京都の寺の手入れの行き届きっぷりとかすごい。観光の目玉になってるところはやっぱりなんだかんだ言ってきれいだ。なんかよくわからんが水辺が好きなもんで、川とか池とか多いのもいい。淀川とちがって鴨川きれいだしな。高台寺の、池に紅葉がライトアップで映りこむってのが実にきれいかった。

さて、目の保養した分、面接の準備がんばろう。

雑感

2008-11-15T21:14:00.002+09:00

大学のパソコンにもInstant Rails入れたんで、Red mine動かしてみている。リンク先に書いてあった設定そのまんまで一応動く。本来はApacheと連携させるべく書いてあるようになんかmod_rails？使うの推奨らしいけど、一応動くもんを下手にいじってぶっとんでもいやだしなぁとInstant railsのデフォルト通りMongrelで動かしてる。自分がローカルでちょろっと使うんなら安定性とかさほど気にしなくていいし。ローカルで動かすと、データベースが自分とこにあるっていう安心感があるな。
たぶんこれ作られた目的とはずれてるんだろうけど、これ物理系研究者にも便利だとおもう。ありがちな、「数値計算とりあえず大量に走らせた」→「数日後、どれが何用のプログラムだったか忘れた」みたいなのもちょこちょこ記録とるので防げるし。日時が入るってのも便利。文章作成して保存する機能もあるので、思いつきとかとりあえず書いとくとかもできる。やらなあかんタスクをチケットとして書き出し→経過を記録ってのはあとあと考えをまとめるのにも重宝。経過の記録をするのはあとで見返した時に「あー今週結構すすんだな」とか妙な満足感を覚えられるので結構好き。レポート用紙にがりがり書くのもありだけど、紙媒体は階層的に扱いにくくて、結局情報がばらばらになってうまく使えない気がするので、日々の研究の経過の記録とかはパソコン使うのが正解だと思う。数式だけはやっぱ手書きが早いので、計算は紙上にするけど。

今日は風邪っぽいのでぐだぐだして過ごすことにした。「非線形科学」っていう新書が面白そうだったので買ってぼちぼち読む。ローレンツアトラクタとか前から気になってた話が数式ナシで割とわかりやすく書いてあるっぽい。相転移とかともこの辺の話はからんでくるっぽいから興味ある。そして物理とかの話をきれいな文章で書けるこの著者はすごいと思う。内容も興味アリなんだが、著者が研究分野変える時の悩んだ様子とかもうっすら書いてあってそっちも興味深々だ。
私らはどっちかというと、ミクロ・量子力学な世界よりなんだが、たぶん量子力学ベースの理屈だけで現実世界の何かを解明するってのはキツイ気がしている。現実世界はマクロな世界でダイナミクスありだし。ダイナミックな状態を量子力学ベースのみで扱おうとするとたぶん時間発展させる必要があって、時間発展させたところで扱えるのはピコ秒単位のスケールになっちゃうのでかなりきついと思う。だからDFT＋KMCとかなんだろな。ミクロとマクロをつなぐには相互作用とかポテンシャルとか量子効果がガチで効いてくるところは量子論ベースで扱って、量子論きっついマクロに近いところは別の手法で補うというのが厳密ではないけど正解に一番近いと思う。数値くりこみ群やってて、補う手法にフラクタルやらくりこみ群やらランダムウォークらへんがつかえるんじゃねーのっていうのが私の勘。なのでいろいろ漁っている

一週間が早い

2008-11-13T21:46:00.003+09:00

今週は週明けからあっという間に木曜日になった。なんか明日が金曜日とか自分の中ではありえない。２日ぐらいどこかに時空のひずみがある気がする。たぶん結構ばたばたしてたせいだろう。
ばたばたしている一番の理由は、面接が２週間後にあるからで、面接のときに紹介する業績をふやそうってことで、それまでに２編論文をなんとかしようっていう流れになってるからだ。一ヶ月で２本なんとかしろって私には無理やろ常識的に（ｒｙ。まぁそうでもないとすごーーーーくキライな論文作成作業をするはずもないのでいい機会なんだけど。そんな状況下なのにやっぱり論文がいやで、「くりこみ固定点について勉強しよーっと」とかいって新しいこと始めたりするから余計やばくなってる。

たぶん水曜日の午前中のTAも一週間が早い原因なんだろう。自分の研究室の演習じゃなくて（４回もやってるのですごく飽きた）わざわざ情報系の演習に出向いてるので結構出先で勉強する。そのぶんおもしろい。二時間半が結構早く過ぎる。普段のFortranじゃなくてJavaなのでこっちもいろいろ調べんとわからん。（てーかなんでJavaなんだろ。進学先の大半の研究室で最も使われてなさそうなんだけど。物理系ならC/C++かFortranやろ）学生のエラーの原因を解明するのに似たようなの組んでみたりしてるうちになんとなく慣れてきた。でも戻り値がクラスとかの例はまだなんとなく気持ち悪い。オブジェクト指向はうまい仕組みなのはわかった、でもファイル入出力とかにいちいちインスタンス生成するのがなんとなく謎。でもなんかこー世の中で広く使われてる言語はやっぱりいろんな機能があって便利だなー、ええなぁと。TAしてて意外だったのは、わからんことがあるときにLinuxとはいえブラウザがあるのにググらない人がおおいこと。自分たちはもっとひどかったけど、今の子はもっとPC慣れしてるかなぁと思ったらそうでもないのか。google先生でだいたい解決すると思うんだけどなぁ。

シェルスクリプトお試し

2008-11-08T21:32:00.002+09:00

シェルスクリプト・・・なんか名前は知ってるけどややこしそうで手を出してなかった代物。
でも、あまりにめんどくさい作業をしなくてはいけなくなって、これを使うとそのへんが自動化できそうなので手を出してみた。

シェルスクリプトはWindowsでいうところのバッチファイルみたいなもんで、基本的にはLinuxのコマンドを羅列していく。それにくわえて、while do~done やif ~fiといったループや条件分岐もかけるし、$1,$2といった記述で実行時に引数を渡すこともできるらしい。なのでおんなじような作業でちょっとだけちがうっていう自力でいちいちコマンド打っていくには超めんどくさい作業を大半自動化することができる。
もっとも基本的ルールは冒頭に#!/bin/shを書いておく、拡張子は.shな模様（bashを使う場合）

今回自動化したかったのはだいたい

・プログラムのあるパラメータだけを変化させたバージョンをそのパラメーターに合わせた名前を持つフォルダに移動して実行。

・フォーマットしたデータをおさめたファイルrev***.txt全部を***.txtって名前にしてnewってフォルダ下にコピーする
の2点。たぶんこういう場面には今後もよく出くわすので解決しとくと以後楽だろう・・・

なので、sedをシェルスクリプトの引数をどう組み合わせて使うかってゆーのと、ファイル名の取得＋一部分を切り出してきて使うのをどうするかってのが問題になった。

sedの実行文のなかでシェルスクリプトの引数とどう組み合わせるかの例
sed 's/a=3.0/a='$1'/g' "aaa.txt" > "newaaa.txt"
aaa中のa=3.0の部分がa=$1に変換されたnewaaa.txtっていうファイルができる。
$1は''で囲っておかないとイケナイ。sed実行文も''で囲っておく。
（$1はシェルスクリプト実行時に./aaa.sh 5.0 みたいにしてシェルスクリプトの名前打った後につづけて入力して与える）

条件に合う全ファイルの検索＋ファイル名を変数に格納、必要な部分のみを取り出して新しいファイルの名前にしてnewってフォルダにコピーていう例
# revから始まるテキストファイルをフォルダ内で検索、filenameにパス付でファイル名を格納
find . -maxdepth 1 -type f -name 'rev*.txt' | while read filename;
do
echo $filename　#filenameは./revaaa.txtみたいにパスまでついてる
# ./rev以降を読み込ませるために、ｃオプション、6文字以降を読むように。
newname=`echo $filename | cut -c 6-`
echo $newname
cp $filename ./new/$newname
done

これくみあわせて自分は./aaa.sh 5.0とか打つだけ一発で作業が終わるようになった。
ってこれ調べてる間に実はこつこつ作業してたらとっくに終わってる気もするんだが、単純＋メンドイことはなるべくしたくないし。なんか得してんか損してんだか。

VMware server + Vine Linux 4.2

2008-11-07T22:55:00.003+09:00

前から使ってみたかったVmware server を自宅機に導入。
お手軽版のVmware Player は大学で使ってるんだけど、VMware server はisoファイルを落としてきたら自分の好きなタイプのLinuxを仮想で動かせたり、Vmware tools使えば仮想マシンの画面上へのマウスポインタの移動がシームレスにできたり（Vmware Playerではこれでイライラしてあまり使わずじまい。Tools導入済みのAppearanceおとせばいいんかもしれんけど）なんだかserverのほうが高機能な予感がしたもんで。
VMware serverのダウンロードはここ。ダウンロード→押すだけ。
あと登録（無料）でプロダクトキーをもらうだけ。

Vmware自体はこんな感じで単純。あとは仮想マシンを作らないと。
VMware server の場合、OSのisoファイルがあると、それからブートしてインストールできるみたいなのでとりあえずisoをダウンロードしてくる。日本語だし、大学の情報教育センターのがVineでTAのときに見慣れてるしな（rpmつかえるしrubyはいってるし）ってことでVine Linux4.2をダウンロード。
VineLinuxをVMware server に入れる記事が結構かいてくれてて、これとこれとこれ統合してやってみる。
ゲストOSの選択肢にVineLinuxがないので代わりに何をつかって仮想マシンの構成を決めるvmxファイルつくるかが問題なんだが、RedHatLinux選択でうまくいった。（書いてる人やVineのバージョンによっても違うっぽい・・・最悪自作するっていう方法もあるらしいが）
あとは参考記事通り。Vineは設定のとことか全部日本語でほっとするわー。
再起動を仮想OSでせず、Vmware consoleのPowerのボタンでしないとだめらしい？
マウスポインタをスムーズに移動させたり、時刻がはげしくずれるのを何とかする用のVmware toolsについてはここに書いてあるとおりにすればなんとかなった。
Linuxの練習用として活用しよう。

面接通知

2008-11-05T22:24:00.002+09:00

さて今のシーズン院生（博士）が心騒がされるのは学振結果通知だろう。
私は公募関係の通知先は基本、大学にしてる。だめだったときも「うあーあかんかったわぁ」とか周りと騒いで「やけ酒じゃぁ～」みたいなノリにもちこんで気分をごまかせるもん。家に夜帰ってだめでした通知があって一人耐え忍ぶとかウチの精神力じゃしばらく浮上できない気がする。弱い子なんで。

さて今回は「面接に来い」って通知が入ってました。わぁー、ハイサヨナラな去年までとくらべりゃすげぇ進歩だ。でもDC2で面接ってまじボーダーライン上らしい、つまり面接突破できるかもしくはまた閉店がらがら～かは半々みたい。3回目にしてやっとそこに到達か、やっぱ研究むいてな（ｒｙ。とはいえ５０％で突破できるなら、がんばればなんとかなるかもしれない。
ポスター作成、気合い入れてやらんとな。

発表ポスターは学会用のは作りなれているんだけれども、学振はなんだか勝手が違うみたいだ。ちょろっと情報収集すべくググってみると学振関連の記事あるわあるわ。申請書の書き方指南とかまであるやん。いままでなーんも見ず我流で書いてたよ。ちゃんと確認しなきゃー。
ポスターもフォントは最低でも５０P以上推奨とかで「うわ、自分のなかのポスターのイメージと全然ちゃうで・・・」って感じで、ほんま調べといてよかったって感じっス。５０Pフォントなんでふつーじゃまず使わん。そんなわけでまた忙しさ倍増。倒れる暇もない。

論文のめんどくささの要因

2008-11-04T23:08:00.003+09:00

風邪だけど論文せにゃいかん。なんか一ヶ月ぐらい論文やだやだ言ってる気がする。書き出したらプログラムばぐってるのみつけで全部やりなおしとか出鼻くじかれたもんでますますやる気ない。基本的になんか新しいことを大まかなとこまでやることは結構うきうきしながらするんだけど、細かいとこをいじる作業はきらいなんだよな。論文とか書いてるととくに「あぁーなんか結局すっきりとは全部説明できないやん」とか「あんまり大したことしてないよなぁ」とか思ってダウナー入るわ。

論文でなにがいやって参考文献の整理がめんどくさい。Bibtexとかそろそろ導入したほうがいいかな。とかおもってZoteroっていうアドオンをFirefoxに導入してみた。Zoteroは対応サイトだったら著者名とかを読み込んで、データベースに記憶してくれる＋PDF本体を添付データとして関連付けてローカルのフォルダに格納してくれる（PDFをブラウザで表示して、「このページのスナップショット」ってやるとPDFがダウソされる。どこがスナップショットやねん）。で、著者名とかはBibtexの形式でエクスポートできる。PDF落とすついでに関連情報も関連づけて記憶できるのがいいよなぁ、タグつけれるし、メモもつけれる。なんだけど対応サイトが少ない。とりあえず物性系がよく使うPRL,PRB,IOP系統、Elsevier系統、は対応してねぇ・・・Arxivは対応してるっていってもな。google scholarの設定で”Endnoteへの文献情報取り込みリンクを表示”にしとけば、そのリンクをクリックするとZoteroが勝手に読み込む機能はアリ。これは便利なんだけどPRBとかが早く対応してくれんことにはあまり意味が・・・。まぁ論文以外にもいろんなWeb上の情報管理につかえるっちゃそーなんだけどそれならgoogle notebookあるしなぁってことで本格的にこれ使うかどうか迷い中。

そういえば某m×ri× claire（なんとなく伏字にしとく）ていう雑誌にとある関係でのったらしい。物理屋でこーゆの載るとかだいぶ想定外だな。ファッション系の雑誌なんだがこの一部分だけには（コメントとか略歴に）物理とか書いてあるという。なもんで研究室ではウケてた。結構でかくのってるとのことなんで買ってみるかとおもへども、ふだん女性誌ほとんど買わんもんで売り場がわからずAmazonでぽちっと買ってみるテスト。友人いわく「清純派ぽく写ってる」らしい。まじか。恥ずかしくて自分ではよく見てない。写真には頭の中身は写らんからのぅ、よかったよかった。とりあえずこの雑誌買う時に一緒に購入したのが（1500円こえないと送料無料にならん）気になってたslipknotの新盤、どこが清純派やねんってことだ。さて論文しあげるか。

暇暇にRailsとやらを試してみる

2008-11-03T21:21:00.008+09:00

風邪である。ひきこもりである。暇である。
ってわけで普段だったらたぶん時間なくてできないことをひたすらやってみる。その一つがInstantRails使ってみよう。
なんで突然Ruby on Railsなのかというと、ひとつにはRedMineっていうサーバーにおいて使う形の管理ツールを使ってみたかったから＋サーバーにアプリケーションおいていじるの一度やってみたかったから。まーたぶん今後はサーバーになんでもおいてそれを共有ってタイプの作業や処理が増えてくるであろうし、全然畑違いの物理ギョーカイのわれわれもその辺知っておいたほうがなにかといいんじゃねーかと。Fortranだけじゃなんもできんしな。たぶんサーバーにデータベースおいてそれいじる簡単な方法は知っといて損ないだろ。
で、流行ってるとことには便利なものがあるもんで、Eclipseを日本語化してそれにInstantRailsやそれようのIDE（Aptana?）やらなんかPythonとかも全部ほりこんでくれているPleiades all in one日本語ディストリビューションってのがあった。それのUltimet版とやらを落として解凍するだけでJavaのみならずからInstantRails、RubyからPythonからC++まで全部使える環境が手に入る。うはーいい時代だなぁ。新しく始める人間にとって実にありがたい。ほっとんどプログラム関係のこと知らんくてもサクッとお試しできるんだもんなー。
とりあえずeclipse起動して、パースペクティブ？をRedRailsのそれにして軽く試す。新規プロジェクトつくるだけで、一通り仕組みに必要なファイルが書き出されて、Mongrelとかゆーサーバーが起動して、初期設定画面が出てくる。そこに書いてあるチュートリアルどおりにいじって動くことは確認。へーこんなあっさり鯖ができてDBができてそれを操作ってRoRすげーとか思う。

さて問題は、ちょっと自分が使えるような簡単なものを作ってみるかって時に大抵発生する。ネットでいろいろやり方をさがして試してみるもなんだかうまくいかん。ネットの例は大抵ROR+mySQLだったんだが調べた結果
・Ruby on Rails 1.x とRuby on Rails 2.x　では動作が全然違う
・2.x ではデータベースがデフォではSQlite3ってやつらしい→プロジェクト作る時にmySQL指定とかしないとmySQL使ってくれないらしい。
・eclipse経由だとInstant Railsを別に起動しとかないとmySQLを使えない（eclipseはDBにmySQL指定してもmySQLを立ち上げるわけじゃないっぽい？）
のがうまくいかない原因だったぽい。
とりあえず自分が「んー学会発表リストをつくってみるか」とおもってそれらしきものをつくった手順をメモっとく。参考にさせてもらったここのサイト見れば、陥りがちな日本語うまくいかねーとかも回避。
1 プロジェクトとしてlistってのを作成。データベースはmysql指定

2 Instant Rails起動

3 ジェネレーターって書いてあるところでscaffoldを選択。　パラメーターは list name:string place:string title:string from:datetime to:datetime coauthor:string　で実行　（ruby script/generate list name:string place:string title:string from:datetime to:datetime coauthor:string ってコマンドラインで打つのと一緒ぽい）
これで、name place title from to coauthor っていう要素をもったテーブルをデータベースに作る用のファイルが001_create_lists.rbとかいう風にdb\migrateのフォルダにできる。同時にそうした要素への処理とかのメソッド？をいろいろ記述したのがapp\views下にいろいろできてるはず。

4 mysqlで必要なデータベースを作る。Instant Rails のI型のボタン押して、Railsアプリケーション→Railsアプリケーションの管理。新規アプリケーションの作成を押すとパスがとおった状態でコマンドプロンプトが開くので、cdとかつかってeclipseのプロジェクト置いてあるとこへ移動。
mysql -u root -p
create database list_development;
create database list_test;
create database list_production;
とやると、上記3つのデータベースができる。これはconfig\database.ymlに書かれている、Railsが使う３つのデータベースの名前に合わせてあって、デフォルトはプロジェクト名＋developmentとかだ。
この段階でたとえば
use list_development;
show table;
とかやると、中身が空なのが分かる。

5 eclipseにもどって、テーブルの中身についてはさっきので指定してあるのでRakeタスクっていうタブ（ペインっちゅーんか？）を開いてdb:migrationを選択して実行する。（ rake db:migration）するとデータベースにテーブルが作られる。
use list_development;
show table;
desclibe lists（←テーブル名）；
でなんか変更が反映されてるのがわかる。

6 サーバーを起動。http://127.0.0.1:サーバーがふり当てられたポート番号/listsで動いていることを確認。

こんなんで動いた。RoRのすげーなところはデフォルトで、show edit destroyなんかのボタンが既にできているところで、単純な仕組みをつくるだけでいいならほんと一瞬でそれらしきものができてしまうあたり。DBのしくみとかほとんど知らなくてもなんか動くものできちゃったもんなぁ。というわけでしばらくこれいじって遊んでそうな気がする。rubyって誕生石だし。明日からは本業もせんとなぁ・・・

法則発動

2008-11-02T20:13:00.008+09:00

VASSCAA-4は良い学会だった。自分が興味ある話もあったし、会場きれいし、そしてやたらと食事がでた。そして、なにがびっくりするってやたら日本酒が提供された。国際会議はビールやワインが結構出るんだけど日本酒はレアケース。松江は地酒と魚がめちゃおいしいなー。で、毎晩晩御飯は酒盛りになってた。この舟盛り4800円とか神すぎる。（＠海鮮問屋　博多）

エクスカーションていう学会主催の観光ツアーには行かず、OBの人と石見銀山へ行ってきた。世界遺産だし、なんか興味があった。島根の観光の目玉らしく駅やパンフでおしていたが、松江からはものすんごく遠かった。山陰はひとつの県がでかいんだと思う。
石見銀山は銀山の跡そのものもすげーけど、周りの景観がよかった。山間の山村というか、のどかな？風景が続く。渓流も水がきれいで、マイナスイオン感抜群だ。通行止めになっている個所があったため、石見銀山エリアを巡回するバスがなく、坑道跡までは駐車場から片道４ｋｍほどを歩く。夜のバンケットに参加するため時間に制限があったので旧市街エリアは観光できず。ちなみにずっと上り坂。正直運動不足がたたってやばい。結構なお年寄りも平然と帰ってきているのでこれでへたばるわけにはいかない。
坑道のあと（龍源寺間歩）はすごくせまくて圧迫感がある。いま通れるところはあとで整備のためいくばくか広げたらしいけど、それでも天井低い・幅狭いで圧迫感あり。当時のままのところをみるとほんとうに人一人が通るのが限界で、こんなとこに一日中閉じこもって作業とか無理だろうと思った。想像して息苦しくなった。それに落盤防ぐ柵が木とか強度的に無理だろ・・・。というか当時のあまり性能がいいとはおもえん工具＋人力であんな奥地まで掘り進むとは。なんか銀とか金って人にそこまで限界な労働やらすとは怖いもんだなと思った。

で、バンケット出て飲んで、明日は大阪か、寝るかーとか思ってたら突然かなりの腹痛が来た。げ、刺身が原因かとかおもって数時間耐えるも治まる気配がない、これまでにないタイプだ。ひどくなったりしても嫌だし大阪帰る間中痛いってのもやばいぞと思って、思い切って夜間外来を使うことにする。フロントに「えー夜中にやってる病院ってこのへんあります？」とか聞いて、タクシー手配してもらって行く。夜中だし不気味に痛いわで貧血りそうになる。遠い。まぁ原因がいまいちわからんし、死ぬほどの症状で無いので鎮痛剤でごまかして大阪に帰ることになった。大阪に帰って検査したがいまいち原因がわからんかった。風邪＋生モノ・油物大量＋連日の酒＋発表で緊張＋毎月の出張＋α＝なにか。出雲大社の神々の思召しかもしれんのでしばらく謹慎してます。というわけで連休はひきこもりな予感。見事に11月の法則が発動だ。

松江初日

2008-10-27T21:46:00.003+09:00

松江来た。
飛行機のゆれがやばすぎて死ぬかとおもった。
チケットとるとき「んーボンバルディアはトラブル多いからいやだなぁ」とおもってそれ以外の機種をえらんだら、定員36人だかのマイクロバスかよってサイズの飛行機だった。サーブとかいう名前の海上保安庁？かでも輸送機としてもつかわれてたりする機体らしい。飴と飛行機の写真入りのポストカードもらったんはいいけどそれ以上にもう揺れとか急降下が激しすぎて涙目すぎる・・・肘掛おもわずつかんだわ。初めて飛行機で酔った、おえー。飛行機のるとき最近よく揺れてたんだけど、小型機はその比じゃねぇわ。

ふらふらしつつ昼について夜のミクサーとかいうタダでごはんが食べられるウェルカムパーティ的なものまで時間があったのでふらっと出雲大社にいく。ノリでふらっといったら予想以上に遠い＆不便だった。出雲と松江って距離あるんだな。運よくバスのれてよかったけど。出雲大社のインパクトはすごい。参道の厳かなかんじとか。神社って寺より若干怖い感じがする気がする。怒らせるとやばそうというか。山がすぐそばにあってなんだかものすごい圧力みたいなものがある。あれは一度は見ておくべきかも。しめ縄にコインをなげて刺さったら願いがかなうというのでやってみたら3回ぐらいでドスッってかんじで刺さったのでなんか嬉しい。たぶんカード投げるのと一緒で手首にスナップかけて回転かけるといいんだと思う。そんなわけで一日中乗り物に乗ってたので初日から疲れた。

だうなー

2008-10-26T20:36:00.003+09:00

なーんかうっすらと憂鬱。風邪気味だからか、いやいや、いろいろ原因を考えてみたところもうすぐ学振の結果が返ってくるシーズンだからかもしんない。学振てのは博士課程の学生を特別研究員だかにして毎月いくらかのお金＋年間なんぼかの研究費をくれる仕組み。あるとないでは大違い。応募書類書くのがとにかくすんごーくめんどくさい。で、結果を見ては毎年落ち込まされてる。周りが結構とってるもんで、あるのが当然ぽくいわれて「え、無いの」とかいわれるとぐさっとくるね。そりゃさ、落ちるほうが多いのはわかってても結果の書類に「研究の将来性」とか「本人の資質」みたいな項目があってそれにCだのDだの書かれると鬱入るってもんです。わざわざ、ただでさえ将来への不安だので激下がりしやすい博士の学生のモチベーションだだ下げるような書き方すんなよなぁと毎年思います。毎年、その書類見る→あぁぁ自分って研究むいてねぇ→就活orなんかの資格修得を本気で考えるの流れ、そしてこの年までのんべんだらりと理論物理のみやってた人間がまともな職に就ける気もせずさらに気分は下降していくわけです。ていうか結果まだ来てすらないのにいまからダウナー入ってどうするんだか。
学振以外にも、とにかく11月てのが鬼門らしく難聴やら肺炎やら体調にもやばいことが起こるのは大抵10月末～11月。今年は平和に済みますように。でも風邪＋明日から遠出ってやばいフラグじゃ。
そういえば新聞記事にのった。別に犯罪をしたわけではない。写真とるときに黒い服＋ななめ立ちはエラの張りとかをごまかせていいことがわかった。文中に数値くりこみ群ってワードがでたので満足。ただ内容の一部に教授陣から不評くらうかもな部分があるので彼らに見せるのがちと怖い。

MinGW+MSYS

2008-10-25T21:24:00.004+09:00

Windowsの何がいやってコマンドプロンプトでlsとか使えないとこ。dirって何よそれ。Linuxのコマンド、CとかC++のコンパイラgcc、fortran用のg95なんぞをWindowsでも使えるようにできる超便利ツールがMinGW。ちなみにg95はMinGWもセットになったものを配布している＋最近Fortran2003のいくつかの機能も追加されたらしい。至れり尽くせりやん。
MinGW単体ではlsとかは使えないので、Linuxライクにコマンドプロンプトを使いたい場合はこれにMSYSも加える必要がある。これ入れとくとtarとかgzipも使えるっぽい。MSYSはvimまで積んでいる模様。
前のパソコンは暇暇にこれら一式インストールして設定してたんだが新しいのにまた設定しないといけない。めんどー。やり方が忘却の彼方だったのでまた検索してやり直し。懲りたので自分用にメモっておく。
1.MinGW編
現在の最新版はこれらしい→http://sourceforge.net/project/showfiles.php?group_id=2435&package_id=240780&release_id=595197
Windows用の.exeをダウンロード
クリックするとほとんど自動でインストール。こっちが選ぶのは構成要素とインストール先ぐらいか。gcc使いたい人はgccにチェックしておく。大量にライブラリとかの圧縮ファイルをダウンロードしてくるので、おわったら.tar.gzとかは消しておく。
おわったら環境変数いじっとく必要あり
システム環境変数に;C\dev\MinGW\bin （今回はdevってフォルダ下にインストールしたので）

2.MSYS編
最新版はこれっぽいな→http://sourceforge.net/project/downloading.php?group_id=2435&use_mirror=nchc&filename=MSYS-1.0.10.exe&57459211
MSYSのインストール先はMinGWとは別フォルダ作って入れろって書いてあった。同じとこに入れるほうが稀有だと思うけど。これもクリックするとほとんど自動進行。途中でコマンドプロンプトが起動して、
・インストールをすすめますか？　→　y
・MinGWはいってますか？ → y
・MinGWどこー？→　C:/dev/MinGW （向こうのフォーマットに合わせてこたえてやらんとエラーでる）
とかきいてくるので適宜答える。
同じく終わったら環境変数に;C\dev\msys\binとか追加。

3.g95編
http://www.g95.org/downloads.shtml
行って、Self-extracting Windows x86版を落としてインストールするだけ。
ほとんど何も考えんでも勝手に進行。

さてこれでLinuxなコマンドがWindowsでも使えるわ、g95でfortranコンパイルできるわ作業が楽だ。
ところで、使っているランチャー、CLaunchからコマンドプロンプトを起動すると環境変数の変更が反映されてないのはなんでなんだろ？OSおかしいんかとおもって焦ったけど、普通にスタート→すべてのプログラム・・・って普通の方法で起動したら反映されてるしなぁ。謎いわ。

奈良行き

2008-10-24T21:49:00.002+09:00

さて、奈良行ってきた。
正倉院展の内覧会？へ。なんか本来は明日かららしいけど開会式とかを今日するらしく、親の代理で。去年も行ったんだが正倉院展は理系がみても結構面白い。同行したのも同じく理系でたぶん専門家が聞いたら憤死しそーなコメントを垂れ流しつつぼけーっと巡回。古文書とかはさぱーりわからんが、装飾品やオリエンタルな感じの鏡や調度品は素直にきれいだなぁと。今回の目玉？は「白瑠璃碗」と「平螺鈿背八角鏡」ぽい。「白瑠璃碗」はガラス製の器で、「平螺鈿背八角鏡」はめちゃめちゃデコラティブな鏡だ。どっちもさすがに綺麗。ガラス器のほうはとくに色が薄い飴色でガラスなのに柔らかい感じがしてよかった。あと帯にじゃらじゃらつける飾りが出てて、魚型だったりカメ形だったりかわゆい小鳥型だったりいい感じなのでストラップとかアイコンで復刻希望。古文書は「お経とかはわかんねーしなー」って大半あっさり流したんだが、ひとつだけ、「奈良時代の休暇届や始末書コーナー」があってこれは面白かった。「家が壊れたのを修理するので休暇ください」「さぼってごめん、布納めます」とか。
で、ついでなので奈良公園をぶらぶらしてきた。奈良公園はたしか遠足で来たことがあるんだけどあんまり記憶にない。この年になってきてみると古くてでかい木がうっそうとしてたり鹿がいたりで何とも癒される。鹿せんべいを買って鹿とまったり戯れようとしたらまったりどころか一瞬で食いつくされた。なくなったら鼻でつつかれまくるし奈良鹿はまったく人を恐れてない。そして浮見堂が実にいい感じだった。池にお堂が浮かんでいるというのはなんともいい。公園の真ん中にあるのでビルや電線や道路が見えないのもよい。基本的に池や水辺がある風景がたぶん好きなんだな。奈良いいな。なんで最寄り駅から大和路快速で40分かからないのにいままでほとんど来てなかったんだか。

さて月曜からは松江だ。宍道湖に期待。準備しなきゃ

gnuplot :: フォントサイズとか

2008-10-23T22:30:00.003+09:00

動的帯磁率のバグもとれてあとはパラメータいじってひたすら計算機にあげるのみ。
今週の目標達成やれやれと思ってたら「論文マダー」がきた。
えー月曜から学会なんですけどとかブルー入りつつとりあえず図を作るかと思ってgnuplotのフォルダみたらbin以下全部なぜか消えてた。なんで？意味がわからん。というわけでこないだ作った図すべて消えた。涙目。bloggerにバッチファイルの内容をあげてたのですぐなんとかなったけど。blogに書いとくのはバックアップとしてもいいのかも知れん。一応文章化してるからあとから意味取りやすいし。

今日もgnuplotのいろんな小技を調べて時間が過ぎる。論文本体からの逃避行動なのかもしれないが。gnuplotは意外に機能満載だ。グラフ描かせるならexcelより高機能。そりゃそうか。

gnuplotつかってて困るのが凡例や軸のタイトルの文字が小さいことだ。gnuplotの出力はA4とかにグラフ一枚ばーんってのをデフォルトで想定しているのか、論文とかのサイズにすると「字がちっちゃくて見えん・・・」とかになる。図を小さくするのも一つの方法らしいが。eps出力してイラレとかでいじってたんだがめんどい。出力の段階でなんとかならんのかとおもってたらなんとかなるらしい。

・postscript eps出力ならset terminalのオプションとして”フォント名”サイズ　で、目盛とか凡例のフォントサイズが決まる。
・x軸y軸のラベルのフォント、サイズはfont "フォント名、サイズ"のオプションで決めることができる。
・Åなんかの記号もencodingつかって出力できてしまう。iso_8859_1のエンコーディングを指定して\305がÅに対応。

例（あとでイラストレーターでいじったりを考えてあえてArialフォントにしてみる。gnuplotのデフォルトのHelveticaだとたまに字化けするorz。Arialは確かWindows版Helvetica）
set terminal postscript eps enhanced color "Arial" 20
set output "dIdV.eps"
set encoding iso_8859_1
set xlabel "bias voltage (V)" font"Arial,30"
set ylabel "dI/dV (a.u.)" font"Arial,30"
set xrange [-0.005:0.005]
plot "a3.0" w lp lt 1 lw 2 ps 0.5 title "a=3.0[{\305}]", \
"a5.0" w lp lt 2 lw 2 ps 0.5 title "a=5.0[{\305}]", \
"a7.0" w lp lt 3 lw 2 ps 0.5 title "a=7.0[{\305}]", \
"a11.0" w lp lt 4 lw 2 ps 0.5 title "a=11.0[{\305}]"

参考文献
http://ryukyu.astr.tohoku.ac.jp/pukiwiki/index.php?Members%2Fchinone%2F%B3%D0%BD%F1%2FGnuplot#i90ad3c5
http://co.maxwell.jp/?page=gnuplot#p4

ギリシア文字も出力できるのか！

ClebschGordanを三行で

2008-10-21T23:07:00.003+09:00

なんか忙しくなってきた。あれー、余裕とか思ってたんだけどなー。新しいパソコン＝玩具が手に入った＋今週も金曜日に出かけるせいなんだろうが。
なんだかんだいって動的帯磁率も計算する部分書けたし仕事は進んでいるんだがなんか忙しい。で、夜に研究室にいると「今日は遅くまでいるんですね」って後輩に言われた。うむすっかり「あまりいない・とっとと帰る人」キャラなわけね。
で、「光のスピンとクレブシュゴルダン係数について」その後輩に聞かれた。ていうか光にスピンがあるんかい。そういえば量子暗号の話で、スピンの増幅に偏光をつかってうんぬんってのを聞いたなぁとか思って少々漁ると、右偏光・直線偏光・左偏光をスピンと同じくパウリマトリックスを用いて書けるんだそうな。光は当然電場を持つから、光を原子にあててやると電気双極子遷移かなんだかが起こるが、そのときの原子の角運動量変化にこの光のスピン的なものがかんできて結局クレブシュゴルダン係数が出てくるらしい。へぇ～角運動量の理論は応用範囲ひろいね。やっててよかったよ。その遷移確率の書き方とかなんか動的帯磁率に似ていて応答関数ちっくなものの形はどんなんでも大抵にてるもんなんかとか。そんなわけでクレブシュゴルダン係数を三行で説明しないといけない。B4にいきなり３ｊシンボル云々いってもしょうがないし、Mathematicaで計算できるよハイ終了も無責任だからなぁ。

というわけで適当にまとめてみた。
・クレブシュゴルダン（Clebsch-Gordan）係数はスピンの合成が絡んでくるときにつかう。
・S=Sa+Sb、Sz=Sza+Szbのような感じで角運動量の合成をおこなうとき、合成後の状態をいきなりひねり出すのはきっついので、合成前の状態|Sa,Saz>,|Sb,Sbz>の積（直積）で可能なものを集めてきてその線形結合で表わすことを考える。
・その時のそれぞれの項にかかってくる係数がクレブシュゴルダン係数で、漸化式とけば求められるけどすげぇややこしい階乗とかがバンバン入った形になるので普通は数表つかうかMathematicaにやらせる。

フォントとか

2008-10-20T21:53:00.002+09:00

そういやWindowsの画面のフォントってあんましきれいじゃないなぁと。なんとかならんのか。というわけでいまフォントがマイブーム。なんじゃそりゃ。
メイリオ＋ClearTypeで滑らかにってのをしてだいぶマシとはいえ、まだLinuxのkterm端末画面のほうがフォントなめらかってどうなん・・・。Linuxとかイラストレーターの画面みてからWindowsみるとなんかがっくし・・・。なんでもWindowsは小さい文字についてはアンチエイリアスしないからどーしてもあまり普通に使うフォントサイズのはギザギザするらしい。LinuxでもがんばってんだからWindowsもそれぐらいしろよと。なんかあとClearTypeTunerとかいうのを落としてきていじったらちょっとはましになるっぽい。M+とかIPAフォントちゅーのも試し中。

そういやWindowsでSSH接続させるのによくつかわれるPuttyもフォントのアンチエイリアスができるらしい。フォントいじってアンチエイリアスするだけでたしかにずいぶん見やすい。なんかLinuxぽくなる。viもつかいやすくなることだろう。

いろいろ周りが不満タラタラのOffice2007のPowerPoint使ってみる。意外に良い。フォントの汚さもマシ。Office２００３系統よりずっと描画とか種類増えたし、クイックスタイルが便利だ。2003で泣かされてた好みの小塚ゴシック系つかっても行間がばかでかくなるっちゅー問題がないのが良い。図形とか効果とかすげー多いのでかなりお手軽にきれいなのを作れてよい。高いだけあるわ。なんであんだけ文句いわれてるのかわからん。あー大学のもこっちにしたいな。そういやKingsoftOfficeてのもなぜかプリインストールされてたので使ってみたがあまりに２００３のパクリすぎて大丈夫かよこれって思った。なにかと怪しい。

あたらしいのが来た

2008-10-19T20:57:00.003+09:00

パソコンが来た。
こないだから調子わるかった（で、なぜか今調子よくなってる意味フ）家の後継機。PCないと仕事できないので壊れる前に思い切って買った。なんつーか普段デカイ買い物をあまりしないので、ものすごくてんぱった。
メーカー製のは某V●IOですごい痛い目にあったので、避けようと。いままで使ってた日立製Prius７７０Dはそこそこ頑張ってくれたし液晶もきれいでよかったんだけど。メーカーのて高いしなぁ。
何買っていいのやらさっぱりわからんかったので相談してみて、Faithで売ってたIntel Core2DuoE8500+Geforce9800GTにメモリ４G、HDD５００Gっていうのにした。グラボってなにさ～レベルな私にはもったいない構成なんかも。たしかにサクサク動く気がする。なんか今は結構CPUもメモリも安くなってて買い時であるらしい。ほぉほぉ。ディスプレイやOS（WindowsXP）やOfficeつけて全部で１６万ぐらいだった。あとで調べたらS●NYとかだったら同じぐらいの性能が３０万弱なのでちょっとニヤついた。

で、今日届いたので、まずWindowsのアップデートしてOffice入れてFirefox入れて・・・ってして一段落したあたりでちょっといらんことしたくなって（見た目をVista風にしたくなった）いじったらOS起動しなくなってめっちゃめんどくさかった。とりあえずXPのCD-Rほりこんでブートして修復すりゃなんとかなったぽい。修復だったのでOfficeも消えなかったしやれやれ。でもそのせいかグラボのドライバがぶっとんだりそれ入れなおしたりめんどくさかった。いまでもこれでちゃんとなおってるのやら不安でしょうがない。素人がいらんことしちゃだめっすね。OS入れなおすとかあんまりやったことないから始終「あぁー買ったばっかしやのにこわしたらどーしよぉ」とか涙目やった。まぁええ経験だけど。ブルースクリーンが定番D●ll様のおかげでパソコンがやばくなることにこの２年でちょっと慣れといたのはよかったかも。
PuttyとWinSCPとUltraVNCとサクラエディタも無事入ったしこれで仕事ができる構成になったかなぁ。
いじっててわかったのをメモっておくと、
１．UltraVNC使うように解放するポートはVNCビューアー用が５９００でブラウザ用が５８００が標準。
２．XPにもメイリオフォント入れれるようになった。
３．フォントがぎざぎざして気持ちわりーのは「コントロールパネル」→「デスクトップの表示とテーマ」→画面→デザインのタブ→効果→フォントの縁を滑らかにする方法→ClearType選択　　でちょっとましになるっぽい。

あとはイラストレーターいれて、その後MinGW+gcc、VMwareとかちょいちょいいじるか。ややこしいの入れるのはやっぱり出張から帰ってきてからかなぁ。またOSお亡くなりとかいやだし。せっかくCore2Duo上位タイプにしたのでVmwareのLinuxでifcいれてLapack+blasコンパイル、数値計算とか重たいことさせてみたいもんだ。研究室から出たときにそういった環境をフリーで作るノウハウもっといたほうがどこいっても強いとおもうし。それにしても新しいの＋なぜか旧機が調子よくなって使えるPCがこれで２台という状況はいろいろできてしまいそうな気がする。

効率を上げたい

2008-10-16T22:50:00.003+09:00

ここんとこ、「美容院いってくる」とか真っ昼間に帰ったりしてたら研究室で何人かに「久しぶり」って言われた。最近来て無い人イメージなんかい。いや、土曜日来てたりもしたんですけど。そして超めずらしく土曜きたら教授に「どないしたんや？」ってものすごい言われたんですけど。焦って土曜日行ったりした効果でサーバーが何日か停止したにもかかわらず余裕を生み出せてたり。

基本、理論系の研究者ってのは裁量労働制ちゅーかあんまり何時間研究室にいるかってのに関係ない気がする。どっちかっちゅうと集中力勝負だと私は思っている。それにうまく行かない点の解決策を思いつくのは実は風呂だったり電車だったりご飯中だったりするし。そして研究の仕方ってのはツールを使えばもっと効率化できる気がする。

昨日テレビみてたらクラウドコンピューティングで業務管理云々って話をしてた。へぇー企業の営業ってどこに電話かけてどんな話してとかをそこまで細かくデータベースに管理するのかぁとか思った。どの企業もベースがおんなじツール使うってのはノウハウとかの均一化生みそうでそれってどうよとかも思うけど、管理ツールが見たとこ出来がよさそうで、いつどんなことをして、スケジュールがどうでってのをブラウザかなんかで一括管理できるようになってる。研究者用のそういうの欲しい。まぁカネにならんしあまりにもニッチだから誰も作ってくれないだろう。なので転用できそうなのを探すかと。

研究においても、「大体何日に何が締め切り（学会とか）」「いまどんなテーマでやってるか」「結果を得るにはどういった工程が必要か」「これまでに得られたデータ」「現状と問題点」とかを整理しなあかんわけで。それをしとかんと学会前に死亡フラグ乱立とかになるわけで。まぁ頭の中でなんとなくやったり、メモつけたりしてはいるけどやっぱりどことなく能率悪い。スケジューリングとタスク管理と作業の分解と問題点とかをまとめてうまいこと記録しておく方法が欲しいわけだ。スケジュール＋タスク管理はそりゃいろいろあるんだが、それだけじゃだめなんだよな。経過とかいろいろ記録する部分もほしい。Ｗｉｋｉ的なのがあれば最強。なんか一番近いのがＩＴ系の人がやってるプロジェクト管理っぽかった。理論系だといやでも計算用のコード作ることがおおいから、あながちはずれではない気も。所詮Ｆｏｒｔｒａｎですがなにか。フリーで使えるの探してみたんだが、自分とこのサーバーに導入するってのが多くて、サーバー立ち上げるなーんてうちにゃ無理っすと諦めかけてたところ9Arrowsっていうのが今すぐ使えるＡＳＰ版あったので試してみることにした。さて効率があがるんだろーか。一ヶ月ほど試してみよう。

readのIOSTAT

2008-10-15T21:32:00.003+09:00

パソコンが不調だ。ガリガリ言う。起動時が最もひどく、しばらく動かしてると安定⇒思い出したようにガリガリガリクソン。卒論・修論・その他もろもろをこいつで作ったのでそれがそろそろ終了っぽいのはなんか寂しい。ＭｉｎＧＷつんだりまた一からかーメンドクサ。まぁ今はなんとか動いているけどイラストレーター起動に5分ぐらいとか処理速度が遅いのも問題だしガリガリ音は精神衛生上よろしくないので買い替えないといけない。無いとは思うけど動かしっぱなしで火を噴かれたりとかいやなもんでなるべく使わないようにしているこの数日。ネット巡回とかしない分、早寝だし本読むしたまにはパソコン使わない方がいいのかも。

ところで、Fortran90のread文には3番目にIOSTATという引数があるのを知った。readでなんのエラーもなくファイルを読み込んでいるときはIOSTAT=0で、ＥＯＦ検出ではIOSTAT=-1になる。まぁほかにも色々あるっぽいが。なので、何個データあるかわからんファイルから読み込んで処理するときには、IOSTATで条件分岐させるのが便利ぽ。ただしFortranは可変長配列をサポートしてないというアレな言語なので（それさえクリアすれば数値計算もっと楽になると思われ、Fortran2003では対応してるとかの噂。Fortran2003はオブジェクト指向もできるとかなんか凄く拡張されてるしSunもＮＡＧもコンパイラ作ったみたいだしつかってみたいやも）90の段階ではいったん行数をカウントして、割付して・・・とかちょっと回りくどくやらなあかんぽいけど。以下つくってみたのでサンプル。data.txtから数値データ読み込んでフォーマットしてformatted.txtに書き込む。

program main
!program for setting the data file formatted form
implicit none
integer ::Ncount,i,j,ISTS
real(8) :: temp
real(8),allocatable :: a(:)
open(11,file="data.txt")
Ncount=1
do while(.true.)
read (11,*,IOSTAT=ISTS) temp
print*,temp
if (ISTS /= 0) then
print*,"IOSTAT=",ISTS
exit
end if
Ncount=Ncount+1
end do
close(11)
print*,"Number of low+1=",Ncount
allocate(a(Ncount-1))
a=0.0d0
open (12,file="data.txt")
do i=1,Ncount-1
read (12,*) a(i)
end do
close (12)
open(13,file="formatted.txt")
do j=1,i-1
write(13,'(E22.14)') a(j)
end do
close (13)
end program

経済と物理

2008-10-10T22:11:00.003+09:00

なんか経済がすごいことになってるみたいだ、恐ろしい。経済はさっぱり分からない。周りにちょくちょく株いじっている人がいたのでいろいろ聞いてみたけどやっぱり分からなかった。

分からなさの原因は定式化できないからだと思う。物理はまぁ複雑系とかもあるけど、相互作用の形とかをそれなりにモデル化できる。固体物理業界では少なくとも原子スケールでの相互作用というのはある程度決まっている。問題は系のでかさに由来する数値的な問題をどうクリアするかってのが大半で、頑張れば何とかなりそうな感じが漂ってるし、基本的な式はすっきりしている。経済はなんかパラメーターが何個あるかすら予想できん。

が、経済物理学というジャンルがある。株式市場の暴落を一種の「臨界現象」とみなしていつ暴落が起こるかとかを予測しようというものらしく、物理学会にもちゃんとセッションがある。なんとこれには「くりこみ群」もつかわれてるらしい。手元にあった「入門・経済物理学」にもくりこみ群の紹介が載ってた。びっくり。ちら見したところによるとほとんどイジングモデルだ。トレーダーが他のトレーダーを模倣する頻度みたいなのを相互作用Ｊみたいに扱い、強気筋・弱気筋をスピンのアップダウンとおなじように扱って実空間くりこみ群ちっくなことをやってる。かなりストレートな方法でなんだかそれっぽい株価の変動グラフが出てきてるんで驚く。人間の挙動もスピン系と似たようなもんなのかいな。経済物理学ではくりこみ群＋αでランダムネスを取り入れようとかしているようで、その手法は逆に物性系でも使えるかもしれない。なんかちょっと興味をもった、調べてみるか。いやそれより締め切りその他諸々が臨界点近傍で明日は休日出勤せなあかんねんけどさ・・・忙しいときほど別のことしたくなるんだよな。

sed+gnuplot

2008-10-08T21:47:00.003+09:00

ノーベル賞日本人4人ってすげー。
南部先生といえば南部演算子とかしか知らんので、てっきり超伝導の人だとおもってたら素粒子の人だったんだ。超伝導も自発的対称性の破れってやつなんかぁ。超伝導は食わず嫌いというか手を出すと危険な香りがぷんぷん漂う恐ろしい領域なので踏み入ったことが無いからよく知らん。それにしても素粒子と物性ってなんだかんだ言って理論のベースは近いんだな。応用になってくると全然違うほうを向いてる気がするけど。見慣れない用語に毛嫌いしてんとマトモに場の理論とかも見といたほうがいいんかもなー。そして物性の人間としては近藤先生のノーベル賞まだー？と言いたい。

gnuplotはsedコマンドと組み合わせて使うと図の大量生産がラクだ。
ちなみにコマンドプロンプト経由でgnuplotを起動する場合、Windowsではwgnuplot.exeがあるディレクトリにパスを通して、wgnuplotで起動する。（Ｌｉｎｕｘだとgnuplotなんだが）
sedはＵｎｉｘコマンドなんだが、ＭｉｎＧＷなりＣｙｇｗｉｎが入ってるとコマンドプロンプトで使える。
sed -e "s/x/y/g" aaa.txt > bbb.txt
のようにするとaaa.txt中のｘをｙで置換したものをbbb.txtとして出力するということになる。
Linuxだと
sed -e "s/x/y/g" aaa.txt | gnuplot
でsedで置換して直でgnuplotにほり込めるみたいなんだが、Windowsではなんかうまく行かない。
たいてい論文とか学会用の図を作るときはある一つのパラメーターだけが違うデータで似たようなグラフを何枚か作るというパターンが多い。ので、sedでパラメーターの部分だけ変換したgnuplot用のバッチつくるバッチファイルを作る＋バッチをまとめてロードするバッチをつくる（文にするとわけわからんな）をつくってやると楽。サンプルを自分用にあげとくと、

固有値の結果をプロットする用（a=5.0って条件）
set terminal gif
set output "eigen_odd_a5.0.gif"
set xrange[1:41]
set yrange[0:2.0]
set xlabel "iteration N+1= odd"
plot "odd_eigen_5.0" w p ps 0.5 pt 7 lt -1 ti "a5.0"

eigen5.0.gpって名前で保存（ちなみに点でプロットするときもその色はlt で決める。lt はラインタイプってだけでは無い模様）
↓
sed でひたすら変換してる部分作る

sed -e "s/5.0/3.0/g" eigen5.0.gp >eigen3.0.gp
sed -e "s/5.0/7.0/g" eigen5.0.gp >eigen7.0.gp
sed -e "s/5.0/9.0/g" eigen5.0.gp >eigen9.0.gp
sed -e "s/5.0/11.0/g" eigen5.0.gp >eigen11.0.gp

これを実行してから

load "eigen3.0.gp"
load "eigen5.0.gp"
load "eigen7.0.gp"
load "eigen9.0.gp"
load "eigen11.0.gp"

って書いたファイルをgnuplotにロードしてやる。と、gif画像が5枚一気に生産される。自分で打つ部分がちょっとでも少なくなるとなんか嬉しい。
ちなみにgnuplotでファイルに出力したグラフはgnuplot起動中には見れなかったりするので、gifを確認したかったらいちどgnuplotを終了しないといけない。

一言ほどに難しい

2008-10-06T21:20:00.004+09:00

たとえば自分の短い発言がなんらかの媒体にのるときに、チェックしてくださいといわれて「え、こんなん言ったっけ、ちゅーか若干言いたいこととちがうぜ」ってときにどう対処するかとかは実に面倒な問題だなと。向こうの事情とかも忖度すれど、自分が言いたいことと違うことを自分のコメントとされるのは嫌だ。でも交渉するのもめんどい。代わりの文かんがえなきゃだめやもん。しかもコンパクトに。しかも相手は文筆業のプロだろ。とはいえこのまま載せられて「他人の高尚なのとくらべておまえはなんちゅーアホなコメントしとんじゃ」と周りにあきれられるのも嫌なのでちゃんと言わんとなぁ。20文字程度で研究のミリョクを語れとか無理ってもんですよ。そりゃ宇宙の階層性うんぬんとか言い出したらアカデミックな薫りただようそれっぽいことはいえるかもしれんけど、カッコつけるのいやだし。高尚な目的があってやってるわけじゃなく、なんかやってみたら面白いっぽくて延々おたっきーなことやってるって感じだからとてもえらそーなことなんておこがましくて言えない。「くりこみ群ヤバイ」、「近藤効果パネェ」が私が訴えたいこととしてもそれを載せてくれるとは思えないしな。女性に向けてのgeneralなコメントなんてねぇっすよ。さて、ちょっと笑いをとれそうな？文言を考えるか。

秋眠夕暮れもおぼえず

2008-10-05T20:21:00.004+09:00

毎月遠方に出かける用事があるというのはこたえるようで眠気とかが半端無い。家で作業してても、「んーわからんな、ちょっと考えるか」ってソファでごろっとすると考える間もなく寝てたりする。仕事で強制的旅行が多いシーズンは大阪にいるときの引きこもり指数が増大する傾向がある。もう家がfixed pointって感じ。そういえば出先ではあんなに観光したおすのに、関西圏の名勝はほとんど行って無い。学生時代の遠足もさぼりまくってたからな。毛越寺までいってるくせにそのモデルの平等院鳳凰堂は見たことないという。なんかもったいない気もするな。

なんか文科省が女性を採用した研究機関には年間600万とかいう案だして叩かれまくってるらしい。女性研究者を増やそうとかいうの多いよな最近。あまりロコツにやられると逆に当事者に風当たりがきつくなったりするんじゃないかね？まぁうちらその恩恵を頂いてるんでなにもいえないですけど。アカポスにいる先輩曰く、工学系で女性をふやそうというのは、フェミ云々以前に、実に単純にこれから人口が減っていくから絶対数を確保するためみたいだ。少子化でこれから絶対数を確保するには未踏のエリアを開拓するしかないってことだな。そういうわけであんまりフェミとかを勘ぐって叩かんといて欲しいもんだ。そして女性が理系に少ないのは、女が理系っていうとドン引きされたりするのが原因だとおもわれる。女性で理系に行くとなんか損しそうな気配が世間に漂ってて、普通は早々に気づくぽい。脳の基本スペックには男女さほど違いは無いかと思うんだけど。研究に必須の、しつこく一個のことを考える性質とかむしろ女性的だと思うんだけど。

まぁそんな流れにのって某大学も女性の支援機構的なものを作ったらしく、そこのイベントで企業ではたらく女性研究者の話を聞くってのがあったんで行って来て見た。話はいろいろ面白かったんだが、なんていうか女性の支援＝育児出産支援っていうのになんかね。質問とかもそれ関係ばっかしだし。重要事項なんだろうが女ってそれしかないみたいにカテゴライズされるのもやだなぁ。うちはむしろそっちにはいまのとこあんまり興味ないっす。女性むけ雑誌とかも仕事と育児の両立とかひたすら前面に押し出してるのとか、一歩引いてしまう。人間そんなに全方向に器用に頑張らなあかんのだろうか・・・どっちかでも十分すぎるやん。もっと気楽に行こうぜとか思う。

parity

2008-09-30T22:42:00.003+09:00

4日ほど悩んでたparityの問題が片付いた感じ。定義とか論文に書いてなくて困ったんだがやっぱり原義通り空間対称性の偶奇ぽい。原義そのまんまだから書いて無いんだろうけどさー知らんかったら分かるかっちゅうの。2不純物アンダーソンモデルだと、不純物については、不純物間の中央を原点とすれば、2不純物の演算子の線形結合でeven（偶関数）とodd（奇関数）の波動関数を定義できる。伝導電子についてもｋ空間でoddとevenのチャンネルを作れるぽいな。対角化時に各基底がoddまたはevenのパリティを持ち、ハミルトニアンも空間対称性を持つことから、同じパリティをもつ基底間にしか行列間要素ができないというわけか。パリティ対称性を使っても行列を2分割するようにしかはたらかないので劇的に計算速度があがるわけじゃないが、それでもいまの丸3日かかる状況からは脱せそう。そしてこれ以外の落とし穴も発見して埋めた。今度こそちゃんとした結果がでるはず・・・。それにしても対称性は奥がふかいな。対称性によって固定点が安定か不安定か、臨界点が存在するかとか変わってくるらしい。やばい、群論とかやってないから用語がわからなすぎる。そろそろ対称性やら固定点とか繰り込みの流れとかの扱いを勉強せねば。いままではひたすら計算方法を理解してプログラム組んでがーっと計算させてるだけだったんだが、ちゃんと考察するにはその辺のことも把握しなあかん気がする。はぁー、計算コードできても一段落ではなかった。近藤効果は底が知れない。

勉強せんとなぁとちょっと切実なのは、レビュートークとやらが近いので、30分しゃべる内容を用意しないといけないせいもある。院生相手に30分かぁ、どうせ皆寝てるだろうけど。というかその前の週は学会で島根いってるんだがいつ準備しろと。

さて明日からＴＡが始まる。Ｆｏｒｔｒａｎ屋がＪａｖａベースの計算機実験のＴＡでいいんかよって感じだが。ＴＡは教科書をかってもらえるとの噂。給料もらって勉強できるんだからだいぶいいかも。

おおまちがい

2008-09-27T19:30:00.003+09:00

なんだかんだいってるうちに9月最後の週末だ。月の最後の週末はいつも家業で線香くさい日々なんだよな。なんか4月以降妙にいそがしくて半年経ったとは信じがたい。今年はやけに濃い一年かもしれない。Ｄ２で忙しいとかいってて来年Ｄ論あるのにどうなるんだ。

とか感慨に浸ってる場合ではない。すごいミスを見つけてしまって、やっべーこないだ学会発表した結果まちがってるかもしんねぇっていうこれはちょっとまずいんではという事態が出来した。木曜の夜中にいろいろ考えてるときに、「あーあのへん怪しいなぁ」とふと嫌な予感がして考察しなおしたら、計算の初期設定にミスを見つけた。あぁぁーヘンな対称性とか制限を入れたようなモデルになってた、まじがっくし。まぁそういうモデルではこうなりますってな結果として発表するのもアリかもしれんけどさ。これも一重に参考にした論文のミスプリが原因だったりする。教科書や論文を信じてはいけないな。1人でやるのはこういうことが怖いのだ。うっすら、「うーんここ自力でチェックしといたほうがいいんじゃ・・・」と思ってて、ま、いっかーで流したとこなのが悔やまれる。いやそこ流すなよ。対称性の間違いはトンデモなく結果をずらすのではないかと恐れていたが、いまんとこ途中結果を見た感じ、定性的な違いは出てなさそうなんだがそれはそれで不安になるな。たぶん近藤温度が変わるぐらいぽい・・・。初期設定でエラーしてて、定性的結果が変わらないとしたら一体なにが系の傾向を決めてるんだって言う・・・。それがくりこみ群なのか？まぁなにが問題って対称性を一個使えなくなったせいで計算が異様に重くなって時間がかかるってことだ。それよりパリティ保存量ってなんやねん。

そしてこないだ取材を受けた。記者のひとって質問うまいねー。理系のことはわかりませんとかいうてはったけど、要点を的確にまとめてはったのが凄い。頭の回転はやいんやな。自分やったら専門外のこと延々いわれたら寝るか、意識どっかに飛んでいくわ。院生してると基本、自分の興味あることしかしないからそういうとこでの忍耐力が減るのがいかんねぇ。さてどんな記事になるんだろう。

岩手ダイジェスト

2008-09-23T21:51:00.010+09:00

岩手すごいわ。
学会で岩手は盛岡に行ってきたわけですが、岩手は見所満載ですな。
8月が異様に忙しくてなんら下調べせずに行ったら、学会の受付に観光ガイドが満載で、見てみると、これは行っとかねばという観光地が多くて大変でした。だって東北に行くチャンスってなかなかないし。
龍泉洞、小岩井牧場、平泉は絶対行きたいと思ったので、気合で予定をねじ込んでまわりました。発表の合間をぬってな弾丸ツアーな日程で若干疲労困憊かも・・・。写真は平泉の毛越寺、龍泉洞、小岩井牧場。
龍泉洞は、地底湖が有名な鍾乳洞。地底湖の透明度が凄いせいで、９０ｍの湖底が透けて見える。そしてその青さが怖いぐらい。鍾乳洞はもともと興味深々で、機会があれば行ってみたかった場所。龍泉洞は地底湖の存在とあいまってとくに神秘的。水の滴りが何万年かけて石灰岩をけずってそれが地底湖になってとか考えると眩暈する。地底湖のひんやりした青さ、奇怪な鍾乳石、もうなんかすごい絶景だ。写真にうまく写らないのが残念。
小岩井牧場はのどかだった。山を背景に牧場ってもう心なごむわ。牛乳にがてなんだがここの出来立てのは美味しかったな。小岩井ってのは創業者3人の頭文字らしい。知らなかった。牧場内には、森なんかもあるのだが、できた当時は、木の一本もないひたすら荒地が広がっていたらしい。そこに一から苗木を植え、牧草をそだて、100年かけて今の風景になったらしい。荒地の状態から牧場を作る苦労とか想像もつかないレベルだと思う。それをはじめようとビジョンを立てて実行する力は凄い。創業者は三菱の当主や日本鉄道の社長だったみたいだが、渋沢栄一といい昔の大会社のトップには凄い人が多いな。たぶんなんか収益とかを越えた理念があったんだろう。財閥という制度は弊害もあったかもしれないけど、普通の企業にできないこういう公共性高いデカイ事業をやったりしてたのは良かったんじゃないかと思う。今の日本の企業でもそういうことやってたりする例はあるんかね、知らんだけで、まぁ見た感じなさそうだけど。それにしても町ひとつが入るほどの牧場。スケールがでかすぎる。
そして平泉。奥州藤原氏。平泉の一帯を世界遺産にしようとかで整備も進んでて、巡回バスあったり、予想以上に観光しやすかった。
日本史でむかーしならったなぁ、前九年の役とか・・・。新幹線のなかに入ってた雑誌にくわしく奥州藤原氏のことが書いてあったので、予習していった。高校時代には「歴史・・・暗記科目つまんねぇ」だったが、テストとかから開放されて、観光とかで触れると結構面白いかも。平安末期～鎌倉の最初のころの平泉はすごかったみたいだ。金山があったのを背景に、北方とも交易し、京都をしのぐ経済力を持っていたらしい。現在の価値にして15億をポンとだして中国から経文を買ったりとか。その凄さの一端が中尊寺金色堂や毛越寺とか。今は首都東京にすべてが集中し東京が全てという感じだが、昔はそれぞれの地域に首都に匹敵するパワーをもった都市があって、オリジナルな文化を持ってたんかな。なんかそっちのほうが断然面白そうだ。首都一極集中はよくないと思う。地方都市が東京ベースのチェーン店の進出とかで均一化するのは面白くない。

まぁそんなこんなで岩手は面白かった（学会発表もちゃんとしたって。）関西にいると東北のことをあまり知らないが、面白いところがいっぱいありそうだ。また行くチャンスないかなぁ。

台風と共に

2008-09-18T22:14:00.005+09:00

明日から岩手・・遠い。そしてすごく風邪。出先で風邪とか熱とかまじ避けたいんやけど。台風もやってきたしお仕事もふってきた、ヒィィ。いろいろ無理っす。台風の影響で飛行機揺れたらやだなぁ。忙しいほうが日々まぎれていいけど、ここまで一点集中しなくても。むしろしばらくグダグダ生活するつもりだったのに。なんかＧＣＯＥがらみのあたらしい企画も始まるそうで、今後どうなるんかなぁとか考えるとさすがにテンパるわ。絶対これ来年一年で終わらんって。ていうか年内にせなあかんこと多杉。どうするんだ。そして私でていったあと後継する人いるんかなぁ？

なぜか岩手に行くのに、直行便がとれなくて仙台経由。東北新幹線に初めてのることになりそうなんだが、学割証を出し忘れたorz。学生の利点って学割ぐらいしかないっちゅうに。まぁそんなこんなで4泊5日を無事乗り切れるのだろうか。それより荷物詰めなければ。