Stability of Planetary Orbital Configurations Klassifier
Let's predict the probaility that a given 3-planet system is stable:
import rebound
from spock import StabilityClassifier
model = StabilityClassifier()
sim = rebound.Simulation()
sim.add(m=1.)
sim.add(m=1.e-5, P=1., e=0.03, l=0.3)
sim.add(m=1.e-5, P=1.2, e=0.03, l=2.8)
sim.add(m=1.e-5, P=1.5, e=0.03, l=-0.5)
sim.move_to_com()
model.predict_stable(sim)
>>> 0.011536411
The best place to start is the example notebooks in jupyter_examples/
pip install spock
SPOCK relies on XGBoost, which has installation issues with OpenMP on Mac OSX. If you have problems (dmlc/xgboost#4477), the easiest way is probably to install homebrew, and
brew install libomp
pip install spock