docs + comments

examples/nbody/README.md

@@ -8,7 +8,7 @@ The lomb-scargle periodogram from our ephemeris fitting code can be used to esti
 
 ![](ttv_prior.png)
 
-The example code has recently been refactored from [Nbody-AI](https://github.com/pearsonkyle/Nbody-ai) and is still in development. Some retrievals not work entirely... please be paitent and report any issues to our slack. 
+The example code has recently been refactored from [Nbody-AI](https://github.com/pearsonkyle/Nbody-ai) and is still in development. Some retrievals not work entirely... please be paitent and report any issues to our slack.
 
 If you make use of this code please cite:
 

examples/nbody/example.py

@@ -139,5 +139,5 @@
 # run the fitter
 nfit = nbody_fitter(data, nbody_prior, nbody_bounds)
 
-print(nfit.parameters)
-print(nfit.errors)
+# print(nfit.parameters)
+# print(nfit.errors)

exotic/api/nbody.py

@@ -50,7 +50,6 @@
 from ultranest import ReactiveNestedSampler
 from astropy.io import fits
 from astropy import units as u
-from astropy import constants as c
 from scipy.interpolate import interp1d
 from astropy.timeseries import LombScargle
 from scipy.signal import find_peaks
@@ -609,7 +608,8 @@ def prior_transform(upars):
         self.quantiles = {}
         self.parameters = copy.deepcopy(self.prior)
 
-        # TODO finish
+        # TODO finish saving final results + posteriors
+
         #for i, key in enumerate(freekeys):
         #    self.parameters[key] = self.results['maximum_likelihood']['point'][i]
         #    self.errors[key] = self.results['posterior']['stdev'][i]
@@ -750,5 +750,5 @@ def prior_transform(upars):
     # run the fitter
     nfit = nbody_fitter(data, nbody_prior, nbody_bounds)
 
-    print(nfit.parameters)
-    print(nfit.errors)
+    # print(nfit.parameters)
+    # print(nfit.errors)