Skip to content

Instantly share code, notes, and snippets.

View hannorein's full-sized avatar

Hanno Rein hannorein

View GitHub Profile
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
<!doctypehtml><html lang=en-us><head><meta charset=utf-8><meta content="text/html; charset=utf-8"http-equiv=Content-Type><title>REBOUND Visualization</title><style>body{font-family:arial;margin:0;padding:none;background:#000}.emscripten{padding-right:0;margin-left:auto;margin-right:auto;display:block}div.emscripten{text-align:center}#output{position:absolute;top:80%;left:0;margin:0;width:100%;height:20%;overflow:hidden;display:block;border:0;padding:0;background-color:#000;color:#fff;font-family:'Lucida Console',Monaco,monospace;outline:0}.reboundlogo{position:relative;top:.15em;float:left;-webkit-filter:invert(1);filter:invert(1);width:1em;height:1em;margin-right:.5em;background-size:contain;background-image:url('data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAACAAAAAgCAYAAABzenr0AAAAAXNSR0IArs4c6QAAAERlWElmTU0AKgAAAAgAAYdpAAQAAAABAAAAGgAAAAAAA6ABAAMAAAABAAEAAKACAAQAAAABAAAAIKADAAQAAAABAAAAIAAAAACshmLzAAAElklEQVRYCZ2We4hVVRSHr5NO4iN7aGY+uj7yCULlf4aRmEVEQ5CSTJFlDmVGmkMw6TQpIv2RWRFFWTFIhZUE5h9R2lgWZopmMlQojGmYQj
@hannorein
hannorein / gist:2a069763cf114f66641c
Created May 14, 2014 22:08
Yes, it's that simple. Querying the Open Exoplanet Catalogue.
#!/usr/bin/python
import xml.etree.ElementTree as ET, urllib, gzip, io
url = "https://github.com/OpenExoplanetCatalogue/oec_gzip/raw/master/systems.xml.gz"
oec = ET.parse(gzip.GzipFile(fileobj=io.BytesIO(urllib.urlopen(url).read())))
# Output mass and radius of all planets
for planet in oec.findall(".//planet"):
print [planet.findtext("mass"), planet.findtext("radius")]
# Find all circumbinary planets
@hannorein
hannorein / sqrt7.ipynb
Last active September 22, 2023 22:09
Sqrt7 IAS15
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@hannorein
hannorein / test.py
Created August 14, 2023 16:38
Garrett
sim = rebound.Simulation()
sim.integrator = "ias15"
sim.add(m=1.)
sim.add(m=1e-3, a=1., r=np.sqrt(1e-3/3.))
sim.add(m=5e-3, a=1.25, r=1.25*np.sqrt(5e-3/3.))
sim.collision="direct"
sim.collision_resolve = "halt"
E0 = sim.energy()
A0 = sim.angular_momentum()[2]
try:
tncmsunscramble = function (sInput) {
var sOutput = '';
for (var i = 0, c = sInput.length; i < c; i++) {
var nChar = sInput.charCodeAt(i);
if (nChar >= 33 && nChar <= 126) {
sTmp = String.fromCharCode(33 + (((nChar - 33) + 47) % 94));
sOutput += sTmp
} else {
sOutput += sInput.charAt(i)
}
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@hannorein
hannorein / 2014 UN271.ipynb
Created June 20, 2021 23:29
2014 UN271 with REBOUND
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.