mabarnes

Re: Bob's theory about the mirror term and finite vpa grid, do you find that increasing vpa_max alleviates the initial phi blowup? If the mirror term is all that's in...

Yes, it’s because in stellarators with multiple wells the velocity space and field line coordinate would get coupled together in a complicated way if using energy rather than vpa

in case it's helpful, I've created a new branch (development/delphi_response), in which one has the option to choose between using phi^n+1 and (phi^n+1-phi^n) when solving for g_inh in parallel streaming....

In case it is helpful, I have started looking at the CBC with ky=0.001 (using the delphi_response branch, but with options set to mimic the master branch). Here is what...

agree with removal and support anyone interested in implementing slab geometry as an option

so should we interpret this as meaning that with the slight tweak to the centering in the EM implementation, the code behaves in the same way as the master branch?...

Happy to go with this auto-formatter. I do have a slight preference for 3 spaces over 2, possibly because that's just what I've been used to with emacs auto-formatting.

Is the proposal that we make the zonal modes periodic at the ends of the flux tube train rather than within each car? I'm not sure which approach is preferable...

the ia dimension is indeed a placeholder for the full-flux-surface version of the code, where alpha-dependence will be present. Possible solutions include the introduction of a dedicated type or the...

I don't have a strong opinion. I decided to separate the stella utils from the gs2 utils at the beginning to eliminate a lot of the complication that wasn't/isn't needed...