wrote formColumnsHB, it compiles but doesn't run. Also a previous commit seems to have changed something about phase 1. Now the iterations converge slower and the first threshold is at D = 0.28, not D = 0.32. Need to go back and fix this.