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

Full amplitude as matrix #36

Merged
merged 27 commits into from
Jun 30, 2024
Merged

Full amplitude as matrix #36

merged 27 commits into from
Jun 30, 2024

Conversation

mmikhasenko
Copy link
Owner

@mmikhasenko mmikhasenko commented Jun 16, 2024

Closes #28
Closes #27

Speed up (x60)

Working with matrices speeds up the code by a factor of x40, 🎉
Another 20% one gains using @tullion for index loops.

Test

The Intervention let do a finding that test coverage on full model functionality was pure.
New tests with 20 chains of different k improve the coverage.

Breaking

  1. Note, the PR might modify how SymbolicThreeBodyDecays.jl work due to
    one_T = one(typeof(two_js[1]))
    d_norm = sqrt(two_j * one_T + 1)

that used to be from two_λs

  1. amplitude with two arguments is used for full_amplitude computation, but also for dalitz-plot point evaluation. The conflict is resolved by specifying the type of sigma for the first call.

It's not ideal, since SymbolicTBDs.jl won't work again

@mmikhasenko mmikhasenko merged commit 9725dcb into main Jun 30, 2024
9 of 10 checks passed
@mmikhasenko mmikhasenko deleted the wigner-d-speed-up branch June 30, 2024 22:32
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.

[FR] Add performance benchmark [FR] Speed up by working with helicity matrices
1 participant