Skip to content

Instantly share code, notes, and snippets.

View hannorein's full-sized avatar

Hanno Rein hannorein

View GitHub Profile
import h5py
import numpy as np
import rebound
def hdf5_to_rebound(row):
row = row.reshape((-1,5))
Pmin = np.min(row[:,1])
TCmin = np.min(row[:,4])
epoch = TCmin - 0.1 * Pmin
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
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 / sqrt7.ipynb
Last active September 22, 2023 22:09
Sqrt7 IAS15
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
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)
}
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.