デスマから解放されて一息つきたかったので、今日は麻雀ぶちましたよ。学食の外のテラスで。 テーブルが丁度いいくらいの大きさだったんで、雀マットを敷いて、そこでぶちました。学食のオープンテラスで麻雀打つなんてマジで頭悪い感じですが、麻雀の発祥の…

今日やっと色々な物から解放されました。といってもやることはこまごまとあるんですけどね。資料のコピーとか研究会の打ち合わせに東京まで行くとか。申請書を書くとか。 そして昼間は健康診断に行ってきました。内科検診と、尿検査と、胸部X線くらいしかや…

最近Web2.0とかいうのを良く見かけますが、一体何のことやらさっぱり分かりません。 何かWeb1.0のビジネスモデルはAmzonみたいに物を売る商売で、Web1.0は広告収入とか、アフィリエイトで稼ぐ商売だとかいうのを見た気がするけど、これもまた定かではないし…

なんつーか東京都民にとっての東京都の地図っていうのは電車の路線図を元に構築されてる訳ですね。っつーことで実際の地図と、自分の方向感覚とか距離感というのは地図上のものとは異なってくる訳です。 で、GoogleMapで実際の地図の上で地下鉄がどこを通っ…

今日は暖かくていい日ですね。 筋トレして、ボスとグダグダ喋ってたら半日終わりました。妙に疲れた。そして学会の参加申し込みを4月の末までにやらないといけないことが判明しました。あとはあれだな、今年は京都に行きたいな。ということを次にあったとき…

宮崎あおいが出てます。マジやっベー。萌えるす。萌え氏ねますよ。なんつーか宮崎あおいってこういうのって当たり役ですね。超嵌ってるんですけど。 因にこのCMを見てからチョコレートは明治のを買うようにしてます。なんつーか嵌められたね、俺。そしてこの…

某K君がSkype(・∀・)イイとかいうので、僕もSkypeをDLして、デバイスも買っちゃいました。 どうやらSkype持ってる人同士の間では無料で通話が出来る模様。そしてSkypeから固定電話、携帯電話に掛ける場合は有料らしい。これで固定電話や携帯にかけられれば最…

なんつーかここ一週間デスマってたので、疲れが溜まりまくりんぐです。昼から朝の5時まで研究室とかマジで頭悪すぎです。最近ずーっとそんな調子です。そして色々な日程の都合から俺が頑張って組んだコードは今回は日の目を見ること無く終わりました。とかい…

hoge重塔っていえば、薬師寺の東塔とかありますが、僕の中では上野動物園の正門から入ってすぐ左手にちんまりとある奴が印象にのこってます。確か小学校の写生会で行った気がする。それでその前に授業で、「五重塔は中に一本柱があって、それが地面から浮い…

行列のかけ算は、テンソル解析の書き方で書くと、 ですね。ここでEinsteinの総和規約を使ってます。 っつーことで、数値的に行列の積を求める場合はこれを書き下すだけでOKです。実際かなりふざけたコードになります。なんつーか上の式そのまんまみたいな。 …

今日は10-10でタッチフットしました。人が多過ぎて訳分かんねー感じでした。 あとは全体的にラインが浅過ぎてどうしょうもない感じでした。もっとためを作らないといけないんだけどなあ。でもあんまし体力が無いのでそこまではやれないのですよね。 っつーか…

今日はこれからラグビーの練習です。その前に学校にきて、俺が書いた壮大なる研究のシナリオを進めました。やっと半分まできました。4段階にフェーズをわけたのですが、最初の第一段がすげー大変でした。新しいことをやるから。でもそれが出来てしまったので…

やったー、やったよままん、自分でコード組んでCLAPACK動かせたよ。固有値求めるところまではうまくいったんだけど、その後の固有ベクトルを求めるところが上手く行かなくて、ずーっと試行錯誤してて、最終的に作業領域の大きさの指定の不整合で上手くいって…

頭の文字でどの変数の型かを選ぶみたいです。zだと複素数で、dだと倍精度実数みたいな感じです。 その次に、行列の種類を選ぶみたいです。一般行列はgeで、エルミート行列なたhe、対称行列はsyでみたいです。 そして最後に求めるものを選ぶ感じです。ぼくの…

ゆで卵ってばくるくる回すとなんとなく飛び跳ねそうな感じがしますが、それについて慶應大学の下村先生が実証して、それがイギリスのRoyal Societyのプロシに掲載されたみたいです。なんでも卵を1800(rpm)*1で回転させると宙に浮くとか。 なんつーかこんなも…

