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

fix vrj reinitialization bug #450

Merged
merged 8 commits into from
Sep 2, 2024
Merged

Conversation

isaacsas
Copy link
Member

@isaacsas isaacsas commented Sep 2, 2024

It seems that the code to reset variable rate jump problems was not correctly broadcasting calls to randexp resulting in the same samples being used for all jumps.

cc: @ChrisRackauckas

Closes #320

@isaacsas
Copy link
Member Author

isaacsas commented Sep 2, 2024

Now that we have a Catalyst VRJ interface I think this is indicative that we really need to add more VRJ correctness tests.

Even with this fix #320 still doesn't seem to work.

edit: Actually I think this fixes that issue. When plotting the error in the mean over 10000 simulations I now see that the jump mean can indeed now sometimes exceed the analytical mean.

@isaacsas
Copy link
Member Author

isaacsas commented Sep 2, 2024

If tests pass this should be good to go (and we should then make an immediate release).

@ChrisRackauckas ChrisRackauckas merged commit e2a7af1 into SciML:master Sep 2, 2024
4 checks passed
@isaacsas isaacsas deleted the vrj_u0_bugfix branch September 2, 2024 02:10
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.

VariableRateJump - Jumps don't match expected distribution for birth-death process
2 participants