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

Time in guiding center #213

Open
Beforerr opened this issue Dec 23, 2024 · 1 comment
Open

Time in guiding center #213

Beforerr opened this issue Dec 23, 2024 · 1 comment

Comments

@Beforerr
Copy link
Contributor

Beforerr commented Dec 23, 2024

What does xu in guiding_center represent for? Because in example we have gc_x0 = gc(stateinit) then xu[end] would be uz not t

function guiding_center(xu, param::FullTPTuple)
   q, m, _, B_field = param
   t = xu[end]
...
end
@henry2004y
Copy link
Owner

In guiding_center, I remember xu is supposed to be a 7-element vector (x, y, z, vx, vy, vz, t). However, in the current implementation the time information is not used at all, so a 6-element vector (x, y, z, vx, vy, vz) also works, but the last element xu[end] actually points to vz as you mentioned. This is indeed wrong.

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

No branches or pull requests

2 participants