Jan/18/2005->
やるよ〜:
GRAPE-7(5) では、粒子分布を包み込む立方体を用意して、その最遠頂点同士を結ぶ距離で最小質量粒子(mmin)を計算したときに力がアンダーフローしないように座標をスケールして計算を行う。この作業は、g5_set_range(xmax,xmin,mmin) で行われ、この際に mmin を素質量とする離散的な対数質量のテーブルの生成も行われる。
内部で固定小数点を使っている関係上、このスケールの元で近接相互作用がオーバーフローしない距離というものが存在し、それは rmin 〜 10^-7 (xmax-xmin) となる。その結果重力ソフトニング (e_i) にも下限が付くことになる。具体的には e_i > rmin sqrt(m_i/mmin) (あるいは次のように書くとなんか意味がありそうな気がする: mmin/rmin^2 > m_i/e_i^2) でなければならないと書いてある (この条件式はなんか変なんじゃないかという気がするが)。
さて、今ツリーを使っているので十分遠方での相互作用は当然質量 mmin よりも十分大きくなっているのでアンダーフローについてはもうちょっと緩い条件で大丈夫のはず。しかしながら、ライブラリ関数の中でスケーリングと対数質量テーブルの生成を同時に行ってしまっているため、mmin として適度に大きい質量を設定した後、mmin としてまっとうな数値を入れるという手続きはそのままでは出来ない。例えばあとから質量テーブルだけこっそり書き換えるとかで対応出来るのであろうか。
今時単精度 GRAPE はツリーと合わせて使うのが大前提なわけだから、ライブラリの方で自発的に対応すべきだと私は思いますがそのつもりはないのかしら。
ついにハローの粒子数を増やすと良い感じに。。。もうちょっと調整したらガスを入れてみよう。
ディスクのスケール長をちょっと大きめにとると割と良い傾向に。ちょっとってでも二倍とか。逆に言うとバルジを入れるとかハローを重くするとかで対応出来るのかな。
いろいろがんばっている、新たに調べたりしない。という感じか。
2-2に対する返答と2-3に対する返答は結局整合性がないんじゃないだろうか。人件費削って、定年延長して、助教人事などが滞らないわけがない気がするが。
精神的余裕があるときにもうちょっとまじめに読んでみよう。
お金も持っていない所詮一万人程度の処遇、なので後回しになりがちなのは解るんだが、もうちょっと実態調査ぐらいしてみたらいいと思われる。
野暮用があってアメリカに電話。「キーボードを見ろ、その数字の上の記号があるだろう?それだ」とか、一生懸命教えてもらう。でも僕のキーボード無刻印。よく分からないけど何となく相づちを打って解ったという。
普段口から出して使わないので簡単な言葉も出てこない。いかんな。
一日に読めるメールの量。
パラメータ変えるとすぐ nan をだしやがる。使えない。なんで nan になっているのか、出てこないんだものなぁ。
そんなブラックボックス使うなという話である。
Disk の打ち切り半径を大きくしたらうまくいった。これはそういうエラーを出すべきなんじゃないだろうか。
二名募集。
こういうネーミングセンスを何とかしないと人生がおぼつかない。
ちなみに、当初思っていたのと微妙に違う振る舞いに見えるのだが、それは相対速度が足りないからだろうか。ということを、調べる前に此処に書いてしまうのはいかがなものかと。
相対速度は十分でござった。正面衝突じゃないからか。ほとんどあり得ん設定で出しても仕方あるまい。
あとは、ひたすら粒子数をつぎ込む。計算機が足りないなぁ。どっかに100台ぐらいあまってないかしら。
あかん、表示領域を間違えていた。つまりこれは単純に二つの銀河のダイナミクスを見ていたんだな。これはまぁそういうことなのかな。いまいち面白味に欠けるけどまぁこれはこれで。
面白いところは最終フェーズなのか?だとバルジ有りにしたいなぁ。でもこれでとりあえず。
ディスクの向きを変えるとかも、いずれは考えるべきだな。いま、完全に同じ面にあるから。x-zでみるときもい。
あ、もういっこ、どうしたらいいのか解らないので適当に決めたパラメータがあった。これどうすっかなぁ。
最後まで見ると結局「私と同じ引き出しを持っていなきゃ嫌だ」といっているだけなんじゃないのかなと。
いやにくわしい。
デスマーチ/
Wolfram|Alpha/
IT戦士 岡田有花リンク集/
Flash Earth/
Garage/
SDSS/
GOLDMine Galaxy On Line Database Milano Network/
Peanuts/
男達が後ろから撃たれるスレ 衛生兵を呼べ
tt-craft/
トムとジェリーの夏休み/
orz/
法治国家つまみぐい/
Happymetal ハッピーメタル/
MARK A. GARLICK'S SPACE-ART.CO.UK
ここはジョークアベニューです
ペリー
(:3)っ⌒つ 幸 [http://p2.chbox.jp/read.php?host=live23.2ch.net&bbs=livetb..]
(:3)っ⌒つ たか [おお、すばらしい。どうもありがとうございます。]