Skip to content

Instantly share code, notes, and snippets.

Avatar
🎯
Focusing

João Faria j-faria

🎯
Focusing
View GitHub Profile
@j-faria
j-faria / matrix_exponential.cpp
Last active Feb 9, 2021
C++ and pybind11 code for matrix exponential
View matrix_exponential.cpp
#include <iostream>
#include <mkl.h>
#include <math.h>
#include <pybind11/pybind11.h>
#include <pybind11/numpy.h>
namespace py = pybind11;
using namespace pybind11::literals;
using namespace std;
@j-faria
j-faria / howto.md
Last active Oct 7, 2020
Auto-reload modules in IPython
View howto.md
  1. Create a file ~/.pythonrc.py
  2. Add
try:
    from IPython import get_ipython
    get_ipython().run_line_magic('load_ext', 'autoreload')
    get_ipython().run_line_magic('autoreload', '2')
except AttributeError:
    pass
View context_chdir.py
import contextlib
import os
@contextlib.contextmanager
def working_directory(path):
"""
A context manager which changes the working directory to the given
path, and then changes it back to its previous value on exit.
"""
prev_cwd = os.getcwd()
@j-faria
j-faria / ExoPlanetCurves
Last active Dec 5, 2019
The Keplerian and transit curves
View ExoPlanetCurves
Two small demonstration notebooks of how the Keplerian and transit curves respond to a change in the parameters.
@j-faria
j-faria / eso_access_phase3.sh
Created Sep 6, 2019
query and download reduced data from the ESO Archive
View eso_access_phase3.sh
#!/bin/sh
#***********************************************************************
#* ESO Science Archive Facility
#* Programmatic Access
#*
#* Script: eso_access_phase3.sh
#* Shell: bash
#* Date: 15-Jul-2015
#* Contact: archive@eso.org
#* Description: Script to query and download reduced data (Phase 3)
@j-faria
j-faria / BayesCorr.py
Created Sep 10, 2015
Bayesian correlation
View BayesCorr.py
#!/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:
@j-faria
j-faria / kepler.cpp
Created Jan 16, 2017
C++ implementation of "A Practical Method for Solving the Kepler Equation"
View kepler.cpp
/**
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.
*/
@j-faria
j-faria / journal_abbreviations.py
Created Jul 2, 2019
Simplified abbreviations of frequently used journals (according to A&A)
View journal_abbreviations.py
# https://www.aanda.org/67-author-information/frequent-abbreviations
name_abbrv = {
'Astronomy and Astrophysics': 'A&A',
'Astronomy & Astrophysics': 'A&A',
'Monthly Notices of the Royal Astronomical Society': 'MNRAS',
'The Astrophysical Journal': 'ApJ',
'The Astronomical Journal': 'AJ',
'Publications of the Astronomical Society of the Pacific': 'PASP',
'Annual Review of Astronomy and Astrophysics': 'ARA&A' ,
@j-faria
j-faria / parameters_torres.py
Last active May 7, 2019
Calculate stellar mass and radius using the Torres calibration
View parameters_torres.py
import numpy as np
try: # numba will provide a ~2x speedup
from numba import jit
except ImportError: # but we can do without it
jit = lambda fn: fn
@jit
def massTorres(teff, erteff, logg, erlogg, feh, erfeh,
ntrials=10000, corrected=True, add_intrinsic=True):
"""