2x acceleration by using last u_factor as initial value for next iteration. 1% speed-up by unifying loops over streamtubes
