Messed up the order of the loop which fills the exparg table.  With this correct, the exparg loops can be in the same energy loop as the actual phi solves.  This prevents the need to store a 3D exparg array for (with energy being the 3rd dimension).
A bunch of random improvements to the energy group being the inner loop.  Removing some arrays that could be replaced by scalars, or were unnecessary altogether with the loop reorder.  This results in some speedup.
Tried to move the final loop where the phis get scaled by xstr and qbar so that it's inside the other group loop.  Commented out, because I can't come up with a way to do it right now that doesn't increase runtime by 2X or more.