I don’t think the code is wrong. Section 3.1, table I:
http://www.plasmafocus.net/IPFS/2010 Papers/saw&leeScalingIJER-doi=10.1002_er.pdf
Check the value for 34KJ. The configuration is somewhat close to the FoFu-1 solution inputed into Lee’s code and indeed the yield is 3,46*10E9, close to what I got here:
https://focusfusion.org/index.php/forums/viewreply/7567/
It seems neutron yields fits very well the actual devices presented, not only on this paper, but on other recent ones too.
So, here is what I think: FoFu-1 achieves much higher densities and ion energy than expected. My bet it is that since the discharge happens by synchrony of several and not by closing a single switch. So, when electrons are free to move from every bank, the speed of light did not send the information to the others to delay the discharge. So, it is like all of banks are firing individualy. By doing simulations, as only 1 bank were discharging individually, the value is about as only 0.25MA, all of them together going to 2MA, at around 6-7Torr.
If this is true, if all of 12 them fire together, the yield will be 50%-100% higher than using tritium-deuterium or the old formula.