You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
To Reproduce
The bug is subtle and will not cause a crash, but may explain some behaviour seen in Fig.7 of Pham et al. (2024).
Moreover, it only features in with a non-default time-step criterion for IAS15.
The text was updated successfully, but these errors were encountered:
Thanks for opening an issue on this. It is a bit subtle indeed. First, just to avoid confusing, note that all the quantities involved in this line are squared (timescale2 is the timescale squared, y3 is the third derivative squared, etc, see here). So I believe we are talking about the difference between
I agree that they are not the same, but the difference is similar to choosing a different norm, which mathematically and for any physical scaling argument shouldn't make a difference. We comment a bit on this in the paper at the end of Section 2. We've tried out various different norms as well as the two different equations above to combine the derivatives into one scalar quantity. We did not notice any systematic difference. Re-running the test shown in Fig 7 for both of the above equations shows the same: both equations give the very similar results. It is not obvious to me why one would be preferred over the other.
Environment
version 4.3.2, C code, MacOS
Describe the bug
src/integrator_ias15.c line 591 reads
but correct is
![image](https://private-user-images.githubusercontent.com/37040636/310546111-b8443747-079b-44de-828b-781234766aac.png?jwt=eyJhbGciOiJIUzI1NiIsInR5cCI6IkpXVCJ9.eyJpc3MiOiJnaXRodWIuY29tIiwiYXVkIjoicmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbSIsImtleSI6ImtleTUiLCJleHAiOjE3MTgwNTI3OTksIm5iZiI6MTcxODA1MjQ5OSwicGF0aCI6Ii8zNzA0MDYzNi8zMTA1NDYxMTEtYjg0NDM3NDctMDc5Yi00NGRlLTgyOGItNzgxMjM0NzY2YWFjLnBuZz9YLUFtei1BbGdvcml0aG09QVdTNC1ITUFDLVNIQTI1NiZYLUFtei1DcmVkZW50aWFsPUFLSUFWQ09EWUxTQTUzUFFLNFpBJTJGMjAyNDA2MTAlMkZ1cy1lYXN0LTElMkZzMyUyRmF3czRfcmVxdWVzdCZYLUFtei1EYXRlPTIwMjQwNjEwVDIwNDgxOVomWC1BbXotRXhwaXJlcz0zMDAmWC1BbXotU2lnbmF0dXJlPWNhZjQ0MTI0ZWRkNWQ4OGI1YmM2ZmE5MTI5OGY4OGRhZmQxNGJiYmVjOTI0NWZiYmQxZjU2YjNhY2I4YjQyNzgmWC1BbXotU2lnbmVkSGVhZGVycz1ob3N0JmFjdG9yX2lkPTAma2V5X2lkPTAmcmVwb19pZD0wIn0.koZzsoY18gJI-1X25dexuTDFb8-QCGG36Pp0WZwamXs)
(equation 6 of Pham, Rein, Spiegel 2024), i.e.
To Reproduce
The bug is subtle and will not cause a crash, but may explain some behaviour seen in Fig.7 of Pham et al. (2024).
Moreover, it only features in with a non-default time-step criterion for IAS15.
The text was updated successfully, but these errors were encountered: