Erol Balkovic

@ludgerpaehler Is this something you're still working on? Would be fantastic to have this feature, even if it's not polished...

Alright, thanks. I am unfortunately not knowledgable enough to help, but I hope someone will take up the task.

Thanks for the work @missing-user. I checked your code and have also successfully ran an axisymmetric free-boundary problem [input_file](https://github.com/PrincetonUniversity/SPEC/blob/master/InputFiles/TestCases/G3V08L3Fr.001.sp). The result matches what we had before and virtual casing calculation...

In the end, I ended up using print inside the njitted function, and contextlib.redirect_stdout before the call to the numba function. This seemed to be the fastest and allowed me...

Hi @Taegeun-Jeong, I would like to point out a couple of things that may help your understanding. When you run SPEC with Nvol=1 you are simply solving for the Beltrami...