fixed one small typo (don't need to zero hb in ComputePsisqCommon) and one MAJOR one: for some reason didn't square the c in the computehb phase!