(1) leapfrog法 - 「Kick-Drift-Kick法」をimplementするために、簡単なKepler問題を解く。 - 重力場：質量1e11 Msunの質点 - 初期条件：(x,y)=(10, 0) kpc, (v_x, v_y)=(0, 100) km/s [eccentricity=0.77] - 積分時間： 1 Gyr - 積分timestep：(a) 1 Myr [1…

mac miniの調子がすこぶる悪い。 (1)セーフモードで起動＝起動時に「Shift」を押しておく>>> 起動に2時間近くかかったりする。しばらくは普通に動くのだが、 (2)普通に起動>>> 起動にx時間かかる(x>1) (3)Apple Diagnostics=起動時に「D」を押しておく>>> 「…

data_0= numpy.loadtxt(inputFileName_0, delimiter=' ', comments='#', usecols=range(13)) データの行の末尾が半角スペースになってしまう場合は、読み込むカラムの個数を与える(この場合は13カラム)。

ネット上の例にあったもの0.1sec x 10 = 1secごとにパラパラアニメ$ convert -delay 10 ./Pic/*.png test.gif ----- [基本] 自分の場合$ convert -delay 1 1with_bar/*.png 1snapshot_with_bar.gif 例えば 1with_bar というフォルダに複数のpngファイルがあ…

googleの検索バーに =(1kpc/km)/(yr/sec)/10^9 と入力すると、0.97781310637を得る。 すなわち、 1kpc/(km/s)=0.978 Gyr

(1) fitsio $pip install fitsio (2) esutil [ https://github.com/esheldon/esutil ] $ pip install esutil (3) isodist [ ] $ python setup.py build$ python setup.py install (4) periodictable-1.4.1 [ https://pypi.python.org/pypi/periodictable ] $…

Pythonの「ラムダ文」が難しいｗw import scipy.optimize as op def lnprob(theta, ell_data, bee_data): lp = lnprior(theta) if not np.isfinite(lp): return -np.inf return lp + lnlike(theta, ell_data, bee_data) theta=[5.05, 250., -200., 1.1] #探…

(1)自分の座標系 太陽：(xGC,yGC,zGC)=(-R0, 0, 0) LSR: (vxGC,vyGC,vzGC)=(0, LSR, 0) (2)galpyの座標系 太陽：(xGC_B,yGC_B,zGC_B)=(R0, 0, 0) LSR: (vxGC_B,vyGC_B,vzGC_B)=(0, -LSR, 0) (3)近いうちに、座標系をもう一度確認しよう x_galpy = -xGC vx_ga…

macbook air (1)emcee (2)emcee test(3)nose (4)Cython(5)triangle_plot(6)astropy(7)scikit-learn(8)astroML(9)healpy(10)pymc (1)emcee install #air の場合、 pip ではうまくいかなかった。#mini はpipでOK $ easy_install emcee Searching for emceeRead…

英国の大御所教授からメールが来た。 教授は３月に自分が投稿した論文の査読者になったのだが、あいにくその教授は俺の論文の仮定の甘さを理由に 「こんなひどい論文を通す気にはなれない」 と真っ向から否定された。 知り合いが以前、こういっていた。 「自…

個人用メモ。持ち物をまとめる： [最低限] ・パスポート ・財布（お金/カード） ・眼鏡 [仕事] ・ポスター ・Macbook Air & 充電器 ・携帯電話 & 充電器 [=目覚まし] ・電源プラグのアダプタ ・筆記用具 ・USB [日常] ・e-ticket ・旅程表 ・腕時計 ・家のカ…

Imagemagickでいろんなことができるのだと学んだ。（１）convert pdf > eps or level-2 epsなどができる。 $ convert test.pdf test_level1.eps$ convert test.pdf eps2:test_level2.eps （例）圧縮率が違うとこんなに違う・pdf: 0.45MB・level-1 eps: 32MB…

紆余曲折を経て yosemite で convert コマンドが使えるようになった。 いろいろやったので何がうまくいったのかわからんが、最後にやった手順を残しておく。 （０）症状 symptoms$ convert test.pdf test.jpg-bash: convert: command not found $ convert -v…

日本語tex文書を作成するために設定。 以前mactex2012を使っていた頃に（設定を無闇にいじりすぎたせいか）日本語がことごとく文字化けするという現象になやまされていた。今回の設定で文字化けも解消。 （１）MacTex2014をインストール 1. MacTeX-2012をア…

c言語のメモリリークを調べるツールとしてcppcheckをmacに入れてみる。 macbook air OSX 10.7.5 Lion mac mini OSX 10.10 Yosemite の両方に入れた。 手順としてはhomebrewをインストールしてhomebrew経由でcoocheckをインストール。 (1)homebrew をインスト…

mac mini で gsl を使うためにやったこと。 OpenMPを使う前提で、gcc47 をMacports 経由で入れた。さらに、gsl のオプション -I/usr/local/include を省略するために.bash_profile に加筆をした。 (1) gcc47 のインストール $ gcc -vConfigured with: --pref…