コード例 #1
0
ファイル: gl_analytic.py プロジェクト: sofiasi/darcoda 
def nrgaiatot_3D_deriv(rad):
    lrho = np.log(rhogaiatot_3D(rad))
    lr   = np.log(rad)
    import gl_helper as gh
    return -gh.derivcoarse(lrho, lr)
コード例 #2
0
ファイル: gl_disc_simple.py プロジェクト: sofiasi/darcoda 
def get_kzpars():
    Kz_zthd = -2.0 * F * zth              # [with pc]
    Sigz_zth = abs(Kz_zthd) / (2.0*np.pi*gp.G1)
    denth = gh.derivcoarse(Sigz_zth,zth)
    kzpars = gh.ipol(zth,abs(denth),gp.xipol)*(2.0*np.pi*gp.G1)
    return kzpars