Rory Conlin

also here's some ideas from last week: ``` """ A winding pack is defined by a width (w) and height (h), centered around a filamentary curve (x), carrying a total...

I agree this is something we should figure out longer term, but as @dpanici points out it's actually a non-trivial amount of work to get it to work well. FWIW,...

Also looks like |B| in the 2nd plot goes negative? and are the big white areas where its NaN or something?

It's probably fine to not mod things in the compute functions (or only mod them when being used in periodic ways) but we do generally need to mod both the...

From looking through STELLOPTs version of Gamma_C [here](https://github.com/PrincetonUniversity/STELLOPT/blob/653cae047b08b7fbfe4ac0166e63247bbe2f6199/STELLOPTV2/Sources/Chisq/chisq_gamma_c.f90#L2) I see a few possible issues with their implementation that might lead to differences: - Ntransits and Nzetapertransit hard coded to 400...

Yeah I looked at that as well but it seems the same? There are more comments explaining some of the coordinate transformations but the bounce points are still rounded to...

Ideally yes, I'd like to keep them separate, but we don't have to do everything at once. The main point was to not remove existing things and replace them with...

I could see something like this working, but I'd vote for moving the extra logic to inside the surface integral function, or a higher order wrapper for it. There may...

I like the idea of partial summation, I'm also wondering if we could figure out how to do it in `map_coordinates` when `inbasis` and `outbasis` share certain coordinates? If we...

I hadn't seen that before, but I'm not sure that completely solves the problem? You would need to evaluate lambda at a generally nonuniform grid in theta which without partial...