Rory Conlin
Rory Conlin
ok, im still a bit worried this is indicative of some larger issue. Can you check to see when it broke? Was it just because we changed the default m...
Please clean up the pr title and description before merging
shouldn't we just need to add `jitable=True` somewhere?
> If you want to get rid of some PRs, I can take this one. It looks interesting and I am bored:D Sure, main thing it needs is some tests,...
> Do you have a symplectic orbit solver to compare with? I have one that I can add to the test. It calculates the the exact particle position instead of...
Several people at the scidac meeting today are using the BMW code: https://github.com/ORNL-Fusion/BMW/tree/master refs: https://scipub.euro-fusion.org/wp-content/uploads/eurofusion/WPS1CPR17_17542_submitted-4.pdf https://bpb-us-e2.wpmucdn.com/wordpress.auburn.edu/dist/5/118/files/2022/05/ISHW2017.pdf basically it does a biot-savart volume integral to compute A (not B), then takes...
I would think that would have the opposite effect? Usually stuff like `jnp.sum(a*b)` is done in full precision, but `jnp.dot(a,b)` can use lower precision for intermediate variables depending on the...
for me I get ``` mulsum error: 8.524274097196166241e-12 dot error: 7.400826775665169309e-12 ``` a bit lower but well within the expected range for floating point rounding error due to different order...
yeah you can take a look at `desc.transform.Transform` to see how we do it for toroidal stuff. I think you'd first do a regular matmul for the radial part, which...