Adrian Price-Whelan

Slicing an Orbit with multiple orbits?

4

I used `evaluate_orbit` with a `PhaseSpacePosition` with shape `(2,)`, then sliced the output orbit: ```python pot = gp.NFWPotential.from_M200_c(Quantity(1e12, "Msun"), 10.0, units="galactic") w0 = gc.PhaseSpacePosition( q=Quantity([[5.0, 0.0, 0.0], [5.05, 0.0, 0.0]],...

Add dispatched function to plot all components of an orbit

4

This adds a dispatch for `plot_components` when no components are specified to plot all pairs of components for an orbit. I had to change the order of arguments to squeeze...

When #727 is merged, do an audit of parameter names to make sure they are consistent with the new `CONVENTIONS.md` file.

Add Zhao potential model

6

This adds the (alpha,beta,gamma) Zhao (1996) double power law model. This replaces #755. Fixes #726. ~~The tests are currently failing but will pass once #760 is merged.~~

Add a Leapfrog solver for diffrax

1

Note that there is an open (stalled) PR to add a similar Verlet integrator to diffrax, but not sure why it hasn't been merged: https://github.com/patrick-kidger/diffrax/pull/303/ In the meantime, since we...

[WIP] Add a default time value

1

I think most potential methods should have a default time value so that you can call, e.g., ```python pot.gradient(xyz) ``` without having to specify a time. I marked this as...

Allow constructing spherical potentials with more constructors

1

All spherical models should allow initializing with: - mass enclosed at a reference radius - circular velocity at a reference radius - ... any others?

This adds a new potential class to represent a spherical Gaussian density profile. cc @mlucey

Make builtin potential list alphabetical (and simplify public API export??)

2

This first sorts the list of potential names alphabetically in the builtin `__init__.py`. I then also was puzzling at why we need another manual list of potential names in `galax/potential/__init__.py`...

For example, to compute for arrays of position and velocity (and time). Failing example: ```python gd.cluster.lagrange_points( pot, orbit_q, orbit_p, u.Quantity(1e5, "Msun"), t_grid ) ``` where `orbit_q.shape == orbit_p.shape == (N,...