そして今日もNHKの朝ドラの再放送を見てから登校です。今日は桜子のとーちゃんが働き過ぎで入院したとかそういうのですた。あんな娘にお見舞されて、りんごむかれた日にゃあんたもうどうしろっていうんですかって感じですよ。もうね宮崎あおいの三つ編みがマ…

今日も数値計算です。 どうやらC++のプログラムで、exten "C"とかやって、そこから何かを書くと、その部分はC言語のコンパイラを呼び出してコンパイルしてくれるらしいです。なんつーかこの一週間で数値計算についての知識がかなり増えました。 そして、どう…

行列の固有値が実数である場合、つまり固有方程式の判別式が正の値をとるときのことを考える。3次元だと面倒なので、2次元で考える。 という風に任意の二次の正方行列の固有値が求まる。 ここでこの固有値が正の値をとる場合について検証する。 上の複号で-…

いい天気です。学内には新入生っぽい人が沢山居ます。なんつーかまだスカートはいてる人が多いですね。そのうちどーでも良くなってジーンズになるんですけどね。そういえばズボンって最近いわないよなあ。パンツとか、ボトムスとかいうよなあ。どーでもいー…

なんつーか、固有値に複素数を許す場合の固有値の問題がすげー楽しいです。 別に行列が複素数でなくても固有値が複素数になるような場合は沢山あります。今、行列Aを、 として、固有値をλとすると、固有方程式は、 になります。 ここで中学とか高校でもやる…

シンガポールの武術の人達が白神山地で達人を探して遭難したっていうネタがありましたが、それの秘伝の書が見付かったらしいです。これって確か継承者探しだかなんだかっていうネタでネット上で話題になってた気が。 今日日山篭りして鍛錬積んでる奴なんてい…

取り敢えずサンプルを動かすのには成功しました。参考にしたのは、 http://www.k.mei.titech.ac.jp/~stamura/NumericalComputation-Tips.html#compile-LAPACK http://www.cazoo.jp/calcsoft/library/lapack/clapack/usage/dgeev/index.html で、上の方はコン…

今日はしけってますね。春のしけってる感じも悪くはないかも。 そしてこないだ家に帰って来たら、鞄にどこかからひらひらと落ちて来た桜の花びらがしけってくっついてた。なんとなく風流なので写真にとってみました。 今朝はボスの授業のTAでした。何かMcrau…

Numerical Recipesをはじめ、色んなところで、C言語では複素数をサポートしてないので、実部と虚部を構造体のメンバにして、そいつらの間の計算の規則を云々とか書いてあるわけだが、そんなダサダサなことしないで、人様の作ったパッケージを使った方が良い…

なんつーかLAPACK++のインスコに挫折しました。だってmakefileどうやって書けば良いか分かんねーんだもんよ。でもってLAPACKの方には既にLinux用のmakefileが入ってるので、それを使いましたよ。 っつーかLAPACK++のマニュアルが超不親切で、「makefileの設…

今日からLinuxで作業です。なんつーかMacのインターフェースが懐かしい。でもLinuxも最近はフォントがアウトライン出てて見やすいので良いとする。Windowsは諸々不便なことが多いのであんまし使う気にならないっす。メーラとかあんまし良いのないし。 そして…

なんつーか、訳あって固有値固有関数を数値的に求める必要が出たのですね。それでMathematica使ったら、N次の正方行列をAとして、 Eigensystem[A];とやったら、それだけで、 {{Eigenvalue1, Eigenvalue2, ..., EigenvalueN}, {1st_Element_of_EigenFunction1…

現在境界値を出す作業をしてます。 特異摂動法をさらにややこしくして、さらに境界層内部をすげー回りくどい解き方で解いてるので、境界層外縁での境界条件がど偉いことになってます。もうね、今シナリオを書いてるんですが、4つくらい手続きを経ないと到達…

何か最近雨が多いっすね。桜の花も落ちちゃう感じですね。 昨日は電力会社とか、ガス会社とか、水道局に電話して、諸々止めてもらいました。そして大家さんちに電話したんですが、大家さんは不在でした。で、今日大家さんから電話がかかってきたんですが、全…

これまで固有値問題って何となく固有値を決めて、それを方程式にぶち込んでってやってたんですが、行列を扱うようになってからちゃんとしたformulationが必要になりまして、ちゃんとやることにしました。 今、解くべき方程式は、 で、ここで、解を、固有ベク…