Rewrite to cut down the stream function evaluations in numerical root finding to a third