Skip to content

Fix constraints missing vperp term - #486

Open
LucasMontoya4 wants to merge 2 commits into
masterfrom
fix_constraints_missing_vperp_term
Open

LucasMontoya4 wants to merge 2 commits into
masterfrom
fix_constraints_missing_vperp_term

Conversation

@LucasMontoya4

Copy link
Copy Markdown
Collaborator

Made to fix Issue #464. Thanks John for noticing this mistake!

…that in moment_kinetic_2V_tests.jl, the changes to the final phi values introduced for the test case of evolve_p = true are roughly 5 times closer on average to the equivalent values of phi in the case where evolve_n,upar = true but evolve_p = false. This suggests that the changes here bring the evolve_p=true case more in line with the other evolve_moments cases, which is good. However, for the trapping tests, the phi values of the evolve_p = true case were not statistically any closer to the equivalent phi values for other evolve_moments cases. This suggests that another error is dominating in this case, perhaps too low a resolution in velocity space of the test, for example. I will write some high resolution 2V tests so that we can ensure that, like this moment_constraints.jl bug, any other code bugs are more likely to be exposed through disparity in test values of phi in the future, rather than being drowned out by the differences between evolve_moments cases that are naturally induced by low resolution.
@johnomotani johnomotani linked an issue Oct 30, 2026 that may be closed by this pull request

This branch has not been deployed

No deployments
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.

2V moment constraints bug?

1 participant