[ExtForces] Add pseudo-Newtonian Paczynsky-Witta potential#1882
[ExtForces] Add pseudo-Newtonian Paczynsky-Witta potential#1882y-lapeyre wants to merge 4 commits into
Conversation
y-lapeyre
commented
Jun 21, 2026
- check it on test problem
|
Thanks @y-lapeyre for opening this PR! You can do multiple things directly here: Once the workflow completes a message will appear displaying informations related to the run. Also the PR gets automatically reviewed by gemini, you can: |
There was a problem hiding this comment.
Code Review
This pull request introduces a new pseudo-Newtonian Paczynsky-Witta (1980) external force potential (AddForcePaczynskyWitta), integrating it into the CMake build, configuration structures, serialization, and Python bindings. The review feedback highlights critical mathematical and performance issues in the acceleration calculation, such as an incorrect sign for the Schwarzschild radius, an incorrect denominator, and loop-invariant computations inside the kernel. Additionally, the new force variant is unhandled in add_ext_forces() and point_mass_accrete_particles(), which would cause runtime crashes and prevent particle accretion. Finally, the LaTeX formula in the documentation needs to be corrected to accurately reflect the potential.
Important
The consumer version of Gemini Code Assist on GitHub is being sunset. Starting June 18, 2026, new organization installations will be blocked, and all code review activity will officially cease on July 17, 2026.
For more details on the timeline and next steps, please review the Help Documentation.
Workflow reportworkflow report corresponding to commit 121ee24 Pre-commit check reportPre-commit check: ✅ Test pipeline can run. Clang-tidy diff reportDoxygen diff with
|
Codecov Report❌ Patch coverage is 📢 Thoughts on this report? Let us know! |