This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#!/usr/env/python | |
# -*- coding: utf-8 -*- | |
# Written by Pedro Figueira. | |
# Original version can be found at https://pedrofigueira@bitbucket.org/pedrofigueira/bayesiancorrelation | |
# Modified by João Faria | |
"""Bayesian Correlation. | |
Usage: |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import numpy as np | |
import cStringIO | |
filename = 'example.rdb' | |
def read_rdb(filename): | |
""" Reads a .rdb file with possible comments '#' and header | |
col1 col2 col3 | |
---- ---- ---- | |
Returns a numpy record array |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import sys | |
from astroquery.eso import Eso | |
eso = Eso() | |
# should use the default username | |
# and the keyring stored password | |
authentication = eso.login() | |
if not authentication: | |
raise RuntimeError('Something went wrong with authentication!') |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#!/usr/bin/env python | |
""" | |
Convert a PDF file to black&white and display it. | |
Run as: python bw.py input.pdf | |
""" | |
viewer = 'evince' | |
from contextlib import contextmanager | |
import tempfile |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
def fun(teff=None, mass=None): | |
if teff and mass: | |
# gave two parameters | |
pass # goto 1 | |
elif teff: | |
# don't have mass, have teff | |
mass = f(teff) | |
elif mass: | |
# don't have teff, have mass | |
# teff = f(mass) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
/** | |
Calculates the eccentric anomaly at time t by solving Kepler's equation. | |
See "A Practical Method for Solving the Kepler Equation", Marc A. Murison, 2006 | |
@param t the time at which to calculate the eccentric anomaly. | |
@param period the orbital period of the planet | |
@param ecc the eccentricity of the orbit | |
@param t_peri time of periastron passage | |
@return eccentric anomaly. | |
*/ |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
ssh -N -f -L port:server2:22 user@server1 | |
# mount the directory locally | |
sshfs -C -p port user@localhost:remote_directory local_directory(must exist) | |
# unmount | |
sudo umount local_directory |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
# matplotlib parameters for publication plots | |
import matplotlib.pyplot as plt | |
import matplotlib.ticker as ticker | |
from matplotlib import gridspec | |
import re | |
import os | |
this_folder = os.path.dirname(__file__) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
import numpy as np | |
from scipy.optimize import bisect | |
from scipy.special import binom as binom_coeff | |
from scipy.integrate import quad | |
from functools import partial | |
def binom_function(N, n, p): | |
c = binom_coeff(N, n) | |
return c * p**n * (1.-p)**(N-n) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
from astropy.constants import G | |
import astropy.units as u | |
# we assume m_planet << m_star | |
def get_planet_mass(P, K, e, star_mass=1.0): | |
""" | |
P in days, K in m/s, star_mass in solar masses | |
output is planet mass in Jupiter masses | |
""" |