From 129d85dd225867797248c94185450882a0beeec7 Mon Sep 17 00:00:00 2001 From: Kyle Pearson Date: Tue, 24 Oct 2023 14:31:30 -0700 Subject: [PATCH] docs + comments --- examples/nbody/README.md | 2 +- examples/nbody/example.py | 4 ++-- exotic/api/nbody.py | 8 ++++---- 3 files changed, 7 insertions(+), 7 deletions(-) diff --git a/examples/nbody/README.md b/examples/nbody/README.md index c1b93395..c45fd839 100644 --- a/examples/nbody/README.md +++ b/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: diff --git a/examples/nbody/example.py b/examples/nbody/example.py index 21fc39b3..64c1dd16 100644 --- a/examples/nbody/example.py +++ b/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) diff --git a/exotic/api/nbody.py b/exotic/api/nbody.py index c3cf6e6d..f55a8ddd 100644 --- a/exotic/api/nbody.py +++ b/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) \ No newline at end of file + # print(nfit.parameters) + # print(nfit.errors) \ No newline at end of file