hacked away phase1 discrepency, now addLasing doesn't call formH if c = 0, (will eventually need to resolve this subset_nonzero_pattern bullshit.

Also, I wrote formPsiMagSq and now formColumnsHB runs (but I'm 99.99999% sure it doesn't give the correct result. If I actually run this right now and the holeburning raph newton works, I'll be a monkey's uncle)