I have used the Surface.cross_section algorithm on the entire QUASR database, and notice that very infrequently, the algorithm fails, see below
I rewrote the algorithm, making it simpler, to fix the bug. Rerunning the new algorithm on the entire dataset, I can't find any incorrect cross sections anymore.
Note that the jagged cross section where varphi=0 is because it uniformly samples the cross section in the Boozer theta angle, which may be insufficient as you can see.
I have used the
Surface.cross_section
algorithm on the entire QUASR database, and notice that very infrequently, the algorithm fails, see belowI rewrote the algorithm, making it simpler, to fix the bug. Rerunning the new algorithm on the entire dataset, I can't find any incorrect cross sections anymore.
Note that the jagged cross section where
varphi=0
is because it uniformly samples the cross section in the Boozer theta angle, which may be insufficient as you can see.