Skip to content

Conversation

@pranav-satheesh
Copy link
Collaborator

Description

Merging the updated cbd hardening model to the dev branch
This branch uses an entirely different evolution class: New_Evolution

Todos

Notable points that this PR has either accomplished or will accomplish.

  • [ ]

Questions

  • Question1

Status

  • Ready to go

mssiwek and others added 30 commits June 21, 2023 14:01
Use time (lookback time) as independent variable instead of binary separation to allow for binary 'softening' (i.e. da/dt > 0).
Integrate each binary, one at a time, from initial conditions until either coalescence or redshift zero.  Evolution arrays thus
are no longer always the same size.  Currently evolution arrays are 1D objects, flattened for all binaries, and the appropriate
indices are stored (Evolution._last_index).  This *seems* to be working (first full implementation).  Lots more changes will be
needed to allow for interpolation, universe-realization, etc.  Its also extremely slow.  There are numerous obvious problematic
places, e.g. interpolation of CBD_Torques and Dynamical Friction halo relations.  Cython could also speed things up a lot, but
the issue will be external function calls.

Really, this should be a second Evolution class, so that both implementations are preserved for the moment... I'll add the old
one back in later.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

enhancement New feature or request

Projects

None yet

Development

Successfully merging this pull request may close these issues.

5 participants