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/euler #325

Merged
merged 3 commits into from
Dec 15, 2023
Merged

Fix/euler #325

merged 3 commits into from
Dec 15, 2023

Conversation

whaeck
Copy link
Member

@whaeck whaeck commented Dec 14, 2023

An error was introduced so many years ago in the value of the Euler-Mascheroni constant when a number of physical constants were put in the phys.f90 file.

Luckily for us, the impact is limited. The constant is used in the calculation of the Coulomb wave functions - I think - so this will only be seen on nuclides that have charged particle channels in their resonance parameters (LRF=7 only). In the tests only test 20 is impacted as this is the only LRF=7 evaluation we use for our non-regression tests (again, I think). The elastic scattering and capture cross section only have a few values where the 6th digit after the decimal point is different. In MT103 and MT600 (the n,p reaction calculated from the resonance parameters), the impact is only really visible in low cross section values where the 4th to 6th digit after the decimal point are different. This is reflected in the few changes in the total cross section.

Copy link
Contributor

@kahlerac kahlerac left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good ... and sadly we had the correct value(s) when this was explicitly typed in multiple code locations! It was only when moving to phys.f90 that the typo occurred, :(.

@whaeck
Copy link
Member Author

whaeck commented Dec 15, 2023

I know. I verified it to see when it was introduced.

@whaeck whaeck merged commit 5c4042a into develop Dec 15, 2023
@whaeck whaeck deleted the fix/euler branch December 15, 2023 07:40
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.

2 participants