Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add covariance propagation / mapping via sampling #82

Merged
merged 5 commits into from
Oct 22, 2023

Conversation

moeyensj
Copy link
Member

This is based on #81. Similar to propagating orbital covariances, this PR allows orbital covariances to be mapped to ephemerides via sampling.

For the time-being I have set the default method to monte-carlo, the sigma point code needs to be investigated for optimal parameter selection. There was an array shape bug that was fixed in d9dab91, that I think impacts the way the random draws were working.

Copy link
Contributor

@akoumjian akoumjian left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

No substantial feedback on the approach. Just some comments on variables/typing and some updated doc strings.

adam_core/orbits/variants.py Outdated Show resolved Hide resolved
adam_core/orbits/variants.py Outdated Show resolved Hide resolved
adam_core/orbits/variants.py Outdated Show resolved Hide resolved
adam_core/orbits/variants.py Show resolved Hide resolved
adam_core/propagator/pyoorb.py Show resolved Hide resolved
orbit_ids_idx
]

ephemeris = Ephemeris.from_kwargs(
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It seems like we are retyping this variable multiple times? First whatever comes back from pyoorb, then a 2D numpy array, and now our Ephemeris / VariantEphemeris classes. Can we split them out into separate variables? Especially from mutable to immutable.

Copy link
Member Author

@moeyensj moeyensj Oct 17, 2023

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Done. Spit it into ephemeris_array_3D, ephemeris_array, ephemeris, and variant_ephemeris (last two being the qv tables) .

@moeyensj moeyensj force-pushed the jm/variant-ephemeris branch 2 times, most recently from 5bd3584 to 9e91de2 Compare October 17, 2023 15:50
@moeyensj moeyensj requested a review from akoumjian October 17, 2023 15:51
@moeyensj moeyensj force-pushed the jm/return-ephemeris branch from 46d1706 to 8f7d9a4 Compare October 17, 2023 16:04
Base automatically changed from jm/return-ephemeris to main October 17, 2023 16:09
@moeyensj moeyensj force-pushed the jm/variant-ephemeris branch from 9e91de2 to 892d076 Compare October 17, 2023 16:10
@moeyensj moeyensj force-pushed the jm/variant-ephemeris branch from 892d076 to 905b550 Compare October 17, 2023 16:15
@moeyensj moeyensj merged commit 44e6642 into main Oct 22, 2023
@moeyensj moeyensj deleted the jm/variant-ephemeris branch October 22, 2023 18:36
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